Protein Structure Analysis Research Articles

The Pfam protein families database: embracing AI/ML.

The Pfam protein families database is a comprehensive collection of protein domains and families used for genome annotation and protein structure and function analysis (https://www.ebi.ac.uk/interpro/). This update describes major developments in Pfam since 2020, including decommissioning the Pfam website and integration with InterPro, harmonization with the ECOD structural classification, and expanded curation of metagenomic, microprotein and repeat-containing families. We highlight how AlphaFold structure predictions are being leveraged to refine domain boundaries and identify new domains. New families discovered through large-scale sequence similarity analysis of AlphaFold models are described. We also detail the development of Pfam-N, which uses deep learning to expand family coverage, achieving an 8.8% increase in UniProtKB coverage compared to standard Pfam. We discuss plans for more frequent Pfam releases integrated with InterPro and the potential for artificial intelligence to further assist curation. Despite recent advances, many protein families remain to be classified, and Pfam continues working toward comprehensive coverage of the protein universe.

Hierarchical Analysis of Protein Structures: From Secondary Structures to Protein Units and Domains.

The three-dimensional structure of proteins is traditionally organized into hierarchical levels, specifically secondary structures and domains. However, different studies suggest the existence of intermediate levels, such as Protein Units (PUs), which provide a refined understanding of protein architecture. PUs, characterized by their compactness and independence, serve as an intermediate organizational level, bridging the gap between secondary structures and domains. This new view not only enhances our comprehension of protein structure, folding, and evolutionary mechanisms but also provides a robust methodology for identifying and categorizing protein domains. Based on the concept of PUs, alternative structural partitioning solutions can be proposed that address the structural ambiguity of proteins, leading to more meaningful domain identification.

Functional characterization of optic photoreception in Lymnaea stagnalis.

Optic photoreception is a critical function for animal survival. Across the evolutionary spectrum, diverse animal models have been used to investigate visual system function and potential mechanisms under physiological or pathophysiological states. However less is known on photoreceptive behaviors and retinal processing in invertebrates, especially molluscs. This study focuses on the freshwater pond snail, Lymnaea stagnalis (L. stagnalis), to explore its visual function and underlying mechanisms. Using anatomical and histological approaches we characterized the L. stagnalis eye structure and demonstrated structural connections and retinal rhodopsin-positive sensory cells potentially critical for phototransduction. To assess the snail phototactic responses, we developed a new neurobehavioral protocol and employed DeepLabCut to track and quantify animal locomotion. We demonstrated that L. stagnalis exhibits a positive locomotory response to intense focal light and has diverse photo-locomotory responses. Further, we conducted phylogenetic and protein structure analyses and demonstrated that L. stagnalis has a unique repertoire of both vertebrate and invertebrate phototransduction genes. Further characterization of a rhodopsin-like gene identified unique characteristics compared to other mollusks and vertebrates, suggesting different mechanisms of phototransduction. Taken together, our work establishes L. stagnalis as a model organism for studying optic photoreception, offering new insights into the evolution and diversity of visual function across animal species.

Characterization of mitochondrial DNA mutations in colorectal cancer progression by in silico approach and use as potential biomarkers for diagnosis and prognosis

BackgroundMitochondrial DNA variants are significant contributors to cancer progression, as evidenced by numerous findings. This study focuses on characterizing mitochondrial DNA mutations in colorectal cancer progression and their potential as biomarkers.MethodologyNext generation sequencing technology was employed to analyze mitochondrial DNA variants in tumor and adjacent normal tissues from 25 patients with colon/rectal cancer. In silico prediction tools (SIFT, Polyphen2, Mutation Assessor, and SNP&GO) were utilized to assess the pathogenicity of these variants. Additionally, homology modeling of mutated protein structures was conducted, and molecular dynamic simulations were performed to assess the impact of mutation on protein function.ResultsEighteen variants were identified across most tumor tissue samples, located in genes from Complex I, IV, and V. Among the identified variants, the V302M and S461 mutations in the MT-ND5 gene and L137F and L220P mutations in the ATP6 gene were predicted to be deleterious, potentially affecting protein function. 3D structural analysis of both wild-type and mutant proteins of MT-ND5 revealed changes in flexibility for the V302M and S461G mutations. The MT-ATP6 mutations L135F and L220P disrupt the interactions with surrounding residues and affect the overall function of protein. Further changes in protein dynamics of the mutated proteins by molecular dynamic simulations also indicate the effects; the mutations have on protein function.ConclusionMT-ND5 and MT-ATP6 variants could serve as potential biomarkers and drug targets in colorectal cancer. This study underscores the significance of mitochondrial DNA variants in cancer progression.

Investigating the role of keratin proteins and microbial associations in hereditary and pathogenic alopecia.

The purpose of this research was to identify the role of keratin proteins in causing inherited as well as pathogenic alopecia, pinpoint deleterious SNPs, and predict structural changes affecting protein-protein interactions in hair disorders. To elucidate the role of keratin proteins and genetic mutations in alopecia by analyzing protein structures through bioinformatics and identifying a mutation in the LPAR6 gene. It sought to identify the microorganisms linked to alopecia and conducted a comprehensive bioinformatics analysis of proteins with unknown experimental structures and molecular simulation analysis. The study identified a genetic mutation (c.188A > T, p.Asp63Val) in the LPAR6 gene associated with hereditary hair loss. Pathogenic alopecia was identified to be associated with S. aureus and two ic keratinophilic fungi namely M. canis, and T. violaceum. Additionally, among 14 proteins lacking prior structural information, four proteins namely Keratin, type II cuticular Hb3 (KR1), Keratin, type II cuticular Hb6 (KR2), Keratin, type II cytoskeletal 74 (KR3) and Keratin, type II cuticular Hb1 (KR4) exhibited common 'K-head' and 'F' domains. Docking analysis revealed five distinct binding sites (C1-C5) for each protein. The 'K-head' displayed the highest predicted binding affinities with Vina scores of -5.6 for KR2 and - 4.7 for KR4 whereas the 'F' domain showed Vina scores of -6.0 for KR3 and - 5.7 for KR2. This research underscores the crucial role of keratin proteins in both hereditary and pathogenic alopecia, emphasizing their significance for future investigations.

Optimal strategy for stabilizing protein folding intermediates.

To manipulate the protein concentration at a certain functional state through chemical stabilizers is crucial for protein-related studies. It not only plays a key role in protein structure analysis and protein folding kinetics, but also affects protein functionality to a large extent and thus has wide applications in medicine, food industry, etc. However, due to concerns about side effects or financial costs of stabilizers, identifying optimal strategies for enhancing protein stability with a minimal amount of stabilizers is of great importance. Here, we prove that either for the fixed terminal time (including both finite and infinite cases) or for the free one, the optimal control strategy for stabilizing the folding intermediates with a linear strategy for stabilizer addition belongs to the class of bang-bang controls. The corresponding optimal switching time is derived analytically, whose phase diagram with respect to several key parameters is explored in detail. The bang-bang control will be broken when nonlinear strategies for stabilizer addition are adopted. Moreover, the above theory is applied to the stabilization of erythropoietin by ten different kinds of chemicals, providing theoretical guidance for the selection and rational usage of stabilizers. Our current study on optimal strategies for protein stabilizers not only offers deep insights into the general picture of protein folding kinetics but also provides valuable theoretical guidance on treatments for protein-related diseases in medicine.

Advancing Native Mass Spectrometry Toward Cellular Biology.

Protein structure, including various post-translational modifications and higher-order structures, regulates diverse biological functions. Native mass spectrometry (native MS) is a powerful analytical technique used to determine the masses of biomolecules, such as proteins and their complexes, while preserving their native folding in solution. This method provides structural information on the composition of monomers or complexes and the stoichiometry of subunits within each complex, significantly contributing to protein structural analysis. Native MS has evolved to incorporate top-down approaches, enabling the characterization of proteoforms and non-covalent interactions between metabolites or proteins and specific targets. This perspective highlights the advancements in native MS for intracellular proteins and protein complexes, and discusses future research directions toward cellular biology.

Clinical report and genetic analysis of a Chinese family with retinitis pigmentosa 79 caused by a novel loss-of-function HK1 variant.

Retinitis pigmentosa (RP) is a genetically heterogeneous disease. RP 79 has been associated with heterozygous variants of hexokinase 1 (HK1). Only two missense HK1 variants have been reported in 11 families. To discover the molecular pathogenic mechanism of RP and validate the biological harm of HK1 through in vitro experiments. We conducted a genetic analysis of a 3-year-old female patient with RP and her family. We also evaluated the ocular phenotypes caused by HK1 (the identified variant). Peripheral blood samples were collected from the patient, her parents, and her brother, and trio whole-exome sequencing was performed. A protein structure analysis was performed to assess the functional impact of the variant, and a mutant plasmid was constructed for the quantitative polymerase chain reaction (qPCR) and western blot (WB) analysis of the effects of the variant on transcription and protein translation. The patient harbored the NM_000188.3: c.613del (p.Ala205Leufs*3) variant, which is a heterozygous variant of HK1. Sanger sequencing confirmed the presence of this variant in the patient; however, the patient's parents and brother had the wild-type variant. The protein structure analysis indicated that the variant resulted in a truncated protein caused by premature termination of amino acid coding. The qPCR results indicated that the variant may not have affected the transcription process. However, the WB analysis demonstrated that the mutant HK-1 protein was not expressed and that the wild-type group exhibited normal expression. Our patient had a loss-of-function (LoF) variant of HK1, which may be the genetic cause of typical features of RP that are observed at an early age. These findings expand the spectrum of HK1 variants and phenotypes and suggest that LoF variants of HK1 may represent a specific pathogenic mechanism of RP.

M, NM-polynomials Based on Reverse, Reduced Reverse Degree and Neighborhood Degree Based Topological Indices with Applications to Bond Energy of Y-Junction Nanotubes.

In graph theory, M polynomials like the matching polynomial are very crucial in examining the matching structures within graphs, while NM polynomials extends this to analyze non-matching edges. These polynomials are important in many fields, including chemistry and network architecture. They support the derivation of topological indices for protein structure analysis, network communication optimization, and drug design in QSAR/QSPR investigations. The aim of this paper is to define novel M and NM polynomials for different topological indices and to derive their closed-form expressions, specifically for Y-junction nanotubes. These new polynomials and indices are employed to create a robust QSPR model to predict bond energy in Y-junction nanotubes, that provide high accuracy and reliability in the model's statistical performance. This paper introduces new forms of M and NM polynomials tailored to specific topological indices related to reverse and neighborhood reverse properties. We derive closed-form expressions for these indices in Y-junction nanotubes. Furthermore, we develop a QSPR model to predict bond energy in Y-junction nanotubes using the newly defined indices. We define novel M and NM polynomials for various topological indices and derive precise expressions for Y-junction nanotubes. Utilizing these indices, we construct a highly accurate QSPR model (R² = 0.999) for predicting bond energy in Y-junction nanotubes, confirming the validity of our polynomial definitions and indices. We have presented new M and NM polynomials for different topological indices and derive their expressions specifically for Y-junction nanotubes. With these newly defined indices, we have developed a highly precise QSPR model to predict bond energy, achieving an R² value of 0.999. This work underscores the effectiveness of our polynomial definitions and indices in predicting material properties.

Determining key residues of engineered scFv antibody variants with improved MMP-9 binding using deep sequencing and machine learning

Given the crucial role of specific matrix metalloproteinases (MMPs) in the extracellular matrix, an imbalance in the regulation of activation of matrix metalloproteinase-9 (MMP-9) zymogen and inhibition of the enzyme can result in various diseases, such as cancer, neurodegenerative, and gynecological diseases. Thus, developing novel therapeutics that target MMP-9 with single-chain antibody fragments (scFvs) is a promising approach. We used fluorescent-activated cell sorting (FACS) to screen a synthetic scFv antibody library displayed on yeast for enhanced binding to MMP-9. The screened scFv mutants demonstrated improved binding to MMP-9 compared to the natural inhibitor of MMPs, tissue inhibitor of metalloproteinases (TIMPs). To identify the molecular determinants of these engineered scFv variants that affect binding to MMP-9, we used next-generation DNA sequencing and computational protein structure analysis. Additionally, a deep-learning language model was trained on the screened scFv library of variants to predict the binding affinities of scFv variants based on their CDR-H3 sequences.

Identification and characterization of ADAR1 mutations and changes in gene expression in human cancers

ADAR1 (Adenosine deaminase action on RNA1) is involved in post-transcriptional RNA editing. ADAR1 mutations have been identified in many cancers but its role tumor formation is still not well understood. Here we used available cancer genomes deposited on CSOMIC and cBioPortal to identify and characterize mutations and changes in ADAR1 expression in cancer cells. We identify several high frequency substitutions including one at R767 which is located in one of the dsRNA interacting domains. In silico protein structure analysis suggest the R767 mutations affect the protein stability and are likely to destabilize interaction with dsRNA. Gene expression analysis shows that in samples with under-expressed ADAR1, there is a statistically significant increase in expression of BLCAP (Bladder Cancer Associated Protein). Although BLCAP was initially identified in bladder cancers, more recent evidence shows that it is a tumor suppressor and BLCAP mutations have been identified in many cancer cells. Epistatic analysis using the cBioPortal mutual exclusivity calculator for the TCGA pan-cancer data shows that co-mutations between ADAR1 and other genes regulated by it are likely in cancer cells except for PTEN, AKT1 and BLCAP. This suggests that when ADAR1 function is impaired, PTEN, AKT1 and BLCAP become essential for survival of cancer cells. We also identified several samples with high mutation burden between ADAR1 and other genes regulated primarily in endometrial cancers. Finally, we show that the deaminase domain is highly conserved in metazoans and mutations within conserved residues do occur in human cancers suggesting that destabilization of the enzyme function is contributing to cancer development.

Structural Prediction and Antigenic Analysis of ROP18, MIC4, and SAG1 Proteins to Improve Vaccine Design against Toxoplasma gondii: An In silico Approach.

Toxoplasmosis is a cosmopolitan infectious disease in warm-blooded mammals that poses a serious worldwide threat due to the lack of effective medications and vaccines. The purpose of this study was to design a multi-epitope vaccine using several bioinfor-matics approaches against the antigens of Toxoplasma gondii (T. gondii). Three proteins of T. gondii, including ROP18, MIC4, and SAG1, were analyzed to predict the most dominant B- and T-cell epitopes. Finally, we designed a chimeric immunogen RMS (ROP18, MIC4, and SAG1) using some domains of ROP18 (N377-E546), MIC4 (D302-G471), and SAG1 (T130-L299) linked by rigid linker A (EAAAK) A. Physicochemical prop-erties, secondary and tertiary structures, antigenicity, and allergenicity of RMS were predicted utilizing immunoinformatic tools and servers. RMS protein had 545 amino acids with a molecular weight (MW) of 58,833.46 Da and a theoretical isoelectric point (IP) of 6.47. The secondary structure of RMS protein con-tained 21.28% alpha-helix, 24.59% extended strand, and 54.13% random coil. In addition, eval-uation of antigenicity and allergenicity showed the protein to be an immunogen and non-aller-gen. The results of the Ramachandran plot indicated that 76.4%, 12.9%, and 10.7% of amino acid residues were incorporated in the favored, allowed, and outlier regions, respectively. ΔG of the best-predicted mRNA secondary structure was -593.80 kcal/mol, which indicated that a stable loop was not formed at the 5' end. Finally, the accuracy and precision of the in silico analysis must be confirmed by successful heterologous expression and experimental studies.

A newly identified Y chromosome gene obp-Y is required for sperm storage in female Zeugodacus tau (Diptera: Tephritidae).

In the organisms with XX/XY sex chromosomes, Y chromosome is unique to males and plays an important role in male reproductive development. The study of Y chromosome genes will contribute to the development of pest genetic prevention and control technology. In this study, we identified 9 Y chromosome genes in Zeugodacus tau (Diptera: Tephritidae), including gene 16222. Protein structure analysis showed that 16222 was highly similar to odorant binding protein, and thus gene 16222 was named obp-Y. Obp-Y knockout (KO) significantly reduced hatching rate of offspring. Sperm detection results showed that obp-Y KO did not affect sperm number in the testes or sperm transfer during mating. We further examined the storage of sperms in females, and found that sperms in females mating with wild-type males began to transfer from spermathecal ducts to the spermathecae at hour 0 after the end of mating (AEM), and at 0-24h AEM, the sperm count in the spermathecae gradually increased. However, no sperms were observed in spermathecae of females mating with mutant males at hours 0, 4, 8, 24 and 48 AEM. In summary, this study revealed that Y chromosome gene obp-Y was necessary for the storage of sperms in females. Our findings not only provide theoretical basis for elucidating the function of the Y chromosome, but also offer a molecular target for the genetic control over Z. tau.

A Novel SPOTTED LEAF1-1 (SPL11-1) Gene Confers Resistance to Rice Blast and Bacterial Leaf Blight Diseases in Rice (Oryza sativa L.)

Programmed cell death (PCD) plays critical roles in plant immunity but must be regulated to prevent excessive damage. In this study, a novel spotted leaf (spl11-1) mutant was identified from an ethyl methane sulfonate (EMS) population. The SPL11-1 gene was genetically mapped to chromosome 12 between the Indel12-37 and Indel12-39 molecular markers, which harbor a genomic region of 27 kb. Annotation of the SPL11-1 genomic region revealed the presence of two candidate genes. Through gene prediction and cDNA sequencing, it was confirmed that the target gene in the spl11-1 mutant is allelic to the rice SPOTTED LEAF (SPL11), hereafter referred to as spl11-1. Sequence analysis of SPL11 revealed a single bp deletion (T) between the spl11-1 mutant and the ‘Shuangkang77009’ wild type. Moreover, protein structure analysis showed that the structural differences between the SPL11-1 and SPL11 proteins might lead to a change in the function of the SPL11 protein. Compared to the ‘Shuangkang77009’ wild type, the spl11-1 mutant showed more disease resistance. The agronomical evaluation showed that the spl11-1 mutant showed more adverse traits. Through further mutagenesis treatment, we obtained the spl11-2 mutant allelic to spl11-1, which has excellent agronomic traits and more improvement and may have certain breeding prospects in future breeding for disease resistance.

Unveiling the Virome of Wild Birds: Exploring CRESS-DNA Viral Dark Matter.

Amid global health concerns and the constant threat of zoonotic diseases, this study delves into the diversity of circular replicase-encoding single-stranded DNA (CRESS-DNA) viruses within Chinese wild bird populations. Employing viral metagenomics to tackle the challenge of "viral dark matter," the research collected and analyzed 3,404 cloacal swab specimens across 26 bird families. Metagenomic analysis uncovered a rich viral landscape, with 67.48% of reads classified as viral dark matter, spanning multiple taxonomic levels. Notably, certain viral families exhibited host-specific abundance patterns, with Galliformes displaying the highest diversity. Diversity analysis categorized samples into distinct groups, revealing significant differences in viral community structure, particularly noting higher diversity in terrestrial birds compared to songbirds and unique diversity in migratory birds versus perching birds. The identification of ten novel Circoviridae viruses, seven Smacoviridae viruses, and 167 Genomoviridae viruses, along with 100 unclassified CRESS-DNA viruses, underscores the expansion of knowledge on avian-associated circular DNA viruses. Phylogenetic and structural analyses of Rep proteins offered insights into evolutionary relationships and potential functional variations among CRESS-DNA viruses. In conclusion, this study significantly enhances our understanding of the avian virome, shedding light on the intricate relationships between viral communities and host characteristics in Chinese wild bird populations. The diverse array of CRESS-DNA viruses discovered opens avenues for future research into viral evolution, spread factors, and potential ecosystem impacts.

Comparative analysis of allosteric rearrangements in FtsZ protein structure induced by benzamide and 4-hydroxycoumarine compounds

Aim. To reveal allosteric rearrangements of FtsZ molecules arising under the influence of benzamide compounds and 4-hydroxycoumarin derivatives. To discover the key molecular mechanisms predetermining the effect of the specified compounds on the cell division in bacteria. Methods. Comparative analysis of FtsZ protein structures and their complexes with ligands. Application of structural bioinformatics software for molecular visualization, measurement of interatomic distances and approximation of intramolecular shifts based on RMSD indicators. Results. Revealed conformational changes in FtsZ protein molecules, induced by allosteric effectors: 4-hydroxycoumarin – 4HC and benzamide – 9PC (PC-190723). Allosteric deformations and their consequences for intact FtsZ protein molecules, there GTPase domains, H7 helixes and C-terminal domains were studied. Conclusions. It was clarified that the binding of benzamides causes more significant shifts in the structure of the FtsZ protein monomer, its C-terminal domain, and H7 helix. At the same time, 4-hydroxycoumarins deform the structure of the GTPase domain almost twofold effectively. Both classes of compounds prosses allosteric action through unique mechanisms that are largely realized through deformations and displacements of the H7 helix. Despite the fact that these compounds demonstrate different allosteric mechanisms of action, their final effect can be summarized to destructions in GTP pocket, protofilament interfaces and the general geometry of molecule.

Quantifying uncertainty of molecular mismatch introduced by mislabeled ancestry using haplotype-based HLA genotype imputation.

Modern histocompatibility algorithms depend on the comparison and analysis of high-resolution HLA protein sequences and structures, especially when considering epitope-based algorithms, which aim to model the interactions involved in antibody or T cell binding. HLA genotype imputation can be performed in the cases where only low/intermediate-resolution HLA genotype is available or if specific loci are missing, and by providing an individuals' race/ethnicity/ancestry information, imputation results can be more accurate. This study assesses the effect of imputing high-resolution genotypes on molecular mismatch scores under a variety of ancestry assumptions. We compared molecular matching scores from "ground-truth" high-resolution genotypes against scores from genotypes which are imputed from low-resolution genotypes. Analysis was focused on a simulated patient-donor dataset and confirmed using two real-world datasets, and deviations were aggregated based on various ancestry assumptions. We observed that using multiple imputation generally results in lower error in molecular matching scores compared to single imputation, and that using the correct ancestry assumptions can reduce error introduced during imputation. We conclude that for epitope analysis, imputation is a valuable and low-risk strategy, as long as care is taken regarding epitope analysis context, ancestry assumptions, and (multiple) imputation strategy.

Selection on the vascular-remodeling BMPER gene is associated with altitudinal adaptation in an insular lizard

Abstract High altitude imposes several extreme constraints on life, such as low oxygen pressure and high levels of ultraviolet radiation, which require specialized adaptations. Many studies have focused on how endothermic vertebrates respond to these challenging environments, but there is still uncertainty on how ectotherms adapt to these conditions. Here, we used whole-genome sequencing of low-altitude (100–600 m) and high-altitude (3,550 m) populations of the wide-ranging Tenerife lizard Gallotia galloti to uncover signatures of selection for altitudinal adaptation. The studied populations show reduced differentiation, sharing similar patterns of genetic variation. Selective sweep mapping suggests that signatures of adaptation to high altitude are not widespread across the genome, clustering in a relatively small number of genomic regions. One of these regions contains BMPER, a gene involved with vascular remodeling, and that has been associated with hypoxia-induced angiogenic response. By genotyping samples across 2 altitudinal transects, we show that allele frequency changes at this locus are not gradual, but rather show a well-defined shift above ca. 1,900 m. Transcript and protein structure analyses on this gene suggest that putative selection likely acts on noncoding variation. These results underline how low oxygen pressure generates the most consistent selective constraint in high-altitude environments, to which vertebrates with vastly contrasting physiological profiles need to adapt in the context of ongoing climate change.

Structural and Evolutionary Analysis of Proteins Endowed with a Nucleotidyltransferase, or Non-canonical Palm, Catalytic Domain

AbstractMany polymerases and other proteins are endowed with a catalytic domain belonging to the nucleotidyltransferase fold, which has also been deemed the non-canonical palm domain, in which three conserved acidic residues coordinate two divalent metal ions. Tertiary structure-based evolutionary analyses provide valuable information when the phylogenetic signal contained in the primary structure is blurry or has been lost, as is the case with these proteins. Pairwise structural comparisons of proteins with a nucleotidyltransferase fold were performed in the PDBefold web server: the RMSD, the number of superimposed residues, and the Qscore were obtained. The structural alignment score (RMSD × 100/number of superimposed residues) and the 1-Qscore were calculated, and distance matrices were constructed, from which a dendogram and a phylogenetic network were drawn for each score. The dendograms and the phylogenetic networks display well-defined clades, reflecting high levels of structural conservation within each clade, not mirrored by primary sequence. The conserved structural core between all these proteins consists of the catalytic nucleotidyltransferase fold, which is surrounded by different functional domains. Hence, many of the clades include proteins that bind different substrates or partake in non-related functions. Enzymes endowed with a nucleotidyltransferase fold are present in all domains of life, and participate in essential cellular and viral functions, which suggests that this domain is very ancient. Despite the loss of evolutionary traces in their primary structure, tertiary structure-based analyses allow us to delve into the evolution and functional diversification of the NT fold.

Enhancing solution structural analysis of large molecular proteins through optimal stereo array isotope labeling of aromatic amino acids

The observation of side-chain peaks of aromatic amino acids is the prerequisite for a high-resolution three-dimensional structure determination of proteins by NMR. However, it becomes difficult with increasing molecular size due to an increased transverse relaxation and the control of the relaxation pathway is needed to achieve the observation. We demonstrated that even for the large molecular size of 82 kDa Malate synthase G (MSG), the aromatic 13C-1H (CH) peaks of Tryptophan (Trp) and Phenylalanine (Phe) residues can be observed with high quality using a systematic stable isotope labeling scheme, Stereo-Array Isotope Labeling (SAIL) method. However, the sequence specific assignments of these peaks relied on the use of amino acid substitutions, employing an inefficient method that required many isotopes labeled samples. In this study, we developed novel SAIL amino acids that allow for the observation of the aromatic ring δ,ζ and the aliphatic β position peak of Phe residues. The application of TROSY-based experiment to the isolated CH moieties resulted in the successful observation of discernible and resolved CH peaks in Phe residues in MSG. In MSG, the sequence-specific assignments of the backbone and Cβ positions have already been confirmed. Therefore, using this labeling method, the δ and β position peaks of Phe residues can be clearly assigned in a sequence-specific and stereospecific manner through experiments based on intra-residue NOE. Furthermore, the NOESY experiment also allows for the acquisition of information pertaining to the conformation of Phe residues, such as the χ1 dihedral angle, providing valuable insights for the determination of accurate protein structures and in dynamic analysis. This new SAIL amino acids open an avenue to achieve a variety of NMR analysis of large molecular proteins, including a high-resolution structure determination and dynamics and interaction analysis.