Single particle cryoelectron tomography characterization of the structure and structural variability of poliovirus-receptor-membrane complex at 30 A resolution

As a long-term goal we want to use cryoelectron tomography to understand how non-enveloped viruses, such as picornaviruses, enter cells and translocate their genomes across membranes. To this end, we developed new image-processing tools using an in vitro system to model viral interactions with membranes. The complex of poliovirus with its membrane-bound receptors was reconstructed by averaging multiple sub-tomograms, thereby producing three-dimensional maps of surprisingly high-resolution (30 A). Recognizable images of the complex could be produced by averaging as few as 20 copies. Additionally, model-free reconstructions of free poliovirus particles, clearly showing the major surface features, could be calculated from 60 virions. All calculations were designed to avoid artifacts caused by missing information typical for tomographic data ("missing wedge"). To investigate structural and conformational variability we applied a principal component analysis classification to specific regions. We show that the missing wedge causes a bias in classification, and that this bias can be minimized by supplementation with data from the Fourier transform of the averaged structure. After classifying images of the receptor into groups with high similarity, we were able to see differences in receptor density consistent with the known variability in receptor glycosylation.     

Prediction of homologous protein structures based on conformational searches and energetics

A "knowledge-based" method of predicting the unknown structure of a protein from a homologous known structure using energetics to determine a sidechain conformation is proposed. The method consists of exchanging the residues in the known structure for the sequence of the unknown protein. Then a conformational search with molecular mechanics energy minimization is done on the exchanged residues. The lowest energy conformer is the one picked to be the predicted structure. In the structure of bovine trypsin, the importance of including a solvation energy term in the search is demonstrated for solvent accessible residues, while molecular mechanics alone is enough to correctly predict the conformation of internal residues. The correctness of the model is assessed by a volume error overlap of the predicted structure compared to the crystal structure. Finally, the structure of rat trypsin is predicted from the crystal structure of bovine trypsin. The sequences of these two proteins are 74% identical and all of the significant changes between them are on external residues. Thus, the inclusion of solvation energy in the conformational search is necessary to accurately predict the structure of the exchanged residues.     

The structure of tumor necrosis factor-alpha at 2.6 A resolution. Implications for receptor binding

The three-dimensional structure of tumor necrosis factor (TNF-alpha), a protein hormone secreted by macrophages, has been determined at 2.6 A resolution by x-ray crystallography. Phases were determined by multiple isomorphous replacement using data collected from five heavy atom derivatives. The multiple isomorphous replacement phases were further improved by real space symmetry averaging, exploiting the noncrystallographic 3-fold symmetry of the TNF-alpha trimer. An atomic model corresponding to the known amino acid sequence of TNF-alpha was readily built into the electron density map calculated with these improved phases. The 17,350-dalton monomer forms an elongated, antiparallel beta-pleated sheet sandwich with a "jelly-roll" topology. Three monomers associate intimately about a 3-fold axis of symmetry to form a compact bell-shaped trimer. Examination of the model and comparison to known protein structures reveals striking structural homology to several viral coat proteins, particularly satellite tobacco necrosis virus. Locations of residues conserved between TNF-alpha and lymphotoxin (TNF-beta, a related cytokine known to bind to the same receptors as TNF-alpha) suggest that lymphotoxin, like TNF-alpha, binds to the receptor as a trimer and that the general site of interaction with the receptor is at the "base" of the trimer.     

How the value of the environment controls persistence in visual search

Classic foraging theory predicts that humans and animals aim to gain maximum reward per unit time. However, in standard instrumental conditioning tasks individuals adopt an apparently suboptimal strategy: they respond slowly when the expected value is low. This reward-related bias is often explained as reduced motivation in response to low rewards. Here we present evidence this behavior is associated with a complementary increased motivation to search the environment for alternatives. We trained monkeys to search for reward-related visual targets in environments with different values. We found that the reward-related bias scaled with environment value, was consistent with persistent searching after the target was already found, and was associated with increased exploratory gaze to objects in the environment. A novel computational model of foraging suggests that this search strategy could be adaptive in naturalistic settings where both environments and the objects within them provide partial information about hidden, uncertain rewards.     

Direct determination of crystallographic phases for diffraction data from phospholipid multilamellar arrays.

Direct determination of crystallographic phases based on probabilistic of sigma 1 and sigma 2 "triplet" structure invariants has been found to be an effective technique for structure analysis with lamellar x-ray or electron diffraction intensity data from phospholipids. In many cases, nearly all phase values are determined, permitting a structure density (electron density for x-ray diffraction; electrostatic potential for electron diffraction) map to be calculated, which is directly interpretable in terms of known bilayer lipid structure. The major source of error is found to be due to the distortion of observed electron diffraction intensity data by incoherent multiple scattering, which can significantly affect the appearance of the electrostatic potential map, but not the success of the phase determination, as long as the observed Patterson function can be interpreted.     

Combined use of freeze-fracture electron microscopy and X-ray diffraction for the structure determination of three-dimensionally ordered specimens

The cubic phases of lipid-water systems have been studied by freeze-fracture electron microscopy. The preservation of the sample structure following cryofixation was verified by low temperature X-ray diffraction. Different types of fracture planes were identified; all display highly ordered two-dimensional domains, each subdivided into sub-domains related to each other by displacements and rotations related to the symmetry of the space group. The images were filtered using cross-correlation averaging techniques and the filtered images were compared to the corresponding planar sections of the electron density maps. Several conclusions were drawn: 1) when properly cryofixed, as assessed by low temperature X-ray diffraction, the structure of the sample was well preserved in the replicas; 2) the symmetry of the space group was faithfully reflected in the electron microscope images; 3) the crystallographic orientations of the most frequently identified'fracture planes coincided with those of the most intense X-ray reflections indicating that the fracture propagates, preferentially, in regions where the electron density variations are the largest; 4) when different structural models are compatible with X-ray diffraction data, it is possible to determine the correct model by comparing the filtered images with sections of the corresponding electron density maps; and 5) this approach constitutes a new and powerful tool of general interest for the low resolution study of three-dimensionally ordered specimens.     

Mischaracterising density dependence biases estimated effects of coloured covariates on population dynamics

Environmental effects on population growth are often quantified by coupling environmental covariates with population time series, using statistical models that make particular assumptions about the shape of density dependence. We hypothesized that faulty assumptions about the shape of density dependence can bias estimated effect sizes of temporally autocorrelated covariates. We investigated the presence of bias using Monte Carlo simulations based on three common per capita growth functions with distinct density dependent forms (θ-Ricker, Ricker and Gompertz), autocorrelated (coloured) ‘known’ environmental covariates and uncorrelated (white) ‘unknown’ noise. Faulty assumptions about the shape of density dependence, combined with overcompensatory intrinsic population dynamics, can lead to strongly biased estimated effects of coloured covariates, associated with lower confidence interval coverage. Effects of negatively autocorrelated (blue) environmental covariates are overestimated, while those of positively autocorrelated (red) covariates can be underestimated, generally to a lesser extent. Prewhitening the focal environmental covariate effectively reduces the bias, at the expense of the estimate precision. Fitting models with flexible shapes of density dependence can also reduce bias, but increases model complexity and potentially introduces other problems of parameter identifiability. Model selection is a good option if an appropriate model is included in the set of candidate models. Under the specific and identifiable circumstances with high risk of bias, we recommend prewhitening or careful modelling of the shape of density dependence.     

X-ray structure of lipoamide dehydrogenase from Azotobacter vinelandii determined by a combination of molecular and isomorphous replacement techniques

The crystal structure of lipoamide dehydrogenase from Azotobacter vinelandii has been determined by a combination of molecular replacement and isomorphous replacement techniques yielding eventually a good-quality 2.8 A electron density map. Initially, the structure determination was attempted by molecular replacement procedures alone using a model of human glutathione reductase, which has 26% sequence identity with this bacterial dehydrogenase. The rotation function yielded the correct orientation of the model structure both when the glutathione reductase dimer and monomer were used as starting model. The translation function could not be solved, however. Consequently, data for two heavy-atom derivatives were collected using the Hamburg synchotron facilities. The derivatives had several sites in common, which was presumably a major reason why the electron density map obtained by isomorphous information alone was of poor quality. Application of solvent flattening procedures cleaned up the map considerably, however, showing clearly the outline of the lipoamide dehydrogenase dimer, which has a molecular weight of 100,000. Application of the "phased translation function", which combines the phase information of both isomorphous and molecular replacement, led to an unambiguous determination of the position of the model structure in the lipoamide dehydrogenase unit cell. The non-crystallographic 2-fold axis of the dimer was optimized by several cycles of constrained-restrained least-squares refinement and subsequently used for phase improvement by 2-fold density averaging. After ten cycles at 3.5 A, the resolution was gradually extended to 2.8 A in another 140 cycles. The 2.8 A electron density distribution obtained in this manner was of much improved quality and allowed building of an atomic model of A. vinelandii lipoamide dehydrogenase. It appears that in the orthorhombic crystals used each dimer is involved in contacts with eight surrounding dimers, leaving unexplained why the crystals are rather fragile. Contacts between subunits within one dimer, which are quite extensive, can be divided into two regions separated by a cavity. In one of the contact regions, the level of sequence identity with glutathione reductase is very low but it is quite high in the other. The folding of the polypeptide chain in each subunit is quite similar to that of glutathione reductase, as is the extended conformation of the co-enzyme FAD.(ABSTRACT TRUNCATED AT 400 WORDS)     

Artificial protein sequences enable recognition of vicinal and distant protein functional relationships

High divergence in protein sequences makes the detection of distant protein relationships through homology-based approaches challenging. Grouping protein sequences into families, through similarities in either sequence or 3-D structure, facilitates in the improved recognition of protein relationships. In addition, strategically designed protein-like sequences have been shown to bridge distant structural domain families by serving as artificial linkers. In this study, we have augmented a search database of known protein domain families with such designed sequences, with the intention of providing functional clues to domain families of unknown structure. When assessed using representative query sequences from each family, we obtain a success rate of 94% in protein domain families of known structure. Further, we demonstrate that the augmented search space enabled fold recognition for 582 families with no structural information available a priori. Additionally, we were able to provide reliable functional relationships for 610 orphan families. We discuss the application of our method in predicting functional roles through select examples for DUF4922, DUF5131, and DUF5085. Our approach also detects new associations between families that were previously not known to be related, as demonstrated through new sub-groups of the RNA polymerase domain among three distinct RNA viruses. Taken together, designed sequences-augmented search databases direct the detection of meaningful relationships between distant protein families. In turn, they enable fold recognition and offer reliable pointers to potential functional sites that may be probed further through direct mutagenesis studies.     

Phase accuracy in high-resolution electron microscopy of trigonal and orthorhombic purple membrane

High-resolution images of orthorhombic purple membrane have been obtained by electron cryomicroscopy with spot-scan illumination, and the projection structure at 3.9 Å resolution calculated after image processing and averaging of the data. Since the phases of the structure factors in the projection down the orthorhombic twofold axis should be either 0 or 180°, this offers the first opportunity to make an independent test of the estimated accuracy of high-resolution phases obtained by electron microscopy. The results show the final phases are less accurate than previously estimated by a small factor (1.3). Careful comparison of the new orthorhombic structure to the known trigonal structure shows only small differences after account is taken of a slight difference in the tilt angle of the molecules in the two crystals. This is consistent with the available kinetic and spectroscopic data which show very small differences in behavior.     

Structure determination of feline panleukopenia virus empty particles

Various crystal forms of the single-stranded DNA, feline panleukopenia virus (FPV), a parvovirus, have been grown of both full virions and empty particles. The structure of empty particles crystallized in an orthorhombic space group P2₁2₁2₁, with unit cell dimensions a = 380.1 Å, b = 379.3 Å, and c = 350.9 Å, has been determined to 3.3 Å resolution. The data were collected using oscillation photography with synchrotron radiation. The orientations of the empty capsids in the unit cell were determined using a self-rotation function and their positions were obtained with an R-factor search using canine parvovirus (CPV) as a model. Phases were then calculated, based on the CPV model, to 6.0 Å resolution and gradually extended to 3.3 Å resolution by molecular replacement electron density averaging. The resultant electron density was readily interpreted in terms of the known amino acid sequence. The structure is contrasted to that of CPV in terms of host range, neutralization by antibodies, hemagglutination properties, and binding of genomic DNA. © Wiley-Liss, Inc.     

Knowledge-based real-space explorations for low-resolution structure determination

The accurate and effective interpretation of low-resolution data in X-ray crystallography is becoming increasingly important as structural initiatives turn toward large multiprotein complexes. Substantial challenges remain due to the poor information content and ambiguity in the interpretation of electron density maps at low resolution. Here, we describe a semiautomated procedure that employs a restraint-based conformational search algorithm, RAPPER, to produce a starting model for the structure determination of ligase interacting factor 1 in complex with a fragment of DNA ligase IV at low resolution. The combined use of experimental data and a priori knowledge of protein structure enabled us not only to generate an all-atom model but also to reaffirm the inferred sequence registry. This approach provides a means to extract quickly from experimental data useful information that would otherwise be discarded and to take into account the uncertainty in the interpretation--an overriding issue for low-resolution data.     

Systematic identification of gene-altering programmed inversions across the bacterial domain

Programmed chromosomal inversions allow bacteria to generate intra-population genotypic and functional heterogeneity, a bet-hedging strategy important in changing environments. Some programmed inversions modify coding sequences, producing different alleles in several gene families, most notably in specificity-determining genes such as Type I restriction-modification systems, where systematic searches revealed cross phylum abundance. Yet, a broad, gene-independent, systematic search for gene-altering programmed inversions has been absent, and little is known about their genomic sequence attributes and prevalence across gene families. Here, identifying intra-species variation in genomes of over 35 000 species, we develop a predictive model of gene-altering inversions, revealing key attributes of their genomic sequence attributes, including gene-pseudogene size asymmetry and orientation bias. The model predicted over 11,000 gene-altering loci covering known targeted gene families, as well as novel targeted families including Type II restriction-modification systems, a protein of unknown function, and a fusion-protein containing conjugative-pilus and phage tail domains. Publicly available long-read sequencing datasets validated representatives of these newly predicted inversion-targeted gene families, confirming intra-population genetic heterogeneity. Together, these results reveal gene-altering programmed inversions as a key strategy adopted across the bacterial domain, and highlight programmed inversions that modify Type II restriction-modification systems as a possible new mechanism for maintaining intra-population heterogeneity.     

The power of two-dimensional dwell-time analysis for model discrimination, temporal resolution, multichannel analysis and level detection

Two-dimensional (2D) dwell-time analysis of time series of single-channel patch-clamp current was improved by employing a Hinkley detector for jump detection, introducing a genetic fit algorithm, replacing maximum likelihood by a least square criterion, averaging over a field of 9 or 25 bins in the 2D plane and normalizing per measuring time, not per events. Using simulated time series for the generation of the "theoretical" 2D histograms from assumed Markov models enabled the incorporation of the measured filter response and noise. The effects of these improvements were tested with respect to the temporal resolution, accuracy of the determination of the rate constants of the Markov model, sensitivity to noise and requirement of open time and length of the time series. The 2D fit was better than the classical hidden Markov model (HMM) fit in all tested fields. The temporal resolution of the two most efficient algorithms, the 2D fit and the subsequent HMM/beta fit, enabled the determination of rate constants 10 times faster than the corner frequency of the low-pass filter. The 2D fit was much less sensitive to noise. The requirement of computing time is a problem of the 2D fit (100 times that of the HMM fit) but can now be handled by personal computers. The studies revealed a fringe benefit of 2D analysis: it can reveal the "true" single-channel current when the filter has reduced the apparent current level by averaging over undetected fast gating.     

A low-resolution study of testicular lactate dehydrogenase using the molecular replacement technique

The orientation of the molecular 2-fold axes of mouse testicular lactate dehydrogenase (LDHase³-C₄) was determined by a rotation function search. These were subsequently identified with the P, Q, and R axes of dogfish LDHase-M₄. Since LDHase-C₄ crystallized with one molecule in a triclinic cell, the origin of the co-ordinate system was arbitrarily fixed at the molecular center. Structure factor phases were derived from an appropriately oriented dogfish apo LDHase-M₄ phasing model and combined with the observed structure amplitudes to produce a hybrid electron density map. Density points related by the molecular 222 point symmetry were averaged so as to remove the bias of the phasing model. At 7.5 Å resolution, the structure of the crystallized mouse LDHase-C₄ was found to be without coenzyme, with a conformation indistinguishable from that of dogfish apo LDHase-M₄.     

2.9 A resolution structure of the N-terminal domain of a variant surface glycoprotein from Trypanosoma brucei

The variant surface glycoprotein (VSG) of Trypanosoma brucei forms a coat on the surface of the parasite; by the expression of a series of antigenically distinct VSGs in the surface coat the parasite escapes the host immune response. The 2.9 A resolution crystal structure of the N-terminal domain of one variant, MITat 1.2, has been determined. The structure was solved using data collected from two crystal forms. Initially a partial model was built into an electron density map based on multiple isomorphous replacement phases and improved by phase combination methods. Subsequently this model was used to obtain the molecular replacement solution for a second crystal form, providing starting phases which were refined using 2-fold non-crystallographic symmetry averaging. The current model includes 362 residues and has been refined using X-PLOR to an R value of 0.22 for data between 7 and 2.9 A. The molecule is a dimer, approximately 100 A long, having an asymmetrical cross section with maximum dimensions of approximately 40 A x 60 A. Two long, approximately 70 A, alpha-helices from each monomer pack together to form, with several other helices, a core helix bundle that extends nearly the full length of the molecule. The "top" of the protein, which in the surface coat may be exposed to the external environment, is formed from the ends of the two long helices, a short three-stranded beta-sheet, and a strand having irregular conformation that packs above these secondary structure elements. Two conserved disulfide bridges are in this part of the molecule. Several elements of the MITat 1.2 sequence, which contribute to the formation of the helix bundle structure, have been identified. These elements can be found in the sequences of several different VSGs, suggesting that to some extent the VSG structure is conserved in those variants.     

Sorting duplicated loci disentangles complexities of polyploid genomes masked by genotyping by sequencing

Many plants and animals of polyploid origin are currently enjoying a genomics explosion enabled by modern sequencing and genotyping technologies. However, routine filtering of duplicated loci in most studies using genotyping by sequencing introduces an unacceptable, but often overlooked, bias when detecting selection. Retained duplicates from ancient whole‐genome duplications (WGDs) may be found throughout genomes, whereas retained duplicates from recent WGDs are concentrated at distal ends of some chromosome arms. Additionally, segmental duplicates can be found at distal ends or nearly anywhere in a genome. Evidence shows that these duplications facilitate adaptation through one of two pathways: neo‐functionalization or increased gene expression. Filtering duplicates removes distal ends of some chromosomes, and distal ends are especially known to harbour adaptively important genes. Thus, filtering of duplicated loci impoverishes the interpretation of genomic data as signals from contiguous duplicated genes are ignored. We review existing strategies to genotype and map duplicated loci; we focus in detail on an overlooked strategy of using gynogenetic haploids (1N) as a part of new genotyping by sequencing studies. We provide guidelines on how to use this haploid strategy for studies on polyploid‐origin vertebrates including how it can be used to screen duplicated loci in natural populations. We conclude by discussing areas of research that will benefit from better inclusion of polyploid loci; we particularly stress the sometimes overlooked fact that basing genomic studies on dense maps provides value added in the form of locating and annotating outlier loci or colocating outliers into islands of divergence.     

Inference from Stratified Samples: Study Design,Bias and Graphical Model Representations

Stratification is a widely used strategy in empirical research to improve efficiency of the sampling design. One concern of stratification is that ignoring it on analysis may bias the relationship between variables. A weighted analysis can only be carried out when sampling weights are known. When these are unknown, valid inference on the relationship between variables then depends on the ignorability of the design, which may be difficult to establish. Here, graphical representations of multivariate dependencies and independencies are used to find necessary conditions for ignorability of stratified sampling designs for inference on conditional and marginal relationships between variables.     

A rapid method for positioning small flexible molecules, nucleic acids, and large protein fragments in experimental electron density maps.

A Monte Carlo procedure, encoded in the program Blob, has been developed and tested for the purpose of positioning large molecular fragments or small flexible molecules in electron density maps. The search performed by the algorithm appears to be sufficiently thorough to accurately position a small flexible ligand in well-defined density while remaining sufficiently random to offer interesting alternate suggestions for density representing disordered binding modes of a ligand. Furthermore, the algorithm is shown to be efficient enough to accurately position large rigid molecular fragments. In the first of the test cases with large molecular fragments, Blob was surprisingly effective in positioning a poly-alanine model of a 53-residue domain in poor electron density resulting from molecular replacement with a partial model. At 3.0 A resolution the domain was positioned consistently within 0.2 A of its experimentally determined position. Even at 6.0 A resolution Blob could consistently position the domain to within 0.75 A of its actual position. A second set of tests with large molecular fragments revealed that Blob could correctly position large molecular fragments with quite significant deviations from the actual structure. In this test case, fragments ranging from a 170-residue protein domain with a 3.8 A rms deviation from the actual structure to a 22-base pair ideal B-form DNA duplex were positioned accurately in a 3.2 A electron density map derived from multiple isomorphous replacement methods. Even when decreasing the quality of the maps, from a figure of merit of 0.57 to as low as 0. 35, Blob could still effectively position the large protein domain and the DNA duplex. Since it is efficient, can handle large molecular fragments, and works in poor and low resolution maps, Blob could be a useful tool for interpreting electron density maps in de novo structure determinations and in molecular replacement studies. Proteins 1999;36:512-525.     

Structural characterization of pyoverdines produced by <Emphasis Type="Italic">Pseudomonas putida</Emphasis> KT2440 and <Emphasis Type="Italic">Pseudomonas taiwanensis</Emphasis> VLB120

The previously unknown sequences of several pyoverdines (PVD) produced by a biotechnologically-relevant bacterium, namely, Pseudomonas taiwanensis VLB120, were characterized by high performance liquid chromatography (HPLC)–high resolution mass spectrometry (HRMS). The same structural characterization scheme was checked before by analysis of Pseudomonas sp. putida KT2440 samples with known PVDs. A new sample preparation strategy based on solid-phase extraction was developed, requiring significantly reduced sample material as compared to existing methods. Chromatographic separation was performed using hydrophilic interaction liquid chromatography with gradient elution. Interestingly, no signals for apoPVDs were detected in these analyses, only the corresponding aluminum(III) and iron(III) complexes were seen. The chromatographic separation readily enabled separation of PVD complexes according to their individual structures. HPLC-HRMS and complementary fragmentation data from collision-induced dissociation and electron capture dissociation enabled the structural characterization of the investigated pyoverdines. In Pseudomonas sp. putida KT2240 samples, the known pyoverdines G4R and G4R A were readily confirmed. No PVDs have been previously described for Pseudomonas sp. taiwanensis VLB120. In our study, we identified three new PVDs, which only differed in their acyl side chains (succinic acid, succinic amide and malic acid). Peptide sequencing by MS/MS provided the sequence Orn-Asp-OHAsn-Thr-AcOHOrn-Ser-cOHOrn. Of particular interest is the presence of OHAsn, which has not been reported as PVD constituent before.