DNA Recombination Research Articles

Appelmans Protocol for in vitro Klebsiella pneumoniae phage host range expansion leads to induction of the novel temperate linear plasmid prophage vB_KpnS-KpLi5.

Adjuvant therapy with bacteriophage (phage) cocktails in combination with antibiotics is a therapeutic approach currently considered for treatment of infections with encapsulated, biofilm forming, and multidrug-resistant Klebsiella pneumoniae (Kp). Klebsiella phage are highly selective in targeting a bacterial capsule type. Considering the numerous Kp capsule types and other host restriction factors, phage treatment could be facilitated when generating phages with a broad host range. A modified 'Appelmans protocol' was used to create phages with an extended host range via invitro forced DNA recombination. Three T7-like Kp phages with highly colinear genomes were subjected to successive propagation on their susceptible host strains representing the capsule types K64, K27, and K23, and five Kp isolates of the same capsule types initially unsusceptible for phage lysis. After 30 propagation cycles, five phages were isolated via plaque assay. Four output phages represented the original input phages, while the fifth lysed a previously non-permissible Kp isolate, which was not lysed by any of the input phages. Surprisingly, sequence analysis revealed a novel N15/phiKO2-like phage genome (vB_KpnS_KpLi5) lacking substantial homologies to any of the used T7-like phages. This phage is not a chimeric recombinant of the applied T7-like phages, but represents a temperate phage that was induced from Kp due to the application of the input phages phages (cocktail), but not by any of them individually. Adapted phages with chimeric genomes and extended host range derived from input phages were not observed. Induction of temperate phages may be a stress response caused by using multiple phages simultaneously (i.e., by destabilization of the cell wall due to an unspecific binding of the phages). Successive use of different phages for therapeutic purposes may be preferable over simultaneous application in cocktail formulations to avoid undesired induction of temperate phages.

RecombiCraft library construction: A novel method for DNA library cloning and expansion using non-enzymatic single-step DNA recombination and liquid culture

In this study, we introduce RecombiCraft, an innovative, rapid, and cost-efficient method for constructing DNA libraries in E. coli. This method uses seamless ligation cloning extract (SLiCE) coupled with liquid culture amplification to effectively minimize sequence biases. The technique capitalizes on the natural homologous recombination capabilities of E. coli cell lysates, eliminating the need for multiple purified enzymes and reducing costs. We first synthesized the library backbone and inserts via PCR, employing high-fidelity polymerase to minimize sequence bias. The SLiCE technique was then used to assemble the DNA fragments introduced into E. coli through electroporation. To ensure the integrity of the library, we optimized culture times based on next-generation sequencing analysis which confirmed the minimal sequence bias. The RecombiCraft method demonstrates that this approach is economical and maintains the library’s uniformity. By using liquid culture, this method can complete DNA library generation in about 12 hours and final extraction is simple, making it a promising tool for genetic research and biotechnology applications.

Ultra-Efficient Integration of Gene Libraries onto Yeast Cytosolic Plasmids.

Efficient methods for diversifying genes of interest (GOIs) are essential in protein engineering. For example, OrthoRep, a yeast-based orthogonal DNA replication system that achieves the rapid in vivo diversification of GOIs encoded on a cytosolic plasmid (p1), has been successfully used to drive numerous protein engineering campaigns. However, OrthoRep-based GOI evolution has almost always started from single GOI sequences, limiting the number of locations on a fitness landscape from where evolutionary search begins. Here, we present a simple approach for the high-efficiency integration of GOI libraries onto OrthoRep. By leveraging integrases, we demonstrate recombination of donor DNA onto the cytosolic p1 plasmid at exceptionally high transformation efficiencies, even surpassing the transformation efficiency of standard circular plasmids into yeast. We demonstrate our method's utility through the straightforward construction of mock nanobody libraries encoded on OrthoRep, from which rare binders were reliably enriched. Overall, integrase-assisted manipulation of yeast cytosolic plasmids should enhance the versatility of OrthoRep in continuous evolution experiments and support the routine construction of large GOI libraries in yeast in general.

HIV Replication Under High-Level Cabotegravir Is Associated with the Appearance of 3′-PPT Mutations, Circular DNA Transcription and Recombination

The HIV integrase inhibitor, dolutegravir (DTG), in the absence of eliciting integrase (int) resistance, has been reported to select mutations in the virus 3′-polypurine tract (3′-PPT) adjacent to the 3′-LTR U3. An analog of DTG, cabotegravir (CAB), has a high genetic barrier to drug resistance and is used in formulations for treatment and long-acting pre-exposure prophylaxis. We examined whether mutations observed for DTG would emerge in vitro with CAB. HIV-1IIIB was cultured in paired experiments of continuous high (300 nM) CAB initiated 2 h or 24 h after infection. After eight months of CAB treatment, no int resistance was detected. Conversely, HIV RNA 3′-PPT mutants were detected within one month and were the majority virus by day 98. The appearance of 3′-PPT variants coincided with a rapid accumulation of HIV 1-LTR and 2-LTR circles. RNA amplification from the 3′-LTR TAR identified transcripts crossing 2-LTR circle junctions, which incorporated the adjacent U5 sequence identical to the 3′-PPT mutants. 3′-PPT variants were only identified in LTR circles and transcripts. Additionally, we found evidence of linear HIV and LTR circle recombination with human DNA at motifs homologous to 3′-PPT sequences. HIV persistence under CAB was associated with transcription and recombination of LTR circle sequences.

Meiotic DNA break resection and recombination rely on chromatin remodeler Fun30.

DNA double-strand breaks (DSBs) are nucleolytically processed to generate single-stranded DNA for homologous recombination. In Saccharomyces cerevisiae meiosis, this resection involves nicking by the Mre11-Rad50-Xrs2 complex (MRX), then exonucleolytic digestion by Exo1. Chromatin remodeling at meiotic DSBs is thought necessary for resection, but the remodeling enzyme was unknown. Here we show that the SWI/SNF-like ATPase Fun30 plays a major, nonredundant role in meiotic resection. A fun30 mutation shortened resection tracts almost as severely as an exo1-nd (nuclease-dead) mutation, and resection was further shortened in a fun30 exo1-nd double mutant. Fun30 associates with chromatin in response to DSBs, and the constitutive positioning of nucleosomes governs resection endpoint locations in the absence of Fun30. We infer that Fun30 promotes both the MRX- and Exo1-dependent steps in resection, possibly by removing nucleosomes from broken chromatids. Moreover, the extremely short resection in fun30 exo1-nd double mutants is accompanied by compromised interhomolog recombination bias, leading to defects in recombination and chromosome segregation. Thus, this study also provides insight about the minimal resection lengths needed for robust recombination.

POLD3 as Controller of Replicative DNA Repair.

Multiple modes of DNA repair need DNA synthesis by DNA polymerase enzymes. The eukaryotic B-family DNA polymerase complexes delta (Polδ) and zeta (Polζ) help to repair DNA strand breaks when primed by homologous recombination or single-strand DNA annealing. DNA synthesis by Polδ and Polζ is mutagenic, but is needed for the survival of cells in the presence of DNA strand breaks. The POLD3 subunit of Polδ and Polζ is at the heart of DNA repair by recombination, by modulating polymerase functions and interacting with other DNA repair proteins. We provide the background to POLD3 discovery, investigate its structure, as well as function in cells. We highlight unexplored structural aspects of POLD3 and new biochemical data that will help to understand the pivotal role of POLD3 in DNA repair and mutagenesis in eukaryotes, and its impact on human health.

DNA damage response inhibitors in cancer therapy: lessons from the past, current status and future implications.

The DNA damage response (DDR) is a network of proteins that coordinate DNA repair and cell-cycle checkpoints to prevent damage being transmitted to daughter cells. DDR defects lead to genomic instability, which enables tumour development, but they also create vulnerabilities that can be used for cancer therapy. Historically, this vulnerability has been taken advantage of using DNA-damaging cytotoxic drugs and radiotherapy, which are more toxic to tumour cells than to normal tissues. However, the discovery of the unique sensitivity of tumours defective in the homologous recombination DNA repair pathway to PARP inhibition led to the approval of six PARP inhibitors worldwide and to a focus on making use of DDR defects through the development of other DDR-targeting drugs. Here, we analyse the lessons learnt from PARP inhibitor development and how these may be applied to new targets to maximize success. We explore why, despite so much research, no other DDR inhibitor class has been approved, and only a handful have advanced to later-stage clinical trials. We discuss why more reliable predictive biomarkers are needed, explore study design from past and current trials, and suggest alternative models for monotherapy and combination studies. Targeting multiple DDR pathways simultaneously and potential combinations with anti-angiogenic agents or immune checkpoint inhibitors are also discussed.

Rare variants of DNA ligase 1 show distinct mechanisms of deficiency

Human DNA ligase 1 (LIG1) performs the final step in DNA repair and recombination pathways by sealing DNA breaks, and it functions as the main replicative ligase. Hypomorphic LIG1 variants R771W and R641L cause immune deficiencies in LIG1 Syndrome patients. In vitro these LIG1 variants have decreased catalytic efficiency and increased abortive ligation and it is not known if either biochemical defect is sufficient on its own to cause immune deficiency. We investigated the enzymatic activity of several new candidate LIG1 Syndrome variants chosen based on their structural proximity to known clinical variants, low minor allele frequency (MAF), high level of conservation, and concurrence in patients with similar symptoms as LIG1 Syndrome patients. The R305Q substitution is in the DNA binding domain, R768W is in the OB-fold domain, and R641S is in the nucleotidyltransferase domain. Biochemical characterization confirmed deficiencies in ligase activity for all three variants, but also revealed marked differences in comparison to the known LIG1 Syndrome variants. Both the R305Q and R768W substitutions increase the KM for DNA and decrease the catalytic efficiency, however, neither exhibit elevated levels of abortive ligation. In contrast, the R641S variant exhibits a greater impairment of activity as well as a more pronounced level of abortive ligation compared to the known LIG1 Syndrome variant, R641L. This work expands the number of LIG1 alleles that are likely candidates for LIG1 Syndrome, and it raises the question of whether distinct enzymatic deficiencies in LIG1 cause unique clinical impacts in patients harboring these alleles.

The Shu complex is an ATPase that regulates Rad51 filaments during homologous recombination in the DNA damage response.

Rad51 filaments are Rad51-coated single-stranded DNA and essential in homologous recombination (HR). The yeast Shu complex (Shu) is a conserved regulator of homologous recombination, working through its modulation on Rad51 filaments to direct HR-associated DNA damage response. However, the biochemical properties of Shu remain unclear, which hinders molecular insights into Shu’s role in HR and the DNA damage response. In this work, we biochemically characterized Shu and analyzed its molecular actions on single-stranded DNA and Rad51 filaments. First, we revealed that Shu preferentially binds fork-shaped DNA with 20nt ssDNA components. Then, we identified and validated, through site-specific mutagenesis, that Shu is an ATPase and hydrolyzes ATP in a DNA-dependent manner. Furthermore, we showed that Shu interacts with ssDNA and Rad51 filaments and alters the properties of ssDNA and the filaments with a 5’-3’ polarity. The alterations depend on the ATP hydrolysis of Shu, suggesting that the ATPase activity of Shu is important in regulating its functions. The preference of Shu for acting on the 5’ end of Rad51 filaments aligns with the observation that Shu promotes lesion bypass at the lagging strand of a replication fork. Our work on Shu, a prototype modulator of Rad51 filaments in eukaryotes, provides a general molecular mechanism for Rad51-mediated error-free DNA lesion bypass.

Pan-genomic characterization and structural variant analysis reveal insights into spore development and species diversity in Ganoderma.

Understanding the genomic diversity and functional implications of Ganoderma species is crucial for elucidating their evolutionary history and biotechnological potential. Here, we present the first pan-genomic analysis of Ganoderma spp., combining five newly sequenced genomes with ten publicly available genomes. Our comprehensive comparative study unveiled a rich genomic landscape, identifying core genes shared among all Ganoderma strains and species-specific gene sets. Additionally, we identified structural variants impacting the expression of key genes, including insights into the MSH4 gene involved in DNA repair and recombination processes, which exhibits a 440 bp insertion in the promoter region and a leucine-to-serine mutation in the gene body, potentially increasing spore production in the S3 strain. Overall, our study provides valuable insights into the genomic architecture and functional diversity of Ganoderma, paving the way for further research on its evolutionary dynamics, biotechnological applications and pharmaceutical potential.

Identification of novel DNA gyrase inhibitor by combined pharmacophore modeling, QSAR analysis, molecular docking, molecular dynamics, ADMET and DFT approaches

DNA gyrase, an ATP-dependent enzyme, plays a critical role in DNA replication, transcription, and recombination in Mycobacterium tuberculosis (MTB). While fluoroquinolones are effective antibacterial agents targeting DNA gyrase, their clinical use is often limited due to side effects and the emergence of bacterial resistance. In this study, we developed a quantitative structure-activity relationship (QSAR) model to predict the anti-tubercular activity of Quinoline-Aminopiperidine derivatives targeting the DNA gyrase enzyme, using a dataset of 48 compounds obtained from the literature. The QSAR model was validated using both internal and external validation metrics. Model 4, the best predictive model, demonstrated a strong fit with an R² of 0.8393, an adjusted R² (R²adj) of 0.8010, and a lack of fit (LOF) parameter of 0.0626. The QSAR results revealed that DNA gyrase inhibition is significantly influenced by factors such as partition coefficient, molecular flexibility, hydrogen bonding potential, and the presence of fluorine atoms. Twelve quinoline-aminopiperidine derivatives were designed, and their anti-tubercular activity was predicted using QSAR model-4. These compounds were further assessed for pharmacokinetic properties, toxicity, and binding affinity to DNA gyrase. Pharmacophore modeling was also performed and validated using MOE software. The final pharmacophore model includes the features of two aromatic hydrophobic features, one hydrogen bond acceptor, and one hydrogen bond donor. The results indicated that designed compounds QA-3 and dataset compounds C-34 exhibit favorable drug-likeness properties. Molecular dynamics simulations confirmed the stable binding of compounds QA-3 and C-34 to the DNA gyrase protein, highlighting their potential as promising anti-tubercular agents. The developed QSAR Model-4 will facilitate the prediction of anti-tubercular activity in Quinoline-Aminopiperidine derivatives.

Tumour-intrinsic PDL1 signals regulate the Chk2 DNA damage response in cancer cells and mediate resistance to Chk1 inhibitors

BackgroundAside from the canonical role of PDL1 as a tumour surface-expressed immune checkpoint molecule, tumour-intrinsic PDL1 signals regulate non-canonical immunopathological pathways mediating treatment resistance whose significance, mechanisms, and therapeutic targeting remain incompletely understood. Recent reports implicate tumour-intrinsic PDL1 signals in the DNA damage response (DDR), including promoting homologous recombination DNA damage repair and mRNA stability of DDR proteins, but many mechanistic details remain undefined.MethodsWe genetically depleted PDL1 from transplantable mouse and human cancer cell lines to understand consequences of tumour-intrinsic PDL1 signals in the DNA damage response. We complemented this work with studies of primary human tumours and inducible mouse tumours. We developed novel approaches to show tumour-intrinsic PDL1 signals in specific subcellular locations. We pharmacologically depleted tumour PDL1 in vivo in mouse models with repurposed FDA-approved drugs for proof-of-concept clinical translation studies.ResultsWe show that tumour-intrinsic PDL1 promotes the checkpoint kinase-2 (Chk2)-mediated DNA damage response. Intracellular but not surface-expressed PDL1 controlled Chk2 protein content post-translationally and independently of PD1 by antagonising PIRH2 E3 ligase-mediated Chk2 polyubiquitination and protein degradation. Genetic tumour PDL1 depletion specifically reduced tumour Chk2 content but not ATM, ATR, or Chk1 DDR proteins, enhanced Chk1 inhibitor (Chk1i) synthetic lethality in vitro in diverse human and murine tumour models, and improved Chk1i efficacy in vivo. Pharmacologic tumour PDL1 depletion with cefepime or ceftazidime replicated genetic tumour PDL1 depletion by reducing tumour Chk2, inducing Chk1i synthetic lethality in a tumour PDL1-dependent manner, and reducing in vivo tumour growth when combined with Chk1i.ConclusionsOur data challenge the prevailing surface PDL1 paradigm, elucidate important and previously unappreciated roles for tumour-intrinsic PDL1 in regulating the ATM/Chk2 DNA damage response axis and E3 ligase-mediated protein degradation, suggest tumour PDL1 as a biomarker for Chk1i efficacy, and support the rapid clinical potential of pharmacologic tumour PDL1 depletion to treat selected cancers.

A history-dependent integrase recorder of plant gene expression with single-cell resolution

During development, most cells experience a progressive restriction of fate that ultimately results in a fully differentiated mature state. Understanding more about the gene expression patterns that underlie developmental programs can inform engineering efforts for new or optimized forms. Here, we present a four-state integrase-based recorder of gene expression history and demonstrate its use in tracking gene expression events in Arabidopsis thaliana in two developmental contexts: lateral root initiation and stomatal differentiation. The recorder uses two serine integrases to mediate sequential DNA recombination events, resulting in step-wise, history-dependent switching between expression of fluorescent reporters. By using promoters that express at different times along each of the two differentiation pathways to drive integrase expression, we tie fluorescent status to an ordered progression of gene expression along the developmental trajectory. In one snapshot of a mature tissue, our recorder is able to reveal past gene expression with single cell resolution. In this way, we are able to capture heterogeneity in stomatal development, confirming the existence of two alternate paths of differentiation.

Hop2-Mnd1 functions as a DNA sequence fidelity switch in Dmc1-mediated DNA recombination

Homologous recombination during meiosis is critical for chromosome segregation and also gives rise to genetic diversity. Genetic exchange between homologous chromosomes during meiosis is mediated by the recombinase Dmc1, which is capable of recombining DNA sequences with mismatches. The Hop2-Mnd1 complex mediates Dmc1 activity. Here, we reveal a regulatory role for Hop2-Mnd1 in restricting substrate selection. Specifically, Hop2-Mnd1 upregulates Dmc1 activity with DNA substrates that are either fully homologous or contain DNA mismatches, and it also acts against DNA strand exchange between substrates solely harboring microhomology. By isolating and examining salient Hop2-Mnd1 separation-of-function variants, we show that suppressing illegitimate DNA recombination requires the Dmc1 filament interaction attributable to Hop2-Mnd1 but not its DNA binding activity. Our study provides mechanistic insights into how Hop2-Mnd1 helps maintain meiotic recombination fidelity.

Cytonuclear evolution in fully heterotrophic plants: lifestyles and gene function determine scenarios

BackgroundEvidence shows that full mycoheterotrophs and holoparasites often have reduced plastid genomes with rampant gene loss, elevated substitution rates, and deeply altered to conventional evolution in mitochondrial genomes, but mechanisms of cytonuclear evolution is unknown. Endoparasitic Sapria himalayana and mycoheterotrophic Gastrodia and Platanthera guangdongensis represent different heterotrophic types, providing a basis to illustrate cytonuclear evolution. Here, we focused on nuclear-encoded plastid / mitochondrial (N-pt / mt) -targeting protein complexes, including caseinolytic protease (ClpP), ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo), oxidative phosphorylation system (OXPHOS), DNA recombination, replication, and repair (DNA-RRR) system, and pentatricopeptide repeat (PPR) proteins, to identify evolutionary drivers for cytonuclear interaction.ResultsThe severity of gene loss of N-pt PPR and pt-RRR genes was positively associated with increased degree of heterotrophy in full mycoheterotrophs and S. himalayana, while N-mt PPR and mt-RRR genes were retained. Substitution rates of organellar and nuclear genes encoding N-pt/mt subunits in protein complexes were evaluated, cytonuclear coevolution was identified in S. himalayana, whereas disproportionate rates of evolution were observed in the OXPHOS complex in full mycoheterotrophs, only slight accelerations in substitution rates were identified in N-mt genes of full mycoheterotrophs.ConclusionsNuclear compensatory evolution was identified in protein complexes encoded by plastid and N-pt genes. Selection shaping codon preferences, functional constraint, mt-RRR gene regulation, and post-transcriptional regulation of PPR genes all facilitate mito-nuclear evolution. Our study enriches our understanding of genomic coevolution scenarios in fully heterotrophic plants.

Isolation and purification of DNA double-strand break repair intermediates for understanding complex molecular mechanisms.

Branched DNA molecules are key intermediates in the molecular pathways of DNA replication, repair and recombination. Understanding their structural details, therefore, helps to envisage the mechanisms underlying these processes. While the configurations of DNA molecules can be effectively analysed in bulk using gel electrophoresis techniques, direct visualization provides a complementary single-molecule approach to investigating branched DNA structures. However, for microscopic examination, the sample needs to be free from impurities that could obscure the molecules of interest, and free from the bulk of unwanted non-specific DNA molecules that would otherwise dominate the field of view. Additionally, in the case of recombination intermediates, the length of the DNA molecules becomes an important factor to consider since the structures can be spread over a large distance on the chromosome in vivo. As a result, apart from sample purity, efficient isolation of large-sized DNA fragments without damaging their branched structures is crucial for further analysis. These factors are illustrated by the example of DNA double-strand break repair in the bacterium E. coli. In E. coli recombination intermediates may be spread over a distance of 40 kb which constitutes less than 1% of the 4.6 Mb genome. This study reveals ways to overcome some of the technical challenges that are associated with the isolation and purification of large and complex branched DNA structures using E. coli DNA double-strand break repair intermediates. High-molecular weight and branched DNA molecules do not run into agarose gels subjected to electrophoresis. However, they can be extracted from the wells of the gels if they are agarose embedded, by using β-agarase digestion, filtration, and concentration. Furthermore, a second round of gel electrophoresis followed by purification is recommended to enhance the purity of the specific DNA samples. These preliminary findings may prove to be pioneering for various single-molecule analyses of large and complex DNA molecules of DNA replication, repair and recombination.

Delivering large genes using adeno-associated virus and the CRE-lox DNA recombination system.

Adeno-associated virus (AAV) is a safe and efficient gene delivery vehicle for gene therapies. However, its relatively small packaging capacity limits its use as a gene transfer vector. Here, we describe a strategy to deliver large genes that exceed the AAV's packaging capacity using up to four AAV vectors and the CRE-lox DNA recombination system. We devised novel lox sites by combining non-compatible and reaction equilibrium-modifying lox site variants. These lox sites facilitate sequence-specific and near-unidirectional recombination of AAV vector genomes, enabling efficient reconstitution of up to 16kb of therapeutic genes in a pre-determined configuration. Using this strategy, we have developed AAV gene therapy vectors to deliver IFT140, PCDH15, CEP290, and CDH23 and demonstrate efficient production of full-length proteins in cultured mammalian cells and mouse retinas. Notably, AAV-IFT140 gene therapy vectors ameliorated retinal degeneration and preserved visual functions in an IFT140-associated retinitis pigmentosa mouse model. The CRE-lox approach described here provides a simple, flexible, and effective platform for generating AAV gene therapy vectors beyond AAV's packaging capacity.

Comprehensive in silico characterization of Arabidopsis thaliana RecQl helicases through structure prediction and molecular dynamics simulations

Helicases are ubiquitous enzymes with specific functions that contribute to almost all nucleic acid metabolic processes. The RecQ helicase family is essential for integrity in all organisms through DNA replication, repair, and recombination. This study investigated five RecQ-like helicases in Arabidopsis thaliana (AtRecQl) that exhibit diverse structural and physiochemical attributes and functions. Cis-regulatory element analysis identified stress, hormone, cell cycle, and development-responsive modules involved in various events in plant growth and development. Gene ontology analysis revealed that the five AtRecQl were associated with various cellular components, molecular functions, and biological processes. Protein-protein interaction analysis also implicated some in various abiotic stress processes. Structural analysis and molecular dynamics (MD) simulations were performed to examine conformational stability through root means square deviation and radius of gyration, showing stable AtRecQl protein structures. Free energy landscape analysis validated thermodynamically stable structures throughout the MD simulation. Principle component analysis and probability density functions from MD simulations provided satisfactory structural variational data for the complexes and limited coordinate movements. These insights might greatly benefit future studies.

Using AI-predicted protein structures as a reference to predict loss-of-function activity in tumor suppressor breast cancer genes

BackgroundThe loss-of-function (LOF) classification of most missense variants in tumor suppressor breast cancer genes BRCA1, BRCA2, PALB2, and RAD51C remains unclassified and confounds clinical actionability. Classifying these variants is challenging due to their rarity, leading clinicians to rely on in silico predictive methods. Protein stability changes are associated with function, making stability predictors valuable. Stability predictions upon missense variant perturbations require high-resolution protein structures. However, the availability of these high-resolution structures is lacking. This study explores using generative AI to predict high-resolution protein structures, which can then be analyzed with in silico protein stability prediction methods to assess LOF activity in ordered regions of the protein. This study also determines the appropriate in silico protein stability and dedicated in silico missense prediction methods in dbNSFP v4.7 database to predict LOF activity in ordered regions of these four genes. Functional classifications from homology recombination DNA repair (HDR) assays and variant classifications from the ClinVar database provide a reliable dataset for evaluating the performance of these in silico prediction methods. ResultsComplex AlphaFold2 structures of the BRCA1-C terminal (BRCT) domain and the DNA-binding (DB) domain of BRCA2, analyzed using protein stability tool FoldX predicts LOF activity from missense variants significantly better than experimentally-derived structures in ordered regions. The BRCT domain achieved an Area Under the Curve (AUC)= 0.861 (95 % CI:0.858–0.863) and AUC= 0.842 (95 % CI:0.840–0.845), while the DB domain achieved an AUC= 0.836 (95 % CI:0.8322–0.841), compared to AUC= 0.847 (95 % CI:0.844–0.850) and AUC= 0.835 (95 % CI:0.832–0.837) from the BRCT domain, and AUC= 0.830 (95 % CI:0.821–0.8320) from the DB domain from experimentally-derived structures. Protein stability does not predict LOF activity from missense variants better than dedicated in silico missense predictors. Overall, we find that AlphaMissense ranks highly, with an average AUC= 0.890 (95 % CI 0.886–0.895) from ordered regions across these four cancer genes, compared to all other in silico missense predictors present in the dbNSFP database. ConclusionsThe study reveals that generative AI protein predicted structures can outperform experimentally-derived structures in evaluating LOF activity from predicted protein stability in ordered regions of genes BRCA1, BRCA2, PALB2 and RAD51C. The study also highlights the predictive performance of AlphaMissense as the premier in silico missense prediction method to predict LOF activity from missense variants in these four tumor suppressor breast cancer genes. The code for this study can be downloaded for free on GitHub (https://github.com/rohandavidg/CarePred)

Context effects on repair of 5'-overhang DNA double-strand breaks induced by Cas12a in Arabidopsis.

Sequence-specific endonucleases have been key to the study of the mechanisms and control of DNA double-strand break (DSB) repair and recombination, and the availability of CRISPR-Cas nucleases over the last decade has driven rapid progress in the understanding and application of targeted recombination in many organisms, including plants. We present here an analysis of recombination at targeted chromosomal 5' overhang DSB generated by the FnCas12a endonuclease in the plant, Arabidopsis thaliana. The much-studied Cas9 nuclease cleaves DNA to generate blunt-ended DSBs, but relatively less is known about the repair of other types of breaks, such as those with 5'-overhanging ends. Sequencing the repaired breaks clearly shows that the majority of repaired DSB carry small deletions and are thus repaired locally by end-joining recombination, confirmed by Nanopore sequencing of larger amplicons. Paired DSBs generate deletions at one or both cut-sites, as well as deletions and reinsertions of the deleted segment between the two cuts, visible as inversions. While differences are seen in the details, overall the deletion patterns are similar between repair at single-cut and double-cut events, notwithstanding the fact that only the former involve cohesive DNA overhangs. A strikingly different repair pattern is however observed at breaks flanked by direct repeats. This change in sequence context results in the presence of a major alternative class of repair events, corresponding to highly efficient repair by single-strand annealing recombination.