Biomedical subjects

Ruth Nussinov

Publications and source records attributed to Ruth Nussinov.

At least 37 records · Page 2Linked to original sources

PubMed2006

From peptide-based material science to protein fibrils: discipline convergence in nanobiology.

This paper illustrates the merits of convergence in nanobiology of two seemingly disparate fields, material science and computational biology. Traditionally, material science has been a discipline involving design and fabrication of synthetic polymers consisting of repeating units. Collaboration with synthetic organic chemists allowed design of new polymers, with a range of altered conformations. Yet, naturally occurring proteins are also materials. Their varied sequences and structures should enrich material science providing more complex shapes, scaffolds and chemical properties. For material scientists, the enhanced coverage of chemical space obtained by integrating proteins and synthetic organic chemistry through the introduction of non-natural residues allows a range of new useful potential applications.

Biocompatible Materials↗

PubMed2006

Comparison of the human and worm p53 structures suggests a way for enhancing stability.

Maintaining the native conformation is essential for the proper function of tumor suppressor protein p53. However, p53 is a low-stability protein that can easily lose its function upon structural perturbations such as those resulting from missense mutations, leading to the development of cancer. Therefore, it is important to develop strategies to design stable p53 which still maintains its normal function. Here, we compare the stabilities of the human and worm p53 core domains using molecular dynamics simulations. We find that the worm p53 is significantly more stable than the human form. Detailed analysis of the structural fluctuations shows that the stability difference lies in the peripheral structural motifs that contrast in their structural features and flexibility. The most dramatic difference in stability originates from loop L1, from the turn between helix H1 and beta-strand S5, and from the turn that connects beta-strands S7 and S8. Structural analysis shows significant differences for these motifs between the two proteins. Loop L1 lacks secondary structure, and the turns between helix H1 and strand S5 and between strands S7 and S8 are much longer in the human form p53. On the basis of these differences, we designed a mutant by shortening the turn between strands S7 and S8 to enhance the stability. Surprisingly, this mutant was very stable when probed by molecular dynamics simulations. In addition, the stabilization was not localized in the turn region. Loop L1 was also significantly stabilized. Our results show that stabilizing peripheral structural motifs can greatly enhance the stability of the p53 core domain and therefore is likely to be a viable alternative in the development of stable p53. In addition, loop- or turn-related mutants with different stabilities may also be used to probe the relationship between function, a particular structural motif, and its flexibility.

Animals↗

PubMed2006

A permissive secondary structure-guided superposition tool for clustering of protein fragments toward protein structure prediction via fragment assembly.

MOTIVATION: Secondary-Structure Guided Superposition tool (SSGS) is a permissive secondary structure-based algorithm for matching of protein structures and in particular their fragments. The algorithm was developed towards protein structure prediction via fragment assembly. RESULTS: In a fragment-based structural prediction scheme, a protein sequence is cut into building blocks (BBs). The BBs are assembled to predict their relative 3D arrangement. Finally, the assemblies are refined. To implement this prediction scheme, a clustered structural library representing sequence patterns for protein fragments is essential. To create a library, BBs generated by cutting proteins from the PDB are compared and structurally similar BBs are clustered. To allow structural comparison and clustering of the BBs, which are often relatively short with flexible loops, we have devised SSGS. SSGS maintains high similarity between cluster members and is highly efficient. When it comes to comparing BBs for clustering purposes, the algorithm obtains better results than other, non-secondary structure guided protein superimposition algorithms.

Algorithms↗

PubMed2006

Paths of lateral gene transfer of lysyl-aminoacyl-tRNA synthetases with a unique evolutionary transition stage of prokaryotes coding for class I and II varieties by the same organisms.

BACKGROUND: While the premise that lateral gene transfer (LGT) is a dominant evolutionary force is still in considerable dispute, the case for widespread LGT in the family of aminoacyl-tRNA synthetases (aaRS) is no longer contentious. aaRSs are ancient enzymes, guarding the fidelity of the genetic code. They are clustered in two structurally unrelated classes. Only lysine aminoacyl-tRNA synthetase (LysRS) is found both as a class 1 and a class 2 enzyme (LysRS1-2). Remarkably, in several extant prokaryotes both classes of the enzyme coexist, a unique phenomenon that has yet to receive its due attention. RESULTS: We applied a phylogenetic approach for determining the extent and origin of LGT in prokaryotic LysRS. Reconstructing species trees for Archaea and Bacteria, and inferring that their last common ancestors encoded LysRS1 and LysRS2, respectively, we studied the gains and losses of both classes. A complex pattern of LGT events emerged. In specific groups of organisms LysRS1 was replaced by LysRS2 (and vice versa). In one occasion, within the alpha proteobacteria, a LysRS2 to LysRS1 LGT was followed by reversal to LysRS2. After establishing the most likely LGT paths, we studied the possible origins of the laterally transferred genes. To this end, we reconstructed LysRS gene trees and evaluated the likely origins of the laterally transferred genes. While the sources of LysRS1 LGTs were readily identified, those for LysRS2 remain, for now, uncertain. The replacement of one LysRS by another apparently transits through a stage simultaneously coding for both synthetases, probably conferring a selective advantage to the affected organisms. CONCLUSION: The family of LysRSs features complex LGT events. The currently available data were sufficient for identifying unambiguously the origins of LysRS1 but not of LysRS2 gene transfers. A selective advantage is suggested to organisms encoding simultaneously LysRS1-2.

Alphaproteobacteria↗

PubMed2006

The stability of monomeric intermediates controls amyloid formation: Abeta25-35 and its N27Q mutant.

The structure and stabilities of the intermediates affect protein folding as well as misfolding and amyloid formation. By applying Kramer's theory of barrier crossing and a Morse-function-like energy landscape, we show that intermediates with medium stability dramatically increase the rate of amyloid formation; on the other hand, very stable and very unstable intermediates sharply decrease amyloid formation. Remarkably, extensive molecular dynamics simulations and conformational energy landscape analysis of Abeta25-35 and its N27Q mutant corroborate the mathematical description. Both experimental and current simulation results indicate that the core of the amyloid structure of Abeta25-35 formed from residues 28-35. A single mutation of N27Q of Abeta25-35 makes the Abeta25-35 N27Q amyloid-free. Energy landscape calculations show that Abeta25-35 has extended intermediates with medium stability that are prone to form amyloids, whereas the extended intermediates for Abeta25-35 N27Q split into stable and very unstable species that are not disposed to form amyloids. The results explain the contribution of both alpha-helical and beta-strand intermediates to amyloid formation. The results also indicate that the structure and stability of the intermediates, as well as of the native folded and the amyloid states can be targeted in drug design. One conceivable approach is to stabilize the intermediates to deter amyloid formation.

Amino Acid Substitution↗

PubMed2006

Atomic-level description of amyloid beta-dimer formation.

All-atom simulations have been carried out on a monomer and dimer of the aggregation-prone fragment (16-22) of amyloid beta peptide, which is implicated in Alzheimer's disease. The replica exchange molecular dynamics method, which has been successfully applied to peptide folding, is utilized as a means to sample the configurational space with proper Boltzmann weighting so that the structural, motional, and thermodynamic description of self-assembly can be obtained. The free energy landscape showing the delicate balance between different monomer and dimer conformations is mapped along carefully chosen reaction coordinates. The canonical ensembles at 38 different temperatures are used to describe the thermodynamics and the relative stabilities of at least six different dimer conformations including that of parallel and antiparallel orientations. We also delineate the nature of the molecular forces that activate and stabilize these different dimer conformations as a function of temperature, especially as related to secondary structural propensity of monomer. We identify parallel loop dimer conformations that are stabilized due to specific interactions with water molecules.

Amyloid beta-Peptides↗

PubMed2006

Concepts and schemes for the re-engineering of physical protein modules: generating nanodevices via targeted replacements with constrained amino acids.

Physically building complex multi-molecular structures from naturally occurring biological macromolecules has aroused a great deal of interest. Here we focus on nanostructures composed of re-engineered, natural 'foldamer' building blocks. Our aim is to provide some of the underlying concepts and schemes for crafting structures utilizing such conformationally relatively stable molecular components. We describe how, via chemical biology strategies, it is further possible to chemically manipulate the foldamer building blocks toward specific shape-driven structures, which in turn could be used toward potential-designed functions. We outline the criteria in choosing candidate foldamers from the vast biological repertoire, and how to enhance their stability through selected targeted replacements by non-proteinogenic conformationally constrained amino acids. These approaches combine bioinformatics, high performance computations and mathematics with synthetic organic chemistry. The resulting artificially engineered self-organizing molecular scale structures take advantage of nature's nanobiology toolkit and at the same time improve on it, since their new targeted function differs from that optimized by evolution. The major challenge facing nanobiology is to be able to exercise fine control over the performance of these target-specific molecular machines.

Combinatorial Chemistry Techniques↗

PubMed2006

Computational study of the fibril organization of polyglutamine repeats reveals a common motif identified in beta-helices.

The formation of fibril aggregates by long polyglutamine sequences is assumed to play a major role in neurodegenerative diseases such as Huntington. Here, we model peptides rich in glutamine, through a series of molecular dynamics simulations. Starting from a rigid nanotube-like conformation, we have obtained a new conformational template that shares structural features of a tubular helix and of a beta-helix conformational organization. Our new model can be described as a super-helical arrangement of flat beta-sheet segments linked by planar turns or bends. Interestingly, our comprehensive analysis of the Protein Data Bank reveals that this is a common motif in beta-helices (termed beta-bend), although it has not been identified so far. The motif is based on the alternation of beta-sheet and helical conformation as the protein sequence is followed from the N to the C termini (beta-alpha(R)-beta-polyPro-beta). We further identify this motif in the ssNMR structure of the protofibril of the amyloidogenic peptide Abeta(1-40). The recurrence of the beta-bend suggests a general mode of connecting long parallel beta-sheet segments that would allow the growth of partially ordered fibril structures. The design allows the peptide backbone to change direction with a minimal loss of main chain hydrogen bonds. The identification of a coherent organization beyond that of the beta-sheet segments in different folds rich in parallel beta-sheets suggests a higher degree of ordered structure in protein fibrils, in agreement with their low solubility and dense molecular packing.

Amino Acid Motifs↗

PubMed2006

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

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

Amino Acid Sequence↗

PubMed2006

Sequence and structure analysis of parallel beta helices: implication for constructing amyloid structural models.

Increasing evidence suggests that amyloids and parallel beta helices may share similar motifs. A systemic analysis of beta helices is performed to examine their sequence and structural characteristics. Ile prefers to occur in beta strands. In contrast, Pro is disfavored, compatible with the underlying assumption in Pro-scanning mutagenesis. Cys, Asn, and Phe form significant homostacking (identical amino acid interactions). Asn is highly conserved in the high-energy, left-handed alpha-helical conformation, where it frequently forms amide stacking. Based on the observed prominent stacking of chemically similar residues in parallel beta helices, we propose that within the "cross-beta" framework, amyloids with longer peptide chains may have common structural features of in-register, parallel alignment, with the side chains forming identical amino acid ladders. The requirement of ladder formation limits the combinations of side chain interactions. Such a limit combined with environmental conditions (e.g., pH, concentration) could be a major reason for the ability of most polypeptides to form amyloids.

Amino Acid Sequence↗

PubMed2006

De novo tubular nanostructure design based on self-assembly of beta-helical protein motifs.

We present an approach for designing self-assembled nanostructures from naturally occurring building block segments obtained from native protein structures. We focus on structural motifs from left-handed beta-helical proteins. We selected 17 motifs. Copies of each of the motifs are stacked one atop the other. The obtained structures were simulated for long periods by using Molecular Dynamics to test their ability to retain their organization over time. We observed that a structural model based on the self-assembly of a motif from E. coli galactoside acetyltransferase produced a very stable tube. We studied the interactions that help maintain the conformational stability of the systems, focusing on the role of specific amino acids at specific positions. Analysis of these systems and a mutational study of selected candidates revealed that the presence of proline and glycine residues in the loops of beta-helical structures greatly enhances the structural stability of the systems.

Acetyltransferases↗

PubMed2006

"Similarity trap" in protein-protein interactions could be carcinogenic: simulations of p53 core domain complexed with 53BP1 and BRCA1 BRCT domains.

Similar binding sites often imply similar protein-protein interactions and similar functions; however, similar binding sites may also constitute traps for nonfunctional associations. How are similar sites distinguished to prevent misassociations? BRCT domains from breast cancer-susceptibility gene product BRCA1 and protein 53BP1 have similar structures yet different binding behaviors with p53 core domain. 53BP1-BRCT domain forms a stable complex with p53. In contrast, BRCA1-p53 interaction is weak or other mechanisms operate. To delineate the difference, we designed 13 BRCA1-BRCT mutants and computationally investigated the structural and stability changes compared to the experimental p53-53BP1 structure. Interestingly, of the 13, the 2 mutations that are cancerous and involve nonconserved residues are those that enforced p53 core domain binding with BRCA1-BRCT in a way similar to p53-53BP1 binding. Hence, falling into the "similarity trap" may disrupt normal BRCA1 and p53 functions. Our results illustrate how this trap is avoided in the native state.

Amino Acid Sequence↗

PubMed2006

The multiple common point set problem and its application to molecule binding pattern detection.

Recognition of binding patterns common to a set of protein structures is important for recognition of function, prediction of binding, and drug design. We consider protein binding sites represented by a set of 3D points with assigned physico-chemical and geometrical properties important for protein-ligand interactions. We formulate the multiple binding site alignment problem as detection of the largest common set of such 3D points. We discuss the computational problem of multiple common point set detection and, particularly, the matching problem in K-partite-epsilon graphs, where K partitions are associated with K structures and edges are defined between epsilon-close points. We show that the K-partite-epsilon matching problem is NP-hard in the Euclidean space with dimension larger than one. Consequently, we show that the largest common point set problem between three point sets is NP-hard. On the practical side, we present a novel computational method, MultiBind, for recognition of binding patterns common to a set of protein structures. It performs a multiple alignment between protein binding sites in the absence of overall sequence, fold, or binding partner similarity. Despite the NP-hardness results, in our applications, we practically overcome the exponential number of multiple alignment combinations by applying an efficient branchand- bound filtering procedure. We show applications of MultiBind to several biological targets. The method recognizes patterns which are responsible for binding small molecules, such as estradiol, ATP/ANP, and transition state analogues.

Adenosine Triphosphate↗

PubMed2006

Association of putative concave protein-binding sites with the fluctuation behavior of residues.

Here, we propose a binding site prediction method based on the high frequency end of the spectrum in the native state of the protein structural dynamics. The spectrum is obtained using an elastic network model (GNM). High frequency vibrating (HFV) residues are determined from the fastest modes dynamics. HFV residue clusters and the associated surface patch residues are tested for their likelihood to locate at the binding interfaces using two different data sets, the Benchmark Set of mainly enzymes and antigen/antibodies and the Cluster Set of more diverse structures. The binding interface is identified to be within 7.5 A of the HFV residue clusters in the Benchmark Set and Cluster Set, for 77% and 70% of the structures, respectively. The success rate increases to 88% and 84%, respectively, by using the surface patches. The results suggest that concave binding interfaces, typically those of enzyme-binding sites, are enriched by the HFV residues. Thus, we expect that the association of HFV residues with the interfaces is mostly for enzymes. If, however, a binding region has invaginations and cavities, as in some of the antigen/antibodies and in cases in the Cluster data set, we expect it would be detected there too. This implies that binding sites possess several (inter-related) properties such as cavities, high packing density, conservation, and disposition for hotspots at binding surfaces. It further suggests that the high frequency vibrating residue-based approach is a potential tool for identification of regions likely to serve as protein-binding sites. The software is available at http://www.prc.boun.edu.tr/PRC/software.html.

Amino Acids↗

PubMed2005

Combinatorial docking approach for structure prediction of large proteins and multi-molecular assemblies.

Protein folding and protein binding are similar processes. In both, structural units combinatorially associate with each other. In the case of folding, we mostly handle relatively small units, building blocks or domains, that are covalently linked. In the case of multi-molecular binding, the subunits are relatively large and are associated only by non-covalent bonds. Experimentally, the difficulty in the determination of the structures of such large assemblies increases with the complex size and the number of components it contains. Computationally, the prediction of the structures of multi-molecular complexes has largely not been addressed, probably owing to the magnitude of the combinatorial complexity of the problem. Current docking algorithms mostly target prediction of pairwise interactions. Here our goal is to predict the structures of multi-unit associations, whether these are chain-connected as in protein folding, or separate disjoint molecules in the assemblies. We assume that the structures of the single units are known, either through experimental determination or modeling. Our aim is to combinatorially assemble these units to predict their structure. To address this problem we have developed CombDock. CombDock is a combinatorial docking algorithm for the structural units assembly problem. Below, we briefly describe the algorithm and present examples of its various applications to folding and to multi-molecular assemblies. To test the robustness of the algorithm, we use inaccurate models of the structural units, derived either from crystal structures of unbound molecules or from modeling of the target sequences. The algorithm has been able to predict near-native arrangements of the input structural units in almost all of the cases, suggesting that a combinatorial approach can overcome the imperfect shape complementarity caused by the inaccuracy of the models. In addition, we further show that through a combinatorial docking strategy it is possible to enhance the predictions of pairwise interactions involved in a multi-molecular assembly.

Algorithms↗

PubMed2005

Generation and analysis of a protein-protein interface data set with similar chemical and spatial patterns of interactions.

Protein-protein interfaces are regions between 2 polypeptide chains that are not covalently connected. Here, we have created a nonredundant interface data set generated from all 2-chain interfaces in the Protein Data Bank. This data set is unique, since it contains clusters of interfaces with similar shapes and spatial organization of chemical functional groups. The data set allows statistical investigation of similar interfaces, as well as the identification and analysis of the chemical forces that account for the protein-protein associations. Toward this goal, we have developed I2I-SiteEngine (Interface-to-Interface SiteEngine) [Data set available at http://bioinfo3d.cs.tau.ac.il/Interfaces; Web server: http://bioinfo3d.cs.tau.ac.il/I2I-SiteEngine]. The algorithm recognizes similarities between protein-protein binding surfaces. I2I-SiteEngine is independent of the sequence or the fold of the proteins that comprise the interfaces. In addition to geometry, the method takes into account both the backbone and the side-chain physicochemical properties of the interacting atom groups. Its high efficiency makes it suitable for large-scale database searches and classifications. Below, we briefly describe the I2I-SiteEngine method. We focus on the classification process and the obtained nonredundant protein-protein interface data set. In particular, we analyze the biological significance of the clusters and present examples which illustrate that given constellations of chemical groups in protein-protein binding sites may be preferred, and are observed in proteins with different structures and different functions. We expect that these would yield further information regarding the forces stabilizing protein-protein interactions.

Chemical Phenomena↗

PubMed2005

ARTS: alignment of RNA tertiary structures.

MOTIVATION: A fast growing number of non-coding RNAs have recently been discovered to play essential roles in many cellular processes. Similar to proteins, understanding the functions of these active RNAs requires methods for analyzing their tertiary structures. However, in contrast to the wide range of structure-based approaches available for proteins, there is still a lack of methods for studying RNA structures. RESULTS: We present a new computational method named ARTS (alignment of RNA tertiary structures). The method compares two nucleic acid structures (RNAs or DNAs) and detects a-priori unknown common substructures. These substructures can be either large global folds containing hundreds and even thousands of nucleotides or small local tertiary motifs with at least two successive base pairs. To the best of our knowledge, this is the first method of this type. The method is highly-efficient and was used to conduct an all-against-all comparison of all the RNA structures currently available in the Protein Data Bank. AVAILABILITY: The program, a web-server and supplementary information are available on http://bioinfo3d.cs.tau.ac.il/ARTS

Algorithms↗

PubMed2005

Fibril modelling by sequence and structure conservation analysis combined with protein docking techniques: beta(2)-microglobulin amyloidosis.

Obtaining atomic resolution structural models of amyloid fibrils is currently impossible, yet crucial for our understanding of the amyloid mechanism. Different pathways in the transformation of a native globular domain to an amyloid fibril invariably involve domain destabilization. Hence, locating the unstable segments of a domain is important for understanding its amyloidogenic transformation and possibly control it. Since relative conservation is suggested to relate to local stability, we performed an extensive, sequence and structure conservation analysis of the beta(2)-microglobulin (beta(2)-m) domain. Our dataset include 51 high resolution structures belonging to the "C1 set domain" family and 132 clustered PSI-BLAST search results. Segments of the beta(2)-m domain corresponding to strands A (residues 12-18), D (45-55) and G (91-95) were found to be less conserved and stable, while the central strands B (residues 22-28), C (36-41), E (62-70) and F (78-83) were found conserved and stable. Our findings are supported by accumulating observations from various experimental methods, including urea denaturation, limited proteolysis, H/D exchange and structure determination by both NMR and X-ray crystallography. We used our conservation findings together with experimental literature information to suggest a structural model for the polymerized unit of beta(2)-m. Pairwise protein docking and subsequent monomer stacking in the same manner suggest a fibril model consistent with the cross-beta structure.

Amyloid↗