Impact of AlphaFold on structure prediction of protein complexes: The CASP15-CAPRI experiment

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homodimers, 3 homo-trimers, 13 heterodimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target. A total of 21 941 models submitted by these groups and by 15 CAPRI scorer groups were evaluated using the CAPRI model quality measures and the DockQ score consolidating these measures. The prediction performance was quantified by a weighted score based on the number of models of acceptable quality or higher submitted by each group among their five best models. Results show substantial progress achieved across a significant fraction of the 60+ participating groups. High-quality models were produced for about 40% of the targets compared to 8% two years earlier. This remarkable improvement is due to the wide use of the AlphaFold2 and AlphaFold2-Multimer software and the confidence metrics they provide. Notably, expanded sampling of candidate solutions by manipulating these deep learning inference engines, enriching multiple sequence alignments, or integration of advanced modeling tools, enabled top performing groups to exceed the performance of a standard AlphaFold2-Multimer version used as a yard stick. This notwithstanding, performance remained poor for complexes with antibodies and nanobodies, where evolutionary relationships between the binding partners are lacking, and for complexes featuring conformational flexibility, clearly indicating that the prediction of protein complexes remains a challenging problem.

Alanine Scanning of the Unstructured Region of Ara h 2 and of a Related Mimotope Reveals Critical Amino Acids for IgE Binding

SCOPE

The unstructured region of Ara h 2, referred to as epitope 3, contains a repeated motif, DYPSh (h = hydroxyproline) that is important for IgE binding.

METHODS AND RESULTS

IgE binding assays to 20mer and shorter peptides of epitope 3, defines a 16mer core sequence containing one copy of the DPYSh motif, DEDSYERDPYShSQDP. This study performs alanine scanning of this and a related 12mer mimotope, LLDPYAhRAWTK. IgE binding, using a pool of 10 sera and with individual sera, is greatly reduced when alanine is substituted for aspartate at position 8 (D8; p < 0.01), tyrosine at position 10 (Y10; p < 0.01), and hydroxyproline at position 12 (h12; p < 0.001). IgE binding to alanine-substituted peptides of a mimotope containing the DPY_h motif confirm the critical importance of Y (p < 0.01) and h (p < 0.01), but not D. Molecular modeling of the core and mimotope suggests an h-dependent conformational basis for the recognition of these sequences by polyclonal IgE.

CONCLUSIONS

IgE from pooled sera and individual sera differentially bound amino acids throughout the sequences of Epitope 3 and its mimotope, with Y10 and h12 being most important for all sera. These results are highly significant for designing hypoallergenic forms of Ara h 2.

The updated Structural Database of Allergenic Proteins (SDAP 2.0) provides 3D models for allergens and incorporated bioinformatics tools

Background

Allergenic proteins can cause IgE-mediated adverse reactions in sensitized individuals. Although the sequences of many allergenic proteins have been identified, bioinformatics data analysis with advanced computational methods and modeling is needed to identify the basis for IgE binding and cross-reactivity.

Objective

We aim to present the features and use of the updated Structural Database of Allergenic Proteins 2.0 (SDAP 2.0) webserver, a unique, publicly available resource to compare allergens using specially designed computational tools and new high-quality 3-D models for most known allergens.

Methods

Previously developed and novel software tools for identifying cross-reactive allergens using sequence and structure similarity are implemented in SDAP 2.0. A comprehensive set of high-quality 3-D models of most allergens was generated with the state-of-the-art AlphaFold 2 software. A graphics tool enables the interactive visualization of IgE epitopes on experimentally determined and modeled 3-D structures.

Results

A user can search for allergens similar to a given input sequence with the FASTA algorithm or the window-based World Health Organization/International Union of Immunological Societies (WHO/IUIS) guidelines on safety concerns of novel food products. Peptides similar to known IgE epitopes can be identified with the property distance tool and conformational epitopes by the Cross-React method. The updated database contains 1657 manually curated sequences including all allergens from the IUIS database, 334 experimentally determined X-ray or NMR structures, and 1565 3-D models. Each allergen/isoallergen is classified according to its protein family.

Conclusions

SDAP provides access to the steadily increasing information on allergenic structures and epitopes with integrated bioinformatics tools to identify and analyze their similarities. In addition to serving the research and regulatory community, it provides clinicians with tools to identify potential coallergies in a sensitive patient and can help companies to design hypoallergenic foods and immunotherapies.

Still SDAPing Along: 20 Years of the Structural Database of Allergenic Proteins

The introduction of plant extracts to mitigate the symptoms of “hay fever”, about a century ago, led to discoveries beginning sixty years ago on determining the sequences and eventually structures of allergenic proteins. As more proteins were cloned, there was a need to rapidly identify and categorize those with significant similarity to known allergens. The Structural Database of Allergenic Proteins (SDAP) was created at the beginning of the 21st century as the first cross-referenced website to allow rapid overview of the structures and sequences of allergenic proteins. SDAP provides a way to identify sequence and functional similarities between these proteins, despite the complex nomenclature system based on the Latin names of their different sources. A rapid FASTA search simplifies grouping allergens from the same structural or functional family. SDAP also provides an overview of the rapidly expanding literature on the sequence, structure and epitopes of allergenic proteins and a way to estimate the potential allergenicity of novel proteins based on rules provided by the IUIS. Twenty years and a pandemic later, the list of allergenic proteins and their attributes continues to grow. SDAP is expanding and improving to allow rapid access to all this information.

Regional and temporal coordinated mutation patterns in SARS-CoV-2 spike protein revealed by a clustering and network analysis

SARS-CoV-2 has steadily mutated during its spread to > 300 million people throughout the world. The WHO has designated strains with certain mutations, "variants of concern" (VOC), as they may have higher infectivity and/or resist neutralization by antibodies in sera of vaccinated individuals and convalescent patients. Methods to detect regionally emerging VOC are needed to guide treatment and vaccine design. Cluster and network analysis was applied to over 1.2 million sequences of the SARS-CoV-2 spike protein from 36 countries in the GISAID database. While some mutations rapidly spread throughout the world, regionally specific groups of variants were identified. Strains circulating in each country contained different sets of high frequency mutations, many of which were known VOCs. Mutations within clusters increased in frequency simultaneously. Low frequency, but highly correlated mutations detected by the method could signal emerging VOCs, especially if they occur at higher frequency in other regions. An automated version of our method to find high frequency mutations in a set of SARS-COV-2 spike sequences is available online at http://curie.utmb.edu/SAR.html .

Design of peptides with high affinity binding to a monoclonal antibody as a basis for immunotherapy

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.

Synthetic proteins for COVID-19 diagnostics

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 E484K or N501Y, can be quickly synthesized in 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.

The importance of the 2S albumins for allergenicity and cross-reactivity of peanuts, tree nuts, and sesame seeds

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.

IgE binding to linear epitopes of Ara h 2 in peanut allergic preschool children undergoing oral Immunotherapy

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 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

Values of PN-sIgE correlated with eliciting dose (P = .001) and with a higher likelihood of achieving SU (P < .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 ≥ 14 kIU/L, binding of IgE to epitopes 5 and 6 of Ara h 2 was associated with a lower eliciting dose at baseline challenge (P < .001; P_c  < .02). In subjects with PN-sIgE ≥ 35 kIU/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 = .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.

Cross-React: a new structural bioinformatics method for predicting allergen cross-reactivity

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.

Contact

ssnegi@utmb.edu.

Prediction of homoprotein and heteroprotein complexes by protein docking and template-based modeling: A CASP-CAPRI experiment

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. Proteins 2016; 84(Suppl 1):323-348. © 2016 Wiley Periodicals, Inc.

Conformational IgE epitopes of peanut allergens Ara h 2 and Ara h 6

BACKGROUND

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.

Engineering proteins with enhanced mechanical stability by force-specific sequence motifs

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 ϕ-value analysis of the transition state.

Validation of a phage display and computational algorithm by mapping a conformational epitope of Bla g 2

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.

Direct detection and serogroup characterization of Neisseria meningitidis from outbreak of meningococcal meningitis in Delhi

BACKGROUND AND OBJECTIVES

Rapid clinical manifestation/progression of the meningococcal meningitis and lacunae in conventional bacteriological test often encourages indiscriminate use of antibiotics much before the etiology is established. Accordingly this study was planned to evaluate ctrA PCR for rapid molecular detection. In addition, multiplex PCR and sequencing was done for serogroup prediction to provide essential epidemiological and laboratory evidence for decision makers of health department of the country for choosing appropriate vaccine and phylogenetic analysis to establish its lineage.

MATERIALS AND METHODS

73 CSF samples, collected from equal number of suspected cases, were investigated by both bacteriological (microscopy, culture, LA and drug sensitivity testing) as well as molecular tests i.e. PCR targeting conserved ctrA gene, multiplex PCR for serogroup characterization and DNA sequencing.

RESULTS

ctrA PCR revealed sensitivity, specificity, positive predictive value and negative predictive values of 93.15%, 100%,100%, and 88.23% respectively. Multiplex PCR based genogrouping followed by DNA sequencing, BLAST and phylogenetic analysis revealed complete homology with earlier submitted Neisseria meningitidis serogroup A strain Z2491 to suggest the sole involvement of only serogroup A in the outbreak. Two strains showed resistance to cefuroxime, ciprofloxacin, nalidixic acid. Only one strain showed resistance to ciprofloxacin, emphasizing the need for a constant surveillance system.

CONCLUSION

These diagnostic molecular tools are of paramount importance in establishing etiology, serogrouping, and epidemiological surveillance especially in developing countries like India.

Characterization of RPO B gene for detection of rifampicin drug resistance by SSCP and sequence analysis

PURPOSE

Because of the emergence of multidrug-resistant tuberculosis in recent times, the rapid detection of resistance to the first-line anti-tuberculosis drug rifampicin was felt worldwide. Accordingly, this study was conducted to evaluate the diagnostic potential of polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) for checking its utility as a rapid screening test for determination of rifampicin drug resistance.

MATERIALS AND METHODS

A total of 34 isolates of Mycobacterium tuberculosis (M. tuberculosis) (22 rifampicin resistant, 11 rifampicin sensitive and one control H37Rv) strains were analysed by PCR-SSCP and DNA sequencing within the 157-bp region of the rpo B gene (Ala 500-Val 550).

RESULTS

Rifampicin resistance was detected successfully by PCR-SSCP in 20/22(90.90%) of rifampicin-resistant strains showing a total of nine different mutations in seven codon positions: codon 513 (CAA-->CCA), 516 (GAC-->GTC), 507 (GGC-->GAC), 526 (CAC-->GAC, TAC), 531 (TCG-->TTG, TGG), 522 (TCG-->TGG) and 533 (GTG-->CCG). Two rifampicin-resistant strains showed an identical PCR-SSCP pattern with the wild type H37Rv; 77.27% rifampicin-resistant strains showed a single point mutation and 9.09% had no mutation. Three rifampicin-resistant strains showed characteristic double mutations at codon positions 526 and 531. Sensitivity and specificity were calculated as 90.90% and 100%.

CONCLUSIONS

Rifampicin-resistant genotypes were mainly found in codon positions 516, 526 and 531. PCR-SSCP seems to be an efficacious method of predicting rifampicin resistance and substantially reduces the time required for susceptibility testing from 4 to 6 weeks to a few weeks.

Automated detection of conformational epitopes using phage display Peptide sequences

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 C₂ 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

Structural analysis of linear and conformational epitopes of allergens

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.

Comprehensive 3D-modeling of allergenic proteins and amino acid composition of potential conformational IgE epitopes

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.

Co-localization of glyceraldehyde-3-phosphate dehydrogenase with ferredoxin-NADP reductase in pea leaf chloroplasts

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.

InterProSurf: a web server for predicting interacting sites on protein surfaces

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/

Statistical analysis of physical-chemical properties and prediction of protein-protein interfaces

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.

Rational engineering of human cytochrome P450 2B6 for enhanced expression and stability: importance of a Leu264->Phe substitution

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.

Diagnostic potential of IS6110, 38kDa, 65kDa and 85B sequence-based polymerase chain reaction in the diagnosis of Mycobacterium tuberculosis in clinical samples

PURPOSE

The correlation between the presence of specific gene sequence of M. tuberculosis and specific diagnosis of clinical tuberculosis is not known. This study compared the results of polymerase chain reaction (PCR) amplification of M. tuberculosis specific DNA sequences (IS6110, 65kDa, 38kDa and mRNA coding for 85B protein) from different clinical samples of pulmonary and extrapulmonary tuberculosis.

METHODS

One hundred and seventy-two clinical samples from suspected tuberculosis patients were tested for smear examination, culture (LJ and rapid BACTEC 460 TB system) and PCR. PCR was performed with specific primers for the targets: IS6110, 65 kDa, 38 kDa and 85 B.

RESULTS

Each PCR test was found to have a much higher positivity than conventional test and BACTEC culture (P < 0.05). Smear positive samples (56) and the samples (36) showing positive results by conventional methods (smear and LJ medium culture) and BACTEC were found to be positive by all PCR protocols. No significant difference was found between the four PCR protocols (P> 0.05). The primer specific for amplifying the 123bp IS6110 fragment gave the highest positivity (83%), followed by 65kDa, 38kDa and 85B RT-PCR in descending order.

CONCLUSIONS

These data suggest that the presence of IS6110 correlates more closely with the diagnosis of clinical tuberculosis than that of 65kDa, 38kDa and 85B proteins.

Dissecting the Thermodynamics and Cooperativity of Ligand Binding in Cytochrome P450eryF

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.

Protein antigen b (Pab) based PCR test in diagnosis of pulmonary and extra-pulmonary tuberculosis

BACKGROUND AND OBJECTIVES

Diagnosis of tuberculosis (TB) is largely based on microscopy and culture examination which are either less sensitive, or time consuming. In the present study a PCR (polymerase chain reaction) test based on DNA sequence coding for a 38-kilodalton protein antigen b (Pab) ,specific for Mycobacterium tuberculosis was compared with Ziehl-Neelsen (ZN) stained AFB (acid fast bacilli) smear examination, culture based on conventional Lowenstein-Jensen (LJ) medium and radiometric BACTEC 460 system for the diagnosis of TB using clinical samples obtained from pulmonary and extra-pulmonary cases of TB.

METHODS

Clinical samples obtained from 168 patients of suspected TB (pulmonary and extrapulmonary) were subjected to ZN smear examination, LJ culture, radiometric BACTEC culture and a PCR test by amplifying 419 bp sequence coding for Pab, a glycoprotein of molecular weight 38 kDa.

RESULTS

A significant difference was seen in the sensitivity of different tests, the figures being 74.2 per cent for PCR test, 53.4 per cent for BACTEC culture, 47.1 per cent for LJ medium based culture and 35.2 per cent for ZN smear examination (P<0.05). However, there was no significant difference between different tests as far as specificity was concerned. PCR test sensitivity in pulmonary and extra-pulmonary clinical samples were 74.3 and 71.5 per cent respectively, being significantly higher (P<0.05) when compared with sensitivity of other tests. The mean detection time for M. tuberculosis was 24.0 days by LJ media culture, 12.8 days by BACTEC culture and less than 1 day by smear examination and PCR test.

INTERPRETATION AND CONCLUSION

PCR test is more sensitive than ZN smear examination, LJ medium culture and BACTEC culture for diagnosing TB in pulmonary and extra-pulmonary clinical samples.

Determining functionally important amino acid residues of the E1 protein of Venezuelan equine encephalitis virus

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.

Comparison of various microbiological tests including polymerase chain reaction for the diagnosis of osteoarticular tuberculosis

PURPOSE

To evaluate the utility of the polymerase chain reaction (PCR) test for diagnosing osteoarticular tuberculosis (TB).

METHODS

Clinical samples (synovial tissue and synovial fluid) obtained from 23 cases of suspected osteoarticular tuberculosis were subjected to Ziehl Neelsen (ZN) smear examination, radiometric BACTEC culture and PCR test for tuberculosis by amplifying 65 kDa antigen coding region of Mycobacterium tuberculosis (M.tb) genome.

RESULTS

PCR test was found to be much sensitive than the ZN smear examination and BACTEC culture (p<0.05) in the diagnosis of osteoarticular TB. In synovial fluid samples, PCR was positive in 73.9%, ZN smear examination in 17.39% and BACTEC culture in 39.13% of cases. The positivities were relatively lower with synovial tissue samples, the corresponding figures being 60.8, 8.6 and 26.08% respectively. Moreover, on combining the results of synovial fluid and tissues, the corresponding figures further increased to 78.2, 21.7 and 43.3% respectively. Further, sensitivity and specificity for PCR employing BACTEC culture as the "gold standard" was 100% respectively. Using BACTEC culture, the earliest positivity was seen in three days using synovial tissue specimen and 13 days with synovial fluid, the average detection times being 23.2 days and 32.6 days respectively. On the other hand, PCR test gave a positive result within 24 hours.

CONCLUSIONS

PCR test was shown to be much more sensitive than ZN smear examination and BACTEC culture test for diagnosing osteoarticular tuberculosis.

Comparison of the conventional diagnostic modalities, bactec culture and polymerase chain reaction test for diagnosis of tuberculosis

PURPOSE

To evaluate the performance of 65 kDa antigen based PCR assay in clinical samples obtained from pulmonary and extrapulmonary cases of tuberculosis.

METHODS

One hundred and fifty six samples were processed for detection of Mycobacterium tuberculosis by ZN smear examination, LJ medium culture, BACTEC radiometric culture and PCR tests.

RESULTS

A significant difference was seen in the sensitivities of different tests, the figures being 74.4% for PCR test, 33.79% for ZN smear examination, 48.9% for LJ culture and 55.8% for BACTEC culture (P< 0.05). However, there was no significant difference (P>0.05) as far as specificity of different tests was concerned. PCR test sensitivity in pulmonary and extrapulmonary clinical samples were 72.7% and 75.9% respectively and found to be significantly higher (P< 0.05) when compared with those of other tests. The mean detection time for M.tuberculosis was 24.03 days by LJ medium culture, 12.89 days by BACTEC culture and less than one day by PCR test.

CONCLUSIONS

PCR is a rapid and sensitive method for the early diagnosis of pulmonary and extrapulmonary tuberculosis.

Benign signet ring cell change with multilayering in the gallbladder mucosa--a case report

We describe a case of benign signet ring cell change in the gallbladder mucosa. On histopathological examination of H&E-stained sections, the gallbladder epithelium showed multilayering. The epithelial cells were large, columnar to polygonal with a small round basal or eccentric nucleus and vacuolated cytoplasm, giving them a signet ring appearance. There was no nuclear atypia, hyperchromatism or mitotic activity. The cells showed uniform positivity with mucicarmine, PAS and Alcian blue stains. The cytoplasmic vacuolations were negative for fat stains (Oil red O and Sudan IV). On immunohistochemistry, the cells showed positivity with antibodies for pancytokeratin (PCK) and epithelial membrane antigen (EMA). A diagnosis of benign signet ring cell change with multilayering in the gall bladder mucosa was made. Thoroughly reviewing the literature, we found only one case of benign signet ring cell aggregates in the gallbladder mucosa documented earlier. The lesion is hereby reported because of the unique histomorphology and the diagnostic dilemma which can occur as a malignant change in situ has to be excluded.

Comparison of magnetic resonance cholangiography and percutaneous transhepatic cholangiography in the evaluation of bile duct strictures after cholecystectomy

BACKGROUND

Percutaneous transhepatic cholangiography (PTC) has been the preferred investigation to delineate the anatomy of the biliary tract in a patient with a bile duct stricture after cholecystectomy. Recently magnetic resonance cholangiography (MRC) has been described to evaluate the obstructed biliary tract. This paper reports a comparison of MRC with PTC in evaluating patients with an iatrogenic bile duct stricture.

METHODS

This was a prospective study of 26 patients who had surgery for a bile duct stricture after cholecystectomy. Before operation all patients underwent both MRC and PTC, the results of which were compared with the intraoperative findings.

RESULTS

Both PTC and MRC were comparable with regard to image quality, detection of intrahepatic bile duct dilatation, assessment of the level of injury and detection of abnormalities such as intraduct calculi, cholangitic liver abscesses and atrophy of liver lobes. MRC provided additional information in four patients, including detection of associated fluid collections (n = 3) and portal hypertension (n = 1). In eight patients more than one puncture had to be performed during PTC to delineate the complete anatomy.

CONCLUSION

MRC is an accurate and non-invasive imaging procedure for preoperative evaluation of patients with a bile duct injury after cholecystectomy, and is capable of providing additional information which may not be available with PTC.

Experience of surgical management of pseudo-aneurysms of branches of the coeliac axis in a North Indian Hospital

BACKGROUND

Bleeding splanchnic artery pseudo-aneurysm is a rare but frequently fatal complication that can be successfully managed by angiographic embolization. However, certain patients because of hemodynamic instability, non-availability of technique or angiographic failure may require primary surgical intervention.

METHOD

Retrospective review of 13 patients presenting with exsanguinating hemorrhage from ruptured pseudo-aneurysm arising from branches of coeliac axis, managed surgically in absence of angiographic embolization.

RESULTS

Splenic artery was most commonly involved (n = 7) followed by hepatic (n = 3), gastroduodenal (n = 2) and left gastric artery (n = 1). The most common underlying aetiology was pancreatitis (n = 8, acute = 2; chronic = 6) followed by iatrogenic (n = 3), liver abscess (n = 1) and gastric ulcer (n = 1). Seven patients presented with upper gastrointestinal (GI) bleed, while 2 each with lower GI bleed, haemobilia and bleeding through tube drains. CT-scan accurately demonstrated the pseudo-aneurysm in 11 (84.6%) patients and additionally demonstrated the underlying pathology. The surgical management chiefly consisted of ligation of offending vessel and additional procedures directed at primary pathology. Overall, 77% patients had a favourable outcome while 23% died consequent to ongoing hemorrhage.

CONCLUSION

Pseudo-aneurysm involving the branches of coeliac axis most commonly arises as a result of pancreatitis and affects splenic artery. CT-scan accurately demonstrates pseudo-aneurysm and associated pathology in majority of cases. Primary surgical management in the presence of hemodynamic instability and non-availability of angiographic embolization is a viable alternative.

Critical states and fractal attractors in fractal tongues: localization in the Harper map

Localized states of Harper's equation correspond to strange nonchaotic attractors in the related Harper mapping. In parameter space, these fractal attractors with nonpositive Lyapunov exponents occur in fractally organized tongue-like regions which emanate from the Cantor set of eigenvalues on the critical line epsilon=1. A topological invariant characterizes wave functions corresponding to energies in the gaps in the spectrum. This permits a unique integer labeling of the gaps and also determines their scaling properties as a function of potential strength.

Paraganglionoma of extrahepatic biliary tract causing obstructive jaundice

We report a young woman with paraganglionoma arising from the extrahepatic bile duct presenting with acute obstructive jaundice. The patient underwent excision of the gall bladder and extrahepatic bile duct with the tumor, and Roux-en-Y hepaticojejunostomy. She is asymptomatic 9 months later, with normal biochemical investigations and imaging.

Mirage formation in the thermal region

Formation of a mirage in the 8-12-µm band has been observed with a high-resolution thermal camera designed and developed by us. Thermal imagery recorded under actual field conditions is presented to illustrate the mirage phenomenon. Further analysis is in progress to estimate the effect of mirages in this waveband for civil and military applications.

Simple determination of microbial protein in rumen liquor

A method for microbial protein determination based on Lowry's assay was devised. Differential centrifugation was used to separate bacterial cells from rumen liquor. Proteins from cells were released by suspending in .25 N sodium hydroxide and heating in boiling water bath for 10 min. Protein was determined by the Folin phenol method. The method is simple, accurate, reproducible, sensitive, specific, brief, and can be performed with small amount of rumen liquor (5 ml). A large number of samples can be handled conveniently.

Changes in blood constituents of guinea pigs in lantana toxicity

Lantana toxicity of guinea pigs elicited an increase in hematocrit, erythrocyte and leukocyte number, hemoglobin, urea-nitrogen and bilirubin contents in the blood of the affected animals. Most of the bilirubin was present in the conjugated form. Enzyme activities of glutamic oxaloacetic transaminase, acid phosphatase, lactate dehydrogenase, glutamate dehydrogenase, sorbitol dehydrogenase in the blood plasma of affect animals exhibited a marked increase. Acid phosphatase activity was inhibited by tartrate. Enzyme activity of alkaline phosphatase remained unchanged while that of glutamic pyruvic transaminase showed a marginal decrease.

Chronic aflatoxicosis in cattle

Pathological, hematological and plasma enzymatic studies were made on milk cattle affected by chronic aflatoxicosis caused by the prolonged feeding of concentrate feed mixtures containing contaminated groundnut cake having aflatoxin B1 (110 ug/kg groundnut cake at the time of sampling), B2, G1 and G2. Clinical and necropsy observations on liver included proliferation of connective tissue along portal triads leaving small group of hepatocytes intact. Liver function tests showed liver damage in three of the four affected animals studied.

A review of the toxicity of Lantana camara (Linn) in animals

Lantana poisoning has been taking a heavy toll of livestock year after year. All aspects of the problem are reviewed. Lantana poisoning in cattle, sheep, buffalo, and guinea pigs caused obstructive jaundice, photosensitization, and rise in serum glutamicoxaloaetic transaminase activity. The symptoms could be reproduced in sheep by administration of purified Lantadene A. Liver and kidneys are the most affected organs during lantana poisoning. Intoxication of guinea pigs with Lantana camara leads to marked alterations in major tissue constituents in liver an kidneys. Hepatic and renal xanthine oxidase activity is also elevated during lantana poisoning. No antidote is available against the toxic section of Lantana camara. Symptomatic treatments have been proposed with limited success. Knowledge of the biochemical mechanism of lantana intoxication at the cellular, subcellular, and molecular levels is essential in order to evolve a successful antidote and more rational therapy during lantana intoxication.

Effect of acetohydroxamic acid on rumen urease activity in vitro

Acetohydroxamic acid at a concentration of 1 X 10(-6) M, 6 X 10(-5) M, and 1 X 10(-3) M inhibited urease of intact rumen microbes in vitro by 11%, 50%, and 74%. Inhibition of rumen urease by the acid reached equilibrium state, unlike jack bean urease. Inhibition was maximum over a broad range of pH (8 to 10) and it did not resemble the pH activity profile of rumen urease. Sulfhydryl compounds did not reverse the inhibition; however, addition of these compounds prior to acetohydroxamic acid addition prevented inhibition. The nature of inhibition was noncompetitive with inhibitor constant 4.8 X 10(-5) M. Acetohydroxamic acid at a concentration that produced 50% urease inhibition did not affect rumen cellulase and proteolytic enzymes in vitro. The complex of acetohydroxamic acid-rumen urease is dissociable on dialysis.

Hepatic and renal toxicity of lantana in the guinea pig

Oral administration of lantana leaf powder to guinea pigs caused a decrease in hepatic and renal tissue dry weight, DNA and protein contents. Total carbohydrate content decreased in liver but was not affected in the kidneys. RNA content (as a ratio of tissue dry weight) showed an increase in both liver and kidneys. Lipid content of liver tissue increased while it decreased in the kidneys. Relative amounts of protein, DNA and RNA showed significant alterations in both the tissues.

Fragility of erythrocytes in animals affected by lantana poisoning

Clinical observations on four male cow calves and three other lantana poisoned animals under field conditions have been made along with hematological studies, including observations on plasma bilirubin content and osmotic fragility of erythrocytes. Four healthy male cow calves served as control. Five clinically advanced cases succumbed within 1 to 10 days of observation. Hematocrit values in four of the seven affected animals were abnormally high. The total plasma bilirubin content in the fatal cases was of the order of 5 to 50 times the normal value while, in calves that recovered, it was not more than 3 times. The osmotic fragility of erythrocytes was measured by the rate of hemolysis in decreasing concentrations of saline. The erythrocytes of normal animals resisted hemolysis till the saline concentration was lowered to 0.60%, and it was complete at the low saline concentration of 0.30%. However, in lantana poisoned animals, hemolysis started even in 0.72% saline and was complete at a saline concentration of 0.44%. The hemolysis curves of animals that died were far removed from the normal, but those of the calves that recovered and the one that died but had the lowest bilirubin content closely approached the normal curve.