PTMdyna: exploring the influence of post-translation modifications on protein conformational dynamics.

Protein post-translational modifications (PTM) play vital roles in cellular regulation, modulating functions by driving changes in protein structure and dynamics. Exploring comprehensively the influence of PTM on conformational dynamics can facilitate the understanding of the related biological function and molecular mechanism. Currently, a series of excellent computation tools have been designed to analyze the time-dependent structural properties of proteins. However, the protocol aimed to explore conformational dynamics of post-translational modified protein is still a blank. To fill this gap, we present PTMdyna to visually predict the conformational dynamics differences between unmodified and modified proteins, thus indicating the influence of specific PTM. PTMdyna exhibits an AUC of 0.884 tested on 220 protein-protein complex structures. The case of heterochromatin protein 1α complexed with lysine 9-methylated histone H3, which is critical for genomic stability and cell differentiation, was used to demonstrate its applicability. PTMdyna provides a reliable platform to predict the influence of PTM on protein dynamics, making it easier to interpret PTM functionality at the structure level. The web server is freely available at http://ccbportal.com/PTMdyna.

HDAC1: An Essential and Conserved Member of the Diverse Zn2+-Dependent HDAC Family Driven by Divergent Selection Pressure.

Zn2+-dependent histone deacetylases (HDACs) are enzymes that regulate gene expression by removing acetyl groups from histone proteins. These enzymes are essential in all living systems, playing key roles in cancer treatment and as potential pesticide targets. Previous phylogenetic analyses of HDAC in certain species have been published. However, their classification and evolutionary origins across biological kingdoms remain unclear, which limits our understanding of them. In this study, we collected the HDAC sequences from 1451 organisms and performed analyses. The HDACs are found to diverge into three classes and seven subclasses under divergent selection pressure. Most subclasses show species specificity, indicating that HDACs have evolved with high plasticity and diversification to adapt to different environmental conditions in different species. In contrast, HDAC1 and HDAC3, belonging to the oldest class, are conserved and crucial in major kingdoms of life, especially HDAC1. These findings lay the groundwork for the future application of HDACs.

LARMD: integration of bioinformatic resources to profile ligand-driven protein dynamics with a case on the activation of estrogen receptor.

Protein dynamics is central to all biological processes, including signal transduction, cellular regulation and biological catalysis. Among them, in-depth exploration of ligand-driven protein dynamics contributes to an optimal understanding of protein function, which is particularly relevant to drug discovery. Hence, a wide range of computational tools have been designed to investigate the important dynamic information in proteins. However, performing and analyzing protein dynamics is still challenging due to the complicated operation steps, giving rise to great difficulty, especially for nonexperts. Moreover, there is a lack of web protocol to provide online facility to investigate and visualize ligand-driven protein dynamics. To this end, in this study, we integrated several bioinformatic tools to develop a protocol, named Ligand and Receptor Molecular Dynamics (LARMD, http://chemyang.ccnu.edu.cn/ccb/server/LARMD/ and http://agroda.gzu.edu.cn:9999/ccb/server/LARMD/), for profiling ligand-driven protein dynamics. To be specific, estrogen receptor (ER) was used as a case to reveal ERß-selective mechanism, which plays a vital role in the treatment of inflammatory diseases and many types of cancers in clinical practice. Two different residues (Ile373/Met421 and Met336/Leu384) in the pocket of ERß/ERα were the significant determinants for selectivity, especially Met336 of ERß. The helix H8, helix H11 and H7-H8 loop influenced the migration of selective agonist (WAY-244). These computational results were consistent with the experimental results. Therefore, LARMD provides a user-friendly online protocol to study the dynamic property of protein and to design new ligand or site-directed mutagenesis.

AIMMS suite: a web server dedicated for prediction of drug resistance on protein mutation.

Drug resistance is one of the most intractable issues for successful treatment in current clinical practice. Although many mutations contributing to drug resistance have been identified, the relationship between the mutations and the related pharmacological profile of drug candidates has yet to be fully elucidated, which is valuable both for the molecular dissection of drug resistance mechanisms and for suggestion of promising treatment strategies to counter resistant. Hence, effective prediction approach for estimating the sensitivity of mutations to agents is a new opportunity that counters drug resistance and creates a high interest in pharmaceutical research. However, this task is always hampered by limited known resistance training samples and accurately estimation of binding affinity. Upon this challenge, we successfully developed Auto In Silico Macromolecular Mutation Scanning (AIMMS), a web server for computer-aided de novo drug resistance prediction for any ligand-protein systems. AIMMS can qualitatively estimate the free energy consequences of any mutations through a fast mutagenesis scanning calculation based on a single molecular dynamics trajectory, which is differentiated with other web services by a statistical learning system. AIMMS suite is available at http://chemyang.ccnu.edu.cn/ccb/server/AIMMS/.

Phylogenetic comparison of 5' splice site determination in central spliceosomal proteins of the U1-70K gene family, in response to developmental cues and stress conditions.

Intron-containing genes have the ability to generate multiple transcript isoforms by splicing, thereby greatly expanding the eukaryotic transcriptome and proteome. In eukaryotic cells, precursor mRNA (pre-mRNA) splicing is performed by a mega-macromolecular complex defined as a spliceosome. Among its splicing components, U1 small nuclear ribonucleoprotein (U1 snRNP) is the smallest subcomplex involved in early spliceosome assembly and 5'-splice site recognition. Its central component, named U1-70K, has been extensively characterized in animals and yeast. Very few investigations on U1-70K genes have been conducted in plants, however. To this end, we performed a comprehensive study to systematically identify 115 U1-70K genes from 67 plant species, ranging from algae to angiosperms. Phylogenetic analysis suggested that the expansion of the plant U1-70K gene family was likely to have been driven by whole-genome duplications. Subsequent comparisons of gene structures, protein domains, promoter regions and conserved splicing patterns indicated that plant U1-70Ks are likely to preserve their conserved molecular function across plant lineages and play an important functional role in response to environmental stresses. Furthermore, genetic analysis using T-DNA insertion mutants suggested that Arabidopsis U1-70K may be involved in response to osmotic stress. Our results provide a general overview of this gene family in Viridiplantae and will act as a reference source for future mechanistic studies on this U1 snRNP-specific splicing factor.

Rational Redesign of Enzyme via the Combination of Quantum Mechanics/Molecular Mechanics, Molecular Dynamics, and Structural Biology Study.

Increasing demands for efficient and versatile chemical reactions have prompted innovations in enzyme engineering. A major challenge in engineering α-ketoglutarate-dependent oxygenases is to develop a rational strategy which can be widely used for directly evolving the desired mutant to generate new products. Herein, we report a strategy for rational redesign of a model enzyme, 4-hydroxyphenylpyruvate dioxygenase (HPPD), based on quantum mechanics/molecular mechanics (QM/MM) calculation and molecular dynamic simulations. This strategy enriched our understanding of the HPPD catalytic reaction pathway and led to the discovery of a series of HPPD mutants producing hydroxyphenylacetate (HPA) as the alternative product other than the native product homogentisate. The predicted HPPD-Fe(IV)═O-HPA intermediate was further confirmed by the crystal structure of Arabidopsis thaliana HPPD/S267W complexed with HPA. These findings not only provide a good understanding of the structure-function relationship of HPPD but also demonstrate a generally applicable platform for the development of biocatalysts.

Structural dynamics and determinants of abscisic acid-receptor binding preference in different aggregation states.

In the 21st century, drought has been the main cause of shortages in world grain production and has created problems with food security. Abscisic acid (ABA) is a key plant hormone involved in the response to abiotic stress, especially drought. The pyrabactin resistance (PYR)/PYR1-like (PYL)/regulatory component of abscisic acid receptor (RCAR) family of proteins (simplified as PYLs) is a well-known ABA receptor family, which can be divided into dimeric and monomeric forms. PYLs can recognize ABA and activate downstream plant drought-resistance signals. However, the difference between monomeric and dimeric receptors in the mechanism of the response to ABA is unclear. Here, we reveal that monomeric receptors have a competitive advantage over dimeric receptors for binding to ABA, driven by the energy penalty resulting from dimer dissociation. ABA also plays different roles with the monomer and the dimer: in the monomer, it acts as a 'conformational stabilizer' for stabilizing the closed gate, whereas for the dimer, it serves as an 'allosteric promoter' for promoting gate closure, which leads to dissociation of the two subunits. This work illustrates how receptor oligomerization could modulate hormonal responses and provides a new concept for novel engineered plants based on ABA binding of monomers.

PIIMS Server: A Web Server for Mutation Hotspot Scanning at the Protein-Protein Interface.

Protein-protein interactions (PPIs) play vital roles in regulating biological processes, such as cellular and signaling pathways. Hotspots are certain residues located at protein-protein interfaces that contribute more in protein-protein binding than other residues. Research on the mutational effects of hotspots is important for understanding basic aspects of protein association. Hence, various computational tools have been developed to explore the impact of mutation hotspots, which will allow a better understanding of the forces that drive PPIs. However, tools that may provide comprehensive substitutions at hotspots are still rare. Hence, there is a strong need for a new free web server to explore mutational effects of hotspots. Herein we introduce a web server named PIIMS that integrates molecular dynamics simulation and one-step free energy perturbation. It contains two main computational functions: (1) computational alanine scanning analysis to identify hotspots and (2) full mutation scanning analysis to evaluate the effects of hotspot mutations. We rigidly validated its ability to predict binding free energy changes by using large and diverse datasets including 1,341 mutations from 50 PPIs with the correlation coefficient R = 0.75. The difference from the existing tools is that PIIMS can perform further evaluation of hotspot residues with regard to their different mutations. The PIIMS web server (accessible at http://chemyang.ccnu.edu.cn/ccb/server/PIIMS/index.php) is free and open to all users without login requirements.

Systematic characterization of the branch point binding protein, splicing factor 1, gene family in plant development and stress responses.

BACKGROUND: Among eukaryotic organisms, alternative splicing is an important process that can generate multiple transcripts from one same precursor messenger RNA, which greatly increase transcriptome and proteome diversity. This process is carried out by a super-protein complex defined as the spliceosome. Specifically, splicing factor 1/branchpoint binding protein (SF1/BBP) is a single protein that can bind to the intronic branchpoint sequence (BPS), connecting the 5' and 3' splice site binding complexes during early spliceosome assembly. The molecular function of this protein has been extensively investigated in yeast, metazoa and mammals. However, its counterpart in plants has been seldomly reported. RESULTS: To this end, we conducted a systematic characterization of the SF1 gene family across plant lineages. In this work, a total of 92 sequences from 59 plant species were identified. Phylogenetic relationships of these sequences were constructed, and subsequent bioinformatic analysis suggested that this family likely originated from an ancient gene transposition duplication event. Most plant species were shown to maintain a single copy of this gene. Furthermore, an additional RNA binding motif (RRM) existed in most members of this gene family in comparison to their animal and yeast counterparts, indicating that their potential role was preserved in the plant lineage. CONCLUSION: Our analysis presents general features of the gene and protein structure of this splicing factor family and will provide fundamental information for further functional studies in plants.

Genome-wide phylogenetic and structural analysis reveals the molecular evolution of the ABA receptor gene family.

The plant hormone abscisic acid (ABA) plays a crucial role during the plant life cycle as well as in adaptive responses to environmental stresses. The core regulatory components of ABA signaling in plants are the pyrabactin resistance1/PYR1-like/regulatory component of ABA receptor family (PYLs), which comprise the largest plant hormone receptor family known. They act as negative regulators of members of the protein phosphatase type 2C family. Due to the biological importance of PYLs, many researchers have focused on their genetic redundancy and consequent functional divergence. However, little is understood of their evolution and its impact on the generation of regulatory diversity. In this study, we identify positive selection and functional divergence in PYLs through phylogenetic reconstruction, gene structure and expression pattern analysis, positive selection analysis, functional divergence analysis, and structure comparison. We found the correlation of desensitization of PYLs under specific modifications in the molecular recognition domain with functional diversification. Hence, an interesting antagonistic co-evolutionary mechanism is proposed for the functional diversification of ABA receptor family proteins. We believe a compensatory evolutionary pathway may have occurred.

Identification, evolution and alternative splicing profile analysis of the splicing factor 30 (SPF30) in plant species.

MAIN CONCLUSION: The work offers a comprehensive evaluation on the phylogenetics and conservation of splicing patterns of the plant SPF30 splicing factor gene family. In eukaryotes, one pre-mRNA can generate multiple mRNA transcripts by alternative splicing (AS), which expands transcriptome and proteome diversity. Splicing factor 30 (SPF30), also known as survival motor neuron domain containing protein 1 (SMNDC1), is a spliceosomal protein that plays an essential role in spliceosomal assembly. Although SPF30 genes have been well characterised in human and yeast, little is known about their homologues in plants. Here, we report the genome-wide identification and phylogenetic analysis of SPF30 genes in the plant kingdom. In total, 82 SPF30 genes were found in 64 plant species from algae to land plants. Alternative transcripts were found in many SPF30 genes and splicing profile analysis revealed that the second intron in SPF30 genome is frequently associated with AS events and contributed to the birth of novel exons in a few SPF30 members. In addition, different conserved sequences were observed at these putative splice sites among moss, monocots and dicots, respectively. Our findings will facilitate further functional characterization of plant SPF30 genes as putative splicing factors.

Genome-wide identification and functional analysis of the splicing component SYF2/NTC31/p29 across different plant species.

MAIN CONCLUSION: This study systematically identifies plant SYF2/NTC31/p29 genes from 62 plant species by a combinatory bioinformatics approach, revealing the importance of this gene family in phylogenetics, duplication, transcriptional, and post-transcriptional regulation. Alternative splicing is a post-transcriptional regulatory mechanism, which is critical for plant development and stress responses. The entire process is strictly attenuated by a complex of splicing-related proteins, designated splicing factors. Human p29, also referred to as synthetic lethal with cdc forty 2 (SYF2) or the NineTeen complex 31 (NTC31), is a core protein found in the NTC complex of humans and yeast. This splicing factor participates in a variety of biological processes, including DNA damage repair, control of the cell cycle, splicing, and tumorigenesis. However, its function in plants has been seldom reported. Thus, we have systematically identified 89 putative plant SYF2s from 62 plant species among the deposited entries in the Phytozome database. The phylogenetic relationships and evolutionary history among these plant SYF2s were carefully examined. The results revealed that plant SYF2s exhibited distinct patterns regarding their gene structure, promoter sequences, and expression levels, suggesting their functional diversity in response to developmental cues or stress treatments. Although local duplication events, such as tandem duplication and retrotransposition, were found among several plant species, most of the plant species contained only one copy of SYF2, suggesting the existence of additional mechanisms to confer duplication resistance. Further investigation using the model dicot and monocot representatives Arabidopsis and rice SYF2s indicated that the splicing pattern and resulting protein isoforms might play an alternative role in the functional diversity.

Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings.

In eukaryotes, mechanisms such as alternative splicing (AS) and alternative translation initiation (ATI) contribute to organismal protein diversity. Specifically, splicing factors play crucial roles in responses to environment and development cues; however, the underlying mechanisms are not well investigated in plants. Here, we report the parallel employment of short-read RNA sequencing, single molecule long-read sequencing and proteomic identification to unravel AS isoforms and previously unannotated proteins in response to abscisic acid (ABA) treatment. Combining the data from the two sequencing methods, approximately 83.4% of intron-containing genes were alternatively spliced. Two AS types, which are referred to as alternative first exon (AFE) and alternative last exon (ALE), were more abundant than intron retention (IR); however, by contrast to AS events detected under normal conditions, differentially expressed AS isoforms were more likely to be translated. ABA extensively affects the AS pattern, indicated by the increasing number of non-conventional splicing sites. This work also identified thousands of unannotated peptides and proteins by ATI based on mass spectrometry and a virtual peptide library deduced from both strands of coding regions within the Arabidopsis genome. The results enhance our understanding of AS and alternative translation mechanisms under normal conditions, and in response to ABA treatment.

PADFrag: A Database Built for the Exploration of Bioactive Fragment Space for Drug Discovery.

Structural analyses of drugs and pesticides can enable the identification of new bioactive compounds with novel and diverse scaffolds as well as improve our understanding of the bioactive fragment space. The Pesticide And Drug Fragments (PADFrag) database is a unique bioinformatic-cheminformatic cross-referencing resource that combines detailed bioactive fragment data and potential targets with a strong focus on quantitative, analytic, and molecular-scale information for the exploration of bioactive fragment space for drug discovery ( http://chemyang.ccnu.edu.cn/ccb/database/PADFrag/ ). The main applications of PADFrag are the analysis of the privileged structures within known bioactive molecules, ab initio molecule library design, and core fragment discovery for fragment-based drug design. Other potential applications include prediction of fragment interactions and general pharmaceutical research.

One-Pot Approach to N-Quinolyl 3'/4'-Biaryl Carboxamides by Microwave-Assisted Suzuki-Miyaura Coupling and N-Boc Deprotection.

N-Quinolyl biaryl carboxamides have received tremendous attention for their notable biological properties. Here we have described a general protocol for the preparation of N-quinolyl 3'/4'-biaryl carboxamides by microwave-assisted Suzuki-Miyaura cross-coupling reaction and N-Boc deprotection in one pot. This method, which did not require acids, was used to produce a series of N-quinolyl 3'/4'-biaryl carboxamides with excellent functional group tolerance and high yields (70% to 95%).

Chromosome-level genome assembly of an oligophagous leaf beetle Ophraella communa (Coleoptera: Chrysomelidae).

The leaf beetle Ophraella communa LeSage (Coleoptera: Chrysomelidae) is an effective biological control agent of the common ragweed. Here, we assembled a chromosome-level genome of the O. communa by combining Illumina, Nanopore, and Hi-C sequencing technologies. The genome size of the final genome assembly is 733.1 Mb, encompassing 17 chromosomes, with an improved contig N50 of 7.05 Mb compared to the original version. Genome annotation reveals 25,873 protein-coding genes, with functional annotations available for 22,084 genes (85.35%). Non-coding sequence annotation identified 204 rRNAs, 626 tRNAs, and 1791 small RNAs. Repetitive elements occupy 414.41 Mb, constituting 57.76% of the genome. This high-quality genome is fundamental for advancing biological control strategies employing O. communa.

CIPDB: A biological structure databank for studying cation and π interactions.

As major forces for modulating protein folding and molecular recognition, cation and π interactions are extensively identified in protein structures. They are even more competitive than hydrogen bonds in molecular recognition, thus, are vital in numerous biological processes. In this review, we introduce the methods for the identification and quantification of cation and π interactions, provide insights into the characteristics of cation and π interactions in the natural state, and reveal their biological function together with our developed database (Cation and π Interaction in Protein Data Bank; CIPDB; http://chemyang.ccnu.edu.cn/ccb/database/CIPDB). This review lays the foundation for the in-depth study of cation and π interactions and will guide the use of molecular design for drug discovery.

Design and Synthesis of Novel Oxathiapiprolin Derivatives as Oxysterol Binding Protein Inhibitors and Their Application in Phytopathogenic Oomycetes.

Oomycetes, particularly those from the genus Phytophthora, are significant threats to global food security and natural ecosystems. Oxathiapiprolin (OXA) is an effective oomycete fungicide that targets an oxysterol binding protein (OSBP), while the binding mechanism of OXA is still unclear, which limits the pesticide design, induced by the low sequence identity of Phytophthora and template models. Herein, we generated the OSBP model of the well-reported Phytophthora capsici using AlphaFold 2 and studied the binding mechanism of OXA. Based on it, a series of OXA analogues were designed. Then, compound 2l, the most potent candidate, was successfully designed and synthesized, showing a control efficiency comparable to that of OXA. Moreover, field trial experiments showed that 2l exhibited nearly the same activity (72.4%) as OXA against cucumber downy mildew at 25 g/ha. The present work indicated that 2l could be used as a leading compound for the discovery of new OSBP fungicides.