Chromosome Replication Research Articles

Chromosomal integrons are genetically and functionally isolated units of genomes.

Integrons are genetic elements that increase the evolvability of bacteria by capturing new genes and stockpiling them in arrays. Sedentary chromosomal integrons (SCIs) can be massive and highly stabilized structures encoding hundreds of genes, whose function remains generally unknown. SCIs have co-evolved with the host for aeons and are highly intertwined with their physiology from a mechanistic point of view. But, paradoxically, other aspects, like their variable content and location within the genome, suggest a high genetic and functional independence. In this work, we have explored the connection of SCIs to their host genome using as a model the Superintegron (SI), a 179-cassette long SCI in the genome of Vibrio cholerae N16961. We have relocated and deleted the SI using SeqDelTA, a novel method that allows to counteract the strong stabilization conferred by toxin-antitoxin systems within the array. We have characterized in depth the impact in V. cholerae's physiology, measuring fitness, chromosome replication dynamics, persistence, transcriptomics, phenomics, natural competence, virulence and resistance against protist grazing. The deletion of the SI did not produce detectable effects in any condition, proving that-despite millions of years of co-evolution-SCIs are genetically and functionally isolated units of genomes.

Pathological R-loops in bacteria from engineered expression of endogenous antisense RNAs whose synthesis is ordinarily terminated by Rho.

In many bacteria, the essential factors Rho and NusG mediate termination of synthesis of nascent transcripts (including antisense RNAs) that are not being simultaneously translated. It has been proposed that in Rho's absence toxic RNA-DNA hybrids (R-loops) may be generated from nascent untranslated transcripts, and genome-wide mapping studies in Escherichia coli have identified putative loci of R-loop formation from more than 100 endogenous antisense transcripts that are synthesized only in a Rho-deficient strain. Here we provide evidence that engineered expression in wild-type E. coli of several such individual antisense regions on a plasmid or the chromosome generates R-loops that, in an RNase H-modulated manner, serve to disrupt genome integrity. Rho inhibition was associated with increased prevalence of antisense R-loops also in Xanthomonas oryzae pv. oryzae and Caulobacter crescentus. Our results confirm the essential role of Rho in several bacterial genera for prevention of toxic R-loops from pervasive yet cryptic endogenous antisense transcripts. Engineered antisense R-looped regions may be useful for studies on both site-specific impediments to bacterial chromosomal replication and the mechanisms of their resolution.

NMR assignment of the conserved bacterial DNA replication protein DnaA domain IV

Chromosomal replication is a ubiquitous and essential cellular process. In bacteria, the master replication initiator DnaA plays a key role in promoting an open complex at the origin (oriC) and recruiting helicase in a tightly regulated process. The C-terminal domain IV specifically recognises consensus sequences of double-stranded DNA in oriC, termed DnaA-boxes, thereby facilitating the initial engagement of DnaA to oriC. Here, we report the 13Cβ and backbone 1H, 15N, and 13C chemical shift assignments of soluble DnaA domain IV from Bacillus subtilis at pH 7.6 and 298 K.

Reconstruction of a robust bacterial replication module.

Chromosomal DNA replication is a fundamental process of life, involving the assembly of complex machinery and dynamic regulation. In this study, we reconstructed a bacterial replication module (pRC) by artificially clustering 23 genes involved in DNA replication and sequentially deleting these genes from their naturally scattered loci on the chromosome of Escherichia coli. The integration of pRC into the chromosome, moving from positions farther away to close to the replication origin, leads to an enhanced efficiency in DNA synthesis, varying from lower to higher. Strains containing replication modules exhibited increased DNA replication by accelerating the replication fork movement and initiating chromosomal replication earlier in the replication cycle. The minimized module pRC16, containing only replisome and elongation encoding genes, exhibited chromosomal DNA replication efficiency comparable to that of pRC. The replication module demonstrated robust and rapid DNA replication, regardless of growth conditions. Moreover, the replication module is plug-and-play, and integrating it into Mb-sized extrachromosomal plasmids improves their genetic stability. Our findings indicate that DNA replication, being a fundamental life process, can be artificially reconstructed into replication functional modules. This suggests potential applications in DNA replication and the construction of synthetic modular genomes.

Structural Insight Into the Function of DnaB Helicase in Bacterial DNA Replication.

In bacteria, chromosome replication is achieved by the coordinations of more than a dozen replisome enzymes. Replication initiation protein DnaA melts DNA duplex at replication origin (oriC) and forms a replication bubble, followed by loading of helicase DnaB with the help of loader protein DnaC. Then the DnaB helicase unwinds the dsDNA and supports the priming of DnaG and the polymerizing of DNA polymerase. The DnaB helicase functions as a platform coupling unwinding, priming, and polymerizing events. The multiple roles of DnaB helicase are underlined by its distinctive architecture and dynamics conformations. In this review, we will discuss the assembling of DnaB hexamer and the conformational changes upon binding of various partners, DnaB in states of closed dilated (CD), closed constricted (CC), closed helical (CH), and open helical (OH) are discussed. These multiple interfaces among DnaB and partners are potential targets for inhibitors design and novel peptide antibiotics development.

Newton's cradle: Cell cycle regulation by two mutually inhibitory oscillators

The cell division cycle is a fundamental physiological process displaying a great degree of plasticity during the course of multicellular development. This plasticity is evident in the transition from rapid and stringently-timed divisions of the early embryo to subsequent size-controlled mitotic cycles. Later in development, cells may pause and restart proliferation in response to myriads of internal or external signals, or permanently exit the cell cycle following terminal differentiation or senescence. Beyond this, cells can undergo modified cell division variants, such as endoreplication, which increases their ploidy, or meiosis, which reduces their ploidy. This wealth of behaviours has led to numerous conceptual analogies intended as frameworks for understanding the proliferative program. Here, we aim to unify these mechanisms under one dynamical paradigm. To this end, we take a control theoretical approach to frame the cell cycle as a pair of arrestable and mutually-inhibiting, doubly amplified, negative feedback oscillators controlling chromosome replication and segregation events, respectively. Under appropriate conditions, this framework can reproduce fixed-period oscillations, checkpoint arrests of variable duration, and endocycles. Subsequently, we use phase plane and bifurcation analysis to explain the dynamical basis of these properties. Then, using a physiologically realistic, biochemical model, we show that the very same regulatory structure underpins the diverse functions of the cell cycle control network. We conclude that Newton's cradle may be a suitable mechanical analogy of how the cell cycle is regulated.

In silico study exploiting chromosomal replication initiator protein as a molecular target for new antibiotics in Shigella dysenteriae

Introduction: DNA chromosomal replication initiator protein, DnaA is a component of bacterial division machinery. Objetives: virtual screening of phytochemicals for new antibiotics at novel molecular targets. Materials and Methods: SWISS-Model tool was used in homology modeling of this protein and validation tools were also used to check for its accuracy. The evaluation tools were ERRAT, PROCHECK and Molprobity servers. The refined model by GalaxyWEB server, then employed in molecular docking experiments. A total of 300 natural products were used in virtual screening for lead compounds targeting the model by AutoDock Vina software. Results and Discussion: Corylin among then best ligands that have binding affinities lower than ADP as control. Interactions in the form of hydrogen bonds and hydrophobic interactions were analyzed by LigPlot. The best compounds were submitted into iGEMDOCK docking tool for consensus scoring approach. admetSAR 2.0 was used to predict for pharmacokinetics, interactions and toxicities of these ligands inside the human body. Conclusions: Therefore, preliminary screening of lead compounds must be accompanied by studying pharmacologic properties.

DnaA modulates the gene expression and morphology of the Lyme disease spirochete.

Lyme disease is the most prevalent tick-borne infection in the Northern Hemisphere. Borrelia burgdorferi , the causative spirochete bacteria, has been maintained in nature for millennia in a consistent enzootic cycle between Ixodes ticks and various small vertebrate hosts. During the tick's blood meal, B. burgdorferi substantially increases its replication rate, alters its repertoire of outer surface proteins, and disseminates into the new vertebrate host. Across eubacteria, DnaA is the master regulatory protein that initiates chromosomal replication and acts as a transcription factor to regulate specific pathways. Here, we describe the roles that B. burgdorferi DnaA has on the physiology and gene expression of this medically important pathogen.

Genomic analysis of the degradation gene clusters and the parABS system in Afipia sp. strain DD3 capable of utilizing 2,4-dichlorophenoxyacetic acid

Here, I report the complete genome of the soil bacterium Afipia strain DD3, the first representative of its genus to utilize the herbicide 2,4-dichlorophenoxyacetic acid (2,4-D) as its sole carbon source. The compact 2,4-D degradation gene cluster of strain DD3 is a set of cadABKC and tfdBFRDEC gene clusters on its chromosome. No tfdA homologs were detected in the genome. The parAB genes responsible for chromosome segregation were found to be close to a candidate chromosomal replication initiation region (oriC) in the chromosome of strain DD3, which is located approximately 180 kb away from dnaA. In 14 strains from five Alphaproteobacterial orders, including strain DD3, parAB was found to be clustered with the genes mnmEG and rsmG coding for tRNA- and rRNA-modifying enzymes. Binding sequences for ParB were found in the promoters of or within these genes, suggesting that parAB genes in Alphaproteobacteria are involved in regulating RNA metabolism during chromosome segregation.

MksB is a novel mycobacterial condensin that orchestrates spatiotemporal positioning of replication machinery

Condensins play important roles in maintaining bacterial chromatin integrity. In mycobacteria, three types of condensins have been characterized: a homolog of SMC and two MksB-like proteins, the recently identified MksB and EptC. Previous studies suggest that EptC contributes to defending against foreign DNA, while SMC and MksB may play roles in chromosome organization. Here, we report for the first time that the condensins, SMC and MksB, are involved in various DNA transactions during the cell cycle of Mycobacterium smegmatis (currently named Mycolicibacterium smegmatis). SMC appears to be required during the last steps of the cell cycle, where it contributes to sister chromosome separation. Intriguingly, in contrast to other bacteria, mycobacterial MksB follows replication forks during chromosome replication and hence may be involved in organizing newly replicated DNA.

The budding yeast Fkh1 Forkhead associated (FHA) domain promotes a G1-chromatin state and the activity of chromosomal DNA replication origins.

In Saccharomyces cerevisiae, the forkhead (Fkh) transcription factor Fkh1 (forkhead homolog) enhances the activity of many DNA replication origins that act in early S-phase (early origins). Current models posit that Fkh1 acts directly to promote these origins' activity by binding to origin-adjacent Fkh1 binding sites (FKH sites). However, the post-DNA binding functions that Fkh1 uses to promote early origin activity are poorly understood. Fkh1 contains a conserved FHA (forkhead associated) domain, a protein-binding module with specificity for phosphothreonine (pT)-containing partner proteins. At a small subset of yeast origins, the Fkh1-FHA domain enhances the ORC (origin recognition complex)-origin binding step, the G1-phase event that initiates the origin cycle. However, the importance of the Fkh1-FHA domain to either chromosomal replication or ORC-origin interactions at genome scale is unclear. Here, S-phase SortSeq experiments were used to compare genome replication in proliferating FKH1 and fkh1-R80A mutant cells. The Fkh1-FHA domain promoted the activity of ≈ 100 origins that act in early to mid- S-phase, including the majority of centromere-associated origins, while simultaneously inhibiting ≈ 100 late origins. Thus, in the absence of a functional Fkh1-FHA domain, the temporal landscape of the yeast genome was flattened. Origins are associated with a positioned nucleosome array that frames a nucleosome depleted region (NDR) over the origin, and ORC-origin binding is necessary but not sufficient for this chromatin organization. To ask whether the Fkh1-FHA domain had an impact on this chromatin architecture at origins, ORC ChIPSeq data generated from proliferating cells and MNaseSeq data generated from G1-arrested and proliferating cell populations were assessed. Origin groups that were differentially regulated by the Fkh1-FHA domain were characterized by distinct effects of this domain on ORC-origin binding and G1-phase chromatin. Thus, the Fkh1-FHA domain controlled the distinct chromatin architecture at early origins in G1-phase and regulated origin activity in S-phase.

The Role of the MCM2-7 Helicase Subunit MCM2 in Epigenetic Inheritance.

Recycling histone proteins from parental chromatin, a process known as parental histone transfer, is an important component in chromosome replication and is essential for epigenetic inheritance. We review recent advances in our understanding of the recycling mechanism of parental histone H3-H4 tetramers (parH3:H4tet), emphasizing the pivotal role of the DNA replisome. In particular, we highlight the function of the MCM2-7 helicase subunit Mcm2 as a histone H3-H4 tetramer chaperone. Disruption of this histone chaperone's functions affects mouse embryonic stem cell differentiation and can lead to embryonic lethality in mice, underscoring the crucial role of the replisome in maintaining epigenomic stability.

Deciphering the Relationship between Cell Growth and Cell Cycle in Individual Escherichia coli Cells by Flow Cytometry.

Accurate coordination of chromosome replication and cell division is essential for cellular processes, yet the regulatory mechanisms governing the bacterial cell cycle remain contentious. The lack of quantitative data connecting key cell cycle players at the single-cell level across large samples hinders consensus. Employing high-throughput flow cytometry, we quantitatively correlated the expression levels of key cell cycle proteins (FtsZ, MreB, and DnaA) with DNA content in individual bacteria. Our findings reveal distinct correlations depending on the chromosome number (CN), specifically whether CN ≤2 or ≥4, unveiling a mixed regulatory scenario in populations where CN of 2 or 4 coexist. We observed function-dependent regulations for these key proteins across nonoverlapping division cycles and various nutrient conditions. Notably, a logarithmic relationship between total protein content and replication origin number across nutrient conditions suggests a unified mechanism governing cell cycle progression, confirming the applicability of Schaechter's growth law to cells with CN ≥4. For the first time, we established a proportional relationship between the synthesis rates of key cell cycle proteins and chromosome dynamics in cells with CN ≥4. Drug experiments highlighted CN 2 and 4 as pivotal turning points influencing cellular resource allocation. This high-throughput, single-cell analysis provides interconnected quantitative insights into key molecular events, facilitating a predictive understanding of the relationship between cell growth and cell cycle.

Dependence of post-segregational killing mediated by Type II restriction-modification systems on the lifetime of restriction endonuclease effective activity.

Plasmid-borne Type II restriction-modification (RM) systems mediate post-segregational killing (PSK). PSK is thought to be caused by the dilution of restriction and modification enzymes during cell division, resulting in accumulation of unmethylated DNA recognition sites and their cleavage by restriction endonucleases. PSK is the likely reason for stabilization of plasmids carrying RM systems in the absence of selection for plasmid maintenance. In this study, we developed a CRISPR interference-based method to eliminate RM-carrying plasmids and study PSK-related phenomena with minimal perturbation to the Escherichia coli host. Plasmids carrying the EcoRV, Eco29kI, and EcoRI RM systems were highly stable, and their loss resulted in SOS response and PSK. In contrast, plasmids carrying the Esp1396I system were poorly stabilized; their loss led to a temporary cessation of growth, followed by full recovery. We demonstrate that this unusual behavior is due to a limited lifetime of the Esp1396I restriction endonuclease activity, which, upon Esp1396I plasmid loss, disappears approximately after two cycles of cell division, i.e., before unmethylated sites appear in significant numbers. Our results indicate that whenever PSK induced by a loss of RM systems, and, possibly, other toxin-antitoxin systems, is considered, the lifetimes of individual system components and the growth rate of host cells shall be taken in account. Mathematical modeling shows, that unlike the situation with classical toxin-antitoxin systems, RM system-mediated PSK is possible when the lifetimes of restriction endonuclease and methyltransferase activities are similar, as long as the toxic restriction endonuclease activity persists for more than two chromosome replication cycles.IMPORTANCEIt is widely accepted that many Type II restriction-modification (RM) systems mediate post-segregational killing (PSK) if plasmids that encode them are lost. In this study, we harnessed an inducible CRISPR-Cas system to remove RM plasmids from Escherichia coli cells to study PSK while minimally perturbing cell physiology. We demonstrate that PSK depends on restriction endonuclease activity lifetime and is not observed when it is less than two replication cycles. We present a mathematical model that explains experimental data and shows that unlike the case of toxin-antitoxin-mediated PSK, the loss of an RM system induced PSK even when the RM enzymes have identical lifetimes.

Influence and mechanism of typical transition metal ions on the denitrification performance of heterotrophic nitrification-aerobic denitrification bacteria

To investigate the inhibitory effects of various transition metal ions on nitrogen removal and their underlying mechanisms, the single and combined effects of Cu2+ Ni2+ and Zn2+ on Heterotrophic nitrification-aerobic denitrification (HN-AD) bacteria Acinetobacter sp. TAC-1 were studied in a batch experiment system. The results revealed that increasing concentrations of Cu2+ and Ni2+ had a detrimental effect on the removal of ammonium nitrogen (NH4+-N) and total nitrogen (TN). Specifically, Cu2+ concentration of 10 mg/L, the TN degradation rate was 55.09%, compared to 77.60% in the control group. Cu2+ exhibited a pronounced inhibitory effect. In contrast, Zn2+ showed no apparent inhibitory effect on NH4+-N removal and even enhanced TN removal at lower concentrations. However, when the mixed ion concentration of Zn2++Ni2+ exceeded 5 mg/L, the removal rates of NH4+-N and TN were significantly reduced. Moreover, transition metal ions did not significantly impact the removal rates of chemical oxygen demand (COD). The inhibition model fitting results indicated that the inhibition sequence was Cu2+ > Zn2+ > Ni2+. Transcriptome analysis demonstrated that metal ions influence TAC-1 activity by modulating the expression of pivotal genes, including zinc ABC transporter substrate binding protein (znuA), ribosomal protein (rpsM), and chromosome replication initiation protein (dnaA) and DNA replication of TAC-1 under metal ion stress, leading to disruptions in transcription, translation, and cell membrane structure. Finally, a conceptual model was proposed by us to summarize the inhibition mechanism and possible response strategies of TAC-1 bacteria under metal ion stress, and to address the lack of understanding regarding the influence mechanism of TAC-1 on nitrogen removal in wastewater co-polluted by metal and ammonia nitrogen. The results provided practical guidance for the management of transition metal and ammonia nitrogen co-polluted water bodies, as well as the removal of high nitrogen.

Stability and gene strand bias of lambda prophages and chromosome organization in Escherichia coli.

Temperate phage-mediated horizontal gene transfer is a potent driver of genetic diversity in the evolution of bacteria. Most lambdoid prophages in Escherichia coli are integrated into the chromosome with the same orientation with respect to the direction of chromosomal replication, and their location on the chromosome is far from homogeneous. To better understand these features, we studied the interplay between lysogenic and lytic states of phage lambda in both native and inverted integration orientations at the wild-type integration site as well as at other sites on the bacterial chromosome. Measurements of free phage released by spontaneous induction showed that the stability of lysogenic states is affected by location and orientation along the chromosome, with stronger effects near the origin of replication. Competition experiments and range expansions between lysogenic strains with opposite orientations and insertion loci indicated that there are no major differences in growth. Moreover, measurements of the level of transcriptional bursts of the cI gene coding for the lambda phage repressor using single-molecule fluorescence in situ hybridization resulted in similar levels of transcription for both orientations and prophage location. We postulate that the preference for a given orientation and location is a result of a balance between the maintenance of lysogeny and the ability to lyse.IMPORTANCEThe integration of genetic material of temperate bacterial viruses (phages) into the chromosomes of bacteria is a potent evolutionary force, allowing bacteria to acquire in one stroke new traits and restructure the information in their chromosomes. Puzzlingly, this genetic material is preferentially integrated in a particular orientation and at non-random sites on the bacterial chromosome. The work described here reveals that the interplay between the maintenance of the stability of the integrated phage, its ability to excise, and its localization along the chromosome plays a key role in setting chromosomal organization in Escherichia coli.

The budding yeast Fkh1 Forkhead associated (FHA) domain promoted a G1-chromatin state and the activity of chromosomal DNA replication origins.

In Saccharomyces cerevisiae, the forkhead (Fkh) transcription factor Fkh1 (forkhead homolog) enhances the activity of many DNA replication origins that act in early S-phase (early origins). Current models posit that Fkh1 acts directly to promote these origins' activity by binding to origin-adjacent Fkh1 binding sites (FKH sites). However, the post-DNA binding functions that Fkh1 uses to promote early origin activity are poorly understood. Fkh1 contains a conserved FHA (forkhead associated) domain, a protein-binding module with specificity for phosphothreonine (pT)-containing partner proteins. At a small subset of yeast origins, the Fkh1-FHA domain enhances the ORC (origin recognition complex)-origin binding step, the G1-phase event that initiates the origin cycle. However, the importance of the Fkh1-FHA domain to either chromosomal replication or ORC-origin interactions at genome scale is unclear. Here, S-phase SortSeq experiments were used to compare genome replication in proliferating FKH1 and fkh1-R80A mutant cells. The Fkh1-FHA domain promoted the activity of 100 origins that act in early to mid- S-phase, including the majority of centromere-associated origins, while simultaneously inhibiting 100 late origins. Thus, in the absence of a functional Fkh1-FHA domain, the temporal landscape of the yeast genome was flattened. Origins are associated with a positioned nucleosome array that frames a nucleosome depleted region (NDR) over the origin, and ORC-origin binding is necessary but not sufficient for this chromatin organization. To ask whether the Fkh1-FHA domain had an impact on this chromatin architecture at origins, ORC ChIPSeq data generated from proliferating cells and MNaseSeq data generated from G1-arrested and proliferating cell populations were assessed. Origin groups that were differentially regulated by the Fkh1-FHA domain were characterized by distinct effects of this domain on ORC-origin binding and G1-phase chromatin. Thus, the Fkh1-FHA domain controlled the distinct chromatin architecture at early origins in G1-phase and regulated origin activity in S-phase.

Genome-wide identification of bZIP transcription factors and their expression analysis in Platycodon grandiflorus under abiotic stress.

The Basic Leucine Zipper (bZIP) transcription factors (TFs) family is among of the largest and most diverse gene families found in plant species, and members of the bZIP TFs family perform important functions in plant developmental processes and stress response. To date, bZIP genes in Platycodon grandiflorus have not been characterized. In this work, a number of 47 PgbZIP genes were identified from the genome of P. grandiflorus, divided into 11 subfamilies. The distribution of these PgbZIP genes on the chromosome and gene replication events were analyzed. The motif, gene structure, cis-elements, and collinearity relationships of the PgbZIP genes were simultaneously analyzed. In addition, gene expression pattern analysis identified ten candidate genes involved in the developmental process of different tissue parts of P. grandiflorus. Among them, Four genes (PgbZIP5, PgbZIP21, PgbZIP25 and PgbZIP28) responded to drought and salt stress, which may have potential biological roles in P. grandiflorus development under salt and drought stress. Four hub genes (PgbZIP13, PgbZIP30, PgbZIP32 and PgbZIP45) mined in correlation network analysis, suggesting that these PgbZIP genes may form a regulatory network with other transcription factors to participate in regulating the growth and development of P. grandiflorus. This study provides new insights regarding the understanding of the comprehensive characterization of the PgbZIP TFs for further exploration of the functions of growth and developmental regulation in P. grandiflorus and the mechanisms for coping with abiotic stress response.

Establishing a biochemical understanding of the initiation of chromosome replication in bacteria

In the mid-1950s, Arthur Kornberg elucidated the enzymatic synthesis of DNA by DNA polymerase, for which he was recognized with the 1959 Nobel Prize in Physiology or Medicine. He then identified many of the proteins that cooperate with DNA polymerase to replicate duplex DNA of small bacteriophages. However, one major unanswered problem was understanding the mechanism and control of the initiation of chromosome replication in bacteria. In a seminal paper in 1981, Fuller, Kaguni, and Kornberg reported the development of a cell-free enzyme system that could replicate DNA that was dependent on the bacterial origin of DNA replication, oriC. This advance opened the door to a flurry of discoveries and important papers that elucidated the process and control of initiation of chromosome replication in bacteria.

RPI-GGCN: Prediction of RNA-Protein Interaction Based on Interpretability Gated Graph Convolution Neural Network and Co-Regularized Variational Autoencoders.

RNA-protein interactions (RPIs) play an important role in several fundamental cellular physiological processes, including cell motility, chromosome replication, transcription and translation, and signaling. Predicting RPI can guide the exploration of cellular biological functions, intervening in diseases, and designing drugs. Given this, this study proposes the RPI-gated graph convolutional network (RPI-GGCN) method for predicting RPI based on the gated graph convolutional neural network (GGCN) and co-regularized variational autoencoder (Co-VAE). First, different types of feature information were extracted from RNA and protein sequences by nine feature extraction methods. Second, Co-VAEs are used to eliminate the redundancy of fused features and generate optimal features. Finally, this study introduces gated cyclic units into graph convolutional networks (GCNs) to construct a model for RPI prediction, which efficiently extracts topological information and improves the model's interpretable feature learning and expression capabilities. In the fivefold cross-validation test, the RPI-GGCN method achieved prediction accuracies of 97.27%, 97.32%, 96.54%, 95.76%, and 94.98% on the RPI369, RPI488, RPI1446, RPI1807, and RPI2241 datasets. To test the generalization performance of the model, we used the model trained on RPI369 to predict the independent NPInter v3.0 dataset and achieved excellent performance in all six independent validation sets. By visualizing the RPI network graph based on the prediction results, we aim to provide a new perspective and reference for studying RPI mechanisms and exploring new RPIs. Extensive experimental results demonstrate that RPI-GGCN can provide an efficient, accurate, and stable RPI prediction method.