Translation-libration-screw Research Articles

A jumping crystal predicted with molecular dynamics and analysed with TLS refinement against powder diffraction data.

By running a temperature series of molecular dynamics (MD) simulations starting from the known low-temperature phase, the experimentally observed phase transition in a 'jumping crystal' was captured, thereby providing a prediction of the unknown crystal structure of the high-temperature phase and clarifying the phase-transition mechanism. The phase transition is accompanied by a discontinuity in two of the unit-cell parameters. The structure of the high-temperature phase is very similar to that of the low-temperature phase. The anisotropic displacement parameters calculated from the MD simulations readily identified libration as the driving force behind the phase transition. Both the predicted crystal structure and the phase-transition mechanism were verified experimentally using TLS (translation, libration, screw) refinement against X-ray powder diffraction data.

Doxycycline hydrate and doxycycline hydrochloride dihydrate – crystal structure and charge density analysis

Abstract High-resolution low-temperature X-ray diffraction experiments for doxycycline monohydrate and hydrochloride dihydrate have been performed. Translation-Libration-Screw (TLS) analysis for both crystal forms as well as the data from neutron diffraction experiment for hydrochloride combined with the Hansen-Coppens formalism resulted in precise charge density distribution models for both the zwitterionic monohydrate and a protonated hydrochloride crystal forms. Their detailed topological analysis suggested that the electron structure of doxycycline’s amide moiety undergoes significant changes during protonation due to formation of a very strong resonance-assisted hydrogen bond. A notably increased participation of amide nitrogen atom and hydrogen-accepting oxygen atom in the resonance upon doxycycline protonation was observed. A comparison of TLS- and neutron data-derived hydrogen parameters confirmed the experimental neutron data to be vital for proper description of intra- and inter-molecular interactions in this compound. Finally, calculated lattice and interaction energies quantified repulsive Dox-Dox interactions in the protonated crystal form of the antibiotic, relating with a good solubility of doxycycline hydrochloride relative to its hydrate.

Yes, one can obtain better quality structures from routine X-ray data collection.

Single-crystal X-ray diffraction structural results for benzidine dihydrochloride, hydrated and protonated N,N,N,N-peri(dimethylamino)naphthalene chloride, triptycene, dichlorodimethyltriptycene and decamethylferrocene have been analysed. A critical discussion of the dependence of structural and thermal parameters on resolution for these compounds is presented. Results of refinements against X-ray data, cut off to different resolutions from the high-resolution data files, are compared to structural models derived from neutron diffraction experiments. The Independent Atom Model (IAM) and the Transferable Aspherical Atom Model (TAAM) are tested. The average differences between the X-ray and neutron structural parameters (with the exception of valence angles defined by H atoms) decrease with the increasing 2θmax angle. The scale of differences between X-ray and neutron geometrical parameters can be significantly reduced when data are collected to the higher, than commonly used, 2θmax diffraction angles (for Mo Kα 2θmax > 65°). The final structural and thermal parameters obtained for the studied compounds using TAAM refinement are in better agreement with the neutron values than the IAM results for all resolutions and all compounds. By using TAAM, it is still possible to obtain accurate results even from low-resolution X-ray data. This is particularly important as TAAM is easy to apply and can routinely be used to improve the quality of structural investigations [Dominiak (2015 ▸). LSDB from UBDB. University of Buffalo, USA]. We can recommend that, in order to obtain more adequate (more accurate and precise) structural and displacement parameters during the IAM model refinement, data should be collected up to the larger diffraction angles, at least, for Mo Kα radiation to 2θmax = 65° (sin θmax/λ < 0.75 Å(-1)). The TAAM approach is a very good option to obtain more adequate results even using data collected to the lower 2θmax angles. Also the results of translation-libration-screw (TLS) analysis and vibrational entropy values are more reliable for 2θmax > 65°.

Estimating temperature-dependent anisotropic hydrogen displacements with the invariom database and a new segmented rigid-body analysis program

Invariom partitioning and notation are used to estimate anisotropic hydrogen displacements for incorporation in crystallographic refinement models. Optimized structures of the generalized invariom database and their frequency computations provide the information required: frequencies are converted to internal atomic displacements and combined with the results of a TLS (translation-libration-screw) fit of experimental non-hydrogen anisotropic displacement parameters to estimate those of H atoms. Comparison with TLS+ONIOM and neutron diffraction results for four example structures where high-resolution X-ray and neutron data are available show that electron density transferability rules established in the invariom approach are also suitable for streamlining the transfer of atomic vibrations. A new segmented-body TLS analysis program called APD-Toolkit has been coded to overcome technical limitations of the established program THMA. The influence of incorporating hydrogen anisotropic displacement parameters on conventional refinement is assessed.

Predicting X-ray diffuse scattering from translation-libration-screw structural ensembles.

Identifying the intramolecular motions of proteins and nucleic acids is a major challenge in macromolecular X-ray crystallography. Because Bragg diffraction describes the average positional distribution of crystalline atoms with imperfect precision, the resulting electron density can be compatible with multiple models of motion. Diffuse X-ray scattering can reduce this degeneracy by reporting on correlated atomic displacements. Although recent technological advances are increasing the potential to accurately measure diffuse scattering, computational modeling and validation tools are still needed to quantify the agreement between experimental data and different parameterizations of crystalline disorder. A new tool, phenix.diffuse, addresses this need by employing Guinier's equation to calculate diffuse scattering from Protein Data Bank (PDB)-formatted structural ensembles. As an example case, phenix.diffuse is applied to translation-libration-screw (TLS) refinement, which models rigid-body displacement for segments of the macromolecule. To enable the calculation of diffuse scattering from TLS-refined structures, phenix.tls_as_xyz builds multi-model PDB files that sample the underlying T, L and S tensors. In the glycerophosphodiesterase GpdQ, alternative TLS-group partitioning and different motional correlations between groups yield markedly dissimilar diffuse scattering maps with distinct implications for molecular mechanism and allostery. These methods demonstrate how, in principle, X-ray diffuse scattering could extend macromolecular structural refinement, validation and analysis.

From deep TLS validation to ensembles of atomic models built from elemental motions.

The translation-libration-screw model first introduced by Cruickshank, Schomaker and Trueblood describes the concerted motions of atomic groups. Using TLS models can improve the agreement between calculated and experimental diffraction data. Because the T, L and S matrices describe a combination of atomic vibrations and librations, TLS models can also potentially shed light on molecular mechanisms involving correlated motions. However, this use of TLS models in mechanistic studies is hampered by the difficulties in translating the results of refinement into molecular movement or a structural ensemble. To convert the matrices into a constituent molecular movement, the matrix elements must satisfy several conditions. Refining the T, L and S matrix elements as independent parameters without taking these conditions into account may result in matrices that do not represent concerted molecular movements. Here, a mathematical framework and the computational tools to analyze TLS matrices, resulting in either explicit decomposition into descriptions of the underlying motions or a report of broken conditions, are described. The description of valid underlying motions can then be output as a structural ensemble. All methods are implemented as part of the PHENIX project.

SHADE3server: a streamlined approach to estimate H-atom anisotropic displacement parameters using periodicab initiocalculations or experimental information

A major update of theSHADEserver (http://shade.ki.ku.dk) is presented. In addition to all of the previous options for estimating H-atom anisotropic displacement parameters (ADPs) that were offered bySHADE2, the newest version offers two new methods. The first method combines the original translation–libration–screw analysis with input from periodicab initiocalculations. The second method allows the user to input experimental information from spectroscopic measurements or from neutron diffraction experiments on related structures and utilize this information to evaluate ADPs of H atoms. Tools are provided to set up theab initiocalculations and to derive the internal motion from the calculations. The new server was tested on a range of compounds where neutron diffraction data were available. In most cases, the results are significantly better than previous estimates, and for strong hydrogen bonds in proton sponges, theab initiocalculations become crucial.

B-factor Analysis and Conformational Rearrangement of Aldose Reductase.

The NADPH-dependent reduction of glucose reaction that is catalyzed by Aldose Reductase (AR) follows a sequential ordered kinetic mechanism in which the co-factor NADPH binds to the enzyme prior to the aldehyde substrate. The kinetic/structural experiments have found a conformational change involving a hinge-like movement of a surface loop (residues 213-224) which is anticipated to take place upon the binding of the diphosphate moiety of NADPH. The reorientation of this loop, expected to permit the release of NADP+, represents the rate-limiting step of the catalytic mechanism. This study reveals: 1) The Translation/Libration/Screw (TLS) analysis of absolute B-factors of apo AR crystal structures indicates that the 212-224 loop might move as a rigid group. 2) Residues that make the flexible loop slide in the AR binary and ternary complexes. 3) The normalized B-factors separate this segment into three different clusters with fewer residues.

Structural Basis for Universal Corrinoid Recognition by the Cobalamin Transport Protein Haptocorrin

Cobalamin (Cbl; vitamin B12) is an essential micronutrient synthesized only by bacteria. Mammals have developed a sophisticated uptake system to capture the vitamin from the diet. Cbl transport is mediated by three transport proteins: transcobalamin, intrinsic factor, and haptocorrin (HC). All three proteins have a similar overall structure but a different selectivity for corrinoids. Here, we present the crystal structures of human HC in complex with cyanocobalamin and cobinamide at 2.35 and 3.0 Å resolution, respectively. The structures reveal that many of the interactions with the corrin ring are conserved among the human Cbl transporters. However, the non-conserved residues Asn-120, Arg-357, and Asn-373 form distinct interactions allowing for stabilization of corrinoids other than Cbl. A central binding motif forms interactions with the e- and f-side chains of the corrin ring and is conserved in corrinoid-binding proteins of other species. In addition, the α- and β-domains of HC form several unique interdomain contacts and have a higher shape complementarity than those of intrinsic factor and transcobalamin. The stabilization of ligands by all of these interactions is reflected in higher melting temperatures of the protein-ligand complexes. Our structural analysis offers fundamental insights into the unique binding behavior of HC and completes the picture of Cbl interaction with its three transport proteins.

On the use of crystal vibrational modes in the estimation of the anisotropic displacement parameters of hydrogen atoms in molecular crystals: para-Nitroaniline as a test case

The use of crystal vibrational modes and frequencies calculated for the periodic system to complement a Translation Libration Screw (TLS) fit in the estimation of anisotropic displacement parameters (ADPs) of hydrogen atoms in molecular crystals is proposed. As a test case we have used the para-nitroaniline crystal, for which a reference set of ADPs has been obtained by performing a single crystal neutron diffraction study at 100K. Although the largest difference between estimated and experimental reference values amounts to 0.06Å2, this value being about six times the experimental uncertainty, the agreement is better than three experimental uncertainties for 33 out of the total of 36 ADPs. The performance of the suggested method, particularly for the amino atoms, is thoroughly analyzed.

Modelling dynamics in protein crystal structures by ensemble refinement

Single-structure models derived from X-ray data do not adequately account for the inherent, functionally important dynamics of protein molecules. We generated ensembles of structures by time-averaged refinement, where local molecular vibrations were sampled by molecular-dynamics (MD) simulation whilst global disorder was partitioned into an underlying overall translation-libration-screw (TLS) model. Modeling of 20 protein datasets at 1.1-3.1 Å resolution reduced cross-validated R(free) values by 0.3-4.9%, indicating that ensemble models fit the X-ray data better than single structures. The ensembles revealed that, while most proteins display a well-ordered core, some proteins exhibit a 'molten core' likely supporting functionally important dynamics in ligand binding, enzyme activity and protomer assembly. Order-disorder changes in HIV protease indicate a mechanism of entropy compensation for ordering the catalytic residues upon ligand binding by disordering specific core residues. Thus, ensemble refinement extracts dynamical details from the X-ray data that allow a more comprehensive understanding of structure-dynamics-function relationships.DOI:http://dx.doi.org/10.7554/eLife.00311.001.

Crystal Structure of C5b-6 Suggests Structural Basis for Priming Assembly of the Membrane Attack Complex

The complement membrane attack complex (MAC) forms transmembrane pores in pathogen membranes. The first step in MAC assembly is cleavage of C5 to generate metastable C5b, which forms a stable complex with C6, termed C5b-6. C5b-6 initiates pore formation via the sequential recruitment of homologous proteins: C7, C8, and 12-18 copies of C9, each of which comprises a central MAC-perforin domain flanked by auxiliary domains. We recently proposed a model of pore assembly, in which the auxiliary domains play key roles, both in stabilizing the closed conformation of the protomers and in driving the sequential opening of the MAC-perforin β-sheet of each new recruit to the growing pore. Here, we describe an atomic model of C5b-6 at 4.2 Å resolution. We show that C5b provides four interfaces for the auxiliary domains of C6. The largest interface is created by the insertion of an interdomain linker from C6 into a hydrophobic groove created by a major reorganization of the α-helical domain of C5b. In combination with the rigid body docking of N-terminal elements of both proteins, C5b becomes locked into a stable conformation. Both C6 auxiliary domains flanking the linker pack tightly against C5b. The net effect is to induce the clockwise rigid body rotation of four auxiliary domains, as well as the opening/twisting of the central β-sheet of C6, in the directions predicted by our model to activate or prime C6 for the subsequent steps in MAC assembly. The complex also suggests novel small molecule strategies for modulating pathological MAC assembly.

To B or not to B: a question of resolution?

In choosing and refining any crystallographic structural model, there is tension between the desire to extract the most detailed information possible and the necessity to describe no more than what is justified by the observed data. A more complex model is not necessarily a better model. Thus, it is important to validate the choice of parameters as well as validating their refined values. One recurring task is to choose the best model for describing the displacement of each atom about its mean position. At atomic resolution one has the option of devoting six model parameters (a 'thermal ellipsoid') to describe the displacement of each atom. At medium resolution one typically devotes at most one model parameter per atom to describe the same thing (a 'B factor'). At very low resolution one cannot justify the use of even one parameter per atom. Furthermore, this aspect of the structure may be described better by an explicit model of bulk displacements, the most common of which is the translation/libration/screw (TLS) formalism, rather than by assigning some number of parameters to each atom individually. One can sidestep this choice between atomic displacement parameters and TLS descriptions by including both treatments in the same model, but this is not always statistically justifiable. The choice of which treatment is best for a particular structure refinement at a particular resolution can be guided by general considerations of the ratio of model parameters to the number of observations and by specific statistics such as the Hamilton R-factor ratio test.

Structural and Functional Relationships between the Lectin and Arm Domains of Calreticulin

Calreticulin and calnexin are key components in maintaining the quality control of glycoprotein folding within the endoplasmic reticulum. Although their lectin function of binding monoglucosylated sugar moieties of glycoproteins is well documented, their chaperone activity in suppressing protein aggregation is less well understood. Here, we use a series of deletion mutants of calreticulin to demonstrate that its aggregation suppression function resides primarily within its lectin domain. Using hydrophobic peptides as substrate mimetics, we show that aggregation suppression is mediated through a single polypeptide binding site that exhibits a K(d) for peptides of 0.5-1 μM. This site is distinct from the oligosaccharide binding site and differs from previously identified sites of binding to thrombospondin and GABARAP (4-aminobutyrate type A receptor-associated protein). Although the arm domain of calreticulin was incapable of suppressing aggregation or binding hydrophobic peptides on its own, it did contribute to aggregation suppression in the context of the whole molecule. The high resolution x-ray crystal structure of calreticulin with a partially truncated arm domain reveals a marked difference in the relative orientations of the arm and lectin domains when compared with calnexin. Furthermore, a hydrophobic patch was detected on the arm domain that mediates crystal packing and may contribute to calreticulin chaperone function.

The Structural Basis of Endocannabinoid Oxygenation by Cyclooxygenase-2

The cyclooxygenases (COX-1 and COX-2) oxygenate arachidonic acid (AA) in the committed step of prostaglandin biogenesis. Substitutions of I434V, H513R, and I523V constitute the only differences in residues lining the cyclooxygenase channel between COX-1 and COX-2. These changes create a hydrophobic pocket in COX-2, with Arg-513 located at the base of the pocket, which has been exploited in the design of COX-2-selective inhibitors. Previous studies have shown that COX-2, but not COX-1, can oxygenate endocannabinoid substrates, including 2-arachidonoyl glycerol (2-AG). To investigate the isoform-specific structural basis of endocannabinoid binding to COX-2, we determined the crystal structure of the 2-AG isomer 1-arachidonoyl glycerol (1-AG) in complex with wild type and R513H murine (mu) COX-2 to 2.2 and 2.35 Å, respectively, and R513H muCOX-2 in complex with AA to 2.45 Å resolution. The 2,3-dihydroxypropyl moiety of 1-AG binds near the opening of the cyclooxygenase channel in the space vacated by the movement of the Leu-531 side chain, validating our previous hypothesis implicating the flexibility of the Leu-531 side chain as a determinant for the ability of COX-2 to oxygenate endocannabinoid substrates. Functional analyses carried out to compliment our structural findings indicated that Y355F and R513H muCOX-2 constructs had no effect on the oxygenation of 1-AG and 2-AG, whereas substitutions that resulted in a shortened side chain for Leu-531 had only modest effects. Both AA and 1-AG bind to R513H muCOX-2 in conformations similar to those observed in the co-crystal structures of these substrates with wild type enzyme.

Processivity and Subcellular Localization of Glycogen Synthase Depend on a Non-catalytic High Affinity Glycogen-binding Site

Glycogen synthase, a central enzyme in glucose metabolism, catalyzes the successive addition of α-1,4-linked glucose residues to the non-reducing end of a growing glycogen molecule. A non-catalytic glycogen-binding site, identified by x-ray crystallography on the surface of the glycogen synthase from the archaeon Pyrococcus abyssi, has been found to be functionally conserved in the eukaryotic enzymes. The disruption of this binding site in both the archaeal and the human muscle glycogen synthases has a large impact when glycogen is the acceptor substrate. Instead, the catalytic efficiency remains essentially unchanged when small oligosaccharides are used as substrates. Mutants of the human muscle enzyme with reduced affinity for glycogen also show an altered intracellular distribution and a marked decrease in their capacity to drive glycogen accumulation in vivo. The presence of a high affinity glycogen-binding site away from the active center explains not only the long-recognized strong binding of glycogen synthase to glycogen but also the processivity and the intracellular localization of the enzyme. These observations demonstrate that the glycogen-binding site is a critical regulatory element responsible for the in vivo catalytic efficiency of GS.

Structure and function of a HECT domain ubiquitin‐binding site

The Rsp5 ubiquitin ligase contains a non-covalent binding site for ubiquitin within the amino-terminal lobe (N-lobe) of the HECT domain, and the X-ray crystal structure of the HECT-ubiquitin complex has been determined. Hydrophobic patch residues of ubiquitin (L8, I44, V70) were crucial for interaction with Rsp5, and amino-acid alterations at the Rsp5-binding interface resulted in defects in polyubiquitination. Our results support a model in which the N-lobe-binding site acts to localize and orient the distal end of the ubiquitin chain to promote conjugation of the next ubiquitin molecule.

Validation of crystallographic models containing TLS or other descriptions of anisotropy.

The use of TLS (translation/libration/screw) models to describe anisotropic displacement of atoms within a protein crystal structure has become increasingly common. These models may be used purely as an improved methodology for crystallographic refinement or as the basis for analyzing inter-domain and other large-scale motions implied by the crystal structure. In either case it is desirable to validate that the crystallographic model, including the TLS description of anisotropy, conforms to our best understanding of protein structures and their modes of flexibility. A set of validation tests has been implemented that can be integrated into ongoing crystallographic refinement or run afterwards to evaluate a previously refined structure. In either case validation can serve to increase confidence that the model is correct, to highlight aspects of the model that may be improved or to strengthen the evidence supporting specific modes of flexibility inferred from the refined TLS model. Automated validation checks have been added to the PARVATI and TLSMD web servers and incorporated into the CCP4i user interface.

Peroxisomal Plant 3-Ketoacyl-CoA Thiolase Structure and Activity Are Regulated by a Sensitive Redox Switch

The breakdown of fatty acids, performed by the beta-oxidation cycle, is crucial for plant germination and sustainability. beta-Oxidation involves four enzymatic reactions. The final step, in which a two-carbon unit is cleaved from the fatty acid, is performed by a 3-ketoacyl-CoA thiolase (KAT). The shortened fatty acid may then pass through the cycle again (until reaching acetoacetyl-CoA) or be directed to a different cellular function. Crystal structures of KAT from Arabidopsis thaliana and Helianthus annuus have been solved to 1.5 and 1.8 A resolution, respectively. Their dimeric structures are very similar and exhibit a typical thiolase-like fold; dimer formation and active site conformation appear in an open, active, reduced state. Using an interdisciplinary approach, we confirmed the potential of plant KATs to be regulated by the redox environment in the peroxisome within a physiological range. In addition, co-immunoprecipitation studies suggest an interaction between KAT and the multifunctional protein that is responsible for the preceding two steps in beta-oxidation, which would allow a route for substrate channeling. We suggest a model for this complex based on the bacterial system.

Comparative Structural Biology of Eubacterial and Archaeal Oligosaccharyltransferases

Oligosaccharyltransferase (OST) catalyzes the transfer of an oligosaccharide from a lipid donor to an asparagine residue in nascent polypeptide chains. In the bacterium Campylobacter jejuni, a single-subunit membrane protein, PglB, catalyzes N-glycosylation. We report the 2.8 A resolution crystal structure of the C-terminal globular domain of PglB and its comparison with the previously determined structure from the archaeon Pyrococcus AglB. The two distantly related oligosaccharyltransferases share unexpected structural similarity beyond that expected from the sequence comparison. The common architecture of the putative catalytic sites revealed a new catalytic motif in PglB. Site-directed mutagenesis analyses confirmed the contribution of this motif to the catalytic function. Bacterial PglB and archaeal AglB constitute a protein family of the catalytic subunit of OST along with STT3 from eukaryotes. A structure-aided multiple sequence alignment of the STT3/PglB/AglB protein family revealed three types of OST catalytic centers. This novel classification will provide a useful framework for understanding the enzymatic properties of the OST enzymes from Eukarya, Archaea, and Bacteria.