Unsupervised classification of single particles by cluster tracking in multi-dimensional space

In cryo-electron microscopy (cryo-EM) single-particle reconstruction, the heterogeneity of two-dimensional projection image data resulting from the co-existence of different conformational or ligand binding states of a macromolecular complex remains a major obstacle as it impairs the validity of reconstructed density maps and limits the progress toward higher resolution. Classification of cryo-EM data according to the different conformations is difficult because of the coexistence of multiple orientations in a single dataset. Here, we present an unsupervised classification method, termed cluster tracking, which utilizes the continuity in multi-dimensional space induced by angular adjacency of projections in large datasets. In a proof of concept, the testing of cluster tracking on simulated projection data, which were generated from multiple conformations and orientations of an existing volume, produced clusters that are consistent with the conformational identity of the data. The application of the method to experimental cryo-EM projection data is found to result in a partition similar to the one generated by supervised classification.     

C-terminal Src kinase-homologous kinase (CHK), a unique inhibitor inactivating multiple active conformations of Src family tyrosine kinases

The Src family of protein kinases (SFKs) mediates mitogenic signal transduction, and constitutive SFK activation is associated with tumorigenesis. To prevent constitutive SFK activation, the catalytic activity of SFKs in normal mammalian cells is suppressed mainly by two inhibitors called C-terminal Src kinase (CSK) and CSK-homologous kinase (CHK), which inactivate SFKs by phosphorylating a consensus tyrosine near the C terminus of SFKs (Y(T)). The phosphorylated Y(T) intramolecularly binds to the SH2 domain of SFKs. This interaction, known as pY(T)/SH2 interaction, together with binding between the SH2 kinase linker and the SH3 domain of SFKs (linker/SH3 interaction) stabilizes SFKs in a "closed" inactive conformation. We previously discovered an alternative mechanism CHK employs to inhibit SFKs. This mechanism, referred to as the non-catalytic inhibitory mechanism, involves tight binding of CHK to SFKs; the binding alone is sufficient to inhibit SFKs. Herein, we constructed multiple active conformations of an SFK member, Hck, by systematically disrupting the two inhibitory interactions. We found that CHK employs the non-catalytic mechanism to inactivate these active conformations of Hck. However, CHK does not bind Hck when it adopts the inactive conformation in which both inhibitory interactions are intact. These data indicate that binding of CHK to SFKs via the non-catalytic mechanism is governed by the conformations of SFKs. Although CSK is also an inhibitor of SFKs, it does not inhibit SFKs by a similar non-catalytic mechanism. Thus, the non-catalytic inhibitory mechanism is a unique property of CHK that allows it to down-regulate multiple active conformations of SFKs.     

Stabilization of Helical Order in the Thick Filaments by Blebbistatin: Further Evidence of Coexisting Multiple Conformations of Myosin

The degree of helical order of the thick filament of mammalian skeletal muscle is highly dependent on temperature and the nature of the ligand. Previously, we showed that there was a close correlation between the conformation of the myosin heads on the surface of the thick filaments and the extent of their helical order. Helical order required the heads to be in the closed conformation. In addition, we showed that, with the same ligand bound at the active site, three conformations of myosin coexisted in equilibrium. Hitherto, however, there was no detectable helical order as measured by x-ray diffraction under the temperatures studied for myosin with MgADP and the nucleotide-free myosin, raising the possibility that the concept of multiple conformations has limited validity. In this study, blebbistatin was used to stabilize the closed conformation of myosin. The degree of helical order is substantially improved with MgATP at low temperature or with MgADP or in the absence of nucleotide. The thermodynamic parameters of the disorder↔order transition and the characteristics of the ordered array were not significantly altered by binding blebbistatin. The simplest explanation is that the binding of blebbistatin increases the proportion of myosin in the closed conformation from being negligible to substantial. These results provide further evidence for the coexistence of multiple conformations of myosin under a wide range of conditions and for the closed conformation being directly coupled to helical order.     

Using least median of squares for structural superposition of flexible proteins

Background

The conventional superposition methods use an ordinary least squares (LS) fit for structural comparison of two different conformations of the same protein. The main problem of the LS fit that it is sensitive to outliers, i.e. large displacements of the original structures superimposed.

Results

To overcome this problem, we present a new algorithm to overlap two protein conformations by their atomic coordinates using a robust statistics technique: least median of squares (LMS). In order to effectively approximate the LMS optimization, the forward search technique is utilized. Our algorithm can automatically detect and superimpose the rigid core regions of two conformations with small or large displacements. In contrast, most existing superposition techniques strongly depend on the initial LS estimating for the entire atom sets of proteins. They may fail on structural superposition of two conformations with large displacements. The presented LMS fit can be considered as an alternative and complementary tool for structural superposition.

Conclusion

The proposed algorithm is robust and does not require any prior knowledge of the flexible regions. Furthermore, we show that the LMS fit can be extended to multiple level superposition between two conformations with several rigid domains. Our fit tool has produced successful superpositions when applied to proteins for which two conformations are known. The binary executable program for Windows platform, tested examples, and database are available from https://engineering.purdue.edu/PRECISE/LMSfit.     

An RNA-binding chameleon

The arginine-rich RNA binding motif is found in a wide variety of proteins, including several viral regulatory proteins. Although related at the primary sequence level, arginine-rich domains from different proteins adopt different conformations depending on the RNA site recognized, and in some cases fold only in the context of RNA. Here we show that the RNA binding domain of the Jembrana disease virus (JDV) Tat protein is able to recognize two different TAR RNA sites, from human and bovine immunodeficiency viruses (HIV and BIV, respectively), adopting different conformations in the two RNA contexts and using different amino acids for recognition. In addition to the conformational differences, the JDV domain requires the cyclin T1 protein for high-affinity binding to HIV TAR, but not to BIV TAR. The "chameleon-like" behavior of the JDV Tat RNA binding domain reinforces the concept that RNA molecules can provide structural scaffolds for protein folding, and suggests mechanisms for evolving distinct RNA binding specificities from a single multifunctional domain.     

A population-based evolutionary search approach to the multiple minima problem in de novo protein structure prediction

Background

Elucidating the native structure of a protein molecule from its sequence of amino acids, a problem known as de novo structure prediction, is a long standing challenge in computational structural biology. Difficulties in silico arise due to the high dimensionality of the protein conformational space and the ruggedness of the associated energy surface. The issue of multiple minima is a particularly troublesome hallmark of energy surfaces probed with current energy functions. In contrast to the true energy surface, these surfaces are weakly-funneled and rich in comparably deep minima populated by non-native structures. For this reason, many algorithms seek to be inclusive and obtain a broad view of the low-energy regions through an ensemble of low-energy (decoy) conformations. Conformational diversity in this ensemble is key to increasing the likelihood that the native structure has been captured.

Methods

We propose an evolutionary search approach to address the multiple-minima problem in decoy sampling for de novo structure prediction. Two population-based evolutionary search algorithms are presented that follow the basic approach of treating conformations as individuals in an evolving population. Coarse graining and molecular fragment replacement are used to efficiently obtain protein-like child conformations from parents. Potential energy is used both to bias parent selection and determine which subset of parents and children will be retained in the evolving population. The effect on the decoy ensemble of sampling minima directly is measured by additionally mapping a conformation to its nearest local minimum before considering it for retainment. The resulting memetic algorithm thus evolves not just a population of conformations but a population of local minima.

Results and conclusions

Results show that both algorithms are effective in terms of sampling conformations in proximity of the known native structure. The additional minimization is shown to be key to enhancing sampling capability and obtaining a diverse ensemble of decoy conformations, circumventing premature convergence to sub-optimal regions in the conformational space, and approaching the native structure with proximity that is comparable to state-of-the-art decoy sampling methods. The results are shown to be robust and valid when using two representative state-of-the-art coarse-grained energy functions.

     

Membrane signalling complexes: implications for development of functionally selective ligands modulating heptahelical receptor signalling

Technological development has considerably changed the way in which we evaluate drug efficacy and has led to a conceptual revolution in pharmacological theory. In particular, molecular resolution assays have revealed that heptahelical receptors may adopt multiple active conformations with unique signalling properties. It is therefore becoming widely accepted that ligand ability to stabilize receptor conformations with distinct signalling profiles may allow to direct the stimulus generated by an activated receptor towards a specific signalling pathway. This capacity to induce only a subset of the ensemble of responses regulated by a given receptor has been termed "functional selectivity" (or "stimulus trafficking"), and provides the bases for a highly specific regulation of receptor signalling. Concomitant with these observations, heptahelical receptors have been shown to associate with G proteins and effectors to form multimeric arrays. These complexes are constitutively formed during protein synthesis and are targeted to the cell surface as integral signalling units. Herein we summarize evidence supporting the existence of such constitutive signalling arrays and analyze the possibility that they may constitute viable targets for developing ligands with "functional selectivity".     

Predicting the preferred conformations of luteolin-4′-O-β-D-glucoside in gas phase: a comparison of two computational approaches

A tree-step computational approach has been applied to determine the lowest-energy conformers of luteolin-4′-O-β-D-glucoside (L4′G). Fifty-seven starting structures of the L4′G have been built, and then by performing with density functional theory (DFT) optimizations and second-order Møller-Plesset (MP2) calculations, the preferred conformations of L4′G are predicted. In order to test the accuracy of the computational approach, a hybrid Monte-Carlo multiple minimum (MCMM)/quantum mechanical (QM) approach is applied to determine the favorable conformers of L4′G. The alternative classification is employed to put similar conformations into the same catalogue according to the dihedral angles among the luteolin rings, glycosidic dihedral angles, and the orientations of hydroxyl and hydroxymethyl groups. The low-energy conformations are located after the optimizations at the HF/6-31G(d) and B3LYP/6-311+G(d) levels. Compared with the hybrid MCMM/QM approach, the tree-step computational approach not only remains accurate but also saves a lot of computing resources.

Mechanism of cross-species prion transmission: an infectious conformation compatible with two highly divergent yeast prion proteins

Efficiency of interspecies prion transmission decreases as the primary structures of the infectious proteins diverge. Yet, a single prion protein can misfold into multiple infectious conformations, and such differences in "strain conformation" also alter infection specificity. Here, we explored the relationship between prion strains and species barriers by creating distinct synthetic prion forms of the yeast prion protein Sup35. We identified a strain conformation of Sup35 that allows transmission from the S. cerevisiae (Sc) Sup35 to the highly divergent C. albicans (Ca) Sup35 both in vivo and in vitro. Remarkably, cross-species transmission leads to a novel Ca strain that in turn can infect the Sc protein. Structural studies reveal strain-specific conformational differences in regions of the prion domain that are involved in intermolecular contacts. Our findings support a model whereby strain conformation is the critical determinant of cross-species prion transmission while primary structure affects transmission specificity by altering the spectrum of preferred amyloid conformations.     

Biased signaling: potential agonist and antagonist of PAR2

Protease activated receptor 2 (PAR2) has emerged as one of the promising therapeutic targets to inhibit rapidly metastasizing breast cancer cells. However, its elusive molecular mechanism of activation and signaling has made it a difficult target for drug development. In this study, in silico methods were used to unfold PAR2 molecular mechanism of signaling based on the concept of GPCR receptor plasticity. Although, there are no conclusive evidences of the presence of specific endogenous ligands for PAR2, the efficacy of synthetic agonist and antagonist in PAR2 signaling has opened up the possibilities of ligand-mediated signaling. Furthermore, it has been proved that ligands specific for one GPCR can induce signaling in GPCRs belonging to other subfamilies. Therefore, the aim of this study was to identify potential agonists and antagonists from the GPCR ligand library (GLL), which may induce biased signaling in PAR2 using the concept of existence of multiple ligand-stabilized receptor conformations. The results of our in silico study suggest that PAR2 may show biased signaling mainly with agonists of serotonin type 1, β-adrenergic type 1,3 and antagonists of substance K (NK1), serotonin type 2, dopamine type 4, and thromboxane receptors. Further, this study also throws light on the putative ligand-specific conformations of PAR2. Thus, the results of this study provide structural insights to putative conformations of PAR2 and also gives initial clues to medicinal chemists for rational drug design targeting this challenging receptor.     

How to choose relevant multiple receptor conformations for virtual screening: a test case of Cdk2 and normal mode analysis

Better treatment of protein flexibility is essential in structure-based drug design projects such as virtual screening and protein-ligand docking. Diversity in ligand-binding mechanisms and receptor conformational changes makes it difficult to treat dynamic features of the receptor during the docking simulation. Thus, the use of pregenerated multiple receptor conformations is applied today in virtual screening studies. However, generation of a small relevant set of receptor conformations remains challenging. To address this problem, we propose a new protocol for the generation of multiple receptor conformations via normal mode analysis and for the selection of several receptor conformations suitable for docking/virtual screening. We validated this protocol on cyclin-dependent kinase 2, which possesses a binding site located at the interface between two subdomains and is known to undergo significant conformational changes in the active site region upon ligand binding. We believe that the suggested rules for the choice of suitable receptor conformations can be applied to other targets when dealing with in silico screening on flexible receptors.     

Cyclic AMP analog blocks kinase activation by stabilizing inactive conformation: conformational selection highlights a new concept in allosteric inhibitor design

The regulatory (R) subunit of protein kinase A serves to modulate the activity of protein kinase A in a cAMP-dependent manner and exists in two distinct and structurally dissimilar, end point cAMP-bound "B" and C-subunit-bound "H"-conformations. Here we report mechanistic details of cAMP action as yet unknown through a unique approach combining x-ray crystallography with structural proteomics approaches, amide hydrogen/deuterium exchange and ion mobility mass spectrometry, applied to the study of a stereospecific cAMP phosphorothioate analog and antagonist((Rp)-cAMPS). X-ray crystallography shows cAMP-bound R-subunit in the B form but surprisingly the antagonist Rp-cAMPS-bound R-subunit crystallized in the H conformation, which was previously assumed to be induced only by C-subunit-binding. Apo R-subunit crystallized in the B form as well but amide exchange mass spectrometry showed large differences between apo, agonist and antagonist-bound states of the R-subunit. Further ion mobility reveals the apo R-subunit as an ensemble of multiple conformations with collisional cross-sectional areas spanning both the agonist and antagonist-bound states. Thus contrary to earlier studies that explained the basis for cAMP action through "induced fit" alone, we report evidence for conformational selection, where the ligand-free apo form of the R-subunit exists as an ensemble of both B and H conformations. Although cAMP preferentially binds the B conformation, Rp-cAMPS interestingly binds the H conformation. This reveals the unique importance of the equatorial oxygen of the cyclic phosphate in mediating conformational transitions from H to B forms highlighting a novel approach for rational structure-based drug design. Ideal inhibitors such as Rp-cAMPS are those that preferentially "select" inactive conformations of target proteins by satisfying all "binding" constraints alone without inducing conformational changes necessary for activation.     

Protein-protein docking with multiple residue conformations and residue substitutions

The protein docking problem has two major aspects: sampling conformations and orientations, and scoring them for fit. To investigate the extent to which the protein docking problem may be attributed to the sampling of ligand side‐chain conformations, multiple conformations of multiple residues were calculated for the uncomplexed (unbound) structures of protein ligands. These ligand conformations were docked into both the complexed (bound) and unbound conformations of the cognate receptors, and their energies were evaluated using an atomistic potential function. The following questions were considered: (1) does the ensemble of precalculated ligand conformations contain a structure similar to the bound form of the ligand? (2) Can the large number of conformations that are calculated be efficiently docked into the receptors? (3) Can near‐native complexes be distinguished from non‐native complexes? Results from seven test systems suggest that the precalculated ensembles do include side‐chain conformations similar to those adopted in the experimental complexes. By assuming additivity among the side chains, the ensemble can be docked in less than 12 h on a desktop computer. These multiconformer dockings produce near‐native complexes and also non‐native complexes. When docked against the bound conformations of the receptors, the near‐native complexes of the unbound ligand were always distinguishable from the non‐native complexes. When docked against the unbound conformations of the receptors, the near‐native dockings could usually, but not always, be distinguished from the non‐native complexes. In every case, docking the unbound ligands with flexible side chains led to better energies and a better distinction between near‐native and non‐native fits. An extension of this algorithm allowed for docking multiple residue substitutions (mutants) in addition to multiple conformations. The rankings of the docked mutant proteins correlated with experimental binding affinities. These results suggest that sampling multiple residue conformations and residue substitutions of the unbound ligand contributes to, but does not fully provide, a solution to the protein docking problem. Conformational sampling allows a classical atomistic scoring function to be used; such a function may contribute to better selectivity between near‐native and non‐native complexes. Allowing for receptor flexibility may further extend these results.     

Human TFIID Binds to Core Promoter DNA in a Reorganized Structural State

1. Download : Download high-res image (197KB)
2. Download : Download full-size image

Highlights? TFIID has distinct conformations with different positions of the TBP-containing lobe ? A rearranged state of TFIID interacts with promoter DNA in a TFIIA-dependent manner ? TFIID induces multiple topological changes on promoter DNA around the start site ? The different conformations of TFIID may serve as targets for regulatory factors     

The role of Sse1 in the de novo formation and variant determination of the [PSI+] prion

Yeast prions are a group of non-Mendelian genetic elements transmitted as altered and self-propagating conformations. Extensive studies in the last decade have provided valuable information on the mechanisms responsible for yeast prion propagation. How yeast prions are formed de novo and what cellular factors are required for determining prion "strains" or variants--a single polypeptide capable of existing in multiple conformations to result in distinct heritable phenotypes--continue to defy our understanding. We report here that Sse1, the yeast ortholog of the mammalian heat-shock protein 110 (Hsp110) and a nucleotide exchange factor for Hsp70 proteins, plays an important role in regulating [PSI+] de novo formation and variant determination. Overproduction of the Sse1 chaperone dramatically enhanced [PSI+] formation whereas deletion of SSE1 severely inhibited it. Only an unstable weak [PSI+] variant was formed in SSE1 disrupted cells whereas [PSI+] variants ranging from very strong to very weak were formed in isogenic wild-type cells under identical conditions. Thus, Sse1 is essential for the generation of multiple [PSI+] variants. Mutational analysis further demonstrated that the physical association of Sse1 with Hsp70 but not the ATP hydrolysis activity of Sse1 is required for the formation of multiple [PSI+] variants. Our findings establish a novel role for Sse1 in [PSI+] de novo formation and variant determination, implying that the mammalian Hsp110 may likewise be involved in the etiology of protein-folding diseases.     

Molecular modelling of glycans: three-dimensional structure and protein fraction interaction. The example of rabbit sero-transferrin

On the basis of experimental data and of computer calculations using the Tripos 5.3 force field in order to examine the three-dimensional structures which are sterically feasible and the conformations which are energetically the most favourable, we have designed a program of molecular modelling of biantennary glycans of the N-acetyllactosaminic type (complex type). We demonstrate that, in absence of any interaction with the protein, a high number of glycan conformations exists which can be classified into five basic conformations, four of which have already been described. In fact, in addition to the Y-, T-, "bird" and "broken-wing" conformations, a "back-folded wing" conformation is energetically feasible. In contrast, the glycan linked to the protein is immobilized into only one conformation: the "broken-wing" conformation. Forming a bridge between the two lobes of the peptide chain, it probably contributes to the maintenance of the protein in a biologically active conformation.     

Reducing the computational complexity of protein folding via fragment folding and assembly

Understanding, and ultimately predicting, how a 1-D protein chain reaches its native 3-D fold has been one of the most challenging problems during the last few decades. Data increasingly indicate that protein folding is a hierarchical process. Hence, the question arises as to whether we can use the hierarchical concept to reduce the practically intractable computational times. For such a scheme to work, the first step is to cut the protein sequence into fragments that form local minima on the polypeptide chain. The conformations of such fragments in solution are likely to be similar to those when the fragments are embedded in the native fold, although alternate conformations may be favored during the mutual stabilization in the combinatorial assembly process. Two elements are needed for such cutting: (1) a library of (clustered) fragments derived from known protein structures and (2) an assignment algorithm that selects optimal combinations to "cover" the protein sequence. The next two steps in hierarchical folding schemes, not addressed here, are the combinatorial assembly of the fragments and finally, optimization of the obtained conformations. Here, we address the first step in a hierarchical protein-folding scheme. The input is a target protein sequence and a library of fragments created by clustering building blocks that were generated by cutting all protein structures. The output is a set of cutout fragments. We briefly outline a graph theoretic algorithm that automatically assigns building blocks to the target sequence, and we describe a sample of the results we have obtained.     

Conformational-functional relationships of tetragastrin analogs]

Sets of low-energy backbone conformations of the active tetragastrin analogue Boc-Trp-Leu-Asp-Phe-NH2 and two competitive antagonists Boc-Trp-Leu psi (CH2NH)-Asp-Phe-NH2 and Boc-Trp-Leu-Asp-O-CH2-CH2-C6H5 were obtained using theoretical conformational analysis methods. Groups of the conformations were selected for the three analogues, allowing a spatial matching of Trp, Asp and Phe residues responsible for the gastrin receptor binding. Three conformations possessing the lowest energies among the geometrically similar structures of these three peptides are suggested as a model for the "receptor-bound" conformations of these analogues. Backbone spatial folding resembling an alpha-helix turn is characteristic of these conformations. The correspondence of the proposed model to the available data on structure--activity relationships for tetragastrin analogues is discussed. Orientations of the putative receptor-bound conformations in a "water--lypophylic medium" two-phase system were investigated.