Optimization of multiple-sequence alignment based on multiple-structure alignment

Routinely used multiple-sequence alignment methods use only sequence information. Consequently, they may produce inaccurate alignments. Multiple-structure alignment methods, on the other hand, optimize structural alignment by ignoring sequence information. Here, we present an optimization method that unifies sequence and structure information. The alignment score is based on standard amino acid substitution probabilities combined with newly computed three-dimensional structure alignment probabilities. The advantage of our alignment scheme is in its ability to produce more accurate multiple alignments. We demonstrate the usefulness of the method in three applications: 1) computing more accurate multiple-sequence alignments, 2) analyzing protein conformational changes, and 3) computation of amino acid structure-sequence conservation with application to protein-protein docking prediction. The method is available at http://bioinfo3d.cs.tau.ac.il/staccato/.     

Human mitotic chromosomes consist predominantly of irregularly folded nucleosome fibres without a 30-nm chromatin structure

How a long strand of genomic DNA is compacted into a mitotic chromosome remains one of the basic questions in biology. The nucleosome fibre, in which DNA is wrapped around core histones, has long been assumed to be folded into a 30-nm chromatin fibre and further hierarchical regular structures to form mitotic chromosomes, although the actual existence of these regular structures is controversial. Here, we show that human mitotic HeLa chromosomes are mainly composed of irregularly folded nucleosome fibres rather than 30-nm chromatin fibres. Our comprehensive and quantitative study using cryo-electron microscopy and synchrotron X-ray scattering resolved the long-standing contradictions regarding the existence of 30-nm chromatin structures and detected no regular structure >11 nm. Our finding suggests that the mitotic chromosome consists of irregularly arranged nucleosome fibres, with a fractal nature, which permits a more dynamic and flexible genome organization than would be allowed by static regular structures.     

1.3-A resolution structure of human glutathione S-transferase with S-hexyl glutathione bound reveals possible extended ligandin binding site

Cytosolic glutathione S-transferases (GSTs) play a critical role in xenobiotic binding and metabolism, as well as in modulation of oxidative stress. Here, the high-resolution X-ray crystal structures of homodimeric human GSTA1-1 in the apo form and in complex with S-hexyl glutathione (two data sets) are reported at 1.8, 1.5, and 1.3A respectively. At this level of resolution, distinct conformations of the alkyl chain of S-hexyl glutathione are observed, reflecting the nonspecific nature of the hydrophobic substrate binding site (H-site). Also, an extensive network of ordered water, including 75 discrete solvent molecules, traverses the open subunit-subunit interface and connects the glutathione binding sites in each subunit. In the highest-resolution structure, three glycerol moieties lie within this network and directly connect the amino termini of the glutathione molecules. A search for ligand binding sites with the docking program Molecular Operating Environment identified the ordered water network binding site, lined mainly with hydrophobic residues, suggesting an extended ligand binding surface for nonsubstrate ligands, the so-called ligandin site. Finally, detailed comparison of the structures reported here with previously published X-ray structures reveal a possible reaction coordinate for ligand-dependent conformational changes in the active site and the C-terminus.     

Improving computational protein design by using structure‐derived sequence profile

Designing a protein sequence that will fold into a predefined structure is of both practical and fundamental interest. Many successful, computational designs in the last decade resulted from improved understanding of hydrophobic and polar interactions between side chains of amino acid residues in stabilizing protein tertiary structures. However, the coupling between main‐chain backbone structure and local sequence has yet to be fully addressed. Here, we attempt to account for such coupling by using a sequence profile derived from the sequences of five residue fragments in a fragment library that are structurally matched to the five‐residue segments contained in a target structure. We further introduced a term to reduce low complexity regions of designed sequences. These two terms together with optimized reference states for amino‐acid residues were implemented in the RosettaDesign program. The new method, called RosettaDesign‐SR, makes a 12% increase (from 34 to 46%) in fraction of proteins whose designed sequences are more than 35% identical to wild‐type sequences. Meanwhile, it reduces 8% (from 22% to 14%) to the number of designed sequences that are not homologous to any known protein sequences according to psi‐blast. More importantly, the sequences designed by RosettaDesign‐SR have 2–3% more polar residues at the surface and core regions of proteins and these surface and core polar residues have about 4% higher sequence identity to wild‐type sequences than by RosettaDesign. Thus, the proteins designed by RosettaDesign‐SR should be less likely to aggregate and more likely to have unique structures due to more specific polar interactions. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Angle-distance image matching techniques for protein structure comparison

Proteins that contain similar structural elements often have analogous functions regardless of the degree of sequence similarity or structure connectivity in space. In general, protein structure comparison (PSC) provides a straightforward methodology for biologists to determine critical aspects of structure and function. Here, we developed a novel PSC technique based on angle-distance image (A-D image) transformation and matching, which is independent of sequence similarity and connectivity of secondary structure elements (SSEs). An A-D image is constructed by utilizing protein secondary structure information. According to various types of SSEs, the mutual SSE pairs of the query protein are classified into three different types of sub-images. Subsequently, corresponding sub-images between query and target protein structures are compared using modified cross-correlation approaches to identify the similarity of various patterns. Structural relationships among proteins are displayed by hierarchical clustering trees, which facilitate the establishment of the evolutionary relationships between structure and function of various proteins.Four standard testing datasets and one newly created dataset were used to evaluate the proposed method. The results demonstrate that proteins from these five datasets can be categorized in conformity with their spatial distribution of SSEs. Moreover, for proteins with low sequence identity that share high structure similarity, the proposed algorithms are an efficient and effective method for structural comparison.     

TASSER-based refinement of NMR structures

The TASSER structure prediction algorithm is employed to investigate whether NMR structures can be moved closer to their corresponding X-ray counterparts by automatic refinement procedures. The benchmark protein dataset includes 61 nonhomologous proteins whose structures have been determined by both NMR and X-ray experiments. Interestingly, by starting from NMR structures, the majority (79%) of TASSER refined models show a structural shift toward their X-ray structures. On average, the TASSER refined models have a root-mean-square-deviation (RMSD) from the X-ray structure of 1.785 A (1.556 A) over the entire chain (aligned region), while the average RMSD between NMR and X-ray structures (RMSD(NMR_X-ray)) is 2.080 A (1.731 A). For all proteins having a RMSD(NMR_X-ray) >2 A, the TASSER refined structures show consistent improvement. However, for the 34 proteins with a RMSD(NMR_X-ray) <2 A, there are only 21 cases (60%) where the TASSER model is closer to the X-ray structure than NMR, which may be due to the inherent resolution of TASSER. We also compare the TASSER models with 12 NMR models in the RECOORD database that have been recalculated recently by Nederveen et al. from original NMR restraints using the newest molecular dynamics tools. In 8 of 12 cases, TASSER models show a smaller RMSD to X-ray structures; in 3 of 12 cases, where RMSD(NMR_X-ray) <1 A, RECOORD does better than TASSER. These results suggest that TASSER can be a useful tool to improve the quality of NMR structures.     

Crystal structure of GlcAT-S, a human glucuronyltransferase, involved in the biosynthesis of the HNK-1 carbohydrate epitope

The HNK-1 carbohydrate epitope is found in various neural cell adhesion molecules. Two glucuronyltransferases (GlcAT-P and GlcAT-S) are involved in the biosynthesis of HNK-1 carbohydrate. Our previous study on the crystal structure of GlcAT-P revealed the reaction and substrate recognition mechanisms of this enzyme. Comparative analyses of the enzymatic activities of GlcAT-S and GlcAT-P showed that there are notable differences in the acceptor substrate specificities of these enzymes. To elucidate differences between their specificities, we now solved the crystal structure of GlcAT-S. Residues interacting with UDP molecule, which is a part of the donor substrate, are highly conserved between GlcAT-P and GlcAT-S. On the other hand, there are some differences between these proteins in the manner they recognize their respective acceptor substrates. Phe245, one of the most important GlcAT-P residues for the recognition of acceptors, is a tryptophan in GlcAT-S. In addition, Val320, which is located on the C-terminal long loop of the neighboring molecule in the dimer and critical in the recognition of the acceptor sugar molecule by the GlcAT-P dimer, is an alanine in GlcAT-S. These differences play key roles in establishing the distinct specificity for the acceptor substrate by GlcAT-S, which is further supported by site-directed mutagenesis of GlcAT-S and a computer-aided model building of GlcAT-S/substrate complexes.     

14.2% Efficiency Dye‐Sensitized Solar Cells by Co‐sensitizing Novel Thieno[3,2‐b]indole‐Based Organic Dyes with a Promising Porphyrin Sensitizer

A new series of 4‐hexyl‐4H‐thieno[3,2‐b]indole (HxTI) based organic chromophores is developed by structural engineering of the electron donor (D) group in the D–HxTI–benzothiadiazole‐phenyl‐acceptor platform with different fluorenyl moieties, such as unsubstituted fluorenyl (SGT‐146) and hexyloxy (SGT‐147), decyloxy (SGT‐148) and hexyloxy‐phenyl substituted (SGT‐149) fluorenyl moieties. In comparison to a reference dye SGT‐137 with a biphenyl‐based donor, the effects of the donating ability and bulkiness of the fluorenyl based donor in this D–π–A‐structured platform on molecular properties and photovoltaic performance are investigated to establish the structure–property relationship. The photovoltaic performance of dye‐sensitized solar cells (DSSCs) is improved according to the bulkiness of the donor groups. As a result, the DSSCs based on SGT‐149 show high power conversion efficiencies (PCEs) of 11.7% and 10.0% with a [Co(bpy)₃]^2+/3+ (bpy = 2,2′‐bipyridine) and an I^?/I₃^? redox electrolyte, respectively. Notably, the co‐sensitization of SGT‐149 with a SGT‐021 porphyrin dye by utilizing a simple “cocktail” method, exhibit state‐of‐the‐art PCEs of 14.2% and 11.6% with a [Co(bpy)₃]^2+/3+ and an I^?/I₃^? redox electrolyte, respectively.     

Determination of protein oligomeric structure from small‐angle X‐ray scattering

Small‐angle X‐ray scattering (SAXS) is useful for determining the oligomeric states and quaternary structures of proteins in solution. The average molecular mass in solution can be calculated directly from a single SAXS curve collected on an arbitrary scale from a sample of unknown protein concentration without the need for beamline calibration or protein standards. The quaternary structure in solution can be deduced by comparing the experimental SAXS curve to theoretical curves calculated from proposed models of the oligomer. This approach is especially robust when the crystal structure of the target protein is known, and the candidate oligomer models are derived from the crystal lattice. When SAXS data are obtained at multiple protein concentrations, this analysis can provide insight into dynamic self‐association equilibria. Herein, we summarize the computational methods that are used to determine protein molecular mass and quaternary structure from SAXS data. These methods are organized into a workflow and demonstrated with four case studies using experimental SAXS data from the published literature.     

Cooperative folding mechanism of a beta-hairpin peptide studied by a multicanonical replica-exchange molecular dynamics simulation

G-peptide is a 16-residue peptide of the C-terminal end of streptococcal protein G B1 domain, which is known to fold into a specific beta-hairpin within 6 micros. Here, we study molecular mechanism on the stability and folding of G-peptide by performing a multicanonical replica-exchange (MUCAREM) molecular dynamics simulation with explicit solvent. Unlike the preceding simulations of the same peptide, the simulation was started from an unfolded conformation without any experimental information on the native conformation. In the 278-ns trajectory, we observed three independent folding events. Thus MUCAREM can be estimated to accelerate the folding reaction more than 60 times than the conventional molecular dynamics simulations. The free-energy landscape of the peptide at room temperature shows that there are three essential subevents in the folding pathway to construct the native-like beta-hairpin conformation: (i) a hydrophobic collapse of the peptide occurs with the side-chain contacts between Tyr45 and Phe52, (ii) then, the native-like turn is formed accompanying with the hydrogen-bonded network around the turn region, and (iii) finally, the rest of the backbone hydrogen bonds are formed. A number of stable native hydrogen bonds are formed cooperatively during the second stage, suggesting the importance of the formation of the specific turn structure. This is also supported by the accumulation of the nonnative conformations only with the hydrophobic cluster around Tyr45 and Phe52. These simulation results are consistent with high phi-values of the turn region observed by experiment.     

Covid‐19.bioreproducibility.org: A web resource for SARS‐CoV‐2‐related structural models

The COVID‐19 pandemic has triggered numerous scientific activities aimed at understanding the SARS‐CoV‐2 virus and ultimately developing treatments. Structural biologists have already determined hundreds of experimental X‐ray, cryo‐EM, and NMR structures of proteins and nucleic acids related to this coronavirus, and this number is still growing. To help biomedical researchers, who may not necessarily be experts in structural biology, navigate through the flood of structural models, we have created an online resource, covid19.bioreproducibility.org, that aggregates expert‐verified information about SARS‐CoV‐2‐related macromolecular models. In this article, we describe this web resource along with the suite of tools and methodologies used for assessing the structures presented therein.     

Compact reduced thioredoxin structure from the thermophilic bacteria Thermus thermophilus

The X-ray crystallographic structure of a thioredoxin from Thermus thermophilus was solved to 1.8 A resolution by molecular replacement. The crystals' space group was C2 with cell dimensions of a = 40.91, b = 95.44, c = 56.68 A, beta =91.41 degrees, with two molecules in the asymmetric unit. Unlike the reported thioredoxin structures, the biological unit of T. thermophilus thioredoxin is a dimer both in solution and in the crystal. The fold conforms to the "thioredoxin fold" that is common over a class of nine protein families including thioredoxin; however, the folded portion of this protein is much more compact than other thioredoxins previously solved by X-ray crystallography being reduced by one alpha-helix and one beta-strand. As with other thioredoxins, the active site is highly conserved even though the variation in sequence can be quite large. The T. thermophilus thioredoxin has some variability at the active site, especially compared with previously solved structures from bacterial sources.     

Cooperativity network of Trp‐cage miniproteins: probing salt‐bridges

Trp‐cage miniprotein was used to investigate the role of a salt‐bridge (Asp⁹–Arg¹⁶) in protein formation, by mutating residues at both sides, we mapped its contribution to overall stability and its role in folding mechanism. We found that both of the above side‐chains are also part of a dense interaction network composed of electrostatic, H‐bonding, hydrophobic, etc. components. To elucidate the fold stabilizing effects, we compared and contrasted electronic circular dichroism and NMR data of miniproteins equipped with a salt‐bridge with those of the salt‐bridge deleted mutants. Data were acquired both in neutral and in acidic aqueous solutions to decipher the pH dependency of both fully and partially charged partners. Our results indicate that the folding of Trp‐cage miniproteins is more complex than a simple two‐state process as we detected an intermediate state that differs significantly from the native fold. The intermediate formation is related to the salt‐bridge stabilization; in the miniprotein variants equipped with salt‐bridge the population of the intermediate state at acidic pH is significantly higher than it is for the salt‐bridge deleted mutants. In this molecular framework Arg¹⁶ stabilizes more than Asp⁹ does, because of its higher degree of 3D‐fold cooperation. In conclusion, the Xxx$^{9} \leftrightarrow$ Yyy¹⁶ salt‐bridge is not an isolated entity of this fold; rather it is an integrated part of a complex interaction network. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

A tetrapeptide-based method for polyproline II-type secondary structure prediction

We describe a new method for polyproline II-type (PPII) secondary structure prediction based on tetrapeptide conformation properties using data obtained from all globular proteins in the Protein Data Bank (PDB). This is the first method for PPII prediction with a relatively high level of accuracy (approximately 60%). Our method uses only frequencies of different conformations among oligopeptides without any additional parameters. We also attempted to predict alpha-helices and beta-strands using the same approach. We find that the application of our method reveals interrelation between sequence and structure even for very short oligopeptides (tetrapeptides).     

Crystal structure of checkpoint kinase 2 in complex with NSC 109555, a potent and selective inhibitor

Checkpoint kinase 2 (Chk2), a ser/thr kinase involved in the ATM‐Chk2 checkpoint pathway, is activated by genomic instability and DNA damage and results in either arrest of the cell cycle to allow DNA repair to occur or apoptosis if the DNA damage is severe. Drugs that specifically target Chk2 could be beneficial when administered in combination with current DNA‐damaging agents used in cancer therapy. Recently, a novel inhibitor of Chk2, NSC 109555, was identified that exhibited high potency (IC₅₀ = 240 nM) and selectivity. This compound represents a new chemotype and lead for the development of novel Chk2 inhibitors that could be used as therapeutic agents for the treatment of cancer. To facilitate the discovery of new analogs of NSC 109555 with even greater potency and selectivity, we have solved the crystal structure of this inhibitor in complex with the catalytic domain of Chk2. The structure confirms that the compound is an ATP‐competitive inhibitor, as the electron density clearly reveals that it occupies the ATP‐binding pocket. However, the mode of inhibition differs from that of the previously studied structure of Chk2 in complex with debromohymenialdisine, a compound that inhibits both Chk1 and Chk2. A unique hydrophobic pocket in Chk2, located very close to the bound inhibitor, presents an opportunity for the rational design of compounds with higher binding affinity and greater selectivity.     

CD studies on films of amyloid proteins and polypeptides: quantitative g-factor analysis indicates a common folding motif

Irrespective of the constituent protein, all amyloid fibrils show similar morphology in the electron microscope and x-ray diffraction patterns characteristic of a "cross-beta" structure, with extended beta-strands perpendicular to the fibril's long axis. Little is known about the amount or type of this structure. I have measured CD spectra of films formed from a number of amyloid proteins and polypeptides, and estimated their contents of extended secondary structure, by analysis of their g-factor spectra, the ratio of the CD and absorbance signals (P. McPhie, Analytical Biochemistry, 2001, Vol. 293, pp. 109-119). Amyloid films of Abeta-(1-40) peptide, beta-2-microglobulin, insulin, and three homopolypeptides show very intense CD spectra, compatible with the presence of a beta-helix-like structure, arranged in a common framework in the fibrils. The extent of this structure was estimated as 45-80% in the protein fibrils and 30-80% in the polypeptide fibrils.     

Toward a structome of Acinetobacter baumannii drug targets

Acinetobacter baumannii is well known for causing hospital‐associated infections due in part to its intrinsic antibiotic resistance as well as its ability to remain viable on surfaces and resist cleaning agents. In a previous publication, A. baumannii strain AB5075 was studied by transposon mutagenesis and 438 essential gene candidates for growth on rich‐medium were identified. The Seattle Structural Genomics Center for Infectious Disease entered 342 of these candidate essential genes into our pipeline for structure determination, in which 306 were successfully cloned into expression vectors, 192 were detectably expressed, 165 screened as soluble, 121 were purified, 52 crystalized, 30 provided diffraction data, and 29 structures were deposited in the Protein Data Bank. Here, we report these structures, compare them with human orthologs where applicable, and discuss their potential as drug targets for antibiotic development against A. baumannii.