Search

The initiation of transcription in Escherichia coli ( E. coli ) is facilitated by promoter specificity factors, also known as σ factors, which may bind a promoter only as part of a complex with RNA polymerase (RNAP). By performing in vitro cross-linking mass spectrometry (CL-MS) of apo-σ 70 , we reveal structural features suggesting a compact conformation compared to the known RNAP-bound extended conformation. Then, we validate the existence of the compact conformation using in vivo CL-MS by identifying cross-links similar to those found in vitro , which deviate from the extended conformation only during the stationary phase of bacterial growth. Conclusively, we provide information in support of a compact conformation of apo-σ 70 that exists in live cells, which might represent a transcriptionally inactive form that can be activated upon binding to RNAP., Competing Interests: All authors declare no competing interests., (© 2024 The Author(s).)

Sepsis is a significant public health challenge. The immune system underlies the pathogenesis of the disease. The liver is both an active player and a target organ in sepsis. Targeting the gut immune system using low-dose colchicine is an attractive method for alleviating systemic inflammation in sepsis without inducing immunosuppression. The present study aimed to determine the use of low-dose colchicine in LPS-induced sepsis in mice. C67B mice were injected intraperitoneal with LPS to induce sepsis. The treatment group received 0.02 mg/kg colchicine daily by gavage. Short and extended models were performed, lasting 3 and 5 days, respectively. We followed the mice for biochemical markers of end-organ injury, blood counts, cytokine levels, and liver pathology and conducted proteomic studies on liver samples. Targeting the gut immune system using low-dose colchicine improved mice's well-being measured by the murine sepsis score. Treatment alleviated the liver injury in septic mice, manifested by a significant decrease in their liver enzyme levels, including ALT, AST, and LDH. Treatment exerted a trend to reduce creatinine levels. Low-dose colchicine improved liver pathology, reduced inflammation, and reduced the pro-inflammatory cytokine TNFα and IL1-β levels. A liver proteomic analysis revealed low-dose colchicine down-regulated sepsis-related proteins, alpha-1 antitrypsin, and serine dehydratase. Targeting the gut immune system using low-dose colchicine attenuated liver injury in LPS-induced sepsis, reducing the pro-inflammatory cytokine levels. Low-dose colchicine provides a safe method for immunomodulation for multiple inflammatory disorders., (© 2023. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Parkinson disease is associated with the aggregation of the protein α-synuclein. While α-synuclein can exist in multiple oligomeric states, the dimer has been a subject of extensive debates. Here, using an array of biophysical approaches, we demonstrate that α-synuclein in vitro exhibits primarily a monomer-dimer equilibrium in nanomolar concentrations and up to a few micromolars. We then use spatial information from hetero-isotopic cross-linking mass spectrometry experiments as restrains in discrete molecular dynamics simulations to obtain the ensemble structure of dimeric species. Out of eight structural sub-populations of dimers, we identify one that is compact, stable, abundant, and exhibits partially exposed β-sheet structures. This compact dimer is the only one where the hydroxyls of tyrosine 39 are in proximity that may promote dityrosine covalent linkage upon hydroxyl radicalization, which is implicated in α-synuclein amyloid fibrils. We propose that this α-synuclein dimer features etiological relevance to Parkinson disease., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Background: We have previously demonstrated that double homeobox 4 centromeric (DUX4C) encoded for a functional DUX4c protein upregulated in dystrophic skeletal muscles. Based on gain- and loss-of-function studies we have proposed DUX4c involvement in muscle regeneration. Here, we provide further evidence for such a role in skeletal muscles from patients affected with facioscapulohumeral muscular dystrophy (FSHD)., Methods: DUX4c was studied at RNA and protein levels in FSHD muscle cell cultures and biopsies. Its protein partners were co-purified and identified by mass spectrometry. Endogenous DUX4c was detected in FSHD muscle sections with either its partners or regeneration markers using co-immunofluorescence or in situ proximity ligation assay., Results: We identified new alternatively spliced DUX4C transcripts and confirmed DUX4c immunodetection in rare FSHD muscle cells in primary culture. DUX4c was detected in nuclei, cytoplasm or at cell-cell contacts between myocytes and interacted sporadically with specific RNA-binding proteins involved, a.o., in muscle differentiation, repair, and mass maintenance. In FSHD muscle sections, DUX4c was found in fibers with unusual shape or central/delocalized nuclei (a regeneration feature) staining for developmental myosin heavy chain, MYOD or presenting intense desmin labeling. Some couples of myocytes/fibers locally exhibited peripheral DUX4c-positive areas that were very close to each other, but in distinct cells. MYOD or intense desmin staining at these locations suggested an imminent muscle cell fusion. We further demonstrated DUX4c interaction with its major protein partner, C1qBP, inside myocytes/myofibers that presented features of regeneration. On adjacent muscle sections, we could unexpectedly detect DUX4 (the FSHD causal protein) and its interaction with C1qBP in fusing myocytes/fibers., Conclusions: DUX4c upregulation in FSHD muscles suggests it contributes not only to the pathology but also, based on its protein partners and specific markers, to attempts at muscle regeneration. The presence of both DUX4 and DUX4c in regenerating FSHD muscle cells suggests DUX4 could compete with normal DUX4c functions, thus explaining why skeletal muscle is particularly sensitive to DUX4 toxicity. Caution should be exerted with therapeutic agents aiming for DUX4 suppression because they might also repress the highly similar DUX4c and interfere with its physiological role., (© 2023. The Author(s).)

Mitochondria play central roles in maintaining cellular metabolic homeostasis, cell survival and cell death, and generate most of the cell's energy. Mitochondria maintain their homeostasis by dynamic (fission and fusion) and quality control mechanisms, including mitophagy, the removal of damaged mitochondria that is mediated mainly by the Pink1/Parkin pathway. Pink1 is a serine/threonine kinase which regulates mitochondrial function, hitherto many molecular mechanisms underlying Pink1 activity in mitochondrial homeostasis and cell fate remain unknown. Peptides are vital biological mediators that demonstrate remarkable potency, selectivity, and low toxicity, yet they have two major limitations, low oral bioavailability and poor stability. Herein, we rationally designed a linear peptide that targets Pink1 and, using straightforward chemistry, we developed molecular probes with drug-like properties to further characterize Pink1. Initially, we conjugated a cell-penetrating peptide and a cross-linker to map Pink1's 3D structure and its interaction sites. Next, we conjugated a fluorescent dye for cell-imaging. Finally, we developed cyclic peptides with improved stability and binding affinity. Overall, we present a facile approach to converting a non-permeable linear peptide into a research tool possessing important properties for therapeutics. This is a general approach using straightforward chemistry that can be tailored for various applications by numerous laboratories.

α-Synuclein plays an important role in synaptic functions by interacting with synaptic vesicle membrane, while its oligomers and fibrils are associated with several neurodegenerative diseases. The specific monomer structures that promote its membrane binding and self-association remain elusive due to its transient nature as an intrinsically disordered protein. Here, we use inter-dye distance distributions from bulk time-resolved Förster resonance energy transfer as restraints in discrete molecular dynamics simulations to map the conformational space of the α-synuclein monomer. We further confirm the generated conformational ensemble in orthogonal experiments utilizing far-UV circular dichroism and cross-linking mass spectrometry. Single-molecule protein-induced fluorescence enhancement measurements show that within this conformational ensemble, some of the conformations of α-synuclein are surprisingly stable, exhibiting conformational transitions slower than milliseconds. Our comprehensive analysis of the conformational ensemble reveals essential structural properties and potential conformations that promote its various functions in membrane interaction or oligomer and fibril formation., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Atomic structures of several proteins from the coronavirus family are still partial or unavailable. A possible reason for this gap is the instability of these proteins outside of the cellular context, thereby prompting the use of in-cell approaches. In situ cross-linking and mass spectrometry (in situ CLMS) can provide information on the structures of such proteins as they occur in the intact cell. Here, we applied targeted in situ CLMS to structurally probe Nsp1, Nsp2, and nucleocapsid (N) proteins from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and obtained cross-link sets with an average density of one cross-link per 20 residues. We then employed integrative modeling that computationally combined the cross-linking data with domain structures to determine full-length atomic models. For the Nsp2, the cross-links report on a complex topology with long-range interactions. Integrative modeling with structural prediction of individual domains by the AlphaFold2 system allowed us to generate a single consistent all-atom model of the full-length Nsp2. The model reveals three putative metal binding sites and suggests a role for Nsp2 in zinc regulation within the replication-transcription complex. For the N protein, we identified multiple intra- and interdomain cross-links. Our integrative model of the N dimer demonstrates that it can accommodate three single RNA strands simultaneously, both stereochemically and electrostatically. For the Nsp1, cross-links with the 40S ribosome were highly consistent with recent cryogenic electron microscopy structures. These results highlight the importance of cellular context for the structural probing of recalcitrant proteins and demonstrate the effectiveness of targeted in situ CLMS and integrative modeling., Competing Interests: The authors declare no competing interest., (Copyright   2021 the Author(s). Published by PNAS.)

Cross-linking of living cells followed by mass spectrometry identification of cross-linked peptides (in situ CLMS) is an emerging technology to study protein structures in their native environment. One of the inherent difficulties of this technology is the high complexity of the samples following cell lysis. Currently, this difficulty largely limits the identification of cross-links to the more abundant proteins in the cell. Here, we describe a targeted approach in which an antibody is used to purify a specific protein-of-interest out of the cell lysate. Mass spectrometry analysis of the protein material that binds to the antibody can then identify considerably more cross-links on the target protein. By using an antibody against the CCT chaperonin, we identified over 200 cross-links that provide in situ evidence for the subunit arrangement of the CCT particle and its interactions with prefoldin. Similar targeting with an antibody against tubulin provided in situ evidence for the structure of the microtubule. Finally, the approach was also successful in identifying cross-links within a protein that expresses at a low level. These results demonstrate the general utility of antibody-based sample simplification for in situ CLMS and greatly expand the scope of protein systems that are amenable to in situ structural studies.

Development of new reagents for protein cross-linking is constantly ongoing. The chemical formulas for the linker adducts formed by these reagents are usually deduced from expert knowledge and then validated by mass spectrometry. Clearly, it would be more rigorous to infer the chemical compositions of the adducts directly from the data without any prior assumptions on their chemistries. Unfortunately, the analysis tools that are currently available to detect chemical modifications on linear peptides are not applicable to the case of two cross-linked peptides. Here, we show that an adaptation of the open search strategy that works on linear peptides can be used to characterize cross-link modifications in pairs of peptides. We benchmark our approach by correctly inferring the linker masses of two well-known reagents, DSS and formaldehyde, to accuracies of a few parts per million. We then investigate the cross-linking chemistries of two poorly characterized reagents: EMCS and glutaraldehyde. In the case of EMCS, we find that the expected cross-linking chemistry is accompanied by a competing chemistry that targets other amino acid types. In the case of glutaraldehyde, we find that the chemical formula of the dominant linker is C 5 H 4 , which indicates a ringed aromatic structure. These results demonstrate how, with very little effort, our approach can yield nontrivial insights to better characterize new cross-linkers.

Whole-cell cross-linking coupled to mass spectrometry is one of the few tools that can probe protein-protein interactions in intact cells. A very attractive reagent for this purpose is formaldehyde, a small molecule which is known to rapidly penetrate into all cellular compartments and to preserve the protein structure. In light of these benefits, it is surprising that identification of formaldehyde cross-links by mass spectrometry has so far been unsuccessful. Here we report mass spectrometry data that reveal formaldehyde cross-links to be the dimerization product of two formaldehyde-induced amino acid modifications. By integrating the revised mechanism into a customized search algorithm, we identify hundreds of cross-links from in situ formaldehyde fixation of human cells. Interestingly, many of the cross-links could not be mapped onto known atomic structures, and thus provide new structural insights. These findings enhance the use of formaldehyde cross-linking and mass spectrometry for structural studies.

Centromeric nucleosomes are at the interface of the chromosome and the kinetochore that connects to spindle microtubules in mitosis. The core centromeric nucleosome complex (CCNC) harbors the histone H3 variant, CENP-A, and its binding proteins, CENP-C (through its central domain; CD) and CENP-N (through its N-terminal domain; NT). CENP-C can engage nucleosomes through two domains: the CD and the CENP-C motif (CM). CENP-C CD is part of the CCNC by virtue of its high specificity for CENP-A nucleosomes and ability to stabilize CENP-A at the centromere. CENP-C CM is thought to engage a neighboring nucleosome, either one containing conventional H3 or CENP-A, and a crystal structure of a nucleosome complex containing two copies of CENP-C CM was reported. Recent structures containing a single copy of CENP-N NT bound to the CENP-A nucleosome in the absence of CENP-C were reported. Here, we find that one copy of CENP-N is lost for every two copies of CENP-C on centromeric chromatin just prior to kinetochore formation. We present the structures of symmetric and asymmetric forms of the CCNC that vary in CENP-N stoichiometry. Our structures explain how the central domain of CENP-C achieves its high specificity for CENP-A nucleosomes and how CENP-C and CENP-N sandwich the histone H4 tail. The natural centromeric DNA path in our structures corresponds to symmetric surfaces for CCNC assembly, deviating from what is observed in prior structures using artificial sequences. At mitosis, we propose that CCNC asymmetry accommodates its asymmetric connections at the chromosome/kinetochore interface. VIDEO ABSTRACT., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Integrative structure modeling computationally combines data from multiple sources of information with the aim of obtaining structural insights that are not revealed by any single approach alone. In the first part of this review, we survey the commonly used sources of structural information and the computational aspects of model building. Throughout the past decade, integrative modeling was applied to various biological systems, with a focus on large protein complexes. Recent progress in the field of cryo-electron microscopy (cryo-EM) has resolved many of these complexes to near-atomic resolution. In the second part of this review, we compare a range of published integrative models with their higher-resolution counterparts with the aim of critically assessing their accuracy. This comparison gives a favorable view of integrative modeling and demonstrates its ability to yield accurate and informative results. We discuss possible roles of integrative modeling in the new era of cryo-EM and highlight future challenges and directions.

The number of publications in the field of chemical cross-linking combined with mass spectrometry (XL-MS) to derive constraints for protein three-dimensional structure modeling and to probe protein-protein interactions has increased during the last years. As the technique is now becoming routine for in vitro and in vivo applications in proteomics and structural biology there is a pressing need to define protocols as well as data analysis and reporting formats. Such consensus formats should become accepted in the field and be shown to lead to reproducible results. This first, community-based harmonization study on XL-MS is based on the results of 32 groups participating worldwide. The aim of this paper is to summarize the status quo of XL-MS and to compare and evaluate existing cross-linking strategies. Our study therefore builds the framework for establishing best practice guidelines to conduct cross-linking experiments, perform data analysis, and define reporting formats with the ultimate goal of assisting scientists to generate accurate and reproducible XL-MS results.

At resolutions worse than 3.5 Å, the electron density is weak or nonexistent at the locations of the side chains. Consequently, the assignment of the protein sequences to their correct positions along the backbone is a difficult problem. In this work, we propose a fully automated computational approach to assign sequence at low resolution. It is based on our surprising observation that standard reciprocal-space indicators, such as the initial unrefined R value, are sensitive enough to detect an erroneous sequence assignment of even a single backbone position. Our approach correctly determines the amino acid type for 15%, 13%, and 9% of the backbone positions in crystallographic datasets with resolutions of 4.0 Å, 4.5 Å, and 5.0 Å, respectively. We implement these findings in an application for threading a sequence onto a backbone structure. For the three resolution ranges, the application threads 83%, 81%, and 64% of the sequences exactly as in the deposited PDB structures., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Cross-linking and mass spectrometry is used more and more for the structural analysis of large proteins and protein complexes. Although essentially a low-resolution method, it avoids the main drawbacks of established structural techniques. Particularly, it is largely insensitive to the inherent flexibility of the studied complexes and is applied under native conditions. It is also applicable to nearly every structural system. Therefore, cross-linking and mass spectrometry is the method of choice for elucidating the general architecture of protein complexes. Advances in instrumentation, techniques, and software now allow every lab that is working with proteins to apply the approach without much difficulty. The most specialized step in the workflow, the mass spectrometry measurement, can be done in most facilities that are performing standard proteomics. We detail here a step-by-step protocol of how to successfully apply the approach in collaboration with the mass spectrometry facility in your institution.

The DNA damage response is an essential process for the survival of living cells. In a subset of stress-responsive genes in humans, Elongin controls transcription in response to multiple stimuli, such as DNA damage, oxidative stress, and heat shock. Yeast Elongin (Ela1-Elc1), along with Def1, is known to facilitate ubiquitylation and degradation of RNA polymerase II (pol II) in response to multiple stimuli, yet transcription activity has not been examined. We have found that Def1 copurifies from yeast whole-cell extract with TFIIH, the largest general transcription factor required for transcription initiation and nucleotide excision repair. The addition of recombinant Def1 and Ela1-Elc1 enhanced transcription initiation in an in vitro reconstituted system including pol II, the general transcription factors, and TFIIS. Def1 also enhanced transcription restart from TFIIS-induced cleavage in a pol II transcribing complex. In the Δdef1 strain, heat shock genes were misregulated, indicating that Def1 is required for induction of some stress-responsive genes in yeast. Taken together, our results extend the understanding of the molecular mechanism of transcription regulation on cellular stress and reveal functional similarities to the mammalian system., Competing Interests: The authors declare no conflict of interest.

The bacterial flagellar apparatus, which involves ∼40 different proteins, has been a model system for understanding motility and chemotaxis. The bacterial flagellar filament, largely composed of a single protein, flagellin, has been a model for understanding protein assembly. This system has no homology to the eukaryotic flagellum, in which the filament alone, composed of a microtubule-based axoneme, contains more than 400 different proteins. The archaeal flagellar system is simpler still, in some cases having ∼13 different proteins with a single flagellar filament protein. The archaeal flagellar system has no homology to the bacterial one and must have arisen by convergent evolution. However, it has been understood that the N-terminal domain of the archaeal flagellin is a homolog of the N-terminal domain of bacterial type IV pilin, showing once again how proteins can be repurposed in evolution for different functions. Using cryo-EM, we have been able to generate a nearly complete atomic model for a flagellar-like filament of the archaeon Ignicoccus hospitalis from a reconstruction at ∼4-Å resolution. We can now show that the archaeal flagellar filament contains a β-sandwich, previously seen in the FlaF protein that forms the anchor for the archaeal flagellar filament. In contrast to the bacterial flagellar filament, where the outer globular domains make no contact with each other and are not necessary for either assembly or motility, the archaeal flagellin outer domains make extensive contacts with each other that largely determine the interesting mechanical properties of these filaments, allowing these filaments to flex., Competing Interests: The authors declare no conflict of interest.

The structure of a 33-protein, 1.5-MDa RNA polymerase II preinitiation complex (PIC) was determined by cryo-EM and image processing at a resolution of 6-11 Å. Atomic structures of over 50% of the mass were fitted into the electron density map in a manner consistent with protein-protein cross-links previously identified by mass spectrometry. The resulting model of the PIC confirmed the main conclusions from previous cryo-EM at lower resolution, including the association of promoter DNA only with general transcription factors and not with the polymerase. Electron density due to DNA was identifiable by the grooves of the double helix and exhibited sharp bends at points downstream of the TATA box, with an important consequence: The DNA at the downstream end coincides with the DNA in a transcribing polymerase. The structure of the PIC is therefore conducive to promoter melting, start-site scanning, and the initiation of transcription.

The protein density and arrangement of subunits of a complete, 32-protein, RNA polymerase II (pol II) transcription pre-initiation complex (PIC) were determined by means of cryogenic electron microscopy and a combination of chemical cross-linking and mass spectrometry. The PIC showed a marked division in two parts, one containing all the general transcription factors (GTFs) and the other pol II. Promoter DNA was associated only with the GTFs, suspended above the pol II cleft and not in contact with pol II. This structural principle of the PIC underlies its conversion to a transcriptionally active state; the PIC is poised for the formation of a transcription bubble and descent of the DNA into the pol II cleft.

In eukaryotes, CCT is essential for the correct and efficient folding of many cytosolic proteins, most notably actin and tubulin. Structural studies of CCT have been hindered by the failure of standard crystallographic analysis to resolve its eight different subunit types at low resolutions. Here, we exhaustively assess the R value fit of all possible CCT models to available crystallographic data of the closed and open forms with resolutions of 3.8 Å and 5.5 Å, respectively. This unbiased analysis finds the native subunit arrangements with overwhelming significance. The resulting structures provide independent crystallographic proof of the subunit arrangement of CCT and map major asymmetrical features of the particle onto specific subunits. The actin and tubulin substrates both bind around subunit CCT6, which shows other structural anomalies. CCT is thus clearly partitioned, both functionally and evolutionary, into a substrate-binding side that is opposite to the ATP-hydrolyzing side., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

The TRiC/CCT chaperonin is a 1-MDa hetero-oligomer of 16 subunits that assists the folding of proteins in eukaryotes. Low-resolution structural studies confirmed the TRiC particle to be composed of two stacked octameric rings enclosing a folding cavity. The exact arrangement of the different proteins in the rings underlies the functionality of TRiC and is likely to be conserved across all eukaryotes. Yet despite its importance it has not been determined conclusively, mainly because the different subunits appear nearly identical under low resolution. This work successfully addresses the arrangement problem by the emerging technique of cross-linking, mass spectrometry, and modeling. We cross-linked TRiC under native conditions with a cross-linker that is primarily reactive toward exposed lysine side chains that are spatially close in the context of the particle. Following digestion and mass spectrometry we were able to identify over 60 lysine pairs that underwent cross-linking, thus providing distance restraints between specific residues in the complex. Independently of the cross-link set, we constructed 40,320 (= 8 factorial) computational models of the TRiC particle, which exhaustively enumerate all the possible arrangements of the different subunits. When we assessed the compatibility of each model with the cross-link set, we discovered that one specific model is significantly more compatible than any other model. Furthermore, bootstrapping analysis confirmed that this model is 10 times more likely to result from this cross-link set than the next best-fitting model. Our subunit arrangement is very different than any of the previously reported models and changes the context of existing and future findings on TRiC.

The identification of catalytic residues is an essential step in functional characterization of enzymes. We present a purely structural approach to this problem, which is motivated by the difficulty of evolution-based methods to annotate structural genomics targets that have few or no homologs in the databases. Our approach combines a state-of-the-art support vector machine (SVM) classifier with novel structural features that augment structural clues by spatial averaging and Z scoring. Special attention is paid to the class imbalance problem that stems from the overwhelming number of non-catalytic residues in enzymes compared to catalytic residues. This problem is tackled by: (1) optimizing the classifier to maximize a performance criterion that considers both Type I and Type II errors in the classification of catalytic and non-catalytic residues; (2) under-sampling non-catalytic residues before SVM training; and (3) during SVM training, penalizing errors in learning catalytic residues more than errors in learning non-catalytic residues. Tested on four enzyme datasets, one specifically designed by us to mimic the structural genomics scenario and three previously evaluated datasets, our structure-based classifier is never inferior to similar structure-based classifiers and comparable to classifiers that use both structural and evolutionary features. In addition to the evaluation of the performance of catalytic residue identification, we also present detailed case studies on three proteins. This analysis suggests that many false positive predictions may correspond to binding sites and other functional residues. A web server that implements the method, our own-designed database, and the source code of the programs are publicly available at http://www.cs.bgu.ac.il/∼meshi/functionPrediction., (Copyright © 2011 Wiley-Liss, Inc.)

The synthesis of Au(102)(p-MBA)(44) nanoparticles on a preparative scale in high yield is described. Various analytical methods are shown to give results consistent with the composition and known structure of the particles, showing the preparation is essentially homogeneous, and attesting to the validity of the methods as well. Derivatization of the particles with proteins and DNA is demonstrated, and conditions are described for imaging individual particles by cryo-EM at low electron dose, close to focus, conditions optimal for recording high-resolution details.

Protein structure refinement is an important but unsolved problem; it must be solved if we are to predict biological function that is very sensitive to structural details. Specifically, critical assessment of techniques for protein structure prediction (CASP) shows that the accuracy of predictions in the comparative modeling category is often worse than that of the template on which the homology model is based. Here we describe a refinement protocol that is able to consistently refine submitted predictions for all categories at CASP7. The protocol uses direct energy minimization of the knowledge-based potential of mean force that is based on the interaction statistics of 167 atom types (Summa and Levitt, Proc Natl Acad Sci USA 2007; 104:3177-3182). Our protocol is thus computationally very efficient; it only takes a few minutes of CPU time to run typical protein models (300 residues). We observe an average structural improvement of 1% in GDT_TS, for predictions that have low and medium homology to known PDB structures (Global Distance Test score or GDT_TS between 50 and 80%). We also observe a marked improvement in the stereochemistry of the models. The level of improvement varies amongst the various participants at CASP, but we see large improvements (>10% increase in GDT_TS) even for models predicted by the best performing groups at CASP7. In addition, our protocol consistently improved the best predicted models in the refinement category at CASP7 and CASP8. These improvements in structure and stereochemistry prove the usefulness of our computationally inexpensive, powerful and automatic refinement protocol., (Copyright 2010 Wiley-Liss, Inc.)

TRiC is an important group II chaperonin that facilitates the folding of many eukaryotic proteins. The TRiC complex consists of two stacked rings, each comprised of eight paralogous subunits with a mutual sequence identity of 30-35%. Each subunit has unique functional roles that are manifested by corresponding sequence conservation. It is generally assumed that the subunit order within each ring is fixed, but this order is still uncertain. Here we address the problem of the intra-ring subunit order by combining two sources of information: evolutionary conservation and a structural hypothesis. Specifically, we identify residues in the TRiC subunits that are likely to be part of the intra-unit interface, based on homology modeling to the solved thermosome structure. Within this set of residues, we search for a subset that shows an evolutionary conservation pattern that is indicative of the subunit order key. This pattern shows considerable conservation across species, but large variation across the eight subunits. By this approach we were able to locate two parts of the interface where complementary interactions seem to favor certain pairing of subunits. This knowledge leads to restrictions on the 5,040 (=7!) possible subunits arrangements in the ring, and limits them to just 72. Although our findings give only partial understanding of the inter-subunit interactions that determine their order, we conclude that they are comprised of complementary charged, polar and hydrophobic interactions that occur in both the equatorial and middle domains of each subunit.

Rotatable torsion angles are the major degrees of freedom in proteins. Adjacent angles are highly correlated and energy terms that rely on these correlations are intensively used in molecular modeling. However, the utility of torsion based terms is not yet fully exploited. Many of these terms do not capture the full scale of the correlations. Other terms, which rely on lookup tables, cannot be used in the context of force-driven algorithms because they are not fully differentiable. This study aims to extend the usability of torsion terms by presenting a set of high-dimensional and fully-differentiable energy terms that are derived from high-resolution structures. The set includes terms that describe backbone conformational probabilities and propensities, side-chain rotamer probabilities, and an elaborate term that couples all the torsion angles within the same residue. The terms are constructed by cubic spline interpolation with periodic boundary conditions that enable full differentiability and high computational efficiency. We show that the spline implementation does not compromise the accuracy of the original database statistics. We further show that the side-chain relevant terms are compatible with established rotamer probabilities. Despite their very local characteristics, the new terms are often able to identify native and native-like structures within decoy sets. Finally, force-based minimization of NMR structures with the new terms improves their torsion angle statistics with minor structural distortion (0.5 A RMSD on average). The new terms are freely available in the MESHI molecular modeling package. The spline coefficients are also available as a documented MATLAB file., (2008 Wiley-Liss, Inc.)

Background: The structural stability of peptides in solution strongly affects their binding affinities and specificities. Thus, in peptide biotechnology, an increase in the structural stability is often desirable. The present work combines two orthogonal computational techniques, Molecular Dynamics and a knowledge-based potential, for the prediction of structural stability of short peptides (< 20 residues) in solution., Results: We tested the new approach on four families of short beta-hairpin peptides: TrpZip, MBH, bhpW and EPO, whose structural stabilities have been experimentally measured in previous studies. For all four families, both computational techniques show considerable correlation (r > 0.65) with the experimentally measured stabilities. The consensus of the two techniques shows higher correlation (r > 0.82)., Conclusion: Our results suggest a prediction scheme that can be used to estimate the relative structural stability within a peptide family. We discuss the applicability of this predictive approach for in-silico screening of combinatorial peptide libraries.

The neocortex has a high capacity for plasticity. To understand the full scope of this capacity, it is essential to know how neurons choose particular partners to form synaptic connections. By using multineuron whole-cell recordings and confocal microscopy we found that axons of layer V neocortical pyramidal neurons do not preferentially project toward the dendrites of particular neighboring pyramidal neurons; instead, axons promiscuously touch all neighboring dendrites without any bias. Functional synaptic coupling of a small fraction of these neurons is, however, correlated with the existence of synaptic boutons at existing touch sites. These data provide the first direct experimental evidence for a tabula rasa-like structural matrix between neocortical pyramidal neurons and suggests that pre- and postsynaptic interactions shape the conversion between touches and synapses to form specific functional microcircuits. These data also indicate that the local neocortical microcircuit has the potential to be differently rewired without the need for remodeling axonal or dendritic arbors.

A model is presented that allows prediction of the probability for the formation of appositions between the axons and dendrites of any two neurons based only on their morphological statistics and relative separation. Statistics of axonal and dendritic morphologies of single neurons are obtained from 3D reconstructions of biocytin-filled cells, and a statistical representation of the same cell type is obtained by averaging across neurons according to the model. A simple mathematical formulation is applied to the axonal and dendritic statistical representations to yield the probability for close appositions. The model is validated by a mathematical proof and by comparison of predicted appositions made by layer 5 pyramidal neurons in the rat somatosensory cortex with real anatomical data. The model could be useful for studying microcircuit connectivity and for designing artificial neural networks.