Mutations in LZTR1 drive human disease by dysregulating RAS ubiquitination.

The leucine zipper-like transcriptional regulator 1 (LZTR1) protein, an adaptor for cullin 3 (CUL3) ubiquitin ligase complex, is implicated in human disease, yet its mechanism of action remains unknown. We found that Lztr1 haploinsufficiency in mice recapitulates Noonan syndrome phenotypes, whereas LZTR1 loss in Schwann cells drives dedifferentiation and proliferation. By trapping LZTR1 complexes from intact mammalian cells, we identified the guanosine triphosphatase RAS as a substrate for the LZTR1-CUL3 complex. Ubiquitome analysis showed that loss of Lztr1 abrogated Ras ubiquitination at lysine-170. LZTR1-mediated ubiquitination inhibited RAS signaling by attenuating its association with the membrane. Disease-associated LZTR1 mutations disrupted either LZTR1-CUL3 complex formation or its interaction with RAS proteins. RAS regulation by LZTR1-mediated ubiquitination provides an explanation for the role of LZTR1 in human disease.

GOSSIP: a method for fast and accurate global alignment of protein structures.

MOTIVATION: The database of known protein structures (PDB) is increasing rapidly. This results in a growing need for methods that can cope with the vast amount of structural data. To analyze the accumulating data, it is important to have a fast tool for identifying similar structures and clustering them by structural resemblance. Several excellent tools have been developed for the comparison of protein structures. These usually address the task of local structure alignment, an important yet computationally intensive problem due to its complexity. It is difficult to use such tools for comparing a large number of structures to each other at a reasonable time. RESULTS: Here we present GOSSIP, a novel method for a global all-against-all alignment of any set of protein structures. The method detects similarities between structures down to a certain cutoff (a parameter of the program), hence allowing it to detect similar structures at a much higher speed than local structure alignment methods. GOSSIP compares many structures in times which are several orders of magnitude faster than well-known available structure alignment servers, and it is also faster than a database scanning method. We evaluate GOSSIP both on a dataset of short structural fragments and on two large sequence-diverse structural benchmarks. Our conclusions are that for a threshold of 0.6 and above, the speed of GOSSIP is obtained with no compromise of the accuracy of the alignments or of the number of detected global similarities. AVAILABILITY: A server, as well as an executable for download, are available at http://bioinfo3d.cs.tau.ac.il/gossip/.

Towards drugs targeting multiple proteins in a systems biology approach.

Protein-protein interactions are increasingly becoming drug targets. This is understandable, since they are crucial at all levels of cellular expression and growth. In practice, targeting specific disease-related interactions has proven difficult, with success varying with specific complexes. Here, we take a Systems Biology approach to targeting protein-protein interactions. Below, we first briefly review drug discovery targeted at protein-protein interactions; we classify protein-protein complexes with respect to their types of interactions and their roles in cellular function and as being targets in drug design; we describe the properties of the interfaces as related to drug design, focusing on hot spots and surface cavities; and finally, in particular, we cast the interactions into the cellular network system, highlighting the challenge of partially targeting multiple interactions in the networks as compared to hitting a specific protein-protein interaction target. The challenge we now face is how to pick the targets and how to improve the efficiency of designed partially-specific multi-target drugs that would block parallel pathways in the network.

MASS: multiple structural alignment by secondary structures.

We present a novel method for multiple alignment of protein structures and detection of structural motifs. To date, only a few methods are available for addressing this task. Most of them are based on a series of pairwise comparisons. In contrast, MASS (Multiple Alignment by Secondary Structures) considers all the given structures at the same time. Exploiting the secondary structure representation aids in filtering out noisy results and in making the method highly efficient and robust. MASS disregards the sequence order of the secondary structure elements. Thus, it can find non-sequential and even non-topological structural motifs. An important novel feature of MASS is subset alignment detection: It does not require that all the input molecules be aligned. Rather, MASS is capable of detecting structural motifs shared only by a subset of the molecules. Given its high efficiency and capability of detecting subset alignments, MASS is suitable for a broad range of challenging applications: It can handle large-scale protein ensembles (on the order of tens) that may be heterogeneous, noisy, topologically unrelated and contain structures of low resolution.

Interdomain interactions in hinge-bending transitions.

BACKGROUND: The mechanisms that allow or constrain protein movement have not been understood. Here we study interdomain interactions in proteins to investigate hinge-bending motions. RESULTS: We find a limited number of salt bridges and hydrogen bonds at the interdomain interface, in both the "closed" and the "open" conformations. Consistently, analysis of 222 salt bridges in an independently selected database indicates that most salt bridges form within rather than between independently folding hydrophobic units. Calculations show that these interdomain salt bridges either destabilize or only marginally stabilize the closed conformation in most proteins. In contrast, the nonpolar buried surface area between the moving parts can be extensive in the closed conformations. However, when the nonpolar buried surface area is large, we find that at the interdomain interface in the open conformation it may be as large or larger than in the closed conformation. Hence, the energetic penalty of opening the closed conformation is overcome. Consistently, a large nonpolar surface area buried in the closed interdomain interface accompanies limited opening of the domains, yielding a larger interface. CONCLUSIONS: Short-range electrostatic interactions are largely absent between moving domains. Interdomain nonpolar buried surface area may be large in the closed conformation, but it is largely offset by the area buried in the open conformation. In such cases the opening of the domains appears to be relatively small. This may allow prediction of the extent of domain opening. Such predictions may have implications for the shape and size of the binding pockets in drug/protein design.

The building block folding model and the kinetics of protein folding.

Here we show that qualitatively, the building blocks folding model accounts for three-state versus the two-state protein folding. Additionally, it is consistent with the faster versus slower folding rates of the two-state proteins. Specifically, we illustrate that the building blocks size, their mode of associations in the native structure, the number of ways they can combinatorially assemble, their population times and the way they are split in the iterative, step-by-step structural dissection which yields the anatomy trees, explain a broad range of folding rates. We further show that proteins with similar general topologies may have different folding pathways, and hence different folding rates. On the other hand, the effect of mutations resembles that of changes in conditions, shifting the population times and hence the energy landscapes. Hence, together with the secondary structure type and the extent of local versus non-local interactions, a coherent, consistent rationale for folding kinetics can be outlined, in agreement with experimental results. Given the native structure of a protein, these guidelines enable a qualitative prediction of the folding kinetics. We further describe these in the context of the protein folding energy landscape. Quantitatively, in principle, the diffusion-collision model for the building block association can be used. However, the folding rates of the building blocks and traps in their formation and association, need to be considered.

Protein folding: binding of conformationally fluctuating building blocks via population selection.

Here we review different aspects of the protein folding literature. We present a broad range of observations, showing them to be consistent with a general hierarchical protein folding model. In such a model, local relatively stable, conformationally fluctuating building blocks bind through population selection, to yield the native state. The model includes several components: (1) the fluctuating building blocks that constitute local minima along the polypeptide chain, which even if unstable still possess higher population times than all alternate conformations; (2) the landscape around the bottom of the funnels; (3) the consideration that protein folding involves intramolecular recognition; (4) similar landscapes are observed for folding and for binding, and that (5) the landscape is dynamic, changing with the conditions. The model considers protein folding to be guided by native interactions. The reviewed literature includes the effects of changing the conditions, intermediates and kinetic traps, mutations, similar topologies, fragment complementation experiments, fragments and pathways, focusing on one specific well-studied example, that of the dihydrofolate reductase, chaperones, and chaperonines, in vivo vs. in vitro folding, still using the dihydrofolate example, amyloid formation, and molecular "disorder". These are consistent with the view that binding and folding are similar events, with the differences stemming from different stabilities and hence population times.

Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.

How do thermophilic proteins deal with heat?

Recent years have witnessed an explosion of sequence and structural information for proteins from hyperthermophilic and thermophilic organisms. Complete genome sequences are available for many hyperthermophilic archaeons. Here, we review some recent studies on protein thermostability along with work from our laboratory. A large number of sequence and structural factors are thought to contribute toward higher intrinsic thermal stability of proteins from these organisms. The most consistent are surface loop deletion, increased occurrence of hydrophobic residues with branched side chains and an increased proportion of charged residues at the expense of uncharged polar residues. The energetic contribution of electrostatic interactions such as salt bridges and their networks toward protein stability can be stabilizing or destabilizing. For hyperthermophilic proteins, the contribution is mostly stabilizing. Macroscopically, improvement in electrostatic interactions and strengthening of hydrophobic cores by branched apolar residues increase the enthalpy change between the folded and unfolded states of a thermophilic protein. At the same time, surface loop deletion contributes to decreased conformational entropy and decreased heat capacity change between the folded and unfolded states of the protein.

Structured disorder and conformational selection.

Traditionally, molecular disorder has been viewed as local or global instability. Molecules or regions displaying disorder have been considered inherently unstructured. The term has been routinely applied to cases for which no atomic coordinates can be derived from crystallized molecules. Yet, even when it appears that the molecules are disordered, prevailing conformations exist, with population times higher than those of all alternate conformations. Disordered molecules are the outcome of rugged energy landscapes away from the native state around the bottom of the funnel. Ruggedness has a biological function, creating a distribution of structured conformers that bind via conformational selection, driving association and multimolecular complex formation, whether chain-linked in folding or unlinked in binding. We classify disordered molecules into two types. The first type possesses a hydrophobic core. Here, even if the native conformation is unstable, it still has a large enough population time, enabling its experimental detection. In the second type, no such hydrophobic core exists. Hence, the native conformations of molecules belonging to this category have shorter population times, hindering their experimental detection. Although there is a continuum of distribution of hydrophobic cores in proteins, an empirical, statistically based hydrophobicity function may be used as a guideline for distinguishing the two disordered molecule types. Furthermore, the two types relate to steps in the protein folding reaction. With respect to protein design, this leads us to propose that engineering-optimized specific electrostatic interactions to avoid electrostatic repulsion would reduce the type I disordered state, driving the molten globule (MG) --> native (N) state. In contrast, for overcoming the type II disordered state, in addition to specific interactions, a stronger hydrophobic core is also indicated, leading to the denatured --> MG --> N state.

MUSTA--a general, efficient, automated method for multiple structure alignment and detection of common motifs: application to proteins.

Here we present an algorithm designed to carry out multiple structure alignment and to detect recurring substructural motifs. So far we have implemented it for comparison of protein structures. However, this general method is applicable to comparisons of RNA structures and to detection of a pharmacophore in a series of drug molecules. Further, its sequence order independence permits its application to detection of motifs on protein surfaces, interfaces, and binding/active sites. While there are many methods designed to carry out pairwise structure comparisons, there are only a handful geared toward the multiple structure alignment task. Most of these tackle multiple structure comparison as a collection of pairwise structure comparison tasks. The multiple structural alignment algorithm presented here automatically finds the largest common substructure (core) of atoms that appears in all the molecules in the ensemble. The detection of the core and the structural alignment are done simultaneously. The algorithm begins by finding small substructures that are common to all the proteins in the ensemble. One of the molecules is considered the reference; the others are the source molecules. The small substructures are stored in special arrays termed combinatorial buckets, which define sets of multistructural alignments from the source molecules that coincide with the same small set of reference atoms (C(alpha)-atoms here). These substructures are initial small fragments that have congruent copies in each of the proteins. The substructures are extended, through the processing of the combinatorial buckets, by clustering the superpositions (transformations). The method is very efficient.

Protein functional epitopes: hot spots, dynamics and combinatorial libraries.

Recent studies increasingly point to the importance of structural flexibility and plasticity in proteins, highlighting the evolutionary advantage. There are an increasing number of cases in which given, presumably specific, binding sites have been shown to bind a range of ligands with different compositions and shapes. These studies have also revealed that evolution tends to find convergent solutions for stable intermolecular associations, largely via conservation of polar residues as hot spots of binding energy. On the other hand, the ability to bind multiple ligands at a given site is largely derived from hinge-based motions. The consideration of these two factors in functional epitopes allows more realism and robustness in the description of protein binding surfaces and, as such, in applications to mutants, modeled structures and design. Efficient multiple structure comparison and hinge-bending structure comparison tools enable the construction of combinatorial binding epitope libraries.

Fluctuations in ion pairs and their stabilities in proteins.

This report investigates the effect of systemic protein conformational flexibility on the contribution of ion pairs to protein stability. Toward this goal, we use all NMR conformer ensembles in the Protein Data Bank (1) that contain at least 40 conformers, (2) whose functional form is monomeric, (3) that are nonredundant, and (4) that are large enough. We find 11 proteins adhering to these criteria. Within these proteins, we identify 22 ion pairs that are close enough to be classified as salt bridges. These are identified in the high-resolution crystal structures of the respective proteins or in the minimized average structures (if the crystal structures are unavailable) or, if both are unavailable, in the "most representative" conformer of each of the ensembles. We next calculate the electrostatic contribution of each such ion pair in each of the conformers in the ensembles. This results in a comprehensive study of 1,201 ion pairs, which allows us to look for consistent trends in their electrostatic contributions to protein stability in large sets of conformers. We find that the contributions of ion pairs vary considerably among the conformers of each protein. The vast majority of the ion pairs interconvert between being stabilizing and destabilizing to the structure at least once in the ensembles. These fluctuations reflect the variabilities in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other, and with respect to other charged groups in the remainder of the protein. The higher crystallographic B-factors for the respective side-chains are consistent with these fluctuations. The major conclusion from this study is that salt bridges observed in crystal structure may break, and new salt bridges may be formed. Hence, the overall stabilizing (or, destabilizing) contribution of an ion pair is conformer population dependent.

Automated multiple structure alignment and detection of a common substructural motif.

While a number of approaches have been geared toward multiple sequence alignments, to date there have been very few approaches to multiple structure alignment and detection of a recurring substructural motif. Among these, none performs both multiple structure comparison and motif detection simultaneously. Further, none considers all structures at the same time, rather than initiating from pairwise molecular comparisons. We present such a multiple structural alignment algorithm. Given an ensemble of protein structures, the algorithm automatically finds the largest common substructure (core) of C(alpha) atoms that appears in all the molecules in the ensemble. The detection of the core and the structural alignment are done simultaneously. Additional structural alignments also are obtained and are ranked by the sizes of the substructural motifs, which are present in the entire ensemble. The method is based on the geometric hashing paradigm. As in our previous structural comparison algorithms, it compares the structures in an amino acid sequence order-independent way, and hence the resulting alignment is unaffected by insertions, deletions and protein chain directionality. As such, it can be applied to protein surfaces, protein-protein interfaces and protein cores to find the optimally, and suboptimally spatially recurring substructural motifs. There is no predefinition of the motif. We describe the algorithm, demonstrating its efficiency. In particular, we present a range of results for several protein ensembles, with different folds and belonging to the same, or to different, families. Since the algorithm treats molecules as collections of points in three-dimensional space, it can also be applied to other molecules, such as RNA, or drugs.

Protein folding and function: the N-terminal fragment in adenylate kinase.

Three-dimensional protein folds range from simple to highly complex architectures. In complex folds, some building block fragments are more important for correct protein folding than others. Such fragments are typically buried in the protein core and mediate interactions between other fragments. Here we present an automated, surface area-based algorithm that is able to indicate which, among all local elements of the structure, is critical for the formation of the native fold, and apply it to structurally well-characterized proteins. In particular, we focus on adenylate kinase. The fragment containing the phosphate binding, P-loop (the "giant anion hole") flanked by a beta-strand and an alpha-helix near the N-terminus, is identified as a critical building block. This building block shows a high degree of sequence and structural conservation in all adenylate kinases. The results of our molecular dynamics simulations are consistent with this identification. In its absence, the protein flips to a stable, non-native state. In this misfolded conformation, the other local elements of the structure are in their native-like conformations; however, their association is non-native. Furthermore, this element is critically important for the function of the enzyme, coupling folding, and function.

A proposed structural model for amyloid fibril elongation: domain swapping forms an interdigitating beta-structure polymer.

We propose a model illustrating how proteins, which differ in their overall sequences and structures, can form the propagating, twisted beta-sheet conformations, characteristic of amyloids. Some cases of amyloid formation can be explained through a "domain swapping" event, where the swapped segment is either a beta-hairpin or an unstable conformation which can partially unfold and assume a beta-hairpin structure. As in domain swapping, here the swapped beta-hairpin is at the edge of the structure, has few (if any) salt bridges and hydrogen bonds connecting it to the remainder of the structure and variable extents of buried non-polar surface areas. Additionally, in both cases the swapped piece constitutes a transient "building block" of the structure, with a high population time. Whereas in domain swapping the swapped fragment has been shown to be an alpha-helix, loop, strand or an entire domain, but so far not a beta-hairpin, despite the large number of cases in which it was already detected, here swapping may involve such a structural motif. We show how the swapping of beta-hairpins would form an interdigitated, twisted beta-sheet conformation, explaining the remarkable high stability of the protofibril in vitro. Such a swapping mechanism is attractive as it involves a universal mechanism in proteins, critical for their function, namely hinge-bending motions. Our proposal is consistent with structural superpositioning of mutational variants. While the overall r.m.s.d.s of the wild-type and mutants are small, the proposed hinge-bending region consistently shows larger deviations. These larger deviations illustrate that this region is more prone to respond to the mutational changes, regardless of their location in the sequence or in the structure. Nevertheless, above all, we stress that this proposition is hypothetical, since it is based on assumptions lacking definitive experimental support.

Molecular dynamics simulation of Escherichia coli dihydrofolate reductase and its protein fragments: relative stabilities in experiment and simulations.

We have carried out molecular dynamics simulations of the native dihydrofolate reductase from Escherichia coli and several of its folded protein fragments at standard temperature. The simulations have shown fragments 1--36, 37--88, and 89--159 to be unstable, with a C(alpha)RMSD (C(alpha) root mean squared deviation) >5 A after 3.0 nsec of simulation. The unfolding of fragment 1--36 was immediate, whereas fragments 37--88 and 89--159 gradually unfolded because of the presence of the beta-sheet core structure. In the absence of residues 1--36, the two distinct domains comprising fragment 39--159 associated with each other, resulting in a stable conformation. This conformation retained most of its native structural elements. We have further simulated fragments derived from computational protein cutting. These were also found to be unstable, with the exception of fragment 104--159. In the absence of alpha(4), the loose loop region of residues 120--127 exhibited a beta-strand-like behavior, associating itself with the beta-sheet core of the protein fragment. The current study suggests that the folding of dihydrofolate reductase involves cooperative folding of distinct domains which otherwise would have been unstable as independent folded units in solution. Finally, the critical role of residues 1--36 in allowing the two distinct domains of fragment 104--159 to fold into the final native conformation is discussed.

Point mutations and sequence variability in proteins: redistributions of preexisting populations.

Here we study the effect of point mutations in proteins on the redistributions of the conformational substates. We show that regardless of the location of a mutation in the protein structure and of its type, the observed movements of the backbone recur largely at the same positions in the structures. Despite the different interactions that are disrupted and formed by the residue substitution, not only are the conformations very similar, but the regions that move are also the same, regardless of their sequential or spatial distance from the mutation. This observation leads us to conclude that, apart from some extreme cases, the details of the interactions are not critically important in determining the protein conformation or in specifying which parts of the protein would be more prone to take on different local conformations in response to changes in the sequence. This finding further illustrates why proteins manifest a robustness toward many mutational events. This nonuniform distribution of the conformer population is consistently observed in a variety of protein structural types. Topology is critically important in determining folding pathways, kinetics, building block cutting, and anatomy trees. Here we show that topology is also very important in determining which regions of the protein structure will respond to sequence changes, regardless of the sequential or spatial location of the mutation.

Transient, highly populated, building blocks folding model.

Protein folding is a hierarchical event, in which transiently formed local structural elements assemble to yield the native conformation. In principle, multiple paths glide down the energy landscape, but, in practice, only a few of the paths are highly traveled. Here, the literature is reviewed in this light, and, particularly, a hierarchical, building block protein-folding model is presented, putting it in the context of a broad range of experimental and theoretical results published over the past few years. The model is based on two premises: First, although the local building block elements may be unstable, they nevertheless have higher population times than all alternate conformations; and, second, protein folding progresses through a combinatorial assembly of these elements. Through the binding of the most favorable building block conformers, there is a redistribution of the conformers in solution, propagating the protein-folding reaction. We describe the algorithm, and illustrate its usefulness, then we focus on its utility in assigning simple vs complex folding pathways, on chaperonin-assisted folding, on its relevance to domain-swapping processes, and on its relevance and relationship to disconnectivity graphs and tree diagrams. Considering protein folding as initiating from local transient structural elements is consistent with available experimental and theoretical results. Here, we have shown that, early in the folding process, sequential interactions are likely to take place, even if the final native fold is a complex, nonsequential one. Such a route is favorable kinetically and entropically. Through the construction of anatomy trees, the model enables derivation of the major folding pathways and their bumps, and qualitatively explains the kinetics of protein folding.

Building blocks, hinge-bending motions and protein topology.

Here we show that the locations of molecular hinges in protein structures fall between building block elements. Building blocks are fragments of the protein chain which constitute local minima. These elements fold first. In the next step they associate through a combinatorial assembly process. While chain-linked building blocks may be expected to trial-associate first, if unstable, alternate more stable associations will take place. Hence, we would expect that molecular hinges will be at such inter-building block locations, or at the less stable, unassigned regions. On the other hand, hinge-bending motions are well known to be critical for protein function. Hence, protein folding and protein function are evolutionarily related. Further, the pathways through which proteins attain their three dimensional folds are determined by protein topology. However, at the same time the locations of the hinges, and hinge-bending motions are also an outcome of protein topology. Thus, protein folding and function appear coupled, and relate to protein topology. Here we provide some results illustrating such a relationship.