Monomer and dimer structures of cytochrome bo3 ubiquinol oxidase from Escherichia coli

The E. coli cytochrome bo₃ ubiquinol oxidase is a four-subunit heme-copper oxidase that serves as a proton pump in the E. coli aerobic respiratory chain. Despite many mechanistic studies, it is unclear whether this ubiquinol oxidase functions as a monomer, or as a dimer in a manner similar to its eukaryotic counterparts - the mitochondrial electron transport complexes. In this study, we determined the monomeric and dimeric structures of the E. coli cytochrome bo₃ ubiquinol oxidase reconstituted in amphipol by cryogenic electron microscopy single particle reconstruction (cryo-EM SPR) to a resolution of 3.15 Å and 3.46 Å, respectively. We have discovered that the protein can form a dimer with C2 symmetry, with the dimerization interface maintained by interactions between the subunit II of one monomer and the subunit IV of the other monomer. Moreover, the dimerization does not induce significant structural changes in the monomers, except the movement of a loop in subunit IV (residues 67-74). This article is protected by copyright. All rights reserved.

CryoEM single particle reconstruction with a complex-valued particle stack

Single particle reconstruction (SPR) in cryoEM is an image processing task with an elaborate hierarchy that starts with many very noisy multi-frame images. Efficient representation of the intermediary image structures is critical for keeping the calculations manageable. One such intermediary structure is called a particle stack and contains cut-out images of particles in square boxes of predefined size. The micrograph that is the source of the boxed images is usually corrected for motion between frames prior to particle stack creation. However, the contrast transfer function (CTF) or its Fourier Transform point spread function (PSF) are not considered at this step. Historically, the particle stack was intended for large particles and for a tighter PSF, which is characteristic of lower resolution data. The field now performs analyses of smaller particles and to higher resolution, and these conditions result in a broader PSF that requires larger padding and slower calculations to integrate information for each particle. Consequently, the approach to handling structures such as the particle stack should be reexamined to optimize data processing. Here we propose to use as a source image for the particle stack a complex-valued image, in which CTF correction is implicitly applied as a real component of the image. We can achieve it by applying an initial CTF correction to the entire micrograph first and perform box cutouts as a subsequent step. The final CTF correction that we refine and apply later has a very narrow PSF, and so cutting out particles from micrographs that were approximately corrected for CTF does not require extended buffering, i.e. the boxes during the analysis only have to be large enough to encompass the particle. The Fourier Transform of an exit-wave reconstruction creates an image that has complex values. This is a complex value image considered in real space, opposed to standard SPR data processing where complex numbers appear only in Fourier space. This extension of the micrograph concept provides multiple advantages because the particle box size can be small and calculations crucial for high resolution reconstruction such as Ewald sphere correction, aberration refinement, and particle-specific defocus refinement can be performed on the small box data.

The His-tag as a decoy modulating preferred orientation in cryoEM

The His-tag is a widely used affinity tag that facilitates purification by means of affinity chromatography of recombinant proteins for functional and structural studies. We show here that His-tag presence affects how coproheme decarboxylase interacts with the air-water interface during grid preparation for cryoEM. Depending on His-tag presence or absence, we observe significant changes in patterns of preferred orientation. Our analysis of particle orientations suggests that His-tag presence can mask the hydrophobic and hydrophilic patches on a protein’s surface that mediate the interactions with the air-water interface, while the hydrophobic linker between a His-tag and the coding sequence of the protein may enhance other interactions with the air-water interface. Our observations suggest that tagging, including rational design of the linkers between an affinity tag and a protein of interest, offer a promising approach to modulating interactions with the air-water interface.

Molecular analysis and essentiality of Aro1 shikimate biosynthesis multi-enzyme in Candida albicans

In the human fungal pathogen Candida albicans, ARO1 encodes an essential multi-enzyme that catalyses consecutive steps in the shikimate pathway for biosynthesis of chorismate, a precursor to folate and the aromatic amino acids. We obtained the first molecular image of C. albicans Aro1 that reveals the architecture of all five enzymatic domains and their arrangement in the context of the full-length protein. Aro1 forms a flexible dimer allowing relative autonomy of enzymatic function of the individual domains. Our activity and in cellulo data suggest that only four of Aro1's enzymatic domains are functional and essential for viability of C. albicans, whereas the 3-dehydroquinate dehydratase (DHQase) domain is inactive because of active site substitutions. We further demonstrate that in C. albicans, the type II DHQase Dqd1 can compensate for the inactive DHQase domain of Aro1, suggesting an unrecognized essential role for this enzyme in shikimate biosynthesis. In contrast, in Candida glabrata and Candida parapsilosis, which do not encode a Dqd1 homolog, Aro1 DHQase domains are enzymatically active, highlighting diversity across Candida species.

Optimal structure determination from sub-optimal diffraction data

Herein we present the newest version of the HKL-3000 system that integrates data collection, data reduction, phasing, model building, refinement, and validation. The system significantly accelerates the process of structure determination and has proven its high value for the determination of very high-quality structures. The heuristic for choosing the best approach for every step of structure determination for various quality samples and diffraction data has been optimized. The latest modifications increase the likelihood of a successful structure determination with challenging data. The HKL-3000 is a successor of HKL and HKL-2000 programs. The use of the HKL family of programs has been reported for over 73,000 PDB deposits, that is, almost 50% of macromolecular structures determined with X-ray diffraction.

Structural basis for IFN antagonism by human respiratory syncytial virus nonstructural protein 2

Human respiratory syncytial virus (RSV) nonstructural protein 2 (NS2) inhibits host interferon (IFN) responses stimulated by RSV infection by targeting early steps in the IFN-signaling pathway. But the molecular mechanisms related to how NS2 regulates these processes remain incompletely understood. To address this gap, here we solved the X-ray crystal structure of NS2. This structure revealed a unique fold that is distinct from other known viral IFN antagonists, including RSV NS1. We also show that NS2 directly interacts with an inactive conformation of the RIG-I-like receptors (RLRs) RIG-I and MDA5. NS2 binding prevents RLR ubiquitination, a process critical for prolonged activation of downstream signaling. Structural analysis, including by hydrogen-deuterium exchange coupled to mass spectrometry, revealed that the N terminus of NS2 is essential for binding to the RIG-I caspase activation and recruitment domains. N-terminal mutations significantly diminish RIG-I interactions and result in increased IFNβ messenger RNA levels. Collectively, our studies uncover a previously unappreciated regulatory mechanism by which NS2 further modulates host responses and define an approach for targeting host responses.

Synchrotron Radiation as a Tool for Macromolecular X-Ray Crystallography: a XXI Century Perspective

Intense X-rays available at powerful synchrotron beamlines provide macromolecular crystallographers with an incomparable tool for investigating biological phenomena on an atomic scale. The resulting insights into the mechanism's underlying biological processes have played an essential role and shaped biomedical sciences during the last 30 years, considered the "golden age" of structural biology. In this review, we analyze selected aspects of the impact of synchrotron radiation on structural biology. Synchrotron beamlines have been used to determine over 70% of all macromolecular structures deposited into the Protein Data Bank (PDB). These structures were deposited by over 13,000 different research groups. Interestingly, despite the impressive advances in synchrotron technologies, the median resolution of macromolecular structures determined using synchrotrons has remained constant throughout the last 30 years, at about 2 Å. Similarly, the median times from the data collection to the deposition and release have not changed significantly. We describe challenges to reproducibility related to recording all relevant data and metadata during the synchrotron experiments, including diffraction images. Finally, we discuss some of the recent opinions suggesting a diminishing importance of X-ray crystallography due to impressive advances in Cryo-EM and theoretical modeling. We believe that synchrotrons of the future will increasingly evolve towards a life science center model, where X-ray crystallography, Cryo-EM, and other experimental and computational resources and knowledge are encompassed within a versatile research facility. The recent response of crystallographers to the COVID-19 pandemic suggests that X-ray crystallography conducted at synchrotron beamlines will continue to play an essential role in structural biology and drug discovery for years to come.

Mutation severity spectrum of rare alleles in the human genome is predictive of disease type

The human genome harbors a variety of genetic variations. Single-nucleotide changes that alter amino acids in protein-coding regions are one of the major causes of human phenotypic variation and diseases. These single-amino acid variations (SAVs) are routinely found in whole genome and exome sequencing. Evaluating the functional impact of such genomic alterations is crucial for diagnosis of genetic disorders. We developed DeepSAV, a deep-learning convolutional neural network to differentiate disease-causing and benign SAVs based on a variety of protein sequence, structural and functional properties. Our method outperforms most stand-alone programs, and the version incorporating population and gene-level information (DeepSAV+PG) has similar predictive power as some of the best available. We transformed DeepSAV scores of rare SAVs in the human population into a quantity termed "mutation severity measure" for each human protein-coding gene. It reflects a gene's tolerance to deleterious missense mutations and serves as a useful tool to study gene-disease associations. Genes implicated in cancer, autism, and viral interaction are found by this measure as intolerant to mutations, while genes associated with a number of other diseases are scored as tolerant. Among known disease-associated genes, those that are mutation-intolerant are likely to function in development and signal transduction pathways, while those that are mutation-tolerant tend to encode metabolic and mitochondrial proteins.

High-resolution cryo-EM reconstructions in the presence of substantial aberrations

Here, an analysis is performed of how uncorrected antisymmetric aberrations, such as coma and trefoil, affect cryo-EM single-particle reconstruction (SPR) results, and an analytical formula quantifying information loss owing to their presence is inferred that explains why Fourier-shell coefficient-based statistics may report significantly overestimated resolution if these aberrations are not fully corrected. The analysis is validated with reference-based aberration refinement for two cryo-EM SPR data sets acquired with a 200 kV microscope in the presence of coma exceeding 40 µm, and 2.3 and 2.7 Å reconstructions for 144 and 173 kDa particles, respectively, were obtained. The results provide a description of an efficient approach for assessing information loss in cryo-EM SPR data acquired in the presence of higher order aberrations, and address inconsistent guidelines regarding the level of aberrations that is acceptable in cryo-EM SPR experiments.

NMR and crystallographic structural studies of the extremely stable monomeric variant of human cystatin C with single amino acid substitution

Human cystatin C (hCC), a member of the superfamily of papain-like cysteine protease inhibitors, is the most widespread cystatin in human body fluids. This small protein, in addition to its physiological function, is involved in various diseases, including cerebral amyloid angiopathy, cerebral hemorrhage, stroke, and dementia. Physiologically active hCC is a monomer. However, all structural studies based on crystallization led to the dimeric structure formed as a result of a three-dimensional exchange of the protein domains (3D domain swapping). The monomeric structure was obtained only for hCC variant V57N and for the protein stabilized by an additional disulfide bridge. With this study, we extend the number of models of monomeric hCC by an additional hCC variant with a single amino acid substitution in the flexible loop L1. The V57G variant was chosen for the X-ray and NMR structural analysis due to its exceptional conformational stability in solution. In this work, we show for the first time the structural and dynamics studies of human cystatin C variant in solution. We were also able to compare these data with the crystal structure of the hCC V57G and with other cystatins. The overall cystatin fold is retained in the solute form. Additionally, structural information concerning the N terminus was obtained during our studies and presented for the first time. DATABASE: Crystallographic structure: structural data are available in PDB databases under the accession number 6ROA. NMR structure: structural data are available in PDB and BMRB databases under the accession numbers 6RPV and 34399, respectively.

Crystal structure of native α-L-rhamnosidase from Aspergillus terreus

α-L-Rhamnosidases cleave terminal nonreducing α-L-rhamnosyl residues from many natural rhamnoglycosides. This makes them catalysts of interest for various biotechnological applications. The X-ray structure of the GH78 family α-L-rhamnosidase from Aspergillus terreus has been determined at 1.38 Å resolution using the sulfur single-wavelength anomalous dispersion phasing method. The protein was isolated from its natural source in the native glycosylated form, and the active site contained a glucose molecule, probably from the growth medium. In addition to its catalytic domain, the α-L-rhamnosidase from A. terreus contains four accessory domains of unknown function. The structural data suggest that two of these accessory domains, E and F, might play a role in stabilizing the aglycon portion of the bound substrate.

The Bear Giant-Skipper genome suggests genetic adaptations to living inside yucca roots

Giant-Skippers (Megathymini) are unusual thick-bodied, moth-like butterflies whose caterpillars feed inside Yucca roots and Agave leaves. Giant-Skippers are attributed to the subfamily Hesperiinae and they are endemic to southern and mostly desert regions of the North American continent. To shed light on the genotypic determinants of their unusual phenotypic traits, we sequenced and annotated a draft genome of the largest Giant-Skipper species, the Bear (Megathymus ursus violae). The Bear skipper genome is the least heterozygous among sequenced Lepidoptera genomes, possibly due to much smaller population size and extensive inbreeding. Their lower heterozygosity helped us to obtain a high-quality genome with an N50 of 4.2 Mbp. The ~ 430 Mb genome encodes about 14000 proteins. Phylogenetic analysis supports placement of Giant-Skippers with Grass-Skippers (Hesperiinae). We find that proteins involved in odorant and taste sensing as well as in oxidative reactions have diverged significantly in Megathymus as compared to Lerema, another Grass-Skipper. In addition, the Giant-Skipper has lost several odorant and gustatory receptors and possesses many fewer (1/3-1/2 of other skippers) anti-oxidative enzymes. Such differences may be related to the unusual life style of Giant-Skippers: they do not feed as adults, and their caterpillars feed inside Yuccas and Agaves, which provide a source of antioxidants such as polyphenols.

Real-space analysis of radiation-induced specific changes with independent component analysis

A method of analysis is presented that allows for the separation of specific radiation-induced changes into distinct components in real space. The method relies on independent component analysis (ICA) and can be effectively applied to electron density maps and other types of maps, provided that they can be represented as sets of numbers on a grid. Here, for glucose isomerase crystals, ICA was used in a proof-of-concept analysis to separate temperature-dependent and temperature-independent components of specific radiation-induced changes for data sets acquired from multiple crystals across multiple temperatures. ICA identified two components, with the temperature-independent component being responsible for the majority of specific radiation-induced changes at temperatures below 130 K. The patterns of specific temperature-independent radiation-induced changes suggest a contribution from the tunnelling of electron holes as a possible explanation. In the second case, where a group of 22 data sets was collected on a single thaumatin crystal, ICA was used in another type of analysis to separate specific radiation-induced effects happening on different exposure-level scales. Here, ICA identified two components of specific radiation-induced changes that likely result from radiation-induced chemical reactions progressing with different rates at different locations in the structure. In addition, ICA unexpectedly identified the radiation-damage state corresponding to reduced disulfide bridges rather than the zero-dose extrapolated state as the highest contrast structure. The application of ICA to the analysis of specific radiation-induced changes in real space and the data pre-processing for ICA that relies on singular value decomposition, which was used previously in data space to validate a two-component physical model of X-ray radiation-induced changes, are discussed in detail. This work lays a foundation for a better understanding of protein-specific radiation chemistries and provides a framework for analysing effects of specific radiation damage in crystallographic and cryo-EM experiments.

When COI barcodes deceive: complete genomes reveal introgression in hairstreaks

Two species of hairstreak butterflies from the genus Calycopis are known in the United States: C. cecrops and C. isobeon Analysis of mitochondrial COI barcodes of Calycopis revealed cecrops-like specimens from the eastern US with atypical barcodes that were 2.6% different from either USA species, but similar to Central American Calycopis species. To address the possibility that the specimens with atypical barcodes represent an undescribed cryptic species, we sequenced complete genomes of 27 Calycopis specimens of four species: C. cecrops, C. isobeon, C. quintana and C. bactra Some of these specimens were collected up to 60 years ago and preserved dry in museum collections, but nonetheless produced genomes as complete as fresh samples. Phylogenetic trees reconstructed using the whole mitochondrial and nuclear genomes were incongruent. While USA Calycopis with atypical barcodes grouped with Central American species C. quintana by mitochondria, nuclear genome trees placed them within typical USA C. cecrops in agreement with morphology, suggesting mitochondrial introgression. Nuclear genomes also show introgression, especially between C. cecrops and C. isobeon About 2.3% of each C. cecrops genome has probably (p-value < 0.01, FDR < 0.1) introgressed from C. isobeon and about 3.4% of each C. isobeon genome may have come from C. cecrops. The introgressed regions are enriched in genes encoding transmembrane proteins, mitochondria-targeting proteins and components of the larval cuticle. This study provides the first example of mitochondrial introgression in Lepidoptera supported by complete genome sequencing. Our results caution about relying solely on COI barcodes and mitochondrial DNA for species identification or discovery.

The first complete genomes of Metalmarks and the classification of butterfly families

Sequencing complete genomes of all major phylogenetic groups of organisms opens unprecedented opportunities to study evolution and genetics. We report draft genomes of Calephelis nemesis and Calephelis virginiensis, representatives of the family Riodinidae. They complete the genomic coverage of butterflies at the family level. At 809 and 855 Mbp, respectively, they become the largest available Lepidoptera genomes. Comparison of butterfly genomes shows that the divergence between Riodinidae and Lycaenidae dates to the time when other families started to diverge into subfamilies. Thus, Riodinidae may be considered a subfamily of Lycaenidae. Calephelis species exhibit unique gene expansions in actin-disassembling factor, cofilin, and chitinase. The functional implications of these gene expansions are not clear, but they may aid molting of caterpillars covered in extensive setae. The two Calephelis species diverged about 5 million years ago and they differ in proteins involved in metabolism, circadian clock, regulation of development, and immune responses.

Complete Genome of Achalarus lyciades, The First Representative of the Eudaminae Subfamily of Skippers

Background:

The Hoary Edge Skipper (Achalarus lyciades) is an eastern North America endemic butterfly from the Eudaminae subfamily of skippers named for an underside whitish patch near the hindwing edge. Its caterpillars feed on legumes, in contrast to Grass skippers (subfamily Hesperiinae) which feed exclusively on monocots.

Results:

To better understand the evolution and phenotypic diversification of Skippers (family Hesperiidae), we sequenced, assembled and annotated a complete genome draft and transcriptome of a wild-caught specimen of A. lyciades and compared it with the available genome of the Clouded Skipper (Lerema accius) from the Grass skipper subfamily. The genome of A. lyciades is nearly twice the size of L. accius (567 Mbp vs. 298 Mbp), however it encodes a smaller number of proteins (15881 vs. 17411). Gene expansions we identified previously in L. accius apparently did not occur in the genome of A. lyciades. For instance, a family of hypothetical cellulases that diverged from an endochitinase (possibly associated with feeding of L. accius caterpillars on nutrient-poor grasses) is absent in A. lyciades. While L. accius underwent gene expansion in pheromone binding proteins, A. lyciades has more opsins. This difference may be related to the mate recognition mechanisms of the two species: visual cues might be more important for the Eudaminae skippers (which have more variable wing patterns), whereas odor might be more important for Grass skippers (that are hardly distinguishable by their wings). Phylogenetically, A. lyciades is a sister species of L. accius, the only other Hesperiidae with a complete genome.

Conclusions:

A new reference genome of a dicot-feeding skippers, the first from the Eudaminae subfamily, reveals its larger size and suggests hypotheses about phenotypic traits and differences from monocot-feeding skippers.

Complete genome of Pieris rapae, a resilient alien, a cabbage pest, and a source of anti-cancer proteins

The Small Cabbage White ( Pieris rapae) is originally a Eurasian butterfly. Being accidentally introduced into North America, Australia, and New Zealand a century or more ago, it spread throughout the continents and rapidly established as one of the most abundant butterfly species. Although it is a serious pest of cabbage and other mustard family plants with its caterpillars reducing crops to stems, it is also a source of pierisin, a protein unique to the Whites that shows cytotoxicity to cancer cells. To better understand the unusual biology of this omnipresent agriculturally and medically important butterfly, we sequenced and annotated the complete genome from USA specimens. At 246 Mbp, it is among the smallest Lepidoptera genomes reported to date. While 1.5% positions in the genome are heterozygous, they are distributed highly non-randomly along the scaffolds, and nearly 20% of longer than 1000 base-pair segments are SNP-free (median length: 38000 bp). Computational simulations of population evolutionary history suggest that American populations started from a very small number of introduced individuals, possibly a single fertilized female, which is in agreement with historical literature. Comparison to other Lepidoptera genomes reveals several unique families of proteins that may contribute to the unusual resilience of Pieris. The nitrile-specifier proteins divert the plant defense chemicals to non-toxic products. The apoptosis-inducing pierisins could offer a defense mechanism against parasitic wasps. While only two pierisins from Pieris rapae were characterized before, the genome sequence revealed eight, offering additional candidates as anti-cancer drugs. The reference genome we obtained lays the foundation for future studies of the Cabbage White and other Pieridae species.

Crystal structure of the cohesin loader Scc2 and insight into cohesinopathy

The ring-shaped cohesin complex topologically entraps chromosomes and regulates chromosome segregation, transcription, and DNA repair. The cohesin core consists of the structural maintenance of chromosomes 1 and 3 (Smc1-Smc3) heterodimeric ATPase, the kleisin subunit sister chromatid cohesion 1 (Scc1) that links the two ATPase heads, and the Scc1-bound adaptor protein Scc3. The sister chromatid cohesion 2 and 4 (Scc2-Scc4) complex loads cohesin onto chromosomes. Mutations of cohesin and its regulators, including Scc2, cause human developmental diseases termed cohesinopathy. Here, we report the crystal structure of Chaetomium thermophilum (Ct) Scc2 and examine its interaction with cohesin. Similar to Scc3 and another Scc1-interacting cohesin regulator, precocious dissociation of sisters 5 (Pds5), Scc2 consists mostly of helical repeats that fold into a hook-shaped structure. Scc2 binds to Scc1 through an N-terminal region of Scc1 that overlaps with its Pds5-binding region. Many cohesinopathy mutations target conserved residues in Scc2 and diminish Ct Scc2 binding to Ct Scc1. Pds5 binding to Scc1 weakens the Scc2-Scc1 interaction. Our study defines a functionally important interaction between the kleisin subunit of cohesin and the hook of Scc2. Through competing with Scc2 for Scc1 binding, Pds5 might contribute to the release of Scc2 from loaded cohesin, freeing Scc2 for additional rounds of loading.

Phylogeny Reconstruction with Alignment-Free Method That Corrects for Horizontal Gene Transfer

Advances in sequencing have generated a large number of complete genomes. Traditionally, phylogenetic analysis relies on alignments of orthologs, but defining orthologs and separating them from paralogs is a complex task that may not always be suited to the large datasets of the future. An alternative to traditional, alignment-based approaches are whole-genome, alignment-free methods. These methods are scalable and require minimal manual intervention. We developed SlopeTree, a new alignment-free method that estimates evolutionary distances by measuring the decay of exact substring matches as a function of match length. SlopeTree corrects for horizontal gene transfer, for composition variation and low complexity sequences, and for branch-length nonlinearity caused by multiple mutations at the same site. We tested SlopeTree on 495 bacteria, 73 archaea, and 72 strains of Escherichia coli and Shigella. We compared our trees to the NCBI taxonomy, to trees based on concatenated alignments, and to trees produced by other alignment-free methods. The results were consistent with current knowledge about prokaryotic evolution. We assessed differences in tree topology over different methods and settings and found that the majority of bacteria and archaea have a core set of proteins that evolves by descent. In trees built from complete genomes rather than sets of core genes, we observed some grouping by phenotype rather than phylogeny, for instance with a cluster of sulfur-reducing thermophilic bacteria coming together irrespective of their phyla. The source-code for SlopeTree is available at: http://prodata.swmed.edu/download/pub/slopetree_v1/slopetree.tar.gz.

Due to their lack of distinct morphological features, bacteria and archaea were extremely difficult to classify until technology was developed to obtain their DNA sequences; these sequences could then be compared to estimate evolutionary relationships. Now, due to technological advances, there is a flood of available sequences from a wide variety of organisms. These advances have spurred the development of algorithms which can estimate evolutionary relationships using whole genomes, in contrast to the more traditional methods which used single genes earlier and now typically use groups of conserved genes. However, there are many challenges when attempting to infer evolutionary relationships, in particular horizontal gene transfer, where DNA is transferred from one organism to another, resulting in an organism’s genome containing DNA that does not reflect its evolution by descent. We developed a new whole-genome method for estimating evolutionary distances which identifies and corrects for horizontal transfer. We found that for SlopeTree and all other whole-genome methods we applied, horizontal transfer causes some evolutionary distances to be grossly underestimated, and that our correction corrects for this.

Estimation of Uncertainties in the Global Distance Test (GDT_TS) for CASP Models

The Critical Assessment of techniques for protein Structure Prediction (or CASP) is a community-wide blind test experiment to reveal the best accomplishments of structure modeling. Assessors have been using the Global Distance Test (GDT_TS) measure to quantify prediction performance since CASP3 in 1998. However, identifying significant score differences between close models is difficult because of the lack of uncertainty estimations for this measure. Here, we utilized the atomic fluctuations caused by structure flexibility to estimate the uncertainty of GDT_TS scores. Structures determined by nuclear magnetic resonance are deposited as ensembles of alternative conformers that reflect the structural flexibility, whereas standard X-ray refinement produces the static structure averaged over time and space for the dynamic ensembles. To recapitulate the structural heterogeneous ensemble in the crystal lattice, we performed time-averaged refinement for X-ray datasets to generate structural ensembles for our GDT_TS uncertainty analysis. Using those generated ensembles, our study demonstrates that the time-averaged refinements produced structure ensembles with better agreement with the experimental datasets than the averaged X-ray structures with B-factors. The uncertainty of the GDT_TS scores, quantified by their standard deviations (SDs), increases for scores lower than 50 and 70, with maximum SDs of 0.3 and 1.23 for X-ray and NMR structures, respectively. We also applied our procedure to the high accuracy version of GDT-based score and produced similar results with slightly higher SDs. To facilitate score comparisons by the community, we developed a user-friendly web server that produces structure ensembles for NMR and X-ray structures and is accessible at http://prodata.swmed.edu/SEnCS. Our work helps to identify the significance of GDT_TS score differences, as well as to provide structure ensembles for estimating SDs of any scores.

Complete genomes of Hairstreak butterflies, their speciation, and nucleo-mitochondrial incongruence

Comparison of complete genomes of closely related species enables research on speciation and how phenotype is determined by genotype. Lepidoptera, an insect order of 150,000 species with diverse phenotypes, is well-suited for such comparative genomics studies if new genomes, which cover additional Lepidoptera families are acquired. We report a 729 Mbp genome assembly of the Calycopis cecrops, the first genome from the family Lycaenidae and the largest available Lepidoptera genome. As detritivore, Calycopis shows expansion in detoxification and digestion enzymes. We further obtained complete genomes of 8 Calycopis specimens: 3 C. cecrops and 5 C. isobeon, including a dry specimen stored in the museum for 30 years. The two species differ subtly in phenotype and cannot be differentiated by mitochondrial DNA. However, nuclear genomes revealed a deep split between them. Genes that can clearly separate the two species (speciation hotspots) mostly pertain to circadian clock, mating behavior, transcription regulation, development and cytoskeleton. The speciation hotspots and their function significantly overlap with those we previously found in Pterourus, suggesting common speciation mechanisms in these butterflies.

Taking a look at the calibration of a CCD detector with a fiber-optic taper

At the Structural Biology Center beamline 19BM, located at the Advanced Photon Source, the operational characteristics of the equipment are routinely checked to ensure they are in proper working order. After performing a partial flat-field calibration for the ADSC Quantum 210r CCD detector, it was confirmed that the detector operates within specifications. However, as a secondary check it was decided to scan a single reflection across one-half of a detector module to validate the accuracy of the calibration. The intensities from this single reflection varied by more than 30% from the module center to the corner of the module. Redistribution of light within bent fibers of the fiber-optic taper was identified to be a source of this variation. The degree to which the diffraction intensities are corrected to account for characteristics of the fiber-optic tapers depends primarily upon the experimental strategy of data collection, approximations made by the data processing software during scaling, and crystal symmetry.

Speciation in Cloudless Sulphurs Gleaned from Complete Genomes

For 200 years, zoologists have relied on phenotypes to learn about the evolution of animals. A glance at the genotype, even through several gene markers, revolutionized our understanding of animal phylogeny. Recent advances in sequencing techniques allow researchers to study speciation mechanisms and the link between genotype and phenotype using complete genomes. We sequenced and assembled a complete genome of the Cloudless Sulphur (Phoebis sennae) from a single wild-caught specimen. This genome was used as reference to compare genomes of six specimens, three from the eastern populations (Oklahoma and north Texas), referred to as a subspeciesPhoebis sennae eubule, and three from the southwestern populations (south Texas) known as a subspeciesPhoebis sennae marcellina While the two subspecies differ only subtly in phenotype and mitochondrial DNA, comparison of their complete genomes revealed consistent and significant differences, which are more prominent than those between tiger swallowtailsPterourus canadensisandPterourus glaucus The two sulphur taxa differed in histone methylation regulators, chromatin-associated proteins, circadian clock, and early development proteins. Despite being well separated on the whole-genome level, the two taxa show introgression, with gene flow mainly fromP. s. marcellinatoP. s. eubule Functional analysis of introgressed genes reveals enrichment in transmembrane transporters. Many transporters are responsible for nutrient uptake, and their introgression may be of selective advantage for caterpillars to feed on more diverse food resources. Phylogenetically, complete genomes place family Pieridae away from Papilionidae, which is consistent with previous analyses based on several gene markers.

Data publication with the structural biology data grid supports live analysis

Access to experimental X-ray diffraction image data is fundamental for validation and reproduction of macromolecular models and indispensable for development of structural biology processing methods. Here, we established a diffraction data publication and dissemination system, Structural Biology Data Grid (SBDG; data.sbgrid.org), to preserve primary experimental data sets that support scientific publications. Data sets are accessible to researchers through a community driven data grid, which facilitates global data access. Our analysis of a pilot collection of crystallographic data sets demonstrates that the information archived by SBDG is sufficient to reprocess data to statistics that meet or exceed the quality of the original published structures. SBDG has extended its services to the entire community and is used to develop support for other types of biomedical data sets. It is anticipated that access to the experimental data sets will enhance the paradigm shift in the community towards a much more dynamic body of continuously improving data analysis.

Skipper genome sheds light on unique phenotypic traits and phylogeny

BACKGROUND

Butterflies and moths are emerging as model organisms in genetics and evolutionary studies. The family Hesperiidae (skippers) was traditionally viewed as a sister to other butterflies based on its moth-like morphology and darting flight habits with fast wing beats. However, DNA studies suggest that the family Papilionidae (swallowtails) may be the sister to other butterflies including skippers. The moth-like features and the controversial position of skippers in Lepidoptera phylogeny make them valuable targets for comparative genomics.

RESULTS

We obtained the 310 Mb draft genome of the Clouded Skipper (Lerema accius) from a wild-caught specimen using a cost-effective strategy that overcomes the high (1.6 %) heterozygosity problem. Comparative analysis of Lerema accius and the highly heterozygous genome of Papilio glaucus revealed differences in patterns of SNP distribution, but similarities in functions of genes that are enriched in non-synonymous SNPs. Comparison of Lepidoptera genomes revealed possible molecular bases for unique traits of skippers: a duplication of electron transport chain components could result in efficient energy supply for their rapid flight; a diversified family of predicted cellulases might allow them to feed on cellulose-enriched grasses; an expansion of pheromone-binding proteins and enzymes for pheromone synthesis implies a more efficient mate-recognition system, which compensates for the lack of clear visual cues due to the similarities in wing colors and patterns of many species of skippers. Phylogenetic analysis of several Lepidoptera genomes suggested that the position of Hesperiidae remains uncertain as the tree topology varied depending on the evolutionary model.

CONCLUSION

Completion of the first genome from the family Hesperiidae allowed comparative analyses with other Lepidoptera that revealed potential genetic bases for the unique phenotypic traits of skippers. This work lays the foundation for future experimental studies of skippers and provides a rich dataset for comparative genomics and phylogenetic studies of Lepidoptera.

An Intrinsically Disordered Peptide from Ebola Virus VP35 Controls Viral RNA Synthesis by Modulating Nucleoprotein-RNA Interactions

During viral RNA synthesis, Ebola virus (EBOV) nucleoprotein (NP) alternates between an RNA-template-bound form and a template-free form to provide the viral polymerase access to the RNA template. In addition, newly synthesized NP must be prevented from indiscriminately binding to noncognate RNAs. Here, we investigate the molecular bases for these critical processes. We identify an intrinsically disordered peptide derived from EBOV VP35 (NPBP, residues 20-48) that binds NP with high affinity and specificity, inhibits NP oligomerization, and releases RNA from NP-RNA complexes in vitro. The structure of the NPBP/ΔNPNTD complex, solved to 3.7 Å resolution, reveals how NPBP peptide occludes a large surface area that is important for NP-NP and NP-RNA interactions and for viral RNA synthesis. Together, our results identify a highly conserved viral interface that is important for EBOV replication and can be targeted for therapeutic development.

Refinement by shifting secondary structure elements improves sequence alignments

Constructing a model of a query protein based on its alignment to a homolog with experimentally determined spatial structure (the template) is still the most reliable approach to structure prediction. Alignment errors are the main bottleneck for homology modeling when the query is distantly related to the template. Alignment methods often misalign secondary structural elements by a few residues. Therefore, better alignment solutions can be found within a limited set of local shifts of secondary structures. We present a refinement method to improve pairwise sequence alignments by evaluating alignment variants generated by local shifts of template-defined secondary structures. Our method SFESA is based on a novel scoring function that combines the profile-based sequence score and the structure score derived from residue contacts in a template. Such a combined score frequently selects a better alignment variant among a set of candidate alignments generated by local shifts and leads to overall increase in alignment accuracy. Evaluation of several benchmarks shows that our refinement method significantly improves alignments made by automatic methods such as PROMALS, HHpred and CNFpred. The web server is available at http://prodata.swmed.edu/sfesa.

A tethered delivery mechanism explains the catalytic action of a microtubule polymerase

Stu2p/XMAP215 proteins are essential microtubule polymerases that use multiple αβ-tubulin-interacting TOG domains to bind microtubule plus ends and catalyze fast microtubule growth. We report here the structure of the TOG2 domain from Stu2p bound to yeast αβ-tubulin. Like TOG1, TOG2 binds selectively to a fully 'curved' conformation of αβ-tubulin, incompatible with a microtubule lattice. We also show that TOG1-TOG2 binds non-cooperatively to two αβ-tubulins. Preferential interactions between TOGs and fully curved αβ-tubulin that cannot exist elsewhere in the microtubule explain how these polymerases localize to the extreme microtubule end. We propose that these polymerases promote elongation because their linked TOG domains concentrate unpolymerized αβ-tubulin near curved subunits already bound at the microtubule end. This tethering model can explain catalyst-like behavior and also predicts that the polymerase action changes the configuration of the microtubule end.

Comparative genome sequencing reveals chemotype-specific gene clusters in the toxigenic black mold Stachybotrys

BACKGROUND

The fungal genus Stachybotrys produces several diverse toxins that affect human health. Its strains comprise two mutually-exclusive toxin chemotypes, one producing satratoxins, which are a subclass of trichothecenes, and the other producing the less-toxic atranones. To determine the genetic basis for chemotype-specific differences in toxin production, the genomes of four Stachybotrys strains were sequenced and assembled de novo. Two of these strains produce atranones and two produce satratoxins.

RESULTS

Comparative analysis of these four 35-Mbp genomes revealed several chemotype-specific gene clusters that are predicted to make secondary metabolites. The largest, which was named the core atranone cluster, encodes 14 proteins that may suffice to produce all observed atranone compounds via reactions that include an unusual Baeyer-Villiger oxidation. Satratoxins are suggested to be made by products of multiple gene clusters that encode 21 proteins in all, including polyketide synthases, acetyltransferases, and other enzymes expected to modify the trichothecene skeleton. One such satratoxin chemotype-specific cluster is adjacent to the core trichothecene cluster, which has diverged from those of other trichothecene producers to contain a unique polyketide synthase.

CONCLUSIONS

The results suggest that chemotype-specific gene clusters are likely the genetic basis for the mutually-exclusive toxin chemotypes of Stachybotrys. A unified biochemical model for Stachybotrys toxin production is presented. Overall, the four genomes described here will be useful for ongoing studies of this mold's diverse toxicity mechanisms.

Structure of the effector-binding domain of deoxyribonucleoside regulator DeoR from Bacillus subtilis

Deoxyribonucleoside regulator (DeoR) from Bacillus subtilis negatively regulates expression of enzymes involved in the catabolism of deoxyribonucleosides and deoxyribose. The DeoR protein is homologous to the sorbitol operon regulator family of metabolic regulators and comprises an N-terminal DNA-binding domain and a C-terminal effector-binding domain. We have determined the crystal structure of the effector-binding domain of DeoR (C-DeoR) in free form and in covalent complex with its effector deoxyribose-5-phosphate (dR5P). This is the first case of a covalently attached effector molecule captured in the structure of a bacterial transcriptional regulator. The dR5P molecule is attached through a Schiff base linkage to residue Lys141. The crucial role of Lys141 in effector binding was confirmed by mutational analysis and mass spectrometry of Schiff base adducts formed in solution. Structural analyses of the free and effector-bound C-DeoR structures provided a structural explanation for the mechanism of DeoR function as a molecular switch.

DATABASES

Atomic coordinates and structure factors for crystal structures of free C-DeoR and the covalent Schiff base complex of C-DeoR with dR5P have been deposited in the Protein Data Bank with accession codes 4OQQ and 4OQP, respectively.

STRUCTURED DIGITAL ABSTRACT

C-DeoR and C-DeoR bind by x-ray crystallography (View interaction) DeoR and DeoR bind by molecular sieving (1, 2).

Extracellular architecture of the SYG-1/SYG-2 adhesion complex instructs synaptogenesis

SYG-1 and SYG-2 are multipurpose cell adhesion molecules (CAMs) that have evolved across all major animal taxa to participate in diverse physiological functions, ranging from synapse formation to formation of the kidney filtration barrier. In the crystal structures of several SYG-1 and SYG-2 orthologs and their complexes, we find that SYG-1 orthologs homodimerize through a common, bispecific interface that similarly mediates an unusual orthogonal docking geometry in the heterophilic SYG-1/SYG-2 complex. C. elegans SYG-1's specification of proper synapse formation in vivo closely correlates with the heterophilic complex affinity, which appears to be tuned for optimal function. Furthermore, replacement of the interacting domains of SYG-1 and SYG-2 with those from CAM complexes that assume alternative docking geometries or the introduction of segmental flexibility compromised synaptic function. These results suggest that SYG extracellular complexes do not simply act as "molecular velcro" and that their distinct structural features are important in instructing synaptogenesis. PAPERFLICK:

In silico derived small molecules bind the filovirus VP35 protein and inhibit its polymerase cofactor activity

The Ebola virus (EBOV) genome only encodes a single viral polypeptide with enzymatic activity, the viral large (L) RNA-dependent RNA polymerase protein. However, currently, there is limited information about the L protein, which has hampered the development of antivirals. Therefore, antifiloviral therapeutic efforts must include additional targets such as protein-protein interfaces. Viral protein 35 (VP35) is multifunctional and plays important roles in viral pathogenesis, including viral mRNA synthesis and replication of the negative-sense RNA viral genome. Previous studies revealed that mutation of key basic residues within the VP35 interferon inhibitory domain (IID) results in significant EBOV attenuation, both in vitro and in vivo. In the current study, we use an experimental pipeline that includes structure-based in silico screening and biochemical and structural characterization, along with medicinal chemistry, to identify and characterize small molecules that target a binding pocket within VP35. NMR mapping experiments and high-resolution x-ray crystal structures show that select small molecules bind to a region of VP35 IID that is important for replication complex formation through interactions with the viral nucleoprotein (NP). We also tested select compounds for their ability to inhibit VP35 IID-NP interactions in vitro as well as VP35 function in a minigenome assay and EBOV replication. These results confirm the ability of compounds identified in this study to inhibit VP35-NP interactions in vitro and to impair viral replication in cell-based assays. These studies provide an initial framework to guide development of antifiloviral compounds against filoviral VP35 proteins.

Novel pyridazinone inhibitors for vascular adhesion protein-1 (VAP-1): old target-new inhibition mode

Vascular adhesion protein-1 (VAP-1) is a primary amine oxidase and a drug target for inflammatory and vascular diseases. Despite extensive attempts to develop potent, specific, and reversible inhibitors of its enzyme activity, the task has proven challenging. Here we report the synthesis, inhibitory activity, and molecular binding mode of novel pyridazinone inhibitors, which show specificity for VAP-1 over monoamine and diamine oxidases. The crystal structures of three inhibitor-VAP-1 complexes show that these compounds bind reversibly into a unique binding site in the active site channel. Although they are good inhibitors of human VAP-1, they do not inhibit rodent VAP-1 well. To investigate this further, we used homology modeling and structural comparison to identify amino acid differences, which explain the species-specific binding properties. Our results prove the potency and specificity of these new inhibitors, and the detailed characterization of their binding mode is of importance for further development of VAP-1 inhibitors.

Structure of the human cohesin inhibitor Wapl

Cohesin, along with positive regulators, establishes sister-chromatid cohesion by forming a ring to circle chromatin. The wings apart-like protein (Wapl) is a key negative regulator of cohesin and forms a complex with precocious dissociation of sisters protein 5 (Pds5) to promote cohesin release from chromatin. Here we report the crystal structure and functional characterization of human Wapl. Wapl contains a flexible, variable N-terminal region (Wapl-N) and a conserved C-terminal domain (Wapl-C) consisting of eight HEAT (Huntingtin, Elongation factor 3, A subunit, and target of rapamycin) repeats. Wapl-C folds into an elongated structure with two lobes. Structure-based mutagenesis maps the functional surface of Wapl-C to two distinct patches (I and II) on the N lobe and a localized patch (III) on the C lobe. Mutating critical patch I residues weaken Wapl binding to cohesin and diminish sister-chromatid resolution and cohesin release from mitotic chromosomes in human cells and Xenopus egg extracts. Surprisingly, patch III on the C lobe does not contribute to Wapl binding to cohesin or its known regulators. Although patch I mutations reduce Wapl binding to intact cohesin, they do not affect Wapl-Pds5 binding to the cohesin subcomplex of sister chromatid cohesion protein 1 (Scc1) and stromal antigen 2 (SA2) in vitro, which is instead mediated by Wapl-N. Thus, Wapl-N forms extensive interactions with Pds5 and Scc1-SA2. Wapl-C interacts with other cohesin subunits and possibly unknown effectors to trigger cohesin release from chromatin.

Structural characterization of V57D and V57P mutants of human cystatin C, an amyloidogenic protein

Wild-type human cystatin C (hCC wt) is a low-molecular-mass protein (120 amino-acid residues, 13,343 Da) that is found in all nucleated cells. Physiologically, it functions as a potent regulator of cysteine protease activity. While the biologically active hCC wt is a monomeric protein, all crystallization efforts to date have resulted in a three-dimensional domain-swapped dimeric structure. In the recently published structure of a mutated hCC, the monomeric fold was preserved by a stabilization of the conformationally constrained loop L1 caused by a single amino-acid substitution: Val57Asn. Additional hCC mutants were obtained in order to elucidate the relationship between the stability of the L1 loop and the propensity of human cystatin C to dimerize. In one mutant Val57 was substituted by an aspartic acid residue, which is favoured in β-turns, and in the second mutant proline, a residue known for broadening turns, was substituted for the same Val57. Here, 2.26 and 3.0 Å resolution crystal structures of the V57D andV57P mutants of hCC are reported and their dimeric architecture is discussed in terms of the stabilization and destabilization effects of the introduced mutations.

Identification of patterns in diffraction intensities affected by radiation exposure

In an X-ray diffraction experiment, the structure of molecules and the crystal lattice changes owing to chemical reactions and physical processes induced by the absorption of X-ray photons. These structural changes alter structure factors, affecting the scaling and merging of data collected at different absorbed doses. Many crystallographic procedures rely on the analysis of consistency between symmetry-equivalent reflections, so failure to account for the drift of their intensities hinders the structure solution and the interpretation of structural results. The building of a conceptual model of radiation-induced changes in macromolecular crystals is the first step in the process of correcting for radiation-induced inconsistencies in diffraction data. Here the complexity of radiation-induced changes in real and reciprocal space is analysed using matrix singular value decomposition applied to multiple complete datasets obtained from single crystals. The model consists of a resolution-dependent decay correction and a uniform-per-unique-reflection term modelling specific radiation-induced changes. This model is typically sufficient to explain radiation-induced effects observed in diffraction intensities. This analysis will guide the parameterization of the model, enabling its use in subsequent crystallographic calculations.

Structure of Escherichia coli RutC, a member of the YjgF family and putative aminoacrylate peracid reductase of the rut operon

RutC is the third enzyme in the Escherichia coli rut pathway of uracil degradation. RutC belongs to the highly conserved YjgF family of proteins. The structure of the RutC protein was determined and refined to 1.95 Å resolution. The crystal belonged to space group P2(1)2(1)2 and contained six molecules in the asymmetric unit. The structure was solved by SAD phasing and was refined to an Rwork of 19.3% (Rfree=21.7%). The final model revealed that this protein has a Bacillus chorismate mutase-like fold and forms a homotrimer with a hydrophobic cavity in the center of the structure and ligand-binding clefts between two subunits. A likely function for RutC is the reduction of peroxy-aminoacrylate to aminoacrylate as a part of a detoxification process.

A multi-faceted analysis of RutD reveals a novel family of α/β hydrolases

The rut pathway of pyrimidine catabolism is a novel pathway that allows pyrimidine bases to serve as the sole nitrogen source in suboptimal temperatures. The rut operon in E. coli evaded detection until 2006, yet consists of seven proteins named RutA, RutB, etc. through RutG. The operon is comprised of a pyrimidine transporter and six enzymes that cleave and further process the uracil ring. Herein, we report the structure of RutD, a member of the α/β hydrolase superfamily, which is proposed to enhance the rate of hydrolysis of aminoacrylate, a toxic side product of uracil degradation, to malonic semialdehyde. Although this reaction will occur spontaneously in water, the toxicity of aminoacrylate necessitates catalysis by RutD for efficient growth with uracil as a nitrogen source. RutD has a novel and conserved arrangement of residues corresponding to the α/β hydrolase active site, where the nucleophile's spatial position occupied by Ser, Cys, or Asp of the canonical catalytic triad is replaced by histidine. We have used a combination of crystallographic structure determination, modeling and bioinformatics, to propose a novel mechanism for this enzyme. This approach also revealed that RutD represents a previously undescribed family within the α/β hydrolases. We compare and contrast RutD with PcaD, which is the closest structural homolog to RutD. PcaD is a 3-oxoadipate-enol-lactonase with a classic arrangement of residues in the active site. We have modeled a substrate in the PcaD active site and proposed a reaction mechanism.

A new generation of crystallographic validation tools for the protein data bank

This report presents the conclusions of the X-ray Validation Task Force of the worldwide Protein Data Bank (PDB). The PDB has expanded massively since current criteria for validation of deposited structures were adopted, allowing a much more sophisticated understanding of all the components of macromolecular crystals. The size of the PDB creates new opportunities to validate structures by comparison with the existing database, and the now-mandatory deposition of structure factors creates new opportunities to validate the underlying diffraction data. These developments highlighted the need for a new assessment of validation criteria. The Task Force recommends that a small set of validation data be presented in an easily understood format, relative to both the full PDB and the applicable resolution class, with greater detail available to interested users. Most importantly, we recommend that referees and editors judging the quality of structural experiments have access to a concise summary of well-established quality indicators.

Structure of the effector-binding domain of the arabinose repressor AraR from Bacillus subtilis

In Bacillus subtilis, the arabinose repressor AraR negatively controls the expression of genes in the metabolic pathway of arabinose-containing polysaccharides. The protein is composed of two domains of different phylogenetic origin and function: an N-terminal DNA-binding domain belonging to the GntR family and a C-terminal effector-binding domain that shows similarity to members of the GalR/LacI family. The crystal structure of the C-terminal effector-binding domain of AraR in complex with the effector L-arabinose has been determined at 2.2 Å resolution. The L-arabinose binding affinity was characterized by isothermal titration calorimetry and differential scanning fluorimetry; the K(d) value was 8.4 ± 0.4 µM. The effect of L-arabinose on the protein oligomeric state was investigated in solution and detailed analysis of the crystal identified a dimer organization which is distinctive from that of other members of the GalR/LacI family.

Crystallization and preliminary X-ray diffraction analysis of Val57 mutants of the amyloidogenic protein human cystatin C

Human cystatin C (hCC) is a low-molecular-mass protein (120 amino-acid residues, 13 343 Da) found in all nucleated cells. Its main physiological role is regulation of the activity of cysteine proteases. Biologically active hCC is a monomeric protein, but all crystallization efforts have resulted in a dimeric domain-swapped structure. Recently, two monomeric structures were reported for cystatin C variants. In one of them stabilization was achieved by abolishing the possibility of domain swapping by the introduction of an additional disulfide bridge connecting the two protein domains (Cys47-Cys69). In the second structure, reported by this group, the monomeric hCC fold was preserved by stabilization of the conformationally constrained loop (L1) by a single-amino-acid substitution (V57N). To further assess the influence of changes in the sequence and properties of loop L1 on the dimerization propensity of cystatin C, two additional hCC mutants were obtained: one with a residue favoured in β-turns (V57D) and another with proline (V57P), a residue that is known to be a structural element that can rigidify but also broaden turns. Here, the expression, purification and crystallization of V57D and V57P variants of recombinant human cystatin C are described. Crystals were grown by the vapour-diffusion method. Several diffraction data sets were collected using a synchrotron source at the Advanced Photon Source, Argonne National Laboratory, Chicago, USA.

Crystal structure of cardiac troponin C regulatory domain in complex with cadmium and deoxycholic acid reveals novel conformation

The amino-terminal regulatory domain of cardiac troponin C (cNTnC) plays an important role as the calcium sensor for the troponin complex. Calcium binding to cNTnC results in conformational changes that trigger a cascade of events that lead to cardiac muscle contraction. The cardiac N-terminal domain of TnC consists of two EF-hand calcium binding motifs, one of which is dysfunctional in binding calcium. Nevertheless, the defunct EF-hand still maintains a role in cNTnC function. For its structural analysis by X-ray crystallography, human cNTnC with the wild-type primary sequence was crystallized under a novel crystallization condition. The crystal structure was solved by the single-wavelength anomalous dispersion method and refined to 2.2 Å resolution. The structure displays several novel features. Firstly, both EF-hand motifs coordinate cadmium ions derived from the crystallization milieu. Secondly, the ion coordination in the defunct EF-hand motif accompanies unusual changes in the protein conformation. Thirdly, deoxycholic acid, also derived from the crystallization milieu, is bound in the central hydrophobic cavity. This is reminiscent of the interactions observed for cardiac calcium sensitizer drugs that bind to the same core region and maintain the "open" conformational state of calcium-bound cNTnC. The cadmium ion coordination in the defunct EF-hand indicates that this vestigial calcium binding site retains the structural and functional elements that allow it to coordinate a cadmium ion. However, it is a result of, or concomitant with, large and unusual structural changes in cNTnC.