Twinning and intertwined microcrystals in an intriguing, yet elusive, mineral.

After a hundred years with no solution, the seemingly simple but actually very complex mineral kaliophilite, KAlSiO4, is finally revealed by electron crystallography.

3D Electron Diffraction: The Nanocrystallography Revolution.

Crystallography of nanocrystalline materials has witnessed a true revolution in the past 10 years, thanks to the introduction of protocols for 3D acquisition and analysis of electron diffraction data. This method provides single-crystal data of structure solution and refinement quality, allowing the atomic structure determination of those materials that remained hitherto unknown because of their limited crystallinity. Several experimental protocols exist, which share the common idea of sampling a sequence of diffraction patterns while the crystal is tilted around a noncrystallographic axis, namely, the goniometer axis of the transmission electron microscope sample stage. This Outlook reviews most important 3D electron diffraction applications for different kinds of samples and problematics, related with both materials and life sciences. Structure refinement including dynamical scattering is also briefly discussed.

A Rare Lysozyme Crystal Form Solved Using Highly Redundant Multiple Electron Diffraction Datasets from Micron-Sized Crystals.

Recent developments of novel electron diffraction techniques have shown to be powerful for determination of atomic resolution structures from micron- and nano-sized crystals, too small to be studied by single-crystal X-ray diffraction. In this work, the structure of a rare lysozyme polymorph is solved and refined using continuous rotation MicroED data and standard X-ray crystallographic software. Data collection was performed on a standard 200 kV transmission electron microscope (TEM) using a highly sensitive detector with a short readout time. The data collection is fast (∼3 min per crystal), allowing multiple datasets to be rapidly collected from a large number of crystals. We show that merging data from 33 crystals significantly improves not only the data completeness, overall I/σ and the data redundancy, but also the quality of the final atomic model. This is extremely useful for electron beam-sensitive crystals of low symmetry or with a preferred orientation on the TEM grid.

Influence of mutations at the proximal histidine position on the Fe-O2 bond in hemoglobin from density functional theory.

Four mutated hemoglobin (Hb) variants and wild type hemoglobin as a reference have been investigated using density functional theory methods focusing on oxygen binding. Dispersion-corrected B3LYP functional is used and found to provide reliable oxygen binding energies. It also correctly reproduces the spin distribution of both bound and free heme groups as well as provides correct geometries at their close vicinity. Mutations in hemoglobin are not only an intrigued biological problem and it is also highly important to understand their effects from a clinical point of view. This study clearly shows how even small structural differences close to the heme group can have a significant effect in reducing the oxygen binding of mutated hemoglobins and consequently affecting the health condition of the patient suffering from the mutations. All of the studied mutated Hb variants did exhibit much weaker binding of molecular oxygen compared to the wild type of hemoglobin.

Three-dimensional electron diffraction as a complementary technique to powder X-ray diffraction for phase identification and structure solution of powders.

Phase identification and structure determination are important and widely used techniques in chemistry, physics and materials science. Recently, two methods for automated three-dimensional electron diffraction (ED) data collection, namely automated diffraction tomography (ADT) and rotation electron diffraction (RED), have been developed. Compared with X-ray diffraction (XRD) and two-dimensional zonal ED, three-dimensional ED methods have many advantages in identifying phases and determining unknown structures. Almost complete three-dimensional ED data can be collected using the ADT and RED methods. Since each ED pattern is usually measured off the zone axes by three-dimensional ED methods, dynamic effects are much reduced compared with zonal ED patterns. Data collection is easy and fast, and can start at any arbitrary orientation of the crystal, which facilitates automation. Three-dimensional ED is a powerful technique for structure identification and structure solution from individual nano- or micron-sized particles, while powder X-ray diffraction (PXRD) provides information from all phases present in a sample. ED suffers from dynamic scattering, while PXRD data are kinematic. Three-dimensional ED methods and PXRD are complementary and their combinations are promising for studying multiphase samples and complicated crystal structures. Here, two three-dimensional ED methods, ADT and RED, are described. Examples are given of combinations of three-dimensional ED methods and PXRD for phase identification and structure determination over a large number of different materials, from Ni-Se-O-Cl crystals, zeolites, germanates, metal-organic frameworks and organic compounds to intermetallics with modulated structures. It is shown that three-dimensional ED is now as feasible as X-ray diffraction for phase identification and structure solution, but still needs further development in order to be as accurate as X-ray diffraction. It is expected that three-dimensional ED methods will become crucially important in the near future.

Influence of antifreeze proteins on the ice/water interface.

Antifreeze proteins (AFP) are responsible for the survival of several species, ranging from bacteria to fish, that encounter subzero temperatures in their living environment. AFPs have been divided into two main families, moderately and hyperactive, depending on their thermal hysteresis activity. We have studied one protein from both families, the AFP from the snow flea (sfAFP) and from the winter flounder (wfAFP), which belong to the hyperactive and moderately active family, respectively. On the basis of molecular dynamics simulations, we have estimated the thickness of the water/ice interface for systems both with and without the AFPs attached onto the ice surface. The calculation of the diffusion profiles along the simulation box allowed us to measure the interface width for different ice planes. The obtained widths clearly show a different influence of the two AFPs on the ice/water interface. The different impact of the AFPs here studied on the interface thickness can be related to two AFPs properties: the protein hydrophobic surface and the number of hydrogen bonds that the two AFPs faces form with water molecules.

Induced ice melting by the snow flea antifreeze protein from molecular dynamics simulations.

Antifreeze proteins (AFP) allow different life forms, insects as well as fish and plants, to survive in subzero environments. AFPs prevent freezing of the physiological fluids. We have studied, through molecular dynamics simulations, the behavior of the small isoform of the AFP found in the snow flea (sfAFP), both in water and at the ice/water interface, of four different ice planes. In water at room temperature, the structure of the sfAFP is found to be slightly unstable. The loop between two polyproline II helices has large fluctuations as well as the C-terminus. Torsional angle analyses show a decrease of the polyproline II helix area in the Ramachandran plots. The protein structure instability, in any case, should not affect its antifreeze activity. At the ice/water interface the sfAFP triggers local melting of the ice surface. Bipyramidal, secondary prism, and prism ice planes melt in the presence of AFP at temperatures below the melting point of ice. Only the basal plane is found to be stable at the same temperatures, indicating an adsorption of the sfAFP on this ice plane as confirmed by experimental evidence.

Glucose oxidase from Penicillium amagasakiense: characterization of the transition state of its denaturation from molecular dynamics simulations.

Glucose oxidase (GOx) is a flavoenzyme having applications in food and medical industries. However, GOx, as many other enzymes when extracted from the cells, has relatively short operational lifetimes. Several recent studies (both experimental and theoretical), carried out on small proteins (or small fractions of large proteins), show that a detailed knowledge of how the breakdown process starts and proceeds on molecular level could be of significant help to artificially improve the stability of fragile proteins. We have performed extended molecular dynamics (MD) simulations to study the denaturation of GOx (a protein dimer containing nearly 1200 amino acids) to identify weak points in its structure and in this way gather information to later make it more stable, for example, by mutations. A denaturation of a protein can be simulated by increasing the temperature far above physiological temperature. We have performed a series of MD simulations at different temperatures (300, 400, 500, and 600 K). The exit from the protein's native state has been successfully identified with the clustering method and supported by other methods used to analyze the simulation data. A common set of amino acids is regularly found to initiate the denaturation, suggesting a moiety where the enzyme could be strengthened by a suitable amino acid based modification.

Three-dimensional rotation electron diffraction: software RED for automated data collection and data processing.

Implementation of a computer program package for automated collection and processing of rotation electron diffraction (RED) data is described. The software package contains two computer programs: RED data collection and RED data processing. The RED data collection program controls the transmission electron microscope and the camera. Electron beam tilts at a fine step (0.05-0.20°) are combined with goniometer tilts at a coarse step (2.0-3.0°) around a common tilt axis, which allows a fine relative tilt to be achieved between the electron beam and the crystal in a large tilt range. An electron diffraction (ED) frame is collected at each combination of beam tilt and goniometer tilt. The RED data processing program processes three-dimensional ED data generated by the RED data collection program or by other approaches. It includes shift correction of the ED frames, peak hunting for diffraction spots in individual ED frames and identification of these diffraction spots as reflections in three dimensions. Unit-cell parameters are determined from the positions of reflections in three-dimensional reciprocal space. All reflections are indexed, and finally a list with hkl indices and intensities is output. The data processing program also includes a visualizer to view and analyse three-dimensional reciprocal lattices reconstructed from the ED frames. Details of the implementation are described. Data collection and data processing with the software RED are demonstrated using a calcined zeolite sample, silicalite-1. The structure of the calcined silicalite-1, with 72 unique atoms, could be solved from the RED data by routine direct methods.

Structures of pseudo-decagonal approximants in Al-Co-Ni.

Quasi-crystals shocked the crystallographic world when they were reported in 1984. We now know that they are not a rare exception, and can be found in many alloy systems. One of the richer systems for quasi-crystals and their approximants is Al-Co-Ni. A large series of pseudo-decagonal (PD) approximants have been found. Only two of them, PD4 and PD8, have been solved by X-ray crystallography. We report here the structures of PD1, PD2, PD3 and PD5, solved from the limited information that is provided by electron diffraction patterns, unit cell dimensions and high-resolution electron microscopy images.

Structure projection reconstruction from through-focus series of high-resolution transmission electron microscopy images.

A structure projection reconstruction method based on contrast transfer function correction of through-focus series of high-resolution transmission electron microscopy images is presented. In this method, defocus values are determined by evaluating phase similarities of the pixels on the Fourier transforms of the images after correction using trial defocus values. Two-fold astigmatism is also determined, by measuring focus variation along different directions. Each image in the series is corrected for the effects of contrast transfer function and then combined into a structure projection image. The method works for both crystalline and non-crystalline objects. Test results with experimental images are presented. Influences of experimental parameters for imaging and effects of crystal thickness on reconstruction are discussed.

Precession electron diffraction using a digital sampling method.

A software-based method for collecting precession electron diffraction (PED) patterns is described. The PED patterns are obtained on a computer controlled transmission electron microscope. A series of electron diffraction (ED) patterns are collected as still ED frames at equal intervals, while the electron beam is precessed by one period (360°) around the optical axis. A PED pattern is obtained by combining the different ED frames, which resembles the sampling of a conventional PED pattern. Since intermediate ED frames are collected, it is possible to perform different post-processing strategies on the ED data. This can be used for geometric corrections to obtain accurate integrated intensities. The alignments and data collection are fully automated and controlled by software. The data quality is comparable to what can be achieved using specialized hardware for precession. The PED data can be used for structure solution and refinement with reasonably good R-values.

A complicated quasicrystal approximant epsilon16 predicted by the strong-reflections approach.

The structure of a complicated quasicrystal approximant epsilon(16) was predicted from a known and related quasicrystal approximant epsilon(6) by the strong-reflections approach. Electron-diffraction studies show that in reciprocal space, the positions of the strongest reflections and their intensity distributions are similar for both approximants. By applying the strong-reflections approach, the structure factors of epsilon(16) were deduced from those of the known epsilon(6) structure. Owing to the different space groups of the two structures, a shift of the phase origin had to be applied in order to obtain the phases of epsilon(16). An electron-density map of epsilon(16) was calculated by inverse Fourier transformation of the structure factors of the 256 strongest reflections. Similar to that of epsilon(6), the predicted structure of epsilon(16) contains eight layers in each unit cell, stacked along the b axis. Along the b axis, epsilon(16) is built by banana-shaped tiles and pentagonal tiles; this structure is confirmed by high-resolution transmission electron microscopy (HRTEM). The simulated precession electron-diffraction (PED) patterns from the structure model are in good agreement with the experimental ones. Epsilon(16) with 153 unique atoms in the unit cell is the most complicated approximant structure ever solved or predicted.

A novel method for accurate one-dimensional protein structure prediction based on fragment matching.

MOTIVATION: The precise prediction of one-dimensional (1D) protein structure as represented by the protein secondary structure and 1D string of discrete state of dihedral angles (i.e. Shape Strings) is a prerequisite for the successful prediction of three-dimensional (3D) structure as well as protein-protein interaction. We have developed a novel 1D structure prediction method, called Frag1D, based on a straightforward fragment matching algorithm and demonstrated its success in the prediction of three sets of 1D structural alphabets, i.e. the classical three-state secondary structure, three- and eight-state Shape Strings. RESULTS: By exploiting the vast protein sequence and protein structure data available, we have brought secondary-structure prediction closer to the expected theoretical limit. When tested by a leave-one-out cross validation on a non-redundant set of PDB cutting at 30% sequence identity containing 5860 protein chains, the overall per-residue accuracy for secondary-structure prediction, i.e. Q3 is 82.9%. The overall per-residue accuracy for three- and eight-state Shape Strings are 85.1 and 71.5%, respectively. We have also benchmarked our program with the latest version of PSIPRED for secondary structure prediction and our program predicted 0.3% better in Q3 when tested on 2241 chains with the same training set. For Shape Strings, we compared our method with a recently published method with the same dataset and definition as used by that method. Our program predicted at 2.2% better in accuracy for three-state Shape Strings. By quantitatively investigating the effect of data base size on 1D structure prediction we show that the accuracy increases by approximately 1% with every doubling of the database size.

Describing and comparing protein structures using shape strings.

Different methods for describing and comparing the structures of the tens of thousands of proteins that have been determined by X-ray crystallography are reviewed. Such comparisons are important for understanding the structures and functions of proteins and facilitating structure prediction, as well as assessing structure prediction methods. We summarize methods in this field emphasizing ways of representing protein structures as one-dimensional geometrical strings. Such strings are based on the shape symbols of clustered regions of phi/Psi dihedral angle pairs of the polypeptide backbones as described by the Ramachandran plot. These one-dimensional expressions are as compact as secondary structure description but contain more information in loop regions. They can be used for fast searching for similar structures in databases and for comparing similarities between proteins and between the predicted and native structures.

Salt-bridge dynamics control substrate-induced conformational change in the membrane transporter GlpT.

Active transport of substrates across cytoplasmic membranes is of great physiological, medical and pharmaceutical importance. The glycerol-3-phosphate (G3P) transporter (GlpT) of the E. coli inner membrane is a secondary active antiporter from the ubiquitous major facilitator superfamily that couples the import of G3P to the efflux of inorganic phosphate (P(i)) down its concentration gradient. Integrating information from a novel combination of structural, molecular dynamics simulations and biochemical studies, we identify the residues involved directly in binding of substrate to the inward-facing conformation of GlpT, thus defining the structural basis for the substrate-specificity of this transporter. The substrate binding mechanism involves protonation of a histidine residue at the binding site. Furthermore, our data suggest that the formation and breaking of inter- and intradomain salt bridges control the conformational change of the transporter that accompanies substrate translocation across the membrane. The mechanism we propose may be a paradigm for organophosphate:phosphate antiporters.

Prediction of zinc-binding sites in proteins from sequence.

MOTIVATION: Motivated by the abundance, importance and unique functionality of zinc, both biologically and physiologically, we have developed an improved method for the prediction of zinc-binding sites in proteins from their amino acid sequences. RESULTS: By combining support vector machine (SVM) and homology-based predictions, our method predicts zinc-binding Cys, His, Asp and Glu with 75% precision (86% for Cys and His only) at 50% recall according to a 5-fold cross-validation on a non-redundant set of protein chains from the Protein Data Bank (PDB) (2727 chains, 235 of which bind zinc). Consequently, our method predicts zinc-binding Cys and His with 10% higher precision at different recall levels compared to a recently published method when tested on the same dataset. AVAILABILITY: The program is available for download at www.fos.su.se/~nanjiang/zincpred/download/

Electron crystallography: imaging and single-crystal diffraction from powders.

The study of crystals at atomic level by electrons - electron crystallography - is an important complement to X-ray crystallography. There are two main advantages of structure determinations by electron crystallography compared to X-ray diffraction: (i) crystals millions of times smaller than those needed for X-ray diffraction can be studied and (ii) the phases of the crystallographic structure factors, which are lost in X-ray diffraction, are present in transmission-electron-microscopy (TEM) images. In this paper, some recent developments of electron crystallography and its applications, mainly on inorganic crystals, are shown. Crystal structures can be solved to atomic resolution in two dimensions as well as in three dimensions from both TEM images and electron diffraction. Different techniques developed for electron crystallography, including three-dimensional reconstruction, the electron precession technique and ultrafast electron crystallography, are reviewed. Examples of electron-crystallography applications are given. There is in principle no limitation to the complexity of the structures that can be solved by electron crystallography.

Docking and homology modeling explain inhibition of the human vesicular glutamate transporters.

As membrane transporter proteins, VGLUT1-3 mediate the uptake of glutamate into synaptic vesicles at presynaptic nerve terminals of excitatory neural cells. This function is crucial for exocytosis and the role of glutamate as the major excitatory neurotransmitter in the central nervous system. The three transporters, sharing 76% amino acid sequence identity in humans, are highly homologous but differ in regional expression in the brain. Although little is known regarding their three-dimensional structures, hydropathy analysis on these proteins predicts 12 transmembrane segments connected by loops, a topology similar to other members in the major facilitator superfamily, where VGLUT1-3 have been phylogenetically classified. In this work, we present a three-dimensional model for the human VGLUT1 protein based on its distant bacterial homolog in the same superfamily, the glycerol-3-phosphate transporter from Escherichia coli. This structural model, stable during molecular dynamics simulations in phospholipid bilayers solvated by water, reveals amino acid residues that face its pore and are likely to affect substrate translocation. Docking of VGLUT1 substrates to this pore localizes two different binding sites, to which inhibitors also bind with an overall trend in binding affinity that is in agreement with previously published experimental data.