Biomedical subjects

Ruth Nussinov

Publications and source records attributed to Ruth Nussinov.

At least 19 recordsLinked to original sources

PubMed2025

Mutations in tumor signaling, metastases, and synthetic lethality establish distinct patterns.

Effective identification of oncogenic mutations is essential for diagnosis, forecasting resistance, and metastasis in remission. It is required for an optimal drug regimen. We develop a framework to discover mutations that co-exist in different oncoproteins, and those that are excluded, likely encoding oncogene-induced senescence. First, mapping the proteins onto pathways assists combinatorial drug selections and helps to detect metastases. Second, it provides the molecular basis for synthetic lethality, to date investigated at the genome level. Our pan-cancer profiles of ~60,000 tumor sequences, detect 3424 co-existing tumor-specific mutations. Mapping them onto pathways indicates that they preferentially promote specific primary tumors. We uncover metastatic mutations and provide metastatic breast-cancer markers. This work not only clarifies the mechanistic basis of intratumor mutational diversity but usefully reveals markers for metastasis in patients' genomes and introduces a novel computational framework for detecting metastasis based on tumor mutational profiles. Mapping the mutations onto pathways provides an invaluable metastasis-targeting resource, guiding drug combinations.

Humans↗

PubMed2025

Deep learning-based cell-specific gene regulatory networks inferred from single-cell multiome data.

Gene regulatory networks (GRNs) provide a global representation of how genetic/genomic information is transferred in living systems and are a key component in understanding genome regulation. Single-cell multiome data provide unprecedented opportunities to reconstruct GRNs at fine-grained resolution. However, the inference of GRNs is hindered by insufficient single omic profiles due to the characteristic high loss rate of single-cell sequencing data. In this study, we developed scMultiomeGRN, a deep learning framework to infer transcription factor (TF) regulatory networks via unique integration of single-cell genomic (single-cell RNA sequencing) and epigenomic (single-cell ATAC sequencing) data. We create scMultiomeGRN to elucidate these networks by conceptualizing TF network graph structures. Specifically, we build modality-specific neighbor aggregators and cross-modal attention modules to learn latent representations of TFs from single-cell multi-omics. We demonstrate that scMultiomeGRN outperforms state-of-the-art models on multiple benchmark datasets involved in diseases and health. Via scMultiomeGRN, we identified Alzheimer's disease-relevant regulatory network of SPI1 and RUNX1 for microglia. In summary, scMultiomeGRN offers a deep learning framework to identify cell type-specific gene regulatory network from single-cell multiome data.

Deep Learning↗

PubMed2025

Pioneer in Molecular Biology: Conformational Ensembles in Molecular Recognition, Allostery, and Cell Function.

In 1978, for my PhD, I developed the efficient O(n3) dynamic programming algorithm for the-then open problem of RNA secondary structure prediction. This algorithm, now dubbed the "Nussinov algorithm", "Nussinov plots", and "Nussinov diagrams", is still taught across Europe and the U.S. As sequences started coming out in the 1980s, I started seeking genome-encoded functional signals, later becoming a bioinformatics trend. In the early 1990s I transited to proteins, co-developing a powerful computer vision-based docking algorithm. In the late 1990s, I proposed the foundational role of conformational ensembles in molecular recognition and allostery. At the time, conformational ensembles and free energy landscapes were viewed as physical properties of proteins but were not associated with function. The classical view of molecular recognition and binding was based on only two conformations captured by crystallography: open and closed. I proposed that all conformational states preexist. Proteins always have not one folded form-nor two-but many folded forms. Thus, rather than inducing fit, binding can work by shifting the ensembles between states, and this shifting, or redistributing the ensembles to maintain equilibrium, is the origin of the allosteric effect and protein, thus cell, function. This transformative paradigm impacted community views in allosteric drug design, catalysis, and regulation. Dynamic conformational ensemble shifts are now acknowledged as the origin of recognition, allostery, and signaling, underscoring that conformational ensembles-not proteins-are the workhorses of the cell, pioneering the fundamental idea that dynamic ensembles are the driving force behind cellular processes. Nussinov was recognized as pioneer in molecular biology by JMB.

Molecular Biology↗

PubMed2007

A hierarchical protein folding scheme based on the building block folding model.

The building block protein folding model states that the native protein structure is the product of a combinatorial assembly of relatively structurally independent contiguous parts of the protein that possess a hydrophobic core, i.e., building blocks (BBs). According to this model, our group proposed a three-stage scheme for a feasible time-wise semi ab-intio protein structure prediction. Given a protein sequence, at the first stage of the prediction scheme, we propose cutting the sequence into structurally assigned BBs. Next, we perform a combinatorial assembly and attempt to predict the relative three-dimensional arrangement of the BBs. In the third stage, we refine and rank the assemblies. The scheme has proven to be very promising in reducing the complexity of the protein folding problem and gaining insight into the protein folding process. In this chapter, we describe the different stages of the scheme and discuss a possible application of the model to protein design.

Algorithms↗

PubMed2006

EMatch: an efficient method for aligning atomic resolution subunits into intermediate-resolution cryo-EM maps of large macromolecular assemblies.

Structural analysis of biological machines is essential for inferring their function and mechanism. Nevertheless, owing to their large size and instability, deciphering the atomic structure of macromolecular assemblies is still considered as a challenging task that cannot keep up with the rapid advances in the protein-identification process. In contrast, structural data at lower resolution is becoming more and more available owing to recent advances in cryo-electron microscopy (cryo-EM) techniques. Once a cryo-EM map is acquired, one of the basic questions asked is what are the folds of the components in the assembly and what is their configuration. Here, a novel knowledge-based computational method, named EMatch, towards tackling this task for cryo-EM maps at 6-10 A resolution is presented. The method recognizes and locates possible atomic resolution structural homologues of protein domains in the assembly. The strengths of EMatch are demonstrated on a cryo-EM map of native GroEL at 6 A resolution.

Algorithms↗

PubMed2006

Structural basis for p53 binding-induced DNA bending.

Specific p53 binding-induced DNA bending has important biological implications such as transcription activation. However, the detailed structures of the bent DNA and the p53-DNA complex are still unavailable, hampering our understanding of the mechanism for p53-induced DNA bending and its consequent biological significance. To gain insight into the p53 binding-induced DNA bending, we performed molecular dynamics simulations on DNA segments with the consensus sequence for p53-specific binding, half site DNA-p53 complexes, and full site DNA-p53 complexes. We show that each DNA-bound p53 core domain caused a local DNA conformational change within the quarter site; upon the binding of the p53 dimer, there was an apparent DNA bending at the center of the half site; when bound with two p53 dimers, the full site DNAs with two different sequences bent 20 and 35 degrees, respectively. These results are in agreement with experimental observations. Our simulations demonstrate that the two p53 dimers favored a staggered conformation in which they make favorable interactions at the interface. This dimer-dimer interface organization necessitated conformational changes in the DNA, leading to the bending at the center of the full site, which in turn is dependent on the DNA sequence. Overall, our results provide the detailed atomic model for the DNA-p53 tetramer complex and delineate the roles of DNA-p53, p53 dimer-dimer interactions, and DNA sequence in specific p53 binding-induced DNA conformational changes.

Binding Sites↗

PubMed2006

Backbone conformational preferences and pseudorotational ring puckering of 1-aminocyclopentane-1-carboxylic acid.

We have used quantum mechanical calculations at the B3LYP/6-311G(d,p) level to determine the conformational preferences of the N-acetyl-N'-methylamide derivative of 1-aminocyclopentane-1-carboxylic acid in the gas phase, chloroform solution, and water solution. The backbone conformation of this dipeptide has been described through the dihedral angles varphi and psi, while the pseudorotational phase angle was used to define the conformation of the cyclopentane ring. Results indicate that the backbone flexibility of this amino acid is restricted by the cyclic nature of the side chain, the relative stability of the different conformations depending on the polarity of the environment. The potential energy of the pseudorotation was also studied as a function of the backbone conformation. Interestingly, the conformation of the cyclic side chain depends on the backbone arrangement. Furthermore, the number of pseudorotational states accessible at room temperature is high in all the investigated environments, especially in aqueous solution. Finally, a set of force-field parameters for classical molecular mechanics calculations was developed for the investigated amino acid. Molecular dynamics simulations in both chloroform and aqueous solutions were performed to demonstrate the reliability of such parameters.

Computer Simulation↗

PubMed2006

Consensus features in amyloid fibrils: sheet-sheet recognition via a (polar or nonpolar) zipper structure.

Amyloid fibrils characterized as highly intractable thread-like species are associated with many neurodegenerative diseases. Although neither the mechanism of amyloid formation nor the origin of amyloid toxicity is currently completely understood, the detailed three-dimensional atomic structures of the yeast protein Sup35 and Abeta amyloid protein determined by recent experiments provide the first and important step towards the comprehension of the pathogenesis and aggregation mechanisms of amyloid diseases. By analyzing these two amyloid peptides which have available crystal structures and other amyloid sequences with proposed structures using computational simulations, we delineate three common features in amyloid organizations and amyloid structures. These could contribute to an improved understanding of the molecular mechanism of amyloid formation, the nature of the aggregation driving forces that stabilize these structures and the development of potential therapeutic agents against amyloid diseases.

Amino Acid Motifs↗

PubMed2006

How different are structurally flexible and rigid binding sites? Sequence and structural features discriminating proteins that do and do not undergo conformational change upon ligand binding.

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue-residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic-hydrophobic, aromatic-aromatic, and hydrophobic-polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (phi, psi) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the L-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of phi, psi changes are in general minimal, although some large changes particularly between right-handed alpha-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein-ligand complexes and for template-based prediction of protein function.

Amino Acids↗

PubMed2006

Simulations as analytical tools to understand protein aggregation and predict amyloid conformation.

Computational tools are increasingly being applied to solve the protein aggregation problem, providing insight into amyloid structures and aggregation mechanisms. The paradigm of Abeta amyloid structure elucidation provides an example of an innovative experimental design and endeavor, echoing the computational testing of possible molecular associations, all reflected in the current Ma-Nussinov-Tycko model of the Abeta amyloid. Simulations have shown that dimer formation can lock some misfolded conformations, and catalyze the shift of the equilibrium away from the native state. In most cases, a stable amyloid seed requires at least two-layered beta-sheets with properly registered side-chains. Under kinetic control, the final protein aggregations are the outcome of maximizing the van der Waals interactions between side chains and backbone hydrogen bonds.

Amyloid↗

PubMed2006

Structure by design: from single proteins and their building blocks to nanostructures.

Nanotechnology realizes the advantages of naturally occurring biological macromolecules and their building-block nature for design. Frequently, assembly starts with the choice of a "good" molecule that is synthetically optimized towards the desired shape. By contrast, we propose starting with a pre-specified nanostructure shape, selecting candidate protein building blocks from a library and mapping them onto the shape and, finally, testing the stability of the construct. Such a shape-based, part-assembly strategy is conceptually similar to protein design through the combinatorial assembly of building blocks. If the conformational preferences of the building blocks are retained and their interactions are favorable, the nanostructure will be stable. The richness of the conformations, shapes and chemistries of the protein building blocks suggests a broad range of potential applications; at the same time, it also highlights their complexity. In this Opinion article, we focus on the first step: validating such a strategy against experimental data.

Drug Design↗

PubMed2006

Residue centrality, functionally important residues, and active site shape: analysis of enzyme and non-enzyme families.

The representation of protein structures as small-world networks facilitates the search for topological determinants, which may relate to functionally important residues. Here, we aimed to investigate the performance of residue centrality, viewed as a family fold characteristic, in identifying functionally important residues in protein families. Our study is based on 46 families, including 29 enzyme and 17 non-enzyme families. A total of 80% of these central positions corresponded to active site residues or residues in direct contact with these sites. For enzyme families, this percentage increased to 91%, while for non-enzyme families the percentage decreased substantially to 48%. A total of 70% of these central positions are located in catalytic sites in the enzyme families, 64% are in hetero-atom binding sites in those families binding hetero-atoms, and only 16% belong to protein-protein interfaces in families with protein-protein interaction data. These differences reflect the active site shape: enzyme active sites locate in surface clefts, hetero-atom binding residues are in deep cavities, while protein-protein interactions involve a more planar configuration. On the other hand, not all surface cavities or clefts are comprised of central residues. Thus, closeness centrality identifies functionally important residues in enzymes. While here we focus on binding sites, we expect to identify key residues for the integration and transmission of the information to the rest of the protein, reflecting the relationship between fold and function. Residue centrality is more conserved than the protein sequence, emphasizing the robustness of protein structures.

Amino Acid Sequence↗

PubMed2006

Interaction of protegrin-1 with lipid bilayers: membrane thinning effect.

Protegrins (PG) are important in defending host tissues, preventing infection via an attack on the membrane surface of invading microorganisms. Protegrins have powerful antibiotic abilities, but the molecular-level mechanisms underlying the interactions of their beta-sheet motifs with the membrane are not known. Protegrin-1 (PG-1) is composed of 18 amino acids with a high content of basic residues and two disulfide bonds. Here we focused on the stability of PG-1 at the amphipathic interface in lipid bilayers and on the details of the peptide-membrane interactions. We simulated all-atom models of the PG-1 monomer with explicit water and lipid bilayers composed of both homogeneous POPC (palmitoyl-oleyl-phosphatidylcholine) lipids and a mixture of POPC/POPG (palmitoyl-oleyl-phosphatidylglycerol) (4:1) lipids. We observed that local thinning of the lipid bilayers mediated by the peptide is enhanced in the lipid bilayer containing POPG, consistent with experimental results of selective membrane targeting. The beta-hairpin motif of PG-1 is conserved in both lipid settings, whereas it is highly bent in aqueous solution. The conformational dynamics of PG-1, especially the highly charged beta-hairpin turn region, are found to be mostly responsible for disturbing the membrane. Even though the eventual membrane disruption requires PG-1 oligomers, our simulations clearly show the first step of the monomeric effects. The thinning effects in the bilayer should relate to pore/channel formation in the lipid bilayer and thus be responsible for further defects in the membrane caused by oligomer.

Antimicrobial Cationic Peptides↗

PubMed2006

The ARTS web server for aligning RNA tertiary structures.

RNA molecules with common structural features may share similar functional properties. Structural comparison of RNAs and detection of common substructures is, thus, a highly important task. Nevertheless, the current available tools in the RNA community provide only a partial solution, since they either work at the 2D level or are suitable for detecting predefined or local contiguous tertiary motifs only. Here, we describe a web server built around ARTS, a method for aligning tertiary structures of nucleic acids (both RNA and DNA). ARTS receives a pair of 3D nucleic acid structures and searches for a priori unknown common substructures. The search is truly 3D and irrespective of the order of the nucleotides on the chain. The identified common substructures can be large global folds with hundreds and even thousands of nucleotides as well as small local motifs with at least two successive base pairs. The method is highly efficient and has been used to conduct an all-against-all comparison of all the RNA structures in the Protein Data Bank. The web server together with a software package for download are freely accessible at http://bioinfo3d.cs.tau.ac.il/ARTS.

Internet↗

PubMed2006

Correlated mutations: advances and limitations. A study on fusion proteins and on the Cohesin-Dockerin families.

Correlated mutations have been repeatedly exploited for intramolecular contact map prediction. Over the last decade these efforts yielded several methods for measuring correlated mutations. Nevertheless, the application of correlated mutations for the prediction of intermolecular interactions has not yet been explored. This gap is due to several obstacles, such as 3D complexes availability, paralog discrimination, and the availability of sequence pairs that are required for inter- but not intramolecular analyses. Here we selected for analysis fusion protein families that bypass some of these obstacles. We find that several correlated mutation measurements yield reasonable accuracy for intramolecular contact map prediction on the fusion dataset. However, the accuracy level drops sharply in intermolecular contacts prediction. This drop in accuracy does not occur always. In the Cohesin-Dockerin family, reasonable accuracy is achieved in the prediction of both intra- and intermolecular contacts. The Cohesin-Dockerin family is well suited for correlated mutation analysis. Because, however, this family constitutes a special case (it has radical mutations, has domain repeats, within each species each Dockerin domain interacts with each Cohesin domain, see below), the successful prediction in this family does not point to a general potential in using correlated mutations for predicting intermolecular contacts. Overall, the results of our study indicate that current methodologies of correlated mutations analysis are not suitable for large-scale intermolecular contact prediction, and thus cannot assist in docking. With current measurements, sequence availability, sequence annotations, and underdeveloped sequence pairing methods, correlated mutations can yield reasonable accuracy only for a handful of families.

Amino Acids↗

PubMed2006

Structural stability and dynamics of an amyloid-forming peptide GNNQQNY from the yeast prion sup-35.

A seven amino acid yeast prion sup-35 fragment (GNNQQNY) forms amyloid fibrils. The availability of its detailed atomic oligomeric structure makes it a good model for studying the early stage of aggregation. Here we perform long all-atom explicit solvent molecular simulations of various sizes and arrangements of oligomer seeds of the wild-type and its mutants to study its stability and dynamics. Previous studies have suggested that the early stage rate-limiting step of oligomer formation occurs in high-order oligomers. Our simulations show that with the increase in the number of strands even from a dimer to a trimer, oligomer stability increases dramatically. This suggests that the minimal nucleus seed for GNNQQNY fibril formation could be small and is likely three or four peptides, in agreement with experiment, and that higher-order oligomers do not dissociate quickly since they have small diffusion coefficients and thus slow kinetics. Further, for the hydrophilic polar GNNQQNY, there are no hydrogen bonds and no hydrophobic interactions between adjacent beta-sheets. Simulations suggest that within the sheet, the driving forces to associate and stabilize are interstrand backbone-backbone and side chain-side chain hydrogen bonds, whereas between the sheets, shape-complementary by the dry polar steric zipper via the side chains of Asn-2, Gln-4, and Asn-6 holds the sheets together, as proposed in an earlier study. Since the polar side chains of Asn-2, Gln-4, and Asn-6 act as a hook to bind two neighboring sheets together, these geometric restraints reduce the conformational search for the correct side chain packing to a two-dimensional problem of intersheet side chain interactions. Mutant simulations show that substitution of Asn-2, Gln-4, or Asn-6 by Ala would disrupt this steric zipper, leading to unstable oligomers.

Amino Acid Sequence↗

PubMed2006

Residues crucial for maintaining short paths in network communication mediate signaling in proteins.

Here, we represent protein structures as residue interacting networks, which are assumed to involve a permanent flow of information between amino acids. By removal of nodes from the protein network, we identify fold centrally conserved residues, which are crucial for sustaining the shortest pathways and thus play key roles in long-range interactions. Analysis of seven protein families (myoglobins, G-protein-coupled receptors, the trypsin class of serine proteases, hemoglobins, oligosaccharide phosphorylases, nuclear receptor ligand-binding domains and retroviral proteases) confirms that experimentally many of these residues are important for allosteric communication. The agreement between the centrally conserved residues, which are key in preserving short path lengths, and residues experimentally suggested to mediate signaling further illustrates that topology plays an important role in network communication. Protein folds have evolved under constraints imposed by function. To maintain function, protein structures need to be robust to mutational events. On the other hand, robustness is accompanied by an extreme sensitivity at some crucial sites. Thus, here we propose that centrally conserved residues, whose removal increases the characteristic path length in protein networks, may relate to the system fragility.

Allosteric Regulation↗

PubMed2006

Designing a nanotube using naturally occurring protein building blocks.

Here our goal is to carry out nanotube design using naturally occurring protein building blocks. Inspection of the protein structural database reveals the richness of the conformations of proteins, their parts, and their chemistry. Given target functional protein nanotube geometry, our strategy involves scanning a library of candidate building blocks, combinatorially assembling them into the shape and testing its stability. Since self-assembly takes place on time scales not affordable for computations, here we propose a strategy for the very first step in protein nanotube design: we map the candidate building blocks onto a planar sheet and wrap the sheet around a cylinder with the target dimensions. We provide examples of three nanotubes, two peptide and one protein, in atomistic model detail for which there are experimental data. The nanotube models can be used to verify a nanostructure observed by low-resolution experiments, and to study the mechanism of tube formation.

HIV-1↗