CCP4 Cloud for structure determination and project management in macromolecular crystallography

Nowadays, progress in the determination of three-dimensional macromolecular structures from diffraction images is achieved partly at the cost of increasing data volumes. This is due to the deployment of modern high-speed, high-resolution detectors, the increased complexity and variety of crystallographic software, the use of extensive databases and high-performance computing. This limits what can be accomplished with personal, offline, computing equipment in terms of both productivity and maintainability. There is also an issue of long-term data maintenance and availability of structure-solution projects as the links between experimental observations and the final results deposited in the PDB. In this article, CCP4 Cloud, a new front-end of the CCP4 software suite, is presented which mitigates these effects by providing an online, cloud-based environment for crystallographic computation. CCP4 Cloud was developed for the efficient delivery of computing power, database services and seamless integration with web resources. It provides a rich graphical user interface that allows project sharing and long-term storage for structure-solution projects, and can be linked to data-producing facilities. The system is distributed with the CCP4 software suite version 7.1 and higher, and an online publicly available instance of CCP4 Cloud is provided by CCP4.

ARP/wARP for crystallographic model building and drug discovery

The ARP/wARP software project combines automated model building and refinement into an unified approach for macromolecular crystal structure determination. The project is based on two decades of extensive research and development in the areas of macromolecular X-ray crystallography, informatics, data mining and statistical pattern recognition. ARP/wARP collects a vast amount of computationally efficient methods and provides easy-to-use pipelines for building models of proteins, nucleotides, ligands, as well as their complexes. All methods are intuitively accessible from the ArpNavigator [1], which grants direct visualisation and real-time interaction with model building results. Structures determined using ARP/wARP include histones, hsp70, viral proteases, an insect antifreeze protein, transferases, deadenylases, synthases, kinases, photolyases and the spliceosome. The novel release of ARP/wARP, version 7.4, comes with notable innovations for determining structures at medium-to-low resolution such as exploitation of non-crystallographic symmetry, improved protocols for model update and estimation of validity of built models. Joint releases with the CCP4 suite improve software development and integration, and make the installation and updates fast and convenient for the user. A novel procedure for the automatic identification of ligands in electron density maps is introduced. It is based on the sparse parameterisation of density clusters and the matching of the pseudo-atomic grids thus created to conformationally variant ligands using mathematical descriptors of molecular shape, size and topology. The integration of the ViCi web-server for in-silico ligand-based drug design and updated stereo-chemical restraints for ligand fitting make ARP/wARP an asset for crystallographic drug discovery pipelines.

Crambin - record in crystallographic resolution for a biological macromolecule

A couple of years ago a study was presented that set the record for the highest X-ray crystallographic resolution for a biological macromolecule [1]. The structure of the small protein crambin was determined to 0.48 Å resolution - this almost doubled the amount of available experimental data. Crambin is a small protein consisting of 46 amino acids belonging to the thionin family. Although the protein and its structure has long been known, it lacks any obvious enzymatic activity and has a hard-to-guess biological function. The protein crystallizes readily and serves as an excellent specimen for exploring the limits of resolution of the diffraction. The results demonstrated the possibilities that can be offered by a high-energy synchrotron source. The structure refined with Refmac, Shelxl and Mopro revealed a wealth of details. Bonding electron density became visible along the main chain. However, no fundamental additional structural features could be detected in comparison to the previously collected data set to 0.54 Å resolution. The availability of extremely high-resolution data is certainly of great help to drive further software development and methods for data interpretation. The question will always remain as to what the true limits are in terms of what can be seen in a biological macromolecule. Here we will present the results of our recent efforts in interpretation of such ultra-high resolution data.

Towards structural characterization of H. polymorpha telomerase components

Telomeres are regions of non-coding DNA that cap the chromosomes, preventing the loss of coding DNA during cell division and contributing to chromosomal stability. In actively dividing cells, such as embryonic stem cells, the telomeres need to elongated by telomerase. The telomerase complex consist of the enzyme telomerase reverse transcriptase (TERT), telomerase RNA (TR) and additional proteins. TERT and TR are required for the telomerase activity in vitro. Telomerase is active in vast majority of the cancer cells ensuring continuous cell division and tumor growth. Syndromes leading to premature aging are often associated with short telomeres. Finding ways to regulate the telomerase activity would help to advance therapies for these conditions. However, the structural information available of the telomerase complex is very limited. We have chosen thermophilic yeast Hansenula polymorpha as a model system due to the stability of its proteins. The N-terminal domain of the TERT is essential for telomerase activity and possibly is involved in binding of TR, telomeric DNA and additional protein components of the telomerase complex. We have crystallised the N-terminal domain of H. polymorpha TERT and, in lack of a homologious structure, produced a seleno-methionine derivative of the protein. MAD data on N-terminal domain has been collected to resolution of 2.0 Å at the PETRA-III beamline P13 (EMBL/DESY) in Hamburg. We will discuss the structure-function relationship of the N-domain and the whole TERT component.

ValiFrag: Evaluating fragment quality during automated protein model building

X-ray diffraction data from flexible macromolecules and their complexes can rarely be measured to a resolution better than 3 Å. Due to a loss of detectable atomic features, the determination of low-resolution structures is beyond the current operational range of crystallographic software and requires a large amount of manual intervention. ARP/wARP [1] v7.4 generates structures that are up to 80% complete at 3.0 Å, but the completeness drops sharply as the resolution gets worse. Reduction of the model completeness is accompanied with an increase in the number of fragments built, which become shorter. Such fragments are applicable for further model building if they are correct. Though, if they are wrong they may cause the formation of incorrectly built regions in the final model. Thus, there is a need to improve fragment quality before automated model completion is applied. We exploit the vast amount of structural information deposited in the Protein Data Bank (PDB) [2], to make use of it for structural validation of built fragments. Precisely, we evaluate the conformation of each fragment. If the conformation is present in several different protein models in the PDB, it is likely to be modelled correctly in the built model and is accepted. If, on the contrary, it cannot be found in any PDB model, it is probably incorrect. Here we present the software implementation of this validation, called ValiFrag, which checks the validity of automatically built protein chain fragments by evaluating their occurrence in the PDB. Protein models from the PDB were broken into dipeptides and conformational parameters for each of these were then stored in a database. For each automatically built fragment, ValiFrag computes the probability of it to be correct according to the conformation of all possible dipeptides. It can, therefore, assess which fragments are likely to be structurally incorrect and should possibly be modified, or even removed, to improve the final model.

Role of κ→λ light-chain constant-domain switch in the structure and functionality of A17 reactibody

The engineering of catalytic function in antibodies requires precise information on their structure. Here, results are presented that show how the antibody domain structure affects its functionality. The previously designed organophosphate-metabolizing reactibody A17 has been re-engineered by replacing its constant κ light chain by the λ chain (A17λ), and the X-ray structure of A17λ has been determined at 1.95 Å resolution. It was found that compared with A17κ the active centre of A17λ is displaced, stabilized and made more rigid owing to interdomain interactions involving the CDR loops from the VL and VH domains. These VL/VH domains also have lower mobility, as deduced from the atomic displacement parameters of the crystal structure. The antibody elbow angle is decreased to 126° compared with 138° in A17κ. These structural differences account for the subtle changes in catalytic efficiency and thermodynamic parameters determined with two organophosphate ligands, as well as in the affinity for peptide substrates selected from a combinatorial cyclic peptide library, between the A17κ and A17λ variants. The data presented will be of interest and relevance to researchers dealing with the design of antibodies with tailor-made functions.

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.

The TB structural genomics consortium: a resource for Mycobacterium tuberculosis biology

The TB Structural Genomics Consortium is an organization devoted to encouraging, coordinating, and facilitating the determination and analysis of structures of proteins from Mycobacterium tuberculosis. The Consortium members hope to work together with other M. tuberculosis researchers to identify M. tuberculosis proteins for which structural information could provide important biological information, to analyze and interpret structures of M. tuberculosis proteins, and to work collaboratively to test ideas about M. tuberculosis protein function that are suggested by structure or related to structural information. This review describes the TB Structural Genomics Consortium and some of the proteins for which the Consortium is in the progress of determining three-dimensional structures.

Atomic (0.94 A) resolution structure of an inverting glycosidase in complex with substrate

The crystal structure of Clostridium thermocellum endoglucanase CelA in complex with cellopentaose has been determined at 0.94 A resolution. The oligosaccharide occupies six D-glucosyl-binding subsites, three on either side of the scissile glycosidic linkage. The substrate and product of the reaction occupy different positions at the reducing end of the cleft, where an extended array of hydrogen-bonding interactions with water molecules fosters the departure of the leaving group. Severe torsional strain upon the bound substrate forces a distorted boat(2,5) B conformation for the glucosyl residue bound at subsite -1, which facilitates the formation of an oxocarbenium ion intermediate and might favor the breakage of the sugar ring concomitant with catalysis.

Detailed analysis of RNA-protein interactions within the ribosomal protein S8-rRNA complex from the archaeon Methanococcus jannaschii

The crystal structure of ribosomal protein S8 bound to its target 16 S rRNA from a hyperthermophilic archaeon Methanococcus jannaschii has been determined at 2.6 A resolution. The protein interacts with the minor groove of helix H21 at two sites located one helical turn apart, with S8 forming a bridge over the RNA major groove. The specificity of binding is essentially provided by the C-terminal domain of S8 and the highly conserved nucleotide core, characterized by two dinucleotide platforms, facing each other. The first platform (A595-A596), which is the less phylogenetically and structurally constrained, does not directly contact the protein but has an important shaping role in inducing cross-strand stacking interactions. The second platform (U641-A642) is specifically recognized by the protein. The universally conserved A642 plays a pivotal role by ensuring the cohesion of the complex organization of the core through an array of hydrogen bonds, including the G597-C643-U641 base triple. In addition, A642 provides the unique base-specific interaction with the conserved Ser105, while the Thr106 - Thr107 peptide link is stacked on its purine ring. Noteworthy, the specific recognition of this tripeptide (Thr-Ser-Thr/Ser) is parallel to the recognition of an RNA tetraloop by a dinucleotide platform in the P4-P6 ribozyme domain of group I intron. This suggests a general dual role of dinucleotide platforms in recognition of RNA or peptide motifs. One prominent feature is that conserved side-chain amino acids, as well as conserved bases, are essentially involved in maintaining tertiary folds. The specificity of binding is mainly driven by shape complementarity, which is increased by the hydrophobic part of side-chains. The remarkable similarity of this complex with its homologue in the T. thermophilus 30 S subunit indicates a conserved interaction mode between Archaea and Bacteria.

Non-covalent interactions in the crystallization of the enantiomers of 1,7-dioxaspiro[5.5]undecane (olive fly sex pheromone) by enantiospecific cyclodextrin hosts, hexakis(2,3,6-tri-O-methyl)-α-cyclodextrin and heptakis(2,3,6-tri-O-methyl)-β-cyclodextrin

The enantiomers of racemic olive fly sex pheromone 1,7-dioxaspiro[5.5]undecane (1) have been isolated by crystallization with enantiospecific cyclodextrin hosts: (S)-(1) crystallizes with heptakis(2,3,6-tri-O-methyl)-β-cyclodextrin (TMβCD) and (R)-(1) with hexakis(2,3,6-tri-O-methyl)-α-cyclodextrin (TMαCD). The crystal structure of TMβCD/(S)-(1) from synchrotron radiation data at 100 K, determined for the first time, proves that TMβCD crystallizes with only the (S)-enantiomer from the racemic mixture. Comparison with the 100 K structure of TMαCD/(R)-(1) redetermined with synchrotron data has provided insight into the interactions between each of the hosts with the corresponding enantiomeric guests. Owing to the high resolution of the data and the unusually high quality of the crystals, localization of the H atoms has been achieved, a rare accomplishment for cyclodextrin X-ray structures. This made possible, apart from the geometry of the complexes, the detailed description of a five-membered-ring water cluster with very well ordered hydrogen bonding. The enantiospecificity exhibited by the described systems reveals the subtle differences of the weak intermolecular forces involved in the selective binding of the two optical antipodes by the two hosts. The binding geometry in the two complexes is different, but it is effected in both by numerous host–guest C—H...O interactions, resulting from induced fit of the hosts toward each of the enantiomeric guests. In TMαCD/(R)-(1) two of these H...O host–guest distances, directed toward the acetal O atoms defining the chirality of the guest, are much shorter than the rest and also shorter than all the H...O distances in TMβCD/(S)-(1). Moreover, (R)-(1) interacts not only with the enclosing host, but with other hosts in the crystal lattice, in contrast to (S)-(1) in the TMβCD/(S)-(1) complex which is isolated inside channels formed by the host molecules. The above differences are reflected in the much higher binding constant of TMαCD/(R)-(1) compared with that of TMβCD/(S)-(1) (∼6800 and ∼935 M−1, respectively), determined by NMR in aqueous solution, and the ability of TMαCD to selectively precipitate (R)-(1) from racemic (1) in much higher yield than TMβCD precipitates (S)-(1).

On the mechanism of biological methane formation: structural evidence for conformational changes in methyl-coenzyme M reductase upon substrate binding

Methyl-coenzyme M reductase (MCR) catalyzes the final reaction of the energy conserving pathway of methanogenic archaea in which methylcoenzyme M and coenzyme B are converted to methane and the heterodisulfide CoM-S-S-CoB. It operates under strictly anaerobic conditions and contains the nickel porphinoid F430 which is present in the nickel (I) oxidation state in the active enzyme. The known crystal structures of the inactive nickel (II) enzyme in complex with coenzyme M and coenzyme B (MCR-ox1-silent) and in complex with the heterodisulfide CoM-S-S-CoB (MCR-silent) were now refined at 1.16 A and 1.8 A resolution, respectively. The atomic resolution structure of MCR-ox1-silent describes the exact geometry of the cofactor F430, of the active site residues and of the modified amino acid residues. Moreover, the observation of 18 Mg2+ and 9 Na+ ions at the protein surface of the 300 kDa enzyme specifies typical constituents of binding sites for either ion. The MCR-silent and MCR-ox1-silent structures differed in the occupancy of bound water molecules near the active site indicating that a water chain is involved in the replenishment of the active site with water molecules. The structure of the novel enzyme state MCR-red1-silent at 1.8 A resolution revealed an active site only partially occupied by coenzyme M and coenzyme B. Increased flexibility and distinct alternate conformations were observed near the active site and the substrate channel. The electron density of the MCR-red1-silent state aerobically co-crystallized with coenzyme M displayed a fully occupied coenzyme M-binding site with no alternate conformations. Therefore, the structure was very similar to the MCR-ox1-silent state. As a consequence, the binding of coenzyme M induced specific conformational changes that postulate a molecular mechanism by which the enzyme ensures that methylcoenzyme M enters the substrate channel prior to coenzyme B as required by the active-site geometry. The three different enzymatically inactive enzyme states are discussed with respect to their enzymatically active precursors and with respect to the catalytic mechanism.

Transthyretin stability as a key factor in amyloidogenesis: X-ray analysis at atomic resolution

Transthyretin (TTR) amyloidosis is a conformational disturbance, which, like other amyloidoses, represents a life threat. Here, we report a TTR variant, TTR Thr119Met, that has been shown to have a protective role in the development of clinical symptoms in carriers of TTR Val30Met, one of the most frequent variants among TTR amyloidosis patients. In order to understand this effect, we have determined the structures of the TTR Val30Met/Thr119Met double mutant isolated from the serum of one patient and of both the native and thyroxine complex of TTR Thr119Met. Major conclusions are: (i) new H-bonds within each monomer and monomer-monomer inter-subunit contacts, e.g. Ser117-Ser117 and Met119-Tyr114, increase protein stability, possibly leading to the protective effect of the TTR Val30Met/Thr119Met variant when compared to the single variant TTR Val30Met. (ii) The mutated residue (Met119) extends across the thyroxine binding channel inducing conformational changes that lead to closer contacts between different dimers within the tetramer. The data, at atomic resolution, were essential to detect, for the first time, the subtle changes in the inter-subunit contacts of TTR, and explain the non-amyloidogenic potential of the TTR Thr119Met variant, improving considerably current research on the TTR amyloid fibril formation pathway.

Structural basis for substrate specificity differences of horse liver alcohol dehydrogenase isozymes

A structure determination in combination with a kinetic study of the steroid converting isozyme of horse liver alcohol dehydrogenase, SS-ADH, is presented. Kinetic parameters for the substrates, 5beta-androstane-3beta,17beta-ol, 5beta-androstane-17beta-ol-3-one, ethanol, and various secondary alcohols and the corresponding ketones are compared for the SS- and EE-isozymes which differ by nine amino acid substitutions and one deletion. Differences in substrate specificity and stereoselectivity are explained on the basis of individual kinetic rate constants for the underlying ordered bi-bi mechanism. SS-ADH was crystallized in complex with 3alpha,7alpha,12alpha-trihydroxy-5beta-cholan -24-acid (cholic acid) and NAD(+), but microspectrophotometric analysis of single crystals proved it to be a mixed complex containing 60-70% NAD(+) and 30-40% NADH. The crystals belong to the space group P2(1) with cell dimensions a = 55.0 A, b = 73.2 A, c = 92.5 A, and beta = 102.5 degrees. A 98% complete data set to 1.54-A resolution was collected at 100 K using synchrotron radiation. The structure was solved by the molecular replacement method utilizing EE-ADH as the search model. The major structural difference between the isozymes is a widening of the substrate channel. The largest shifts in C(alpha) carbon positions (about 5 A) are observed in the loop region, in which a deletion of Asp115 is found in the SS isozyme. SS-ADH easily accommodates cholic acid, whereas steroid substrates of similar bulkiness would not fit into the EE-ADH substrate site. In the ternary complex with NAD(+)/NADH, we find that the carboxyl group of cholic acid ligates to the active site zinc ion, which probably contributes to the strong binding in the ternary NAD(+) complex.

Accurate protein crystallography at ultra-high resolution: valence electron distribution in crambin

The charge density distribution of a protein has been refined experimentally. Diffraction data for a crambin crystal were measured to ultra-high resolution (0.54 A) at low temperature by using short-wavelength synchrotron radiation. The crystal structure was refined with a model for charged, nonspherical, multipolar atoms to accurately describe the molecular electron density distribution. The refined parameters agree within 25% with our transferable electron density library derived from accurate single crystal diffraction analyses of several amino acids and small peptides. The resulting electron density maps of redistributed valence electrons (deformation maps) compare quantitatively well with a high-level quantum mechanical calculation performed on a monopeptide. This study provides validation for experimentally derived parameters and a window into charge density analysis of biological macromolecules.

Ab initio solution and refinement of two high-potential iron protein structures at atomic resolution

The crystal structure of the reduced high-potential iron protein (HiPIP) from Chromatium vinosum has been redetermined in a new orthorhombic crystal modification, and the structure of its H42Q mutant has been determined in orthorhombic (H42Q-1) and cubic (H42Q-2) modifications. The first two were solved by ab initio direct methods using data collected to atomic resolution (1.20 and 0. 93 A, respectively). The recombinant wild type (rc-WT) with two HiPIP molecules in the asymmetric unit has 1264 protein atoms and 335 solvent sites, and is the second largest structure reported so far that has been solved by pure direct methods. The solutions were obtained in a fully automated way and included more than 80% of the protein atoms. Restrained anisotropic refinement for rc-WT and H42Q-1 converged to R(1) = summation operator||F(o)| - |F(c)|| / summation operator|F(o)| of 12.0 and 13.6%, respectively [data with I > 2sigma(I)], and 12.8 and 15.5% (all data). H42Q-2 contains two molecules in the asymmetric unit and diffracted only to 2.6 A. In both molecules of rc-WT and in the single unique molecule of H42Q-1 the [Fe(4)S(4)](2+) cluster dimensions are very similar and show a characteristic tetragonal distortion with four short Fe-S bonds along four approximately parallel cube edges, and eight long Fe-S bonds. The unique protein molecules in H42Q-2 and rc-WT are also very similar in other respects, except for the hydrogen bonding around the mutated residue that is at the surface of the protein, supporting the hypothesis that the difference in redox potentials at lower pH values is caused primarily by differences in the charge distribution near the surface of the protein rather than by structural differences in the cluster region.