Structure-based energetics of protein interfaces guides foot-and-mouth disease virus vaccine design

Virus capsids are primed for disassembly, yet capsid integrity is key to generating a protective immune response. Foot-and-mouth disease virus (FMDV) capsids comprise identical pentameric protein subunits held together by tenuous noncovalent interactions and are often unstable. Chemically inactivated or recombinant empty capsids, which could form the basis of future vaccines, are even less stable than live virus. Here we devised a computational method to assess the relative stability of protein-protein interfaces and used it to design improved candidate vaccines for two poorly stable, but globally important, serotypes of FMDV: O and SAT2. We used a restrained molecular dynamics strategy to rank mutations predicted to strengthen the pentamer interfaces and applied the results to produce stabilized capsids. Structural analyses and stability assays confirmed the predictions, and vaccinated animals generated improved neutralizing-antibody responses to stabilized particles compared to parental viruses and wild-type capsids.

Evaluation and use of in-silico structure-based epitope prediction with foot-and-mouth disease virus

Understanding virus antigenicity is of fundamental importance for the development of better, more cross-reactive vaccines. However, as far as we are aware, no systematic work has yet been conducted using the 3D structure of a virus to identify novel epitopes. Therefore we have extended several existing structural prediction algorithms to build a method for identifying epitopes on the appropriate outer surface of intact virus capsids (which are structurally different from globular proteins in both shape and arrangement of multiple repeated elements) and applied it here as a proof of principle concept to the capsid of foot-and-mouth disease virus (FMDV). We have analysed how reliably several freely available structure-based B cell epitope prediction programs can identify already known viral epitopes of FMDV in the context of the viral capsid. To do this we constructed a simple objective metric to measure the sensitivity and discrimination of such algorithms. After optimising the parameters for five methods using an independent training set we used this measure to evaluate the methods. Individually any one algorithm performed rather poorly (three performing better than the other two) suggesting that there may be value in developing virus-specific software. Taking a very conservative approach requiring a consensus between all three top methods predicts a number of previously described antigenic residues as potential epitopes on more than one serotype of FMDV, consistent with experimental results. The consensus results identified novel residues as potential epitopes on more than one serotype. These include residues 190-192 of VP2 (not previously determined to be antigenic), residues 69-71 and 193-197 of VP3 spanning the pentamer-pentamer interface, and another region incorporating residues 83, 84 and 169-174 of VP1 (all only previously experimentally defined on serotype A). The computer programs needed to create a semi-automated procedure for carrying out this epitope prediction method are presented.

Nearest-cell: a fast and easy tool for locating crystal matches in the PDB

When embarking upon X-ray diffraction data collection from a potentially novel macromolecular crystal form, it can be useful to ascertain whether the measured data reflect a crystal form that is already recorded in the Protein Data Bank and, if so, whether it is part of a large family of related structures. Providing such information to crystallographers conveniently and quickly, as soon as the first images have been recorded and the unit cell characterized at an X-ray beamline, has the potential to save time and effort as well as pointing to possible search models for molecular replacement. Given an input unit cell, and optionally a space group, Nearest-cell rapidly scans the Protein Data Bank and retrieves near-matches.

xtalPiMS: a PiMS-based web application for the management and monitoring of crystallization trials

A major advance in protein structure determination has been the advent of nanolitre-scale crystallization and (in a high-throughput environment) the development of robotic systems for storing and imaging crystallization trials. Most of these trials are carried out in 96-well (or higher density) plates and managing them is a significant information management challenge. We describe xtalPiMS, a web-based application for the management and monitoring of crystallization trials. xtalPiMS has a user-interface layer based on the standards of the Protein Information Management System (PiMS) and a database layer which links the crystallization trial images to the meta-data associated with a particular crystallization trial. The user interface has been optimized for the efficient monitoring of high-throughput environments with three different automated imagers and work to support a fourth imager is in progress, but it can even be of use without robotics. The database can either be a PiMS database or a legacy database for which a suitable mapping layer has been developed.

Recording information on protein complexes in an information management system

The Protein Information Management System (PiMS) is a laboratory information management system (LIMS) designed for use with the production of proteins in a research environment. The software is distributed under the CCP4 licence, and so is available free of charge to academic laboratories. Like most LIMS, the underlying PiMS data model originally had no support for protein–protein complexes. To support the SPINE2-Complexes project the developers have extended PiMS to meet these requirements. The modifications to PiMS, described here, include data model changes, additional protocols, some user interface changes and functionality to detect when an experiment may have formed a complex. Example data are shown for the production of a crystal of a protein complex. Integration with SPINE2-Complexes Target Tracker application is also described.

The Protein Information Management System (PiMS): a generic tool for any structural biology research laboratory

The Protein Information Management System (PiMS) is described together with a discussion of how its features make it well suited to laboratories of all sizes.

The techniques used in protein production and structural biology have been developing rapidly, but techniques for recording the laboratory information produced have not kept pace. One approach is the development of laboratory information-management systems (LIMS), which typically use a relational database schema to model and store results from a laboratory workflow. The underlying philosophy and implementation of the Protein Information Management System (PiMS), a LIMS development specifically targeted at the flexible and unpredictable workflows of protein-production research laboratories of all scales, is described. PiMS is a web-based Java application that uses either Postgres or Oracle as the underlying relational database-management system. PiMS is available under a free licence to all academic laboratories either for local installation or for use as a managed service.

An ion-channel modulator from the saliva of the brown ear tick has a highly modified Kunitz/BPTI structure

Ra-KLP, a 75 amino acid protein secreted by the salivary gland of the brown ear tick Rhipicephalus appendiculatus has a sequence resembling those of Kunitz/BPTI proteins. We report the detection, purification and characterization of the function of Ra-KLP. In addition, determination of the three-dimensional crystal structure of Ra-KLP at 1.6 A resolution using sulphur single-wavelength anomalous dispersion reveals that much of the loop structure of classical Kunitz domains, including the protruding protease-binding loop, has been replaced by beta-strands. Even more unusually, the N-terminal portion of the polypeptide chain is pinned to the "Kunitz head" by two disulphide bridges not found in classical Kunitz/BPTI proteins. The disulphide bond pattern has been further altered by the loss of the bridge that normally stabilizes the protease-binding loop. Consistent with the conversion of this loop into a beta-strand, Ra-KLP shows no significant anti-protease activity; however, it activates maxiK channels in an in vitro system, suggesting a potential mechanism for regulating host blood supply during feeding.

The crystal structure of UMP kinase from Bacillus anthracis (BA1797) reveals an allosteric nucleotide-binding site

Uridine monophosphate (UMP) kinase is a conserved enzyme that catalyzes the ATP-driven conversion of uridylate monophosphate into uridylate diphosphate, an essential metabolic step. In prokaryotes, the enzyme exists as a homohexamer that is regulated by various metabolites. Whereas the enzymatic mechanism of UMP kinase (UK) is well-characterized, the molecular basis of its regulation remains poorly understood. Here we report the crystal structure of UK from Bacillus anthracis (BA1797) in complex with ATP at 2.82 A resolution. It reveals that the cofactor, in addition to binding in the active sites, also interacts with separate binding pockets located near the center of the hexameric structure. The existence of such an allosteric binding site had been predicted by biochemical studies, but it was not identified in previous crystal structures of prokaryotic UKs. We show that this putative allosteric pocket is conserved across different bacterial species, suggesting that it is a feature common to bacterial UKs, and we present a structural model for the allosteric regulation of this enzyme.

Structures of an alanine racemase from Bacillus anthracis (BA0252) in the presence and absence of (R)-1-aminoethylphosphonic acid (L-Ala-P)

Bacillus anthracis, the causative agent of anthrax, has been targeted by the Oxford Protein Production Facility to validate high-throughput protocols within the Structural Proteomics in Europe project. As part of this work, the structures of an alanine racemase (BA0252) in the presence and absence of the inhibitor (R)-1-aminoethylphosphonic acid (L-Ala-P) have determined by X-ray crystallography to resolutions of 2.1 and 1.47 A, respectively. Difficulties in crystallizing this protein were overcome by the use of reductive methylation. Alanine racemase has attracted much interest as a possible target for anti-anthrax drugs: not only is D-alanine a vital component of the bacterial cell wall, but recent studies also indicate that alanine racemase, which is accessible in the exosporium, plays a key role in inhibition of germination in B. anthracis. These structures confirm the binding mode of L-Ala-P but suggest an unexpected mechanism of inhibition of alanine racemase by this compound and could provide a basis for the design of improved alanine racemase inhibitors with potential as anti-anthrax therapies.

Structure of 5-formyltetrahydrofolate cyclo-ligase from Bacillus anthracis (BA4489)

Bacillus anthracis is a spore-forming bacterium and the causative agent of the disease anthrax. The Oxford Protein Production Facility has been targeting proteins from B. anthracis in order to develop high-throughput technologies within the Structural Proteomics in Europe project. As part of this work, the structure of 5-formyltetrahydrofolate cyclo-ligase (BA4489) has been determined by X-ray crystallography to 1.6 A resolution. The structure, solved in complex with magnesium-ion-bound ADP and phosphate, gives a detailed picture of the proposed catalytic mechanism of the enzyme. Chemical differences from other cyclo-ligase structures close to the active site that could be exploited to design specific inhibitors are also highlighted.

SPINE high-throughput crystallization, crystal imaging and recognition techniques: current state, performance analysis, new technologies and future aspects

This paper reviews the developments in high-throughput and nanolitre-scale protein crystallography technologies within the remit of workpackage 4 of the Structural Proteomics In Europe (SPINE) project since the project's inception in October 2002. By surveying the uptake, use and experience of new technologies by SPINE partners across Europe, a picture emerges of highly successful adoption of novel working methods revolutionizing this area of structural biology. Finally, a forward view is taken of how crystallization methodologies may develop in the future.

Application of high-throughput technologies to a structural proteomics-type analysis ofBacillus anthracis

A collaborative project between two Structural Proteomics In Europe (SPINE) partner laboratories, York and Oxford, aimed at high-throughput (HTP) structure determination of proteins from Bacillus anthracis, the aetiological agent of anthrax and a biomedically important target, is described. Based upon a target-selection strategy combining ;low-hanging fruit' and more challenging targets, this work has contributed to the body of knowledge of B. anthracis, established and developed HTP cloning and expression technologies and tested HTP pipelines. Both centres developed ligation-independent cloning (LIC) and expression systems, employing custom LIC-PCR, Gateway and In-Fusion technologies, used in combination with parallel protein purification and robotic nanolitre crystallization screening. Overall, 42 structures have been solved by X-ray crystallography, plus two by NMR through collaboration between York and the SPINE partner in Utrecht. Three biologically important protein structures, BA4899, BA1655 and BA3998, involved in tRNA modification, sporulation control and carbohydrate metabolism, respectively, are highlighted. Target analysis by biophysical clustering based on pI and hydropathy has provided useful information for future target-selection strategies. The technological developments and lessons learned from this project are discussed. The success rate of protein expression and structure solution is at least in keeping with that achieved in structural genomics programs.

Honing thein silicotoolkit for detecting protein disorder

Not all proteins form well defined three-dimensional structures in their native states. Some amino-acid sequences appear to strongly favour the disordered state, whereas some can apparently transition between disordered and ordered states under the influence of changes in the biological environment, thereby playing an important role in processes such as signalling. Although important biologically, for the structural biologist disordered regions of proteins can be disastrous even preventing successful structure determination. The accurate prediction of disorder is therefore important, not least for directing the design of expression constructs so as to maximize the chances of successful structure determination. Such design criteria have become integral to the construct-design strategies of laboratories within the Structural Proteomics In Europe (SPINE) consortium. This paper assesses the current state of the art in disorder prediction in terms of prediction reliability and considers how best to use these methods to guide construct design. Finally, it presents a brief discussion as to how methods of prediction might be improved in the future.

SPINE bioinformatics and data-management aspects of high-throughput structural biology

SPINE (Structural Proteomics In Europe) was established in 2002 as an integrated research project to develop new methods and technologies for high‐throughput structural biology. Development areas were broken down into workpackages and this article gives an overview of ongoing activity in the bioinformatics workpackage. Developments cover target selection, target registration, wet and dry laboratory data management and structure annotation as they pertain to high‐throughput studies. Some individual projects and developments are discussed in detail, while those that are covered elsewhere in this issue are treated more briefly. In particular, this overview focuses on the infrastructure of the software that allows the experimentalist to move projects through different areas that are crucial to high‐throughput studies, leading to the collation of large data sets which are managed and eventually archived and/or deposited.

SPINE workshop on automated X-ray analysis: a progress report

The Structural Proteomics In Europe (SPINE) consortium contained a workpackage to address the automated X-ray analysis of macromolecules. The aim of this workpackage was to increase the throughput of three-dimensional structures while maintaining the high quality of conventional analyses. SPINE was able to bring together developers of software with users from the partner laboratories. Here, the results of a workshop organized by the consortium to evaluate software developed in the member laboratories against a set of bacterial targets are described. The major emphasis was on molecular-replacement suites, where automation was most advanced. Data processing and analysis, use of experimental phases and model construction were also addressed, albeit at a lower level.

Computational model for the IGF-II/IGF2r complex that is predictive of mutational and surface plasmon resonance data

Insulin-like growth factors (IGFs) are key regulators of cell proliferation, differentiation, and transformation, and are thus pivotal in cancer, especially breast, prostate, and colon neoplasm. Their potent mitogenic and anti-apoptotic actions depend primarily on their availability to bind to the signaling IGF cell surface receptors. One mechanism by which IGF-II availability is thought to be modulated is by binding to the nonsignaling IGF-II receptor (IGF2R). This binding is essentially mediated by domain 11 in the multidomain IGF2R extracellular region. The crystal structure of domain 11 of the human IGF-II receptor (IGF2R-d11) has identified a putative IGF-II binding site, and a nuclear magnetic resonance (NMR) solution structure for the IGF-II ligand has also been characterized. These structures have now been used to model in silico the protein-protein interaction between IGF-II and IGF2R-d11 using the program 3D-Dock. Because the IGF-II data comprise an ensemble of 20 structures, all of which satisfy the NMR constraints, the docking procedure was applied to each member of the ensemble. Only those models in which residue Ile1572 of IGF2R-d11, known to be essential for the binding of IGF-II, was at the interface were considered further. These plausible complexes were then critically assessed using an array of analysis techniques including consideration of additional mutagenesis data. One model was strongly supported by these analyses and is discussed here in detail. Furthermore, we demonstrate in vitro experimental support for this model by studying the binding of chimeras of IGF-I and IGF-II to IGF2R fragments.

Structure determination of human semaphorin 4D as an example of the use of MAD in non-optimal cases

Semaphorins are an important class of signalling molecules involved in axon guidance, immune function and angiogenesis. They are characterized by having an extracellular sema domain of about 500 residues. The steps involved in the determination of the structure of human semaphorin 4D are described here as a case study of selenium MAD phasing in a difficult case with low symmetry, moderate diffraction and low selenium content. A particular feature of this study was the large number of diffraction images required to give data of sufficient quality for structure determination and these data are re-analyzed here to investigate the effects of radiation damage on eventual data quality and to suggest strategies for successful MAD phasing in similar difficult cases.

Benefits of automated crystallization plate tracking, imaging, and analysis

We describe the design of a database and software for managing and organizing protein crystallization data. We also outline the considerations behind the design of a fast web interface linking protein production data, crystallization images, and automated image analysis. The database and associated interfaces underpin the Oxford Protein Production Facility (OPPF) crystallization laboratory, collecting, in a routine and automatic manner, up to 100,000 images per day. Over 17 million separate images are currently held in this database. We discuss the substantial scientific benefits automated tracking, imaging, and analysis of crystallizations offers to the structural biologist: analysis of the time course of the trial and easy analysis of trials with related crystallization conditions. Features of this system address requirements common to many crystallographic laboratories that are currently setting up (semi-)automated crystallization imaging systems.

Structural and kinetic basis for heightened immunogenicity of T cell vaccines

Analogue peptides with enhanced binding affinity to major histocompatibility class (MHC) I molecules are currently being used in cancer patients to elicit stronger T cell responses. However, it remains unclear as to how alterations of anchor residues may affect T cell receptor (TCR) recognition. We correlate functional, thermodynamic, and structural parameters of TCR–peptide–MHC binding and demonstrate the effect of anchor residue modifications of the human histocompatibility leukocyte antigens (HLA)–A2 tumor epitope NY–ESO-1_157–165–SLLMWITQC on TCR recognition. The crystal structure of the wild-type peptide complexed with a specific TCR shows that TCR binding centers on two prominent, sequential, peptide sidechains, methionine–tryptophan. Cysteine-to-valine substitution at peptide position 9, while optimizing peptide binding to the MHC, repositions the peptide main chain and generates subtly enhanced interactions between the analogue peptide and the TCR. Binding analyses confirm tighter binding of the analogue peptide to HLA–A2 and improved soluble TCR binding. Recognition of analogue peptide stimulates faster polarization of lytic granules to the immunological synapse, reduces dependence on CD8 binding, and induces greater numbers of cross-reactive cytotoxic T lymphocyte to SLLMWITQC. These results provide important insights into heightened immunogenicity of analogue peptides and highlight the importance of incorporating structural data into the process of rational optimization of superagonist peptides for clinical trials.

RONN: the bio-basis function neural network technique applied to the detection of natively disordered regions in proteins

MOTIVATION

Recent studies have found many proteins containing regions that do not form well-defined three-dimensional structures in their native states. The study and detection of such disordered regions is important both for understanding protein function and for facilitating structural analysis since disordered regions may affect solubility and/or crystallizability.

RESULTS

We have developed the regional order neural network (RONN) software as an application of our recently developed 'bio-basis function neural network' pattern recognition algorithm to the detection of natively disordered regions in proteins. The results of blind-testing a panel of nine disorder prediction tools (including RONN) against 80 protein sequences derived from the Protein Data Bank shows that, based on the probability excess measure, RONN performed the best.

A procedure for setting up high-throughput nanolitre crystallization experiments. Crystallization workflow for initial screening, automated storage, imaging and optimization

Crystallization trials at the Division of Structural Biology in Oxford are now almost exclusively carried out using a high‐throughput workflow implemented in the Oxford Protein Production Facility. Initial crystallization screening is based on nanolitre‐scale sitting‐drop vapour‐diffusion experiments (typically 100 nl of protein plus 100 nl of reservoir solution per droplet) which use standard crystallization screening kits and 96‐well crystallization plates. For 294 K crystallization trials the barcoded crystallization plates are entered into an automated storage system with a fully integrated imaging system. These plates are imaged in accordance with a pre‐programmed schedule and the resulting digital data for each droplet are harvested into a laboratory information‐management system (LIMS), scored by crystal recognition software and displayed for user analysis via a web‐based interface. Currently, storage for trials at 277 K is not automated and for imaging the crystallization plates are fed by hand into an imaging system from which the data enter the LIMS. The workflow includes two procedures for nanolitre‐scale optimization of crystallization conditions: (i) a protocol for variation of pH, reservoir dilution and protein:reservoir ratio and (ii) an additive screen. Experience based on 592 crystallization projects is reported.

The sema domain

The sema domain was first defined from sequence by Kolodkin and colleagues in the early 1990s, and constitutes the distinctive structural and functional element of semaphorins, their plexin receptors and the receptor tyrosine kinases MET and RON, three protein families with major roles in development, tissue regeneration and cancer. Recently determined crystal structures of two semaphorins (SEMA3A and SEMA4D) and the MET receptor have shown that the sema domain consists of a highly conserved variant form of the seven-blade beta-propeller fold. The structures, however, also suggest differences between these families with respect to the mode of dimerisation and the regions of the domain involved in ligand-receptor interactions. This reflects the considerable plasticity and adaptation of the sema domain in order to meet different binding requirements, properties that may underlie the vast array of ligand-receptor specificities and functions of the semaphorin superfamily.

Polyalanine reconstruction from Calpha positions using the program CALPHA can aid initial phasing of data by molecular replacement procedures

The Calpha positions for a protein can provide a scaffold for the reconstruction of more complete models. Reconstructions can be by manual rebuilding, from geometric solutions to the constraints on main-chain torsion angles or from databases of known protein structures. The last method is usually the most convenient and reliable. This paper describes a database reconstruction program, CALPHA, and assesses its accuracy and reliability by test reconstructions of well refined structures. Typically, backbone atoms are repositioned to within 0.3 A of their original positions. This corresponds to regenerating main-chain torsion angles to within 15 degrees. Uses of CALPHA for automating refinement procedures are discussed. In particular, the uses of Calpha-only and rcconstructed polyalanine models of HIV-1 reverse transcriptase for cross-rotation and translation function searches are compared. The CALPHA polyalanine model is found to provide more selectivity for approximately correct orientations. The effect on the translation function is dependent on the resolution shell employed. It is expected that these observations will be applicable in other cases.

Crystal structure of a soluble CD28-Fab complex

Naive T cell activation requires signaling by the T cell receptor and by nonclonotypic cell surface receptors. The most important costimulatory protein is the monovalent homodimer CD28, which interacts with CD80 and CD86 expressed on antigen-presenting cells. Here we present the crystal structure of a soluble form of CD28 in complex with the Fab fragment of a mitogenic antibody. Structural comparisons redefine the evolutionary relationships of CD28-related proteins, antigen receptors and adhesion molecules and account for the distinct ligand-binding and stoichiometric properties of CD28 and the related, inhibitory homodimer CTLA-4. Cryo-electron microscopy-based comparisons of complexes of CD28 with mitogenic and nonmitogenic antibodies place new constraints on models of antibody-induced receptor triggering. This work completes the initial structural characterization of the CD28-CTLA-4-CD80-CD86 signaling system.

The ligand-binding face of the semaphorins revealed by the high-resolution crystal structure of SEMA4D

Semaphorins, proteins characterized by an extracellular sema domain, regulate axon guidance, immune function and angiogenesis. The crystal structure of SEMA4D (residues 1-657) shows the sema topology to be a seven-bladed beta-propeller, revealing an unexpected homology with integrins. The sema beta-propeller contains a distinctive 77-residue insertion between beta-strands C and D of blade 5. Blade 7 is followed by a domain common to plexins, semaphorins and integrins (PSI domain), which forms a compact cysteine knot abutting the side of the propeller, and an Ig-like domain. The top face of the beta-propeller presents prominent loops characteristic of semaphorins. In addition to limited contact between the Ig-like domains, the homodimer is stabilized through extensive interactions between the top faces in a sector of the beta-propeller used for heterodimerization in integrins. This face of the propeller also mediates ligand binding in integrins, and functional data for semaphorin-receptor interactions map to the equivalent surface.

Missense mutations in the homeodomain of HOXD13 are associated with brachydactyly types D and E

HOXD13, the most 5' gene of the HOXD cluster, encodes a homeodomain transcription factor with important functions in limb patterning and growth. Heterozygous mutations of human HOXD13, encoding polyalanine expansions or frameshifts, are believed to act by dominant negative or haploinsufficiency mechanisms and are predominantly associated with synpolydactyly phenotypes. Here, we describe two mutations of HOXD13 (923C-->G encoding Ser308Cys and 940A-->C encoding Ile314Leu) that cause missense substitutions within the homeodomain. Both are associated with distinctive limb phenotypes in which brachydactyly of specific metacarpals, metatarsals, and phalangeal bones is the most constant feature, exhibiting overlap with brachydactyly types D and E. We investigated the binding of synthetic mutant proteins to double-stranded DNA targets in vitro. No consistent differences were found for the Ser308Cys mutation compared with the wild type, but the Ile314Leu mutation (which resides at the 47th position of the homeodomain) exhibited increased affinity for a target containing the core recognition sequence 5'-TTAC-3' but decreased affinity for a 5'-TTAT-3' target. Molecular modeling of the Ile314Leu mutation indicates that this mixed gain and loss of affinity may be accounted for by the relative positions of methyl groups in the amino acid side chain and target base.

Modeling study of human renal chloride channel (hCLC-5) mutations suggests a structural-functional relationship

BACKGROUND

Dent's disease, a renal tubular disorder characterized by low-molecular-weight proteinuria, hypercalciuria, and nephrolithiasis, is due to inactivating mutations in the X-linked renal-specific chloride channel, hCLC-5. The x-ray crystal structures of two bacterial chloride channels (CLCs) have recently been established, thereby allowing us to construct a model for hCLC-5 and further examine the role of its mutations.

METHODS

The data regarding 49 hCLC-5 mutations were reviewed. Thirty-four mutations that predicted absent or truncated channels were excluded. The remaining 15 mutations (one in-frame insertion and 14 missense mutations), 12 of which have been studied electrophysiologically, were assessed. The hCLC-5 sequence was aligned with the Salmonella typhimurium and Escherichia coli sequences and used to map the hCLC-5 mutations onto a three-dimensional model.

RESULTS

hCLC-5 is a homodimeric protein, with each subunit consisting of 18 helices. None of the missense mutations involved the chloride (Cl-) selectivity filter, but 12 of the 15 mutations were found to be clustered at the interface of the two subunits. Six of these mutations occurred in two of the helices that either form part of the interface or lie in close proximity to the interface, and three other mutations that did not lead to complete loss of Cl- conductance were at the edge of the interface.

CONCLUSION

These results demonstrate a crucial role for the interaction between the two subunits at the interface of the homodimeric hCLC-5.

Structure of a functional IGF2R fragment determined from the anomalous scattering of sulfur

Insulin-like growth factor II receptor (IGF2R) is a multifunctional cell surface receptor implicated in tumour suppression. Its growth inhibitory activity has been associated with an ability to bind IGF-II. IGF2R contains 15 homologous extracellular domains, with domain 11 primarily responsible for IGF-II binding. We report a 1.4 A resolution crystal structure of domain 11, solved using the anomalous scattering signal of sulfur. The structure consists of two crossed beta-sheets forming a flattened beta-barrel. Structural analysis identifies the putative IGF-II binding site at one end of the beta-barrel whilst crystal lattice contacts suggest a model for the full-length IGF2R extracellular region. The structure factors and coordinates of IGF2R domain 11 have been deposited in the Protein Data Bank (accession codes 1GP0 and 1GP3).

Gene polymorphism in Netherton and common atopic disease

Atopic dermatitis (AD) and asthma are characterized by IgE-mediated atopic (allergic) responses to common proteins (allergens), many of which are proteinases. Loci influencing atopy have been localized to a number of chromosomal regions, including the chromosome 5q31 cytokine cluster. Netherton disease is a rare recessive skin disorder in which atopy is a universal accompaniment. The gene underlying Netherton disease (SPINK5) encodes a 15-domain serine proteinase inhibitor (LEKTI) which is expressed in epithelial and mucosal surfaces and in the thymus. We have identified six coding polymorphisms in SPINK5 (Table 1) and found that a Glu420-->Lys variant shows significant association with atopy and AD in two independent panels of families. Our results implicate a previously unrecognized pathway for the development of common allergic illnesses.

Presence of 2',5'-Bis-O-(tert-butyldimethylsilyl)-3'-spiro-5"-(4"-amino-1",2"-oxath iole-2",2"-dioxide) (TSAO)-resistant virus strains in TSAO-inexperienced HIV patients

HIV-1 samples from six patients undergoing diverse anti-HIV therapies possessed the E138A mutation in their reverse transcriptase (RT) genome. Patients were receiving the following therapies: TIBO monotherapy (one patient); zidovudine plus didanosine combination therapy (one); zidovudine monotherapy (one); sequential therapy with zidovudine, then stavudine and finally zalcitabine plus didanosine (one); and two were drug naive. E138K, not E138A, is a known TSAO-specific resistance mutation, emerging under selective pressure in vitro. Our phenotypic data on the patient isolates, confirmed by data on an E138A mutant acquired through in vitro mutagenesis, indicated that an alanine substitution for glutamate at codon 138 of the HIV-1 RT renders the virus TSAO resistant, confirming the importance of this amino acid residue in the activity of TSAO derivatives. In addition, we have demonstrated through phenotypic analysis of the E138A and A98S mutants (after in vitro mutagenesis) that the mutation A98S, found in one of these patients, could be partially responsible for the phenotypic reversal of TSAO resistance. This reversal could be explained by the restoration of a hydrogen bond between 98S and the main-chain residue L349, which compensates for the loss of the E138-G99 main-chain hydrogen bond. As TSAO derivatives have not been used in the clinical setting, the presence of the E138A mutation at a frequency of 6.7% in our study of 90 TSAO-inexperienced HIV-seropositive individuals implies that 138A of the RT must be a natural variant and that the mutant virus is replication competent. Our observations suggest that the E138A mutation may likely arise in patients under the selective pressure of TSAO or related compounds that show a decreased antiviral potency toward the E138A variant.

Phenylethylthiazolylthiourea (PETT) non-nucleoside inhibitors of HIV-1 and HIV-2 reverse transcriptases. Structural and biochemical analyses

Most non-nucleoside reverse transcriptase (RT) inhibitors are specific for HIV-1 RT and demonstrate minimal inhibition of HIV-2 RT. However, we report that members of the phenylethylthiazolylthiourea (PETT) series of non-nucleoside reverse transcriptase inhibitors showing high potency against HIV-1 RT have varying abilities to inhibit HIV-2 RT. Thus, PETT-1 inhibits HIV-1 RT with an IC(50) of 6 nM but shows only weak inhibition of HIV-2 RT, whereas PETT-2 retains similar potency against HIV-1 RT (IC(50) of 5 nM) and also inhibits HIV-2 RT (IC(50) of 2.2 microM). X-ray crystallographic structure determinations of PETT-1 and PETT-2 in complexes with HIV-1 RT reveal the compounds bind in an overall similar conformation albeit with some differences in their interactions with the protein. To investigate whether PETT-2 could be acting at a different site on HIV-2 RT (e.g. the dNTP or template primer binding site), we compared modes of inhibition for PETT-2 against HIV-1 and HIV-2 RT. PETT-2 was a noncompetitive inhibitor with respect to the dGTP substrate for both HIV-1 and HIV-2 RTs. PETT-2 was also a noncompetitive inhibitor with respect to a poly(rC).(dG) template primer for HIV-2 RT. These results are consistent with PETT-2 binding in corresponding pockets in both HIV-1 and HIV-2 RT with amino acid sequence differences in HIV-2 RT affecting the binding of PETT-2 compared with PETT-1.

Crystallographic analysis of the binding modes of thiazoloisoindolinone non-nucleoside inhibitors to HIV-1 reverse transcriptase and comparison with modeling studies

We have determined the crystal structures of thiazoloisoindolinone non-nucleoside inhibitors in complex with HIV-1 reverse transcriptase to high-resolution limits of 2.7 A (BM +21.1326) and 2. 52 A (BM +50.0934). We find that the binding modes of this series of inhibitors closely resemble that of "two-ring" non-nucleoside reverse transcriptase inhibitors. The structures allow rationalization of stereochemical requirements, structure-activity data, and drug resistance data. Comparisons with our previous structures suggest modifications to the inhibitors that might improve resilience to drug-resistant mutant forms of reverse transcriptase. Comparison with earlier modeling studies reveals that the predicted overlap of thiazoloisoindolinones with TIBO was largely correct, while that with nevirapine was significantly different.

Further additions to MolScript version 1.4, including reading and contouring of electron-density maps

MolScript is one of the most popular programs for the generation of publication-quality figures and the recent re-working of the program should ensure its continued popularity. However, some functionality of particular interest to crystallographers is not part of the standard program. A modified MolScript version 1.4 has been described previously, with more flexible colouring schemes among its new features. This report describes further enhancements to MolScript version 1.4, including facilities for drawing rods for helices and ribbons for oligonucleotides and allowing several formats of electron-density maps to be read and contoured using either lines or smoothed triangulated surfaces.

Mutational analysis of Tyr-318 within the non-nucleoside reverse transcriptase inhibitor binding pocket of human immunodeficiency virus type I reverse transcriptase

The highly conserved Tyr-318 is part of the non-nucleoside reverse transcriptase inhibitor (NNRTI)-specific lipophilic pocket of human immunodeficiency virus type I reverse transcriptase (RT) and makes contact within 4 A with the NNRTIs in all reported RT/NNRTI complexes. Using site-directed mutagenesis, six mutant RTs were constructed bearing the mutations Y318H, Y318K, Y318L, Y318C, Y318W, and Y318F. We found that only the Y318W and Y318F mutant RTs retained substantial RT activity, whereas the catalytic activities of the Y318K, Y318C, Y318H, and Y318L RT mutants were less than 5% of the wild-type activity. The Y318F mutant RT retained substantial sensitivity to the majority of NNRTIs tested, whereas the Y318W mutant RT showed varying degrees of resistance to NNRTIs. Subunit-specific site-directed mutagenesis revealed that there was no difference in the catalytic activity or resistance/sensitivity spectrum toward NNRTIs regardless of whether the Tyr-318 mutation was introduced in both subunits or only in the p66 subunit of RT. Recombinant viruses harboring the Y318F or Y318W mutation in the RT showed a similar resistance/sensitivity pattern to NNRTIs as their corresponding 318 mutant recombinant RTs. Our findings stress a functional or structural role for Tyr-318 in wild-type RT and argue for the design of novel NNRTIs that interact more closely with this amino acid in the NNRTI-specific pocket of human immunodeficiency virus type I RT.

Crystal structures of HIV-1 reverse transcriptase in complex with carboxanilide derivatives

The carboxanilides are nonnucleoside inhibitors (NNIs) of HIV-1 reverse transcriptase (RT), of potential clinical importance. The compounds differ in potency and in their retention of potency in the face of drug resistance mutations. Whereas UC-84, the prototype compound, only weakly inhibits many RTs bearing single point resistance mutations, inhibition by UC-781 is little affected. It has been proposed that UC-38 and UC-781 may form quaternary complexes with RT at a site other than the known binding pocket of other NNIs. X-ray crystal structures of four HIV-1 RT-carboxanilide complexes (UC-10, UC-38, UC-84, and UC-781) reported here reveal that all four inhibitors bind in the usual NNI site, forming binary 1:1 complexes with RT in the absence of substrates with the amide/thioamide bond in cis conformations. For all four complexes the anilide rings of the inhibitors overlap aromatic rings of many other NNIs bound to RT. In contrast, the second rings of UC-10, UC-84, and UC-781 do not bind in equivalent positions to those of other "two-ring" NNIs such as alpha-APA or HEPT derivatives. The binding modes most closely resemble that of the structurally dissimilar NNI, Cl-TIBO, with a common hydrogen bond between each carboxanilide NH- group and the main-chain carbonyl oxygen of Lys101. The binding modes differ slightly between the UC-10/UC-781 and UC-38/UC-84 pairs of compounds, apparently related to the shorter isopropylmethanoyl substituents of the anilide rings of UC-38/UC-84, which draws these rings closer to residues Tyr181 and Tyr188. This in turn explains the differences in the effect of mutated residues on the binding of these compounds.

Continuous and discontinuous changes in the unit cell of HIV-1 reverse transcriptase crystals on dehydration

A crystal form of HIV-1 reverse transcriptase (RT) complexed with inhibitors showed diffraction to a high-resolution limit of 3.7 A. Instability in the unit-cell dimensions of these crystals was observed during soaking experiments, but the range of this variability and consequent change in lattice order was revealed by a chance observation of dehydration. Deliberately induced dehydration results in crystals having a variety of unit cells, the best-ordered of which show diffraction to a minimum Bragg spacing of 2.2 A. In order to understand the molecular basis for this phenomenon, the initial observation of dehydration, the data sets from dehydrated crystals, the crystal packing and the domain conformation of RT are analysed in detail here. This analysis reveals that the crystals undergo remarkable changes following a variety of possible dehydration pathways: some changes occur gradually whilst others are abrupt and require significant domain rearrangements. Comparison of domain arrangements in different crystal forms gives insight into the flexibility of RT which, in turn, may reflect the internal motions allowing this therapeutically important enzyme to fulfill its biological function.

3'-Azido-3'-deoxythymidine drug resistance mutations in HIV-1 reverse transcriptase can induce long range conformational changes

HIV reverse transcriptase (RT) is one of the main targets for the action of anti-AIDS drugs. Many of these drugs [e.g., 3'-azido-3'-deoxythymidine (AZT) and 2',3'-dideoxyinosine (ddI)] are analogues of the nucleoside substrates used by the HIV RT. One of the main problems in anti-HIV therapy is the selection of a mutant virus with reduced drug sensitivity. Drug resistance in HIV is generated for nucleoside analogue inhibitors by mutations in HIV RT. However, most of these mutations are situated some distance from the polymerase active site, giving rise to questions concerning the mechanism of resistance. To understand the possible structural bases for this, the crystal structures of AZT- and ddI-resistant RTs have been determined. For the ddI-resistant RT with a mutation at residue 74, no significant conformational changes were observed for the p66 subunit. In contrast, for the AZT-resistant RT (RTMC) bearing four mutations, two of these (at 215 and 219) give rise to a conformational change that propagates to the active site aspartate residues. Thus, these drug resistance mutations produce an effect at the RT polymerase site mediated simply by the protein. It is likely that such long-range effects could represent a common mechanism for generating drug resistance in other systems.

A proline-to-histidine substitution at position 225 of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) sensitizes HIV-1 RT to BHAP U-90152

Two mutant virus strains in which the novel P225H mutation appeared in a V106A reverse transcriptase (RT)-mutated genetic background upon treatment of human immunodeficiency virus type 1 (HIV-1) with quinoxaline S-2720 were isolated. Surprisingly, the addition of the P225H mutation to the V106A RT mutant genetic background resensitized the V106A RT mutant virus to the non-nucleoside RT inhibitor (NNRTI) BHAP U-90152, but not to other NNRTIs. Construction of both recombinant viruses and recombinant RTs containing the V106A, P225H and V106A+P225H mutations revealed that P225H was indeed responsible for the marked potentiation of the antiviral activity of BHAP against the P225H single-mutant virus and the V106A+P225H double-mutant virus when compared to wild-type and V106A single-mutant viruses, respectively. An explanation for the markedly increased sensitivity of the P225H mutant HIV-1 RT to BHAP and not to the other NNRTIs was provided by the unique features of the X-ray structure of the RT-BHAP complex.

Models which explain the inhibition of reverse transcriptase by HIV-1-specific (thio)carboxanilide derivatives

The (thio)carboxanilide derivatives are potent and selective inhibitors of HIV-1 reverse transcriptase (RT) and have a favourable antiviral activity spectrum. To understand better their mode of action, and to provide a structural basis for further improvement, models of RT complexed with four (thio)carboxanilide inhibitors (UC781, UC10, UC38 and UC84) have been constructed based on the X-ray structure of RT complexed with 9-chloro-TIBO. In the models, the protein conformation is similar to that of the RT-TIBO complex and the complexes are stabilised by hydrogen bonding between the inhibitors and the main chain oxygen of Lys101. Significant hydrophobic interactions include those with Leu100, Val106, Val179, Tyr188, Phe227, Leu234, and His235. The thiocarboxanilides UC781 and UC10 also make important hydrophobic interactions with Trp229. The models are consistent with the inhibitors' relative antiviral potencies and the observed resistance data. They further predict that mutations to Phe227, Trp229, or Leu234 might confer resistance. Since these are not observed, some constraining structural or functional role for these residues in the active enzyme is suggested.

Unique features in the structure of the complex between HIV-1 reverse transcriptase and the bis(heteroaryl)piperazine (BHAP) U-90152 explain resistance mutations for this nonnucleoside inhibitor

The viral reverse transcriptase (RT) provides an attractive target in the search for anti-HIV therapies. The nonnucleoside inhibitors (NNIs) are a diverse set of compounds (usually HIV-1 specific) that function by distorting the polymerase active site upon binding in a nearby pocket. Despite being potent and of generally low toxicity, their clinical use has been limited by rapid selection for resistant viral populations. The 2.65-A resolution structure of the complex between HIV-1 RT and the bis(heteroaryl)piperazine (BHAP) NNI, 1-(5-methanesulfonamido-1H-indol-2-yl-carbonyl)-4- [3-(1-methyl-ethylamino) pyridinyl] piperazine (U-90152), reveals the inhibitor conformation and bound water molecules. The bulky U-90152 molecule occupies the same pocket as other NNIs, but the complex is stabilized quite differently, in particular by hydrogen bonding to the main chain of Lys-103 and extensive hydrophobic contacts with Pro-236. These interactions rationalize observed resistance mutations, notably Pro-236-Leu, which occurs characteristically for BHAPs. When bound, part of U-90152 protrudes into the solvent creating a channel between Pro-236 and the polypeptide segments 225-226 and 105-106, giving the first clear evidence of the entry mode for NNIs. The structure allows prediction of binding modes for related inhibitors [(altrylamino)piperidine-BHAPs] and suggests changes to U-90152, such as the addition of a 6 amino group to the pyridine ring, which may make binding more resilient to mutations in the RT. The observation of novel hydrogen bonding to the protein main chain may provide lessons for the improvement of quite different inhibitors.

An extensively modified version of MolScript that includes greatly enhanced coloring capabilities

Owing to its flexibility, MolScript has become one of the most widely used programs for generating publication-quality molecular graphics. Integration with the Raster3D package, to allow the production of photorealistic rendered images, has increased its popularity still further. However, this intensive use has shown the need for enhancement of some areas of the program, especially for controlling the coloring of atoms, bonds, and molecules. This work describes a heavily modified version of MolScript that has added syntax for describing complicated coloring schemes and also has new graphics commands. Enhancements include drawing split-bond ball-and-stick models, smoothly varying the color of molecules (color ramping), abrupt color changes within secondary structural units, and the creation of dashed bonds. Making use of these added features is simple because all MolScript syntax is still supported and one typically needs only to add a few control commands. The final section of this article suggests some uses for this modified MolScript and provides illustrative examples.

Complexes of HIV-1 reverse transcriptase with inhibitors of the HEPT series reveal conformational changes relevant to the design of potent non-nucleoside inhibitors

Crystal structures of HIV-1 reverse transcriptase (RT) complexed with a range of chemically diverse non-nucleoside inhibitors (NNIs) have shown a single pocket in which the inhibitors bind and details of the inhibitor-protein interactions. To delineate the structural requirements for an effective inhibitor, we have determined the structures of three closely related NNIs which vary widely in their potencies. Crystal structures of HIV-1 RT complexed with two very potent inhibitors, MKC-442 and TNK-651, at 2.55 angstroms resolution complement our previous analysis of the complex with the less effective inhibitor, HEPT. These structures reveal conformational changes which correlate with changes in potency. We suggest that a major determinant of increased potency in the analogues of HEPT is an improved interaction between residue Tyr181 in the protein and the 6-benzyl ring of the inhibitors which stabilizes the structure of the complex. This arises through a conformational switching of the protein structure triggered by the steric bulk of the 5-substituent of the inhibitor pyrimidine ring.

Crystals of HIV-1 reverse transcriptase diffracting to 2.2 A resolution

Reverse transcriptase (RT) from the human immunodeficiency virus type 1 has been crystallized in four closely related forms, the best of which diffract X-rays to 2.2 A resolution. The RT was crystallized as a complex with a non-nucleoside inhibitor, either nevirapine or a nevirapine analogue. Crystals grew from 6% PEG 3400 buffered at pH 5. These were of space group P2(1)2(1)2(1) with unit cell parameters a = 147 A, b = 112 A, c = 79 A (form A), with one RT heterodimer in the asymmetric unit. Changes in unit cell parameters and degree of crystalline order were observed on soaking pregrown crystals in various solutions, giving three further sets of unit cells. These were a = 143 A, b = 112, A, c = 79 A (form B), a = 141 A, b = 111 A, c = 73 A (form C), a = 143 A, b = 117 A, c = 66.5 A (form D). The last two forms diffract X-rays to 2.2 A resolution. Structure determinations of these latter crystal forms of RT should give a detailed atomic model for this therapeutically important drug target.