Search

BACKGROUND: For patients with peanut allergy, there are currently no methods to predict who will develop sustained unresponsiveness (SU) after oral immunotherapy (OIT). OBJECTIVE: Assess IgE binding to peanut (PN), Ara h 2 and to specific linear epitopes of Ara h 2 as predictors of the important clinical parameters: eliciting dose threshold and attainment of SU following OIT. METHODS: Samples and clinical data were collected from children undergoing OIT. PN- and Ara h 2-sIgE were quantified by ImmunoCap(®). IgE binding to linear peptides of Ara h 2 and Ara h 6 was measured with peptide microarrays. RESULTS: Lower values of PN-sIgE correlated with eliciting dose (P = 0.001) and with a higher likelihood of achieving SU (P < 0.0001) but these relationships were lost at higher values for PN-sIgE (≥ 14 kIU for eliciting dose and ≥ 35 kIU/L for SU). In subjects with PN-sIgE ≥ 14kIU/L, binding of IgE to epitopes 5 & 6 of Ara h 2 was associated with a lower eliciting dose at baseline challenge (P < 0.001; Pc < 0.02). In subjects with PN-sIgE ≥ 35kIU/L, a combined model of IgE binding to epitopes 1, 5 and 6 with PN-sIgE was highly predictive of attainment of SU (AUC of 0.86; P = 0.0067). CONCLUSION: In young patients with peanut allergy, measurement of PN-sIgE and IgE binding to specific linear epitopes of Ara h 2 in baseline samples may allow stratification of patients regarding sensitivity to challenge and outcome of OIT.

About half of the US population is sensitized to one or more allergens, as found by a National Health and Nutrition Examination Survey (NHANES). The most common treatment for seasonal allergic responses is the daily use of oral antihistamines, which can control some of the symptoms, but are not effective for nasal congestion, and can be debilitating in many patients. Peptide immunotherapy is a promising new approach to treat allergic airway diseases. The small size of the immunogens cannot lead to an unwanted allergic reaction in sensitized patients, and the production of peptides with sufficient amounts for immunotherapy is time- and cost-effective. However, it is not known what peptides are the most effective for an immunotherapy of allergens. We previously produced a unique monoclonal antibody (mAb) E58, which can inhibit the binding of multiple groups of mAbs and human IgEs from patients affected by the major group 1 allergens of ragweed (Amb a 1) and conifer pollens (Jun a 1, Cup s 1, and Cry j 1). Here, we demonstrated that a combined approach, starting from two linear E58 epitopes of the tree pollen allergen Jun a 1 and the ragweed pollen allergen Amb a 1, and residue modifications suggested by molecular docking calculations and peptide design could identify a large number of high affinity binding peptides. We propose that this combined experimental and computational approach by structural analysis of linear IgE epitopes and peptide design, can lead to potential new candidates for peptide immunotherapy.

There is an urgent need for inexpensive, rapid and specific antigen-based assays to test for infection with SARS-CoV-2 and distinguish variants arising as the COVID-19 pandemic spreads. We have identified a small, synthetic protein (JS7), representing a region of maximum variability within the receptor binding domain (RBD), which binds antibodies in sera from nine patients with PCR-verified COVID-19 of varying severity. Antibodies binding to either JS7 or the SARS-CoV-2 recombinant RBD, as well as those that disrupt binding between a fragment of the ACE2 receptor and the RBD, are proportional to disease severity and clinical outcome. Binding to JS7 was inhibited by linear peptides from the RBD interface with ACE2. Variants of JS7, such as N501Y, can be quickly synthesized in a pure form in large quantities by automated methods. JS7 and related synthetic antigens can provide a basis for specific diagnostics for SARS-CoV-2 infections.

Allergies to peanuts, tree nuts, and sesame seeds are among the most important food-related causes of anaphylaxis. Important clinical questions include: Why is there a variable occurrence of coallergy among these foods and Is this immunologically mediated? The clinical and immunologic data summarized here suggest an immunologic basis for these coallergies that is based on similarities among the 2S albumins. Data from component resolved diagnostics have highlighted the relationship between IgE binding to these allergens and the presence of IgE-mediated food allergy. Furthermore, in vitro and in vivo experiments provide strong evidence that the 2S albumins are the most important allergens in peanuts for inducing an allergic effector response. Although the 2S albumins are diverse, they have a common disulfide-linked core with similar physicochemical properties that make them prime candidates to explain much of the observed coallergy among peanuts, tree nuts, and sesame seeds. The well-established frequency of cashew and pistachio nut coallergy (64%-100%) highlights how the structural similarities among their 2S albumins may account for observed clinical cross-reactivity. A complete understanding of the physicochemical properties of the 2S albumins in peanuts, tree nuts, and sesame seeds will enhance our ability to diagnose, treat, and ultimately prevent these allergies.

Proteins are fundamental to life and exhibit a wide diversity of activities, some of which are toxic. Therefore, assessing whether a specific protein is safe for consumption in foods and feeds is critical. Simple BLAST searches may reveal homology to a known toxin, when in fact the protein may pose no real danger. Another challenge to answer this question is the lack of curated databases with a representative set of experimentally validated toxins. Here we have systematically analyzed over 10,000 manually curated toxin sequences using sequence clustering, network analysis, and protein domain classification. We also developed a functional sequence signature method to distinguish toxic from non-toxic proteins. The current database, combined with motif analysis, can be used by researchers and regulators in a hazard screening capacity to assess the potential of a protein to be toxic at early stages of development. Identifying key signatures of toxicity can also aid in redesigning proteins, so as to maintain their desirable functions while reducing the risk of potential health hazards.

There is a pressing need for new vaccines against alphaviruses, which can cause fatal encephalitis (Venezuelan equine encephalitis virus (VEEV) and others) and severe arthralgia (e.g. Chikungunya virus, CHIKV). These positive-strand RNA viruses are diverse and evolve rapidly, meaning that the sequence of any vaccine should cover multiple strains that may be quite different from any previous isolate. Here, consensus proteins were produced to represent the common physicochemical properties (PCPs) of the epitope rich, B domain of the E2 envelope protein. PCP-consensus proteins were based on multiple strains of VEEV (VEEVcon) and CHIKV (CHIKVcon) or the conserved PCPs of 24 different alphaviruses (AllAVcon). The AllAVcon was altered to include binding sites for neutralizing antibodies of both VEEV and CHIKV strains (Mosaikcon). All four designed proteins were produced solubly in E. coli and purified. They formed the β-strand core expected from experimental structures of this region of the wild type E2 proteins as indicated by circular dichroism (CD) spectra. Furthermore, the CHIKVcon protein bound to a structure dependent, CHIKV neutralizing monoclonal antibody. The AllAVcon and Mosaikcon proteins bound to polyclonal antibodies generated during natural infection with either VEEV or CHIKV, indicating they contained epitopes of both serotypes. The Mosaikcon antigen induced antibodies in rabbit sera that recognized both the VEEVcon and CHIKVcon spike proteins. These PCP-consensus antigens are promising starting points for novel, broad-spectrum alphavirus vaccines.

There is an urgent need for inexpensive, rapid and specific antigen-based assays to test for vaccine efficacy and detect infection with SARS-CoV-2 and its variants. We have identified a small, synthetic protein (JS7), representing a region of maximum variability within the receptor binding domain (RBD), which binds antibodies in sera from nine patients with PCR-verified COVID-19 of varying severity. Antibodies binding to either JS7 or the SARS-CoV-2 recombinant RBD, as well as those that disrupt binding between a fragment of the ACE2 receptor and the RBD, are proportional to disease severity and clinical outcome. Binding to JS7 was inhibited by linear peptides from the RBD interface with ACE2. Variants of JS7, such as E484 K or N501Y, can be quickly synthesized in a pure form in large quantities by automated methods. JS7 and related synthetic antigens can provide a basis for specific diagnostics for SARS-CoV-2 infections.

SummaryBackground Cross-linking of IgE antibody by specific epitopes on the surface of mast cells is a prerequisite for triggering symptoms of peanut allergy. IgE epitopes are frequently categorized as linear or conformational epitopes. Although linear IgE-binding epitopes of peanut allergens have been defined, little is known about conformational IgE-binding epitopes. Objective To identify clinically relevant conformational IgE epitopes of the two most important peanut allergens, Ara h 2 and Ara h 6, using phage peptide library. Methods A phage 12mer peptide library was screened with allergen-specific IgE from 4 peanut-allergic patients. Binding of the mimotopes to IgE from a total of 29 peanut-allergic subjects was measured by ELISA. The mimotope sequences were mapped on the surface areas of Ara h 2 and Ara h 6 using EpiSearch. Results Forty-one individual mimotopes were identified that specifically bind anti- Ara h 2/Ara h 6 IgE as well as rabbit anti-Ara h 2 and anti-Ara h 6 IgG. Sequence alignment showed that none of the mimotope sequences match a linear segment of the Ara h 2 or Ara h 6 sequences. EpiSearch analysis showed that all the mimotopes mapped to surface patches of Ara h 2 and Ara h 6. Eight of the mimotopes were recognized by more than 90% of the patients, suggesting immunodominance. Each patient had distinct IgE recognition patterns but the recognition frequency was not correlated to the concentration of peanut specific IgE or to clinical history. Conclusions The mimotopes identified in this study represent conformational epitopes. Identification of similar surface patches on Ara h 2 and Ara h 6 further underscores the similarities between these two potent allergens.

Ebolaviruses Zaire (EBOV), Bundibugyo (BDBV), and Sudan (SUDV) cause human disease with high case fatality rates. Experimental monovalent vaccines, which all utilize the sole envelope glycoprotein (GP), do not protect against heterologous ebolaviruses. Human parainfluenza virus type 3-vectored vaccines offer benefits, including needle-free administration and induction of mucosal responses in the respiratory tract. Multiple approaches were taken to induce broad protection against the three ebolaviruses. While GP consensus-based antigens failed to elicit neutralizing antibodies, polyvalent vaccine immunization induced neutralizing responses to all three ebolaviruses and protected animals from death and disease caused by EBOV, SUDV, and BDBV. As immunization with a cocktail of antigenically related antigens can skew the responses and change the epitope hierarchy, we performed comparative analysis of antibody repertoire and Fc-mediated protective mechanisms in animals immunized with monovalent versus polyvalent vaccines. Compared to sera from guinea pigs receiving the monovalent vaccines, sera from guinea pigs receiving the trivalent vaccine bound and neutralized EBOV and SUDV at equivalent levels and BDBV at only a slightly reduced level. Peptide microarrays revealed a preponderance of binding to amino acids 389 to 403, 397 to 415, and 477 to 493, representing three linear epitopes in the mucin-like domain known to induce a protective antibody response. Competition binding assays with monoclonal antibodies isolated from human ebolavirus infection survivors demonstrated that the immune sera block the binding of antibodies specific for the GP glycan cap, the GP1-GP2 interface, the mucin-like domain, and the membrane-proximal external region. Thus, administration of a cocktail of three ebolavirus vaccines induces a desirable broad antibody response, without skewing of the response toward preferential recognition of a single virus. IMPORTANCE The symptoms of the disease caused by the ebolaviruses Ebola, Bundibugyo, and Sudan are similar, and their areas of endemicity overlap. However, because of the limited antigenic relatedness of the ebolavirus glycoprotein (GP) used in all candidate vaccines against these viruses, they protect only against homologous and not against heterologous ebolaviruses. Therefore, a broadly specific pan-ebolavirus vaccine is required, and this might be achieved by administration of a cocktail of vaccines. The effects of cocktail administration of ebolavirus vaccines on the antibody repertoire remain unknown. Here, an in-depth analysis of the antibody responses to administration of a cocktail of human parainfluenza virus type 3-vectored vaccines against individual ebolaviruses was performed, which included analysis of binding to GP, neutralization of individual ebolaviruses, epitope specificity, Fc-mediated functions, and protection against the three ebolaviruses. The results demonstrated potent and balanced responses against individual ebolaviruses and no significant reduction of the responses compared to that induced by individual vaccines.

Quantitative guidelines to distinguish allergenic proteins from related, but non-allergenic ones are urgently needed for regulatory agencies, biotech companies and physicians. In a previous study, we found that allergenic proteins populate a relatively small number of protein families, as characterized by the Pfam database. However, these families also contain non-allergenic proteins, meaning that allergenic determinants must lie within more discrete regions of the sequence. Thus, new methods are needed to discriminate allergenic proteins within those families. Physical-Chemical Properties (PCP)-motifs specific for allergens within a Pfam class were determined for 17 highly populated protein domains. A novel scoring method based on PCP-motifs that characterize known allergenic proteins within these families was developed, and validated for those domains. The motif scores distinguished sequences of allergens from a large selection of 80,000 randomly selected non-allergenic sequences. The motif scores for the birch pollen allergen (Bet v 1) family, which also contains related fruit and nut allergens, correlated better than global sequence similarities with clinically observed cross-reactivities among those allergens. Further, we demonstrated that the average scores of allergen specific motifs for allergenic profilins are significantly different from the scores of non-allergenic profilins. Several of the selective motifs coincide with experimentally determined IgE epitopes of allergenic profilins. The motifs also discriminated allergenic pectate lyases, including Jun a 1 from mountain cedar pollen, from similar proteins in the human microbiome, which can be assumed to be non-allergens. The latter lacked key motifs characteristic of the known allergens, some of which correlate with known IgE binding sites.

Western equine encephalitis virus (WEEV) is an arbovirus from the genus Alphavirus , family Togaviridae , which circulates in North America between birds and mosquitoes, occasionally causing disease in humans and equids. In recent decades, human infection has decreased dramatically; the last documented human case in North America occurred in 1994, and the virus has not been detected in mosquito pools since 2008. Because limited information exists regarding the evolution of WEEV, we analyzed the genomic sequences of 33 low-passage-number strains with diverse geographic and temporal distributions and performed comprehensive phylogenetic analyses. Our results indicated that WEEV is a highly conserved alphavirus with only approximately 5% divergence in its most variable genes. We confirmed the presence of the previously determined group A and B lineages and further resolved group B into three sublineages. We also observed an increase in relative genetic diversity during the mid-20th century, which correlates with the emergence and cocirculation of several group B sublineages. The estimated WEEV population size dropped in the 1990s, with only the group B3 lineage being sampled in the past 20 years. Structural mapping showed that the majority of substitutions in the envelope glycoproteins occurred at the E2-E2 interface. We hypothesize that an event occurred in the mid-20th century that resulted in the increased genetic diversity of WEEV in North America, followed by genetic constriction due to either competitive displacement by the B3 sublineage or stochastic events resulting from a population decline. IMPORTANCE Western equine encephalitis virus (WEEV) has caused several epidemics that resulted in the deaths of thousands of humans and hundreds of thousands of equids during the past century. During recent decades, human infection decreased drastically and the virus has not been found in mosquito pools since 2008. Because limited information exists regarding the evolution of WEEV, we analyzed 33 complete genome sequences and conducted comprehensive phylogenetic analyses. We confirmed the presence of two major lineages, one of which diverged into three sublineages. Currently, only one of those sublineages is found circulating in nature. Understanding the evolution of WEEV over the past century provides a unique opportunity to observe an arbovirus that is in decline and to better understand what factors can cause said decline.

Filoviruses are the causative agents of lethal hemorrhagic fever in human and non-human primates (NHP). The family of Filoviridae is composed of three genera, Ebolavirus, Marburgvirus and Cuevavirus. There are currently no approved vaccines or antiviral therapeutics for the treatment of filovirus infections in humans. Passive transfer of neutralizing antibodies targeting the Ebola virus (EBOV) glycoprotein (GP) has proven effective in protecting mice, guinea pigs and NHP from lethal challenges with EBOV. In this study, we generated two neutralizing monoclonal antibodies (MAbs), termed S9 and M4 that recognize the GP of EBOV or multiple strains of Marburg virus (MARV), respectively. We characterized the putative binding site of S9 as a linear epitope on the glycan cap of the GP1 subunit of the EBOV-GP. The M4 antibody recognizes an unknown conformational epitope on MARV-GP. Additionally, we demonstrated the post-exposure protection potential of these antibodies in both the mouse and guinea pig models of filovirus infection. These data indicate that MAbs S9 and M4 would be good candidates for inclusion in an antibody cocktail for the treatment of filovirus infections.

We present the results for CAPRI Round 30, the first joint CASP-CAPRI experiment, which brought together experts from the protein structure prediction and protein–protein docking communities. The Round comprised 25 targets from amongst those submitted for the CASP11 prediction experiment of 2014. The targets included mostly homodimers, a few homotetramers, and two heterodimers, and comprised protein chains that could readily be modeled using templates from the Protein Data Bank. On average 24 CAPRI groups and 7 CASP groups submitted docking predictions for each target, and 12 CAPRI groups per target participated in the CAPRI scoring experiment. In total more than 9500 models were assessed against the 3D structures of the corresponding target complexes. Results show that the prediction of homodimer assemblies by homology modeling techniques and docking calculations is quite successful for targets featuring large enough subunit interfaces to represent stable associations. Targets with ambiguous or inaccurate oligomeric state assignments, often featuring crystal contact-sized interfaces, represented a confounding factor. For those, a much poorer prediction performance was achieved, while nonetheless often providing helpful clues on the correct oligomeric state of the protein. The prediction performance was very poor for genuine tetrameric targets, where the inaccuracy of the homology-built subunit models and the smaller pair-wise interfaces severely limited the ability to derive the correct assembly mode. Our analysis also shows that docking procedures tend to perform better than standard homology modeling techniques and that highly accurate models of the protein components are not always required to identify their association modes with acceptable accuracy. We are most grateful to the PDBe at the European Bioinformatics Institute in Hinxton, UK, for hosting the CAPRI website. Our deepest thanks go to all the structural biologists and to the following structural genomics initiatives: Northeast Structural Genomics Consortium, Joint Center for Structural Genomics, NatPro PSI:Biology, New York Structural Genomics Research Center, Midwest Center for Structural Genomics, Structural Genomics Consortium, for contributing the targets for this joint CASP-CAPRI experiment. MFL acknowledges support from the FRABio FR3688 Research Federation “Structural & Functional Biochemistry of Biomolecular Assemblies.”

The phenomenon of cross-reactivity between allergenic proteins plays an important role to understand how the immune system recognizes different antigen proteins. Allergen proteins are known to cross-react if their sequence comparison shows a high sequence identity which also implies that the proteins have a similar 3D fold. In such cases, linear sequence alignment methods are frequently used to predict cross-reactivity between allergenic proteins. However, the prediction of cross-reactivity between distantly related allergens continues to be a challenging task. To overcome this problem, we developed a new structure-based computational method, Cross-React, to predict cross-reactivity between allergenic proteins available in the Structural Database of Allergens (SDAP). Our method is based on the hypothesis that we can find surface patches on 3D structures of potential allergens with amino acid compositions similar to an epitope in a known allergen. We applied the Cross-React method to a diverse set of seven allergens, and successfully identified several cross-reactive allergens with high to moderate sequence identity which have also been experimentally shown to cross-react. Based on these findings, we suggest that Cross-React can be used as a predictive tool to assess protein allergenicity and cross-reactivity. Availability and Implementation: Cross-React is available at: http://curie.utmb.edu/Cross-React.html

The homeobox encodes a DNA-binding domain found in transcription factors regulating key developmental processes. The most notable examples of homeobox containing genes are the Hox genes, arranged on chromosomes in the same order as their expression domains along the body axis. The mechanisms responsible for the synchronous regulation of Hox genes and the molecular function of their colinearity remain unknown. Here we report the discovery of a conserved structural signature of the 180-base pair DNA fragment comprising the homeobox. We demonstrate that the homeobox DNA has a characteristic 3-base-pair periodicity in the hydroxyl radical cleavage pattern. This periodic pattern is significant in most of the 39 mammalian Hox genes and in other homeobox-containing transcription factors. The signature is present in segmented bilaterian animals as evolutionarily distant as humans and flies. It remains conserved despite the fact that it would be disrupted by synonymous mutations, which raises the possibility of evolutionary selective pressure acting on the structure of the coding DNA. The homeobox coding DNA may therefore have a secondary function, possibly as a regulatory element. The existence of such element may have important consequences for understanding how these genes are regulated.

The measles virus (MV) fusion (F) protein trimer executes membrane fusion after receiving a signal elicited by receptor binding to the hemagglutinin (H) tetramer. Where and how this signal is received is understood neither for MV nor for other paramyxoviruses. Because only the prefusion structure of the parainfluenza virus 5 (PIV5) F-trimer is available, to study signal receipt by the MV F-trimer, we generated and energy-refined a homology model. We used two approaches to predict surface residues of the model interacting with other proteins. Both approaches measured interface propensity values for patches of residues. The second approach identified, in addition, individual residues based on the conservation of physical chemical properties among F-proteins. Altogether, about 50 candidate interactive residues were identified. Through iterative cycles of mutagenesis and functional analysis, we characterized six residues that are required specifically for signal transmission; their mutation interferes with fusion, although still allowing efficient F-protein processing and cell surface transport. One residue is located adjacent to the fusion peptide, four line a cavity in the base of the F-trimer head, while the sixth residue is located near this cavity. Hydrophobic interactions in the cavity sustain the fusion process and contacts with H. The cavity is flanked by two different subunits of the F-trimer. Tetrameric H-stalks may be lodged in apposed cavities of two F-trimers. Because these insights are based on a PIV5 homology model, the signal receipt mechanism may be conserved among paramyxoviruses.

In many countries regulatory agencies have adopted safety guidelines, based on bioinformatics rules from the WHO/FAO and EFSA recommendations, to prevent potentially allergenic novel foods or agricultural products from reaching consumers. We created the Structural Database of Allergenic Proteins (SDAP, http://fermi.utmb.edu/SDAP/) to combine data that had previously been available only as flat files on Web pages or in the literature. SDAP was designed to be user friendly, to be of maximum use to regulatory agencies, clinicians, as well as to scientists interested in assessing the potential allergenic risk of a protein. We developed methods, unique to SDAP, to compare the physicochemical properties of discrete areas of allergenic proteins to known IgE epitopes. We developed a new similarity measure, the property distance (PD) value that can be used to detect related segments in allergens with clinical observed cross-reactivity. We have now expanded this work to obtain experimental validation of the PD index as a quantitative predictor of IgE cross-reactivity, by designing peptide variants with predetermined PD scores relative to known IgE epitopes. In complementary work we show how sequence motifs characteristic of allergenic proteins in protein families can be used as fingerprints for allergenicity.

Background Precise determination of conformational epitopes of neutralizing antibodies represents a key step in the rational design of novel vaccines. A powerful experimental method to gain insights on the physical chemical nature of conformational epitopes is the selection of linear peptides that bind with high affinities to a monoclonal antibody of interest by phage display technology. However, the structural characterization of conformational epitopes from these mimotopes is not straightforward, and in the past the interpretation of peptide sequences from phage display experiments focused on linear sequence analysis to find a consensus sequence or common sequence motifs. Results We present a fully automated search method, EpiSearch that predicts the possible location of conformational epitopes on the surface of an antigen. The algorithm uses peptide sequences from phage display experiments as input, and ranks all surface exposed patches according to the frequency distribution of similar residues in the peptides and in the patch. We have tested the performance of the EpiSearch algorithm for six experimental data sets of phage display experiments, the human epidermal growth factor receptor-2 (HER-2/neu), the antibody mAb Bo2C11 targeting the C2 domain of FVIII, antibodies mAb 17b and mAb b12 of the HIV envelope protein gp120, mAb 13b5 targeting HIV-1 capsid protein and 80R of the SARS coronavirus spike protein. In all these examples the conformational epitopes as determined by the X-ray crystal structures of the antibody-antigen complexes, were found within the highest scoring patches of EpiSearch, covering in most cases more than 50% residues of experimental observed conformational epitopes. Input options of the program include mapping of a single peptide or a set of peptides on the antigen structure, and the results of the calculation can be visualized on our interactive web server. Availability Users can access the EpiSearch from our web server http://curie.utmb.edu/episearch.html Contact webraun@utmb.edu

Structural plasticity of mammalian cytochromes P450 (CYP) has recently been explored in our laboratory and elsewhere to understand the ligand-binding promiscuity. CYP2B4 exhibits very different conformations and thermodynamic signatures in binding the small inhibitor 4-(4-chlorophenyl)imidazole (4-CPI) versus the large bifonazole. Using four key active-site mutants (F296A, T302A, I363A, and V367L) that are involved in binding one or both inhibitors, we dissected the thermodynamic basis for the ability of CYP2B4 to bind substrates and inhibitors of different sizes and chemistry. In all cases, 1:1 binding stoichiometry was observed. The inhibitors 4-CPI, 1-(4-chlorophenyl)imidazole, and 1-(2-(benzyloxy)ethyl)imidazole bind to the mutants with a free energy difference (Delta Delta G) of approximately 0.5 to 1 kcal/mol compared with the wild type but with a large entropy-enthalpy compensation of up to 50 kcal/mol. The substrate testosterone binds to all four mutants with a Delta Delta G of approximately 0.5 kcal/mol but with as much as 40 kcal/mol of entropy-enthalpy compensation. In contrast, benzphetamine binding to V367L and F296A is accompanied by a Delta Delta G of approximately 1.5 and 3 kcal/mol, respectively. F296A, I363A, and V367L exhibit very different benzphetamine metabolite profiles, indicating the different substrate-binding orientations in the active site of each mutant. Overall, the findings indicate that malleability of the active site allows mammalian P450s to exhibit a high degree of thermodynamic fidelity in ligand binding.

We have developed a fully automated method, InterProSurf, to predict interacting amino acid residues on protein surfaces of monomeric 3D structures. Potential interacting residues are predicted based on solvent accessible surface areas, a new scale for interface propensities, and a cluster algorithm to locate surface exposed areas with high interface propensities. Previous studies have shown the importance of hydrophobic residues and specific charge distribution as characteristics for interfaces. Here we show differences in interface and surface regions of all physical chemical properties of residues as represented by five quantitative descriptors. In the current study a set of 72 protein complexes with known 3D structures were analyzed to obtain interface propensities of residues, and to find differences in the distribution of five quantitative descriptors for amino acid residues. We also investigated spatial pair correlations of solvent accessible residues in interface and surface areas, and compared log-odds ratios for interface and surface areas. A new scoring method to predict potential functional sites on the protein surface was developed and tested for a new dataset of 21 protein complexes, which were not included in the original training dataset. Empirically we found that the algorithm achieves a good balance in the accuracy of precision and sensitivity by selecting the top eight highest scoring clusters as interface regions. The performance of the method is illustrated for a dimeric ATPase of the hyperthermophile, Methanococcus jannaschii, and the capsid protein of Human Hepatitis B virus. An automated version of the method can be accessed from our web server at http://curie.utmb.edu/prosurf.html.

Despite the emerging importance of human P450 2B6 in xenobiotic metabolism, thorough biochemical and biophysical characterization has been impeded as a result of low expression in Escherichia coli. Comparison with similar N-terminal truncated and C-terminal His-tagged constructs (rat P450 2B1dH, rabbit 2B4dH, and dog 2B11dH) revealed that P450 2B6dH showed the lowest thermal stability, catalytic tolerance to temperature, and chemical stability against guanidinium chloride-induced denaturation. Eleven P450 2B6dH mutants were rationally engineered based on sequence comparison with the three other P450 2B enzymes and the solvent accessibility of residues in the ligand-free crystal structure of P450 2B4dH. L198M, L264F, and L390P showed approximately 3-fold higher expression than P450 2B6dH. L264F alone showed enhanced stability against thermal and chemical denaturation compared with P450 2B6dH and was characterized further functionally. L264F showed similar preferential inhibition by pyridine over imidazole derivatives as P450 2B6dH. The Leu(264)-->Phe substitution did not alter the K(s) for inhibitors or the substrate benzphetamine, the K(m) for 7-ethoxy-4-(trifluoromethyl)coumarin, or the benzphetamine metabolite profiles. The enhanced stability and monodisperse nature of L264F made it suitable for isothermal titration calorimetry studies. Interaction of 1-benzylimidazole with L264F yielded a clear binding isotherm with a distinctly different thermodynamic signature from P450 2B4dH. The inhibitor docked differently in the binding pocket of a P450 2B6 homology model than in 2B4, highlighting the different chemistry of the active site of these two enzymes. Thus, L264F is a good candidate to further explore the unique structure-function relationships of P450 2B6 using X-ray crystallography and solution thermodynamics.

A new method for predicting interacting residues in protein complexes, InterProSurf, was applied to the E1 envelope protein of Venezuelan equine encephalitis (VEEV). Monomeric and trimeric models of VEEV-E1 were constructed with our MPACK program, using the crystal structure of the E1 protein of Semliki forest virus as a template. An alignment of the E1 sequences from representative alphavirus sequences was used to determine physical chemical property motifs (likely functional areas) with our PCPMer program. Information on residue variability, propensity to be in protein interfaces, and surface exposure on the model was combined to predict surface clusters likely to interact with other viral or cellular proteins. Mutagenesis of these clusters indicated that the predictions accurately detected areas crucial for virus infection. In addition to the fusion peptide area in domain 2, at least two other surface areas play an important role in virus infection. We propose that these may be sites of interaction between the E1-E1 and E1-E2 subdomains of the envelope proteins that are required to assemble the functional unit. The InterProSurf method is, thus, an important new tool for predicting viral protein interactions. These results can aid in the design of new vaccines against alphaviruses and other viruses.

Recent x-ray structures of cytochrome P450 2B4 (CYP2B4) reveal an open form that undergoes a large-scale structural transition to a closed form upon binding to 4-(4-chlorophenyl)imidazole (4-CPI). Here, we report for the first time a complete solution thermodynamic study using isothermal titration calorimetry supported by spectroscopic studies to elucidate the conformational flexibility of CYP2B4 in binding imidazole inhibitors with different ring chemistry and side chains: 4-CPI, 1-benzylimidazole (1-BI), 1-CPI, 4-phenylimidazole (4-PI), 1-(2-(benzyloxy)ethyl)imidazole (BEI), and 1-PI. Each of the inhibitors induced type II spectral changes, and IC50 values for enzyme inhibition ranged from 0.1 to 2.4 microM, following the order 1-BI < 4-CPI < 1-CPI < 4-PI < BEI < 1-PI. Calorimetric titrations using monomeric enzyme yielded a 1:1 binding stoichiometry, with the associated KD values ranging from 0.3 to 4.8 microM and following the same rank order as the IC50 values. Changes in enthalpy at 25 degrees C ranged from -6.5 to -8.8 kcal mol(-1). The largest difference in binding entropy (+5.9 versus -4.1 cal mol(-1) K(-1)) was observed between 4-CPI and BEI, respectively, with a 2-fold difference in heat capacity changes (-604 versus -331 cal mol(-1) K(-1)), which is inferred to result from the reduction of apolar surface area of the enzyme ensuing from a conformational change upon 4-CPI binding. Accessibility to acrylamide of the only tryptophan (Trp121), which is located in helix C, was greatly decreased only in protein bound to 4-CPI. Steric restrictions hindered the perfect docking of only BEI to the closed conformation of the enzyme. The thermodynamic signature obtained for structurally similar inhibitors suggests remarkable plasticity of CYP2B4.

We study a driven quasiperiodic skew-product dynamical mapping in which orbits with all Lyapunov exponents equal to zero lie on fractal attractors. These form a special category of strange nonchaotic attractors (SNAs), and we describe the scenario for their formation as well as methods for their characterization.

Aperiodic dynamics which is nonchaotic is realized on Strange Nonchaotic attractors (SNAs). Such attractors are generic in quasiperiodically driven nonlinear systems, and like strange attractors, are geometrically fractal. The largest Lyapunov exponent is zero or negative: trajectories do not show exponential sensitivity to initial conditions. In recent years, SNAs have been seen in a number of diverse experimental situations ranging from quasiperiodically driven mechanical or electronic systems to plasma discharges. An important connection is the equivalence between a quasiperiodically driven system and the Schr\"odinger equation for a particle in a related quasiperiodic potential, giving a correspondence between the localized states of the quantum problem with SNAs in the related dynamical system. In this review we discuss the main conceptual issues in the study of SNAs, including the different bifurcations or routes for the creation of such attractors, the methods of characterization, and the nature of dynamical transitions in quasiperiodically forced systems. The variation of the Lyapunov exponent, and the qualitative and quantitative aspects of its local fluctuation properties, has emerged as an important means of studying fractal attractors, and this analysis finds useful application here. The ubiquity of such attractors, in conjunction with their several unusual properties, suggest novel applications., Comment: 34 pages, 9 figures(5 figures are in ps format and four figures are in gif format)

We show that it is possible to devise a large class of skew-product dynamical systems which have strange nonchaotic attractors (SNAs): the dynamics is asymptotically on fractal attractors and the largest Lyapunov exponent is non-positive. Furthermore, we show that quasiperiodic forcing, which has been a hallmark of essentially all hitherto known examples of such dynamics is not necessary for the creation of SNAs.

Use of atomic force microscopy (AFM) has recently led to a better understanding of the molecular mechanisms of the unfolding process by mechanical forces; however, the rational design of novel proteins with specific mechanical strength remains challenging. We have approached this problem from a new perspective that generates linear physical-chemical properties (PCP) motifs from a limited AFM data set. Guided by our linear sequence analysis we designed and analyzed four new mutants of the titin I1 domain with the goal of increasing the domain's mechanical strength. All four mutants could be cloned and expressed as soluble proteins. AFM data indicate that at least two of the mutants have increased molecular mechanical strength. This observation suggests that the PCP method is useful to graft sequences specific for high mechanical stability to weak proteins to increase their mechanical stability, and represents an additional tool in the design of novel proteins besides steered molecular dynamics calculations, coarse grained simulations and phi-value analysis of the transition state.

Background: Bla g 2, one of the major cockroach allergens, induces a strong IgE response against conformational epitopes, and on reexposure, sensitized individuals often display symptoms of allergic rhinitis and asthma. The aim of the current study was to perform a test of the efficacy of a modified phage display screening, characterization of selected phages and an automated algorithm, EpiSearch, in locating an important conformational epitope. Methods: The monoclonal antibody 7C11, which partially inhibits the binding of patient IgE antibodies to Bla g 2, was used to screen a random peptide phage library. After 3 rounds of panning, 32 phage clones were isolated and the amino acid sequences of their peptides were determined. The relative affinity and specificity of the binding of these peptides to 7C11 were tested in ELISAs. The amino acid composition of these peptides was then matched with clusters of residues on the surface of the 3-dimensional (3D) structure of Bla g 2, using our EpiSearch algorithm. Results: The amino acid sequences of the peptides on selected phages differed at only one position, occupied by 1 of 2 negatively charged residues. The two 12-mer sequences bound to 7C11 with similar avidity and specificity. There was good concordance between the residues in the 3D clusters identified from our phage display/computational method with the co-crystal structural analysis. Conclusion: Conformational epitopes may be mapped through screening of clones from random peptide phage display libraries and EpiSearch.

Similarities in sequences and 3D structures of allergenic proteins provide vital clues to identify clinically relevant immunoglobulin E (IgE) cross-reactivities. However, experimental 3D structures are available in the Protein Data Bank for only 5% (45/829) of all allergens catalogued in the Structural Database of Allergenic Proteins (SDAP, http://fermi.utmb.edu/SDAP ). Here, an automated procedure was used to prepare 3D-models of all allergens where there was no experimentally determined 3D structure or high identity (95%) to another protein of known 3D structure. After a final selection by quality criteria, 433 reliable 3D models were retained and are available from our SDAP Website. The new 3D models extensively enhance our knowledge of allergen structures. As an example of their use, experimentally derived “continuous IgE epitopes” were mapped on 3 experimentally determined structures and 13 of our 3D-models of allergenic proteins. Large portions of these continuous sequences are not entirely on the surface and therefore cannot interact with IgE or other proteins. Only the surface exposed residues are constituents of “conformational IgE epitopes” which are not in all cases continuous in sequence. The surface exposed parts of the experimental determined continuous IgE epitopes showed a distinct statistical distribution as compared to their presence in typical protein–protein interfaces. The amino acids Ala, Ser, Asn, Gly and particularly Lys have a high propensity to occur in IgE binding sites. The 3D-models will facilitate further analysis of the common properties of IgE binding sites of allergenic proteins.

Summary: A new web server, InterProSurf, predicts interacting amino acid residues in proteins that are most likely to interact with other proteins, given the 3D structures of subunits of a protein complex. The prediction method is based on solvent accessible surface area of residues in the isolated subunits, a propensity scale for interface residues and a clustering algorithm to identify surface regions with residues of high interface propensities. Here we illustrate the application of InterProSurf to determine which areas of Bacillus anthracis toxins and measles virus hemagglutinin protein interact with their respective cell surface receptors. The computationally predicted regions overlap with those regions previously identified as interface regions by sequence analysis and mutagenesis experiments.Availability: The InterProSurf web server is available at http://curie.utmb.edu/Contact: webraun@utmb.eduSupplementary information: Other test examples are available as Supplementary Material at Bioinformatics online.

In immunogold double-labeling of pea leaf thin sections with antibodies raised against ferredoxin-NADP reductase (EC 1.18.1.2, FNR) and antibodies directed against the A or B subunits of the NADP-linked glyceraldehyde-3-P dehydrogenase (GAPD) (EC 1.2.1.13), many small and large gold particles were found together over the chloroplasts. Nearest neighbor analysis of the distribution of the gold particles indicates that FNR and the NADP-linked GAPD are co-localized, in situ. This suggests that FNR might carry FADH2 or NADPH from the thylakoid membrane to GAPD, or that ferredoxin might carry electrons to FNR co-localized with GAPD in the stroma. Crystal structures of the spinach enzymes are available. When they are docked computationally, the proteins appear, as modeled, to be able to form at least two different complexes. One involves a single GAPD monomer and an FNR monomer (or dimer). The amino acid residues located at the putative interface are highly conserved on the chloroplastic forms of both enzymes. The other potential complex involves the GAPD A2B2 tetramer and an FNR monomer (or dimer). The interface residues are conserved in this model as well. Ferredoxin is able to interact with FNR in either complex.

Conformational flexibility and cooperativity in ligand recognition are two key aspects of the catalytic diversity of cytochrome P450 enzymes. In this study, we dissect the ligand binding stoichiometry and energetics of the soluble bacterial P450eryF by isothermal titration calorimetry (ITC) using three allosteric and two non-allosteric ligands of diverse chemistry. Complementary spectral binding studies and sequential, two-ligand docking simulations were performed to help assign the binding sites. Binding of 4-phenylpyridine (4-PP) or 4-(4-chlorophenyl)imidazole (4-CPI) showed 1:1 stoichiometry in ITC, consistent with the lack of cooperativity observed in spectral binding studies. The larger ligands 9-aminophenanthrene (9-AP), 1-pyrenebutanol (1-PB), and alpha-naphthoflavone (ANF) show cooperative spectral binding and yielded 2:1 stoichiometry. The associated thermodynamic parameters for the sites were calculated using a sequential binding mechanism. The binding constant (KD) for the first site was almost two times lower than that of the second site for all three compounds. Ligand binding at site 1 was entropically favored, whereas binding at site 2 was weak and entropically unfavorable. Simulations showed that two molecules of 9-AP, ANF or 1-PB can be adequately docked to two individual sub-sites within a large binding pocket. The absence of hydrophobic tethering and ligand stacking are consistent with the single low affinity binding site observed for 4-CPI and 4-PP. Competitive binding studies with P450eryF preloaded with either 1-PB or ANF showed a decrease in the affinities for 9-AP at both the sites, with large entropy-enthalpy compensation, indicating the ability of the binding pocket to accommodate two ligands of diverse chemistry and enable cooperativity.

We discuss several bifurcation phenomena that occur in the quasiperiodically driven logistic map. This system can have strange nonchaotic attractors (SNAs) in addition to chaotic and regular attractors; on SNAs the dynamics is aperiodic, but the largest Lyapunov exponent is nonpositive. There are a number of different transitions that occur here, from periodic attractors to SNAs, from SNAs to chaotic attractors, etc. We describe some of these transitions by examining the behavior of the largest Lyapunov exponent, distributions of finite time Lyapunov exponents and the invariant densities in the phase space., Comment: 17 Pages, 8 Figures(four figures are in ps format and four figures are in gif format

Protein−protein interactions mediated by the carboxy terminus domains of breast cancer gene 1 (BRCT) play a critical role in the regulation of cell cycle checkpoints in response to DNA damage. DNA damage-induced checkpoint regulation bestows radiation resistance to cancers, and inhibition of the proteins involved has been suggested as potential drug targets. Structural studies show that BRCT binding peptides, namely BACH1 and CtIP that are important for cell cycle regulation, occupy the exact same site on BRCT through a pSXXF consensus recognition sequence. In the present study we use two complementary techniques, isothermal titration calorimetry (ITC) and fluorescence polarization (FP), to define structural minima of a BACH1 peptide. Remarkably the minimal tetrapeptide retains low-micromolar affinity for BRCT with a thermodynamic signature comparable to the decapeptide. The study also identifies a potential hydrophobic pocket adjacent to the pSXXF motif and intramolecular hydrogen bonds that can be explo...