Structure-function Analysis Research Articles

Tomosyns attenuate SNARE assembly and synaptic depression by binding to VAMP2-containing template complexes

Tomosyns are widely thought to attenuate membrane fusion by competing with synaptobrevin-2/VAMP2 for SNARE-complex assembly. Here, we present evidence against this scenario. In a novel mouse model, tomosyn-1/2 deficiency lowered the fusion barrier and enhanced the probability that synaptic vesicles fuse, resulting in stronger synapses with faster depression and slower recovery. While wild-type tomosyn-1m rescued these phenotypes, substitution of its SNARE motif with that of synaptobrevin-2/VAMP2 did not. Single-molecule force measurements indeed revealed that tomosyn’s SNARE motif cannot substitute synaptobrevin-2/VAMP2 to form template complexes with Munc18-1 and syntaxin-1, an essential intermediate for SNARE assembly. Instead, tomosyns extensively bind synaptobrevin-2/VAMP2-containing template complexes and prevent SNAP-25 association. Structure-function analyses indicate that the C-terminal polybasic region contributes to tomosyn’s inhibitory function. These results reveal that tomosyns regulate synaptic transmission by cooperating with synaptobrevin-2/VAMP2 to prevent SNAP-25 binding during SNARE assembly, thereby limiting initial synaptic strength and equalizing it during repetitive stimulation.

Multi-Gene Recombinant Baculovirus Expression Systems: From Inception to Contemporary Applications.

Many protein expression systems are primarily utilised to produce a single, specific recombinant protein. In contrast, most biological processes such as virus assembly rely upon a complex of several interacting proteins rather than the activity of a sole protein. The high complexity of the baculovirus genome, coupled with a multiphase replication cycle incorporating distinct transcriptional steps, made it the ideal system to manipulate for high-level expression of a single, or co-expression of multiple, foreign proteins within a single cell. We have developed and utilised a series of recombinant baculovirus systems to unravel the sequential assembly process of a complex non-enveloped model virus, bluetongue virus (BTV). The high protein yields expressed by the baculovirus system not only facilitated structure-function analysis of each viral protein but were also advantageous to crystallography studies and supported the first atomic-level resolution of a recombinant viral protein, the major BTV capsid protein. Further, the formation of recombinant double-shelled virus-like particles (VLPs) provided insights into the structure-function relationships among the four major structural proteins of the BTV whilst also representing a potential candidate for a viral vaccine. The baculovirus multi-gene expression system facilitated the study of structurally complex viruses (both non-enveloped and enveloped viruses) and heralded a new generation of viral vaccines.

Abstract 1718: The lysine demethylase KDM2A regulates tumor cell clustering to potentiate metastasis

Abstract Approximately 70% of patients diagnosed with non-small cell lung cancer (NSCLC) are not eligible for curative resection, often due to disseminated disease at diagnosis. Thus, the advent of novel therapeutics, especially those efficacious for metastatic disease, is of primary importance. Metastasis correlates with epigenetic alterations that drive tumor plasticity and adaptive heterogeneity in response to the multiple microenvironments metastatic cells see, yet the molecular mechanisms underlying these epigenetic alterations remain poorly understood. We performed a functional genomic screen and identified the lysine 36 on histone H3 (H3K36) demethylase KDM2A as a mediator of in vivo and in vitro metastatic competence in multiple models of NSCLC. KDM2A depletion in highly metastatic NSCLC cell lines reduces their capacity to form multi-tumor cell clusters. In vivo, KDM2A depletion does not impact lung orthotopic tumor cell growth or the total number of circulating tumor cells (CTCs) but does reduce their ability to colonize distant metastatic sites. This is consistent with the fact that human CTCs found in clusters are more metastatic than CTCs which disseminate as single-cells. Because the molecular and biological consequences of H3K36 regulation are context dependent, we performed for the first time an integrated transcriptomic-epigenomic analysis of target gene regulation, KDM2A genomic binding, H3K36 methylation, and H3K27 acetylation, in metastatic cell clusters. KDM2A enhances the expression of cell adhesion and anti-oxidant genes, while suppressing inflammatory gene responses, and this regulation preferentially occurs in tumor cell clusters. Transcriptomic analysis of CTCs from patients confirms that KDM2A regulated transcriptional responses correlate with metastatic relapse across multiple cancer types. Moreover, KDM2A may function both as a transcriptional activator and repressor in metastatic clusters. As KDM2A encodes for several domains which can independently regulate E3-ubiqutin ligation, H3K36me2 demethylation, and DNA methylation, we have also performed a structure function analysis of K2MDA in the context of metastatic phenotypes. Thus, the distinct biochemical functions of KDM2A in relation to its regulation of CTC gene transcription and histone modification will also be presented. In conclusion, we have identified KDM2A as a novel regulator of NSCLC metastatic spread. The ability to inhibit KDM2A with small molecule inhibitors may be translated into effective therapeutic interventions for NSCLC patients with metastasis disease. Citation Format: Carolyn Kravitz, Kiran D. Patel, Wesley L. Cai, Anna Arnal-Estapé, Dejian Zhao, Emily Wingrove, Laura Stevens, Thomas F. Westbrook, Qin Yan, Don X. Nguyen. The lysine demethylase KDM2A regulates tumor cell clustering to potentiate metastasis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1718.

Structure–function analysis of PorXFj, the PorX homolog from Flavobacterium johnsioniae, suggests a role of the CheY-like domain in type IX secretion motor activity

The type IX secretion system (T9SS) is a large multi-protein transenvelope complex distributed into the Bacteroidetes phylum and responsible for the secretion of proteins involved in pathogenesis, carbohydrate utilization or gliding motility. In Porphyromonas gingivalis, the two-component system PorY sensor and response regulator PorX participate to T9SS gene regulation. Here, we present the crystal structure of PorXFj, the Flavobacterium johnsoniae PorX homolog. As for PorX, the PorXFj structure is comprised of a CheY-like N-terminal domain and an alkaline phosphatase-like C-terminal domain separated by a three-helix bundle central domain. While not activated and monomeric in solution, PorXFj crystallized as a dimer identical to active PorX. The CheY-like domain of PorXFj is in an active-like conformation, and PorXFj possesses phosphodiesterase activity, in agreement with the observation that the active site of its phosphatase-like domain is highly conserved with PorX.

Tetramerization is essential for the enzymatic function of the Pseudomonas aeruginosa virulence factor UDP-glucose pyrophosphorylase.

Multidrug-resistant bacteria such as the opportunistic pathogen Pseudomonas aeruginosa, which causes life-threatening infections especially in immunocompromised individuals and cystic fibrosis patients, pose an increasing threat to public health. In the search for new treatment options, P. aeruginosa uridine diphosphate-glucose pyrophosphorylase (PaUGP) has been proposed as a novel drug target because it is required for the biosynthesis of important virulence factors and linked to pathogenicity in animal models. Here, we show that UGP-deficient P. aeruginosa exhibits severely reduced virulence against human lung tissue and cells, emphasizing the enzyme's suitability as a drug target. To establish a basis for the development of selective PaUGP inhibitors, we solved the product-bound crystal structure of tetrameric PaUGP and conducted a comprehensive structure-function analysis, identifying key residues at two different molecular interfaces that are essential for tetramer integrity and catalytic activity and demonstrating that tetramerization is pivotal for PaUGP function. Importantly, we show that part of the PaUGP oligomerization interface is uniquely conserved across bacterial UGPs but does not exist in the human enzyme, therefore representing an allosteric site that may be targeted to selectively inhibit bacterial UGPs.IMPORTANCEInfections with the opportunistic bacterial pathogen Pseudomonas aeruginosa are becoming increasingly difficult to treat due to multidrug resistance. Here, we show that the enzyme uridine diphosphate-glucose pyrophosphorylase (UGP) is involved in P. aeruginosa virulence toward human lung tissue and cells, making it a potential target for the development of new antibacterial drugs. Our exploration of P. aeruginosa (Pa)UGP structure-function relationships reveals that the activity of PaUGP depends on the formation of a tetrameric enzyme complex. We found that a molecular interface involved in tetramer formation is conserved in all bacterial UGPs but not in the human enzyme, and therefore hypothesize that it provides an ideal point of attack to selectively inhibit bacterial UGPs and exploit them as drug targets.

Structure-Function Analysis of the A1 Subunit of Shiga Toxin 2 with Peptides That Target the P-Stalk Binding Site and Inhibit Activity.

Shiga toxin 2a (Stx2a) is the virulence factor of Escherichia coli (STEC), which is associated with hemolytic uremic syndrome, the leading cause of pediatric kidney failure. The A1 subunit of Stx2a (Stx2A1) binds to the conserved C-terminal domain (CTD) of the ribosomal P-stalk proteins to remove an adenine from the sarcin-ricin loop (SRL) in the 28S rRNA, inhibiting protein synthesis. There are no antidotes against Stx2a or any other ribosome-inactivating protein (RIP). The structural and functional details of the binding of Stx2A1 to the P-stalk CTD are not known. Here, we carry out a deletion analysis of the conserved P-stalk CTD and show that the last eight amino acids (P8) of the P-stalk proteins are the minimal sequence required for optimal affinity and maximal inhibitory activity against Stx2A1. We determined the first X-ray crystal structure of Stx2A1 alone and in complex with P8 and identified the exact binding site. The C-terminal aspartic acid of the P-stalk CTD serves as an anchor, forming key contacts with the conserved arginine residues at the P-stalk binding pocket of Stx2A1. Although the ricin A subunit (RTA) binds to the P-stalk CTD, the last aspartic acid is more critical for the interaction with Stx2A1, indicating that RIPs differ in their requirements for the P-stalk. These results demonstrate that the catalytic activity of Stx2A1 is inhibited by blocking its interactions with the P-stalk, providing evidence that P-stalk binding is an essential first step in the recruitment of Stx2A1 to the SRL for depurination.

Increased LL37 in psoriasis and other inflammatory disorders promotes LDL uptake and atherosclerosis.

Patients with chronic inflammatory disorders such as psoriasis have an increased risk of cardiovascular disease and elevated levels of LL37, a cathelicidin host defense peptide that has both antimicrobial and proinflammatory properties. To explore whether LL37 could contribute to the risk of heart disease, we examined its effects on lipoprotein metabolism and show that LL37 enhanced LDL uptake in macrophages through the LDL receptor (LDLR), scavenger receptor class B member 1 (SR-B1), and CD36. This interaction led to increased cytosolic cholesterol in macrophages and changes in expression of lipid metabolism genes consistent with increased cholesterol uptake. Structure-function analysis and synchrotron small-angle x-ray scattering showed structural determinants of the LL37-LDL complex that underlie its ability to bind its receptors and promote uptake. This function of LDL uptake is unique to cathelicidins from humans and some primates and was not observed with cathelicidins from mice or rabbits. Notably, Apoe-/- mice expressing LL37 developed larger atheroma plaques than did control mice, and a positive correlation between plasma LL37 and oxidized phospholipid on apolipoprotein B (OxPL-apoB) levels was observed in individuals with cardiovascular disease. These findings provide evidence that LDL uptake can be increased via interaction with LL37 and may explain the increased risk of cardiovascular disease associated with chronic inflammatory disorders.

Structure-function analysis of the cyclic β-1,2-glucan synthase from Agrobacterium tumefaciens

The synthesis of complex sugars is a key aspect of microbial biology. Cyclic β-1,2-glucan (CβG) is a circular polysaccharide critical for host interactions of many bacteria, including major pathogens of humans (Brucella) and plants (Agrobacterium). CβG is produced by the cyclic glucan synthase (Cgs), a multi-domain membrane protein. So far, its structure as well as the mechanism underlining the synthesis have not been clarified. Here we use cryo-electron microscopy (cryo-EM) and functional approaches to study Cgs from A. tumefaciens. We determine the structure of this complex protein machinery and clarify key aspects of CβG synthesis, revealing a distinct mechanism that uses a tyrosine-linked oligosaccharide intermediate in cycles of polymerization and processing of the glucan chain. Our research opens possibilities for combating pathogens that rely on polysaccharide virulence factors and may lead to synthetic biology approaches for producing complex cyclic sugars.

Structure-function analysis of the ATPase domain of African swine fever virus topoisomerase.

The ATPase domain of type II topoisomerase provides energy by hydrolyzing ATP and coordinates with the DNA-binding/cleavage domain to drive and control DNA transport. The precise molecular mechanisms of how these domains respond to DNA binding and ATP hydrolysis signals and communicate with each other remain elusive. We determine the first high-resolution crystal structure of the ATPase domain of African swine fever virus (ASFV) topo II in complex with AMPPNP and biochemically investigate its function in ATPase and DNA relaxation activities. Importantly, we find that mutations at three characteristic regions of the ASFV ATPase domain produce parallel effects on the basal/DNA-stimulated ATPase and relaxation activities, implying the tight coupling of the ATP hydrolysis and strand passage process. Therefore, our data provide important implications for understanding the strand passage mechanism of the type II topoisomerase and the structural basis for developing ATPase domain-targeting antivirals against ASFV.

Deciphering the mechanism of glutaredoxin-catalyzed roGFP2 redox sensing reveals a ternary complex with glutathione for protein disulfide reduction

Glutaredoxins catalyze the reduction of disulfides and are key players in redox metabolism and regulation. While important insights were gained regarding the reduction of glutathione disulfide substrates, the mechanism of non-glutathione disulfide reduction remains highly debated. Here we determined the rate constants for the individual redox reactions between PfGrx, a model glutaredoxin from Plasmodium falciparum, and redox-sensitive green fluorescent protein 2 (roGFP2), a model substrate and versatile tool for intracellular redox measurements. We show that the PfGrx-catalyzed oxidation of roGFP2 occurs via a monothiol mechanism and is up to three orders of magnitude faster when roGFP2 and PfGrx are fused. The oxidation kinetics of roGFP2-PfGrx fusion constructs reflect at physiological GSSG concentrations the glutathionylation kinetics of the glutaredoxin moiety, thus allowing intracellular structure-function analysis. Reduction of the roGFP2 disulfide occurs via a monothiol mechanism and involves a ternary complex with GSH and PfGrx. Our study provides the mechanistic basis for understanding roGFP2 redox sensing and challenges previous mechanisms for protein disulfide reduction.

Pathogenic CANVAS (AAGGG)n repeats stall DNA replication due to the formation of alternative DNA structures.

CANVAS is a recently characterized repeat expansion disease, most commonly caused by homozygous expansions of an intronic (A2G3)n repeat in the RFC1 gene. There are a multitude of repeat motifs found in the human population at this locus, some of which are pathogenic and others benign. In this study, we conducted structure-functional analyses of the pathogenic (A2G3)n and nonpathogenic (A4G)n repeats. We found that the pathogenic, but not the nonpathogenic, repeat presents a potent, orientation-dependent impediment to DNA polymerization in vitro. The pattern of the polymerization blockage is consistent with triplex or quadruplex formation in the presence of magnesium or potassium ions, respectively. Chemical probing of both repeats in vitro reveals triplex H-DNA formation by only the pathogenic repeat. Consistently, bioinformatic analysis of S1-END-seq data from human cell lines shows preferential H-DNA formation genome-wide by (A2G3)n motifs over (A4G)n motifs. Finally, the pathogenic, but not the nonpathogenic, repeat stalls replication fork progression in yeast and human cells. We hypothesize that the CANVAS-causing (A2G3)n repeat represents a challenge to genome stability by folding into alternative DNA structures that stall DNA replication.

Pion distribution functions from low-order Mellin moments

Exploiting an evolution scheme for parton distribution functions (DFs) that is all-orders exact, contemporary lattice-QCD (lQCD) results for low-order Mellin moments of the pion valence quark DF are shown to be mutually consistent. The analysis introduces a means by which key odd moments can be obtained from the even moments in circumstances where only the latter are available. Combining these elements, one arrives at parameter-free lQCD-based predictions for the pointwise behaviour of pion valence, glue, and sea DFs, with sound uncertainty estimates. The behaviour of the pion DFs at large light-front momentum fraction, x≳0.85, is found to be consistent with QCD expectations and continuum analyses of pion structure functions, i.e., damping like (1−x)βparton, with βvalence≈2.4, βglue≈3.6, βsea≈4.6. It may be possible to test these predictions using data from forthcoming experiments.

Structure functional analysis of Ltc1 and Ltc2a: Antimicrobial peptides from spider venom

Structure functional analysis of Ltc1 and Ltc2a: Antimicrobial peptides from spider venom

Fluorescent human RPA to track assembly dynamics on DNA

DNA metabolic processes including replication, repair, recombination, and telomere maintenance occur on single-stranded DNA (ssDNA). In each of these complex processes, dozens of proteins function together on the ssDNA template. However, when double-stranded DNA is unwound, the transiently open ssDNA is protected and coated by the high affinity heterotrimeric ssDNA binding Replication Protein A (RPA). Almost all downstream DNA processes must first remodel/remove RPA or function alongside to access the ssDNA occluded under RPA. Formation of RPA-ssDNA complexes trigger the DNA damage checkpoint response and is a key step in activating most DNA repair and recombination pathways. Thus, in addition to protecting the exposed ssDNA, RPA functions as a gatekeeper to define functional specificity in DNA maintenance and genomic integrity. RPA achieves functional dexterity through a multi-domain architecture utilizing several DNA binding and protein-interaction domains connected by flexible linkers. This flexible and modular architecture enables RPA to adopt a myriad of configurations tailored for specific DNA metabolic roles. To experimentally capture the dynamics of the domains of RPA upon binding to ssDNA and interacting proteins we here describe the generation of active site-specific fluorescent versions of human RPA (RPA) using 4-azido-L-phenylalanine (4AZP) incorporation and click chemistry. This approach can also be applied to site-specific modifications of other multi-domain proteins. Fluorescence-enhancement through non-canonical amino acids (FEncAA) and Förster Resonance Energy Transfer (FRET) assays for measuring dynamics of RPA on DNA are also described. The fluorescent human RPA described here will enable high-resolution structure-function analysis of RPA-ssDNA interactions.

Identification and function characterization of BnaBOR4 genes reveal their potential for Brassica napus cultivation under high boron stress

Boron (B) is essential for plant growth, but toxic in excess. In several countries, soil toxic B levels are always a severe agricultural problem in arid and semi-arid regions. Phytoremediation of excess B containing soil is still in its infancy, while high B tolerant plants with elevated protein abundance of B efflux transporter were successfully established or explored. Brassica napus (B. napus) is one of the most important oil crops. However, B efflux transporters underlying excess B tolerance in B. napus remain unknown. Here, we reported that in Brassicaceae species, B. napus had four homologous genes of Arabidopsis AtBOR4 , which were renamed BnaBOR4.1, BnaBOR4.2, BnaBOR4.3 and BnaBOR4.4. BnaBOR4.1, BnaBOR4.2 and BnaBOR4.3 showed constitutive expression and BnaBOR4.4 appeared to be a pseudogene. BnaBOR4.2 and BnaBOR4.3 were expressed in inner cell layers and BnaBOR4.1 in the outer cell layer in root tip, and all were expressed in vascular tissue in the mature zone. B efflux activity assays in yeast demonstrated that BnaBOR4.1, BnaBOR4.2 and AtBOR4 but not BnaBOR4.3 had comparable levels of B transport activity. Structure-functional analysis between BnaBOR4.3 and BnaBOR4.2 demonstrated that amino acid residue substitution at position 297 (Ala vs Pro) and 427 (Met vs Leu) is critical for the B transport activity. Mutant BnaBOR4.3M427L partially restored the B efflux activity, and both mutants BnaBOR4.3A297P and BnaBOR4.3A297P&M427L fully restored B efflux activity, indicating that the Pro297 residue is critical for their function. Further validation of BnaBOR4 was accomplished by growing transgenic Arabidopsis plants under high B conditions. Taken together, our study identified two functional B efflux genes BnaBOR4.1 and BnaBOR4.2 in B. napus, and a key amino acid residue proline 297 associated with B efflux activity. This study highlights the potential of BanBOR4 genes for B. napus cultivation under high B stress.

Exploring biodiversity patterns at the landscape scale by linking landscape energy and land use/land cover heterogeneity

Agricultural Biodiversity dynamics has been evaluated by social metabolism or by landscape structure-function analysis. In this study, by using ELIA modeling, we used both methods in combination to understand how the interplay between social metabolism and landscape structure-function can affect biodiversity pattern distribution. We used energy reinvestment (E) as an indicator of social metabolism and landscape heterogeneity (Le) as an indicator of landscape structure-function. We propose a research hypothesis to analyze biodiversity patterns considering four different clusters identified based on high or low E or Le. As cluster 1, we defined E as high and Le as low and associated natural ecosystems to it. These ecosystems are expected to contain high species abundance but low richness. As cluster 2, both E and Le were defined as high and semi-natural ecosystems were associated to it, where nature friendly farm system developed. In these ecosystems, high species abundance and richness are expected. Cluster 3 with low E and Le was associated intensive farmland, which is due to the simplification of the landscape. Here, low energy reinvestment and landscape heterogeneity confirm that ecosystem services related to biodiversity have been drastically reduced. Lastly, cluster 4 with low E but high Le refers to intensive mosaics of farmland and pasture. In this cluster, the biodiversity richness index is high due to spatial landscape diversity, but the biodiversity abundance index is low due to the lack of energy reinvestment. We evaluate the proposed hypothesis for biodiversity analysis in the Qazvin province, emphasizing the interplay between energy availability and landscape heterogeneity in shaping ecological communities. This study highlights the importance of understanding biodiversity patterns at spatial scale and emphasizes the need for interdisciplinary research to address conservation and sustainability challenges. Our approach would be very useful where there is lack of biodiversity data.

A Versatile Aldehyde: Ferredoxin Oxidoreductase from the Organic Acid Reducing Thermoanaerobacter sp. Strain X514.

Aldehyde:ferredoxin oxidoreductases (AORs) have been isolated and biochemically-characterized from a handful of anaerobic or facultative aerobic archaea and bacteria. They catalyze the ferredoxin (Fd)-dependent oxidation of aldehydes to acids. Recently, the involvement of AOR in the reduction of organic acids to alcohols with electrons derived from sugar or synthesis gas was demonstrated, with alcohol dehydrogenases (ADHs) carrying out the reduction of the aldehyde to the alcohol (AOR-ADH pathway). Here, we describe the biochemical characterization of an AOR of the thermophilic fermentative bacterium Thermoanaerobacter sp. strain X514 (AORX514). The putative aor gene (Teth514_1380) including a 6x-His-tag was introduced into the genome of the genetically-accessible, related species Thermoanaerobacter kivui. The protein was purified to apparent homogeneity, and indeed revealed AOR activity, as measured by acetaldehyde-dependent ferredoxin reduction. AORX514 was active over a wide temperature (10 to 95 °C) and pH (5.5 to 11.5) range, utilized a wide variety of aldehydes (short and branched-chained, aliphatic, aromatic) and resembles archaeal sensu stricto AORs, as the protein is active in a homodimeric form. The successful, recombinant production of AORX514 in a related, well-characterized and likewise strict anaerobe paves the road towards structure-function analyses of this enzyme and possibly similar oxygen-sensitive or W/Mo-dependent proteins in the future.

Binary Mixtures of Choline Acetate and Tetrabutylammonium Acetate with Natural Organic Acids by Vibrational Spectroscopy and Molecular Dynamics Simulations.

We present a study of several mixtures obtained by the mixing of two organic acetate-based salts (choline acetate, ChAc, or tetrabutylammonium acetate, TBAAc) with three different natural organic acids (ascorbic acid, AA, citric acid, CA, and maleic acid, MA). The structures of the starting materials and of the mixtures were characterized by infrared spectroscopy (FT-IR) and classic molecular dynamics simulations (MD). The thermal behavior was characterized by differential scanning calorimetry (DSC) and thermogravimetry analysis (TGA). The obtained mixtures, especially the ChAc-based ones, strongly tend to vitrify at low temperatures and are stable up to 100-150 °C. The FTIR measurements suggest the formation of a strong H-bond network: the coordination between acids and ChAc or TBAAc takes place by the donation of the H-bond by the acids to the oxygen of the acetate anion, which acts as an acceptor (HBA). The comparison with MD analysis indicates that acids predominantly exploit their more acidic hydrogens. In particular, we observe the progressive shift of νC═O and νOH when the ratios of acids increase. The structural differences between the two studied cations influence the spatial distribution of the components in the mixture bulk phases. In particular, the analysis of the theoretical structure function I(q) of the TBAAc-based systems shows the presence of important prepeaks at low q, a sign of the formation of apolar domain, due to the nanosegregation of the alkyl chains.

Comparison of the structure-function of five newly members of the calcin family

Calcins are a group of scorpion toxin peptides specifically binding to ryanodine receptors (RyRs) with high affinity, and have the ability to activate and stabilize RyR in a long-lasting subconductance state. Five newly calcins synthesized compounds exhibit typical structural characteristics of a specific family through chemical synthesis and virtual analysis. As the calcins from the same species, Petersiicalcin1 and Petersiicalcin2, Jendekicalcin2 and Jendekicalcin3, have only one residue difference. Both Petersiicalcin1 and Petersiicalcin2 exhibited different affinities in stimulating [3H]ryanodine binding, but the residue mutation resulted in a 2.7 folds difference. Other calcins also exhibited a stimulatory effect on [3H]ryanodine binding to RyR1, however, their affinities were significantly lower than that of Petersiiicalcin1 and Petersiiicalcin2. The channel domain of RyR1 was found to be capable of binding with the basic residues of these calcins, which also exhibited interactions with the S6 helices on RyR1. Dynamic simulations were conducted for Petersiicalcin1 and Petersiicalcin2, which demonstrated their ability to form a highly stable conformation and resulting in an asymmetric tetramer structure of RyR1. The discovery of five newly calcins further enriches the diversity of the natural calcin family, which provides more native peptides for the structure-function analysis between calcin and RyRs.

Taxonomic Distribution, Phylogenetic Relationship, and Domain Conservation of CRISPR-Associated Cas Proteins

CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) is a naturally occurring genetic defense system in bacteria and archaea. It is comprised of a series of DNA sequence repeats with spacers derived from previous exposures to plasmid or phage. Further understanding and applications of CRISPR system have revolutionized our capacity for gene or genome editing of prokaryotes and eukaryotes. The CRISPR systems are classified into 3 distinct types: type I, type II, and type III, each of which possesses an associated signature protein, Cas3, Cas9, and Cas10, respectively. As the CRISPR loci originated from earlier independent exposures of foreign genetic elements, it is likely that their associated signature proteins may have evolved rapidly. Also, their functional domain structures might have experienced different selective pressures, and therefore, they have differentially diverged in their amino acid sequences. We employed genomic, phylogenetic, and structure-function constraint analyses to reveal the evolutionary distribution, phylogenetic relationship, and structure-function constraints of Cas3, Cas9, and Cas10 proteins. Results reveal that all 3 Cas-associated proteins are highly represented in the phyla Bacteroidetes, Firmicutes, and Proteobacteria, including both pathogenic and non-pathogenic species. Genomic analysis of homologous proteins demonstrates that the proteins share 30% to 50% amino acid identity; therefore, they are low to moderately conserved and evolved rapidly. Phylogenetic analysis shows that 3 proteins originated monophyletically; however, the evolution rates were different among different branches of the clades. Furthermore, structure-function constraint analysis reveals that both Cas3 and Cas9 proteins experiences low to moderate levels of negative selection, and several protein domains of Cas3 and Cas9 proteins are highly conserved. To the contrary, most protein domains of Cas10 proteins experience neutral or positive selection, which supports rapid genetic divergence and less structure-function constraints.