An official website of the United States government
Abstract Graph theory, a branch of mathematics that focuses on the study of graphs (networks of nodes and edges), provides a robust framework for analysing the structural and functional properties of biomolecules. By leveraging molecular dynamics (MD) simulations, atoms or groups of atoms can be represented as nodes, while their dynamic interactions are depicted as edges. This network-based approach facilitates the characterization of properties such as connectivity, centrality, and modularity, which are essential for understanding the behaviour of molecular systems. This review details the application and development of graph theory-based models in studying biomolecular systems. We introduce key concepts in graph theory and demonstrate their practical applications, illustrating how innovative graph theory approaches can be employed to design biomolecular systems with enhanced functionality. Specifically, we explore the integration of graph theoretical methods with MD simulations to gain deeper insights into complex biological phenomena, such as allosteric regulation, conformational dynamics, and catalytic functions. Ultimately, graph theory has proven to be a powerful tool in the field of molecular dynamics, offering valuable insights into the structural properties, dynamics, and interactions of molecular systems. This review establishes a foundation for using graph theory in molecular design and engineering, highlighting its potential to transform the field and drive advancements in the understanding and manipulation of biomolecular systems.
more »
« less
Network Theory Analysis of Allosteric Drug-Rescue Mechanisms in the Tumor Suppressor Protein p53 Y220C Mutant
Network theory analysis has emerged as a powerful approach for investigating the complex behavior of dynamic and interactive systems, including proteomic systems. One key application of these methods is the study of long-range signaling dynamics in proteins, a phenomenon known as allostery. In this study, we applied computational models using network theory analysis to explore long-range electrostatic interactions and allosteric drug rescue mechanisms in the DNA-binding domain (DBD) of the p53 protein, a critical tumor suppressor whose dysfunction, often caused by missense mutations, is implicated in over 50% of human cancers. Using heat kernel and Wasserstein distance-based analyses, we explored the allosteric behavior of p53-DBD constructs with the Y220C mutation in the presence or absence of allosteric effector drugs. Our results demonstrated that these network theory-based protocols effectively detected the differential efficacies of small molecule allosteric effector drug compounds in restoring long-range electrostatic dynamics in the Y220C mutant. Furthermore, our approach identified key long-range electrostatic interactions critical to both the nominal and drug-rescued functionality of the p53-DBD, providing valuable insights into allosteric modulation and its therapeutic potential.
more »
« less
- Award ID(s):
- 2320718
- PAR ID:
- 10634276
- Publisher / Repository:
- MDPI
- Date Published:
- Journal Name:
- International Journal of Molecular Sciences
- Volume:
- 26
- Issue:
- 14
- ISSN:
- 1422-0067
- Page Range / eLocation ID:
- 6884
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Allosteric regulation is common in protein–protein interactions and is thus promising in drug design. Previous experimental and simulation work supported the presence of allosteric regulation in the SARS-CoV-2 spike protein. Here the route of allosteric regulation in SARS-CoV-2 spike protein is examined by all-atom explicit solvent molecular dynamics simulations, contrastive machine learning, and the Ohm approach. It was found that peptide binding to the polybasic cleavage sites, especially the one at the first subunit of the trimeric spike protein, activates the fluctuation of the spike protein's backbone, which eventually propagates to the receptor-binding domain on the third subunit that binds to ACE2. Remarkably, the allosteric regulation routes starting from the polybasic cleavage sites share a high fraction (39–67%) of the critical amino acids with the routes starting from the nitrogen-terminal domains, suggesting the presence of an allosteric regulation network in the spike protein. Our study paves the way for the rational design of allosteric antibody inhibitors.more » « less
-
Abstract Cas13a is a recent addition to the CRISPR-Cas toolkit that exclusively targets RNA, which makes it a promising tool for RNA detection. It utilizes a CRISPR RNA (crRNA) to target RNA sequences and trigger a composite active site formed by two ‘Higher Eukaryotes and Prokaryotes Nucleotide’ (HEPN) domains, cleaving any solvent-exposed RNA. In this system, an intriguing form of allosteric communication controls the RNA cleavage activity, yet its molecular details are unknown. Here, multiple-microsecond molecular dynamics simulations are combined with graph theory to decipher this intricate activation mechanism. We show that the binding of a target RNA acts as an allosteric effector, by amplifying the communication signals over the dynamical noise through interactions of the crRNA at the buried HEPN1-2 interface. By introducing a novel Signal-to-Noise Ratio (SNR) of communication efficiency, we reveal critical allosteric residues—R377, N378, and R973—that rearrange their interactions upon target RNA binding. Alanine mutation of these residues is shown to select target RNA over an extended complementary sequence beyond guide-target duplex for RNA cleavage, establishing the functional significance of these hotspots. Collectively our findings offer a fundamental understanding of the Cas13a mechanism of action and pave new avenues for the development of highly selective RNA-based cleavage and detection tools.more » « less
-
Complex mechanisms regulate the cellular distribution of cholesterol, a critical component of eukaryote membranes involved in regulation of membrane protein functions directly and through the physiochemical properties of membranes. StarD4, a member of the steroidogenic acute regulator-related lipid-transfer (StART) domain (StARD)-containing protein family, is a highly efficient sterol-specific transfer protein involved in cholesterol homeostasis. Its mechanism of cargo loading and release remains unknown despite recent insights into the key role of phosphatidylinositol phosphates in modulating its interactions with target membranes. We have used large-scale atomistic Molecular dynamics (MD) simulations to study how the dynamics of cholesterol bound to the StarD4 protein can affect interaction with target membranes, and cargo delivery. We identify the two major cholesterol (CHL) binding modes in the hydrophobic pocket of StarD4, one near S136&S147 (the Ser-mode), and another closer to the putative release gate located near W171, R92&Y117 (the Trp-mode). We show that conformational changes of StarD4 associated directly with the transition between these binding modes facilitate the opening of the gate. To understand the dynamics of this connection we apply a machine-learning algorithm for the detection of rare events in MD trajectories (RED), which reveals the structural motifs involved in the opening of a front gate and a back corridor in the StarD4 structure occurring together with the spontaneous transition of CHL from the Ser-mode of binding to the Trp-mode. Further analysis of MD trajectory data with the information-theory based NbIT method reveals the allosteric network connecting the CHL binding site to the functionally important structural components of the gate and corridor. Mutations of residues in the allosteric network are shown to affect the performance of the allosteric connection. These findings outline an allosteric mechanism which prepares the CHL-bound StarD4 to release and deliver the cargo when it is bound to the target membrane.more » « less
-
Skolnick, Jeffrey (Ed.)The link between p53 tumor suppressive functions and organismal lifespan is multifaceted. Its DNA-repair mechanism is longevity-enhancing while its role in cellular senescence pathways induces pro-aging phenotypes. To understand how p53 may regulate organismal lifespan, cross-species genotype-phenotype (GP) studies of the p53 DNA-binding domain (DBD) have been used to assess the correlation of amino acid changes to lifespan. Amino acid changes in non-DNA-binding regions such as the transactivation (TAD), proline-rich (PRD), regulatory (REG), and tetramerization (TET) are largely unexplored. In addition, existing GP correlation tools such as SigniSite do not account for phylogenetic relationships between aligned sequences in correlating genotypic differences to phenotypes such as lifespan. To identify phylogenetically significant, longevity-correlated residues in full-length p53 alignments, we developed a Python- and R-based workflow, Relative Evolutionary Scoring (RES). While RES-predicted longevity-associated residues (RPLARs) are concentrated primarily in the DBD, the PRD, TET, and REG domains also house RPLARs. While yeast functional assay enrichment reveals that RPLARs may be dispensable for p53-mediated transactivation, PEPPI and Rosetta-based protein-protein interaction prediction suggests a role for RPLARs in p53 stability and interaction interfaces of tumor suppressive protein-protein complexes. With experimental validation of the RPLARs’ roles in p53 stability, transactivation, and involvement in senescence-regulatory pathways, we can gain crucial insights into mechanisms underlying dysregulated tumor suppression and accelerated aging.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.3390/ijms26146884
