Automated protein structure modeling with SWISS-MODEL Workspace and the Protein Model Portal

Characterization of a novel multifunctional glycoside hydrolase family in the metagenome-assembled genomes of horse gut

The gut microbiota is a treasure trove of carbohydrate-active enzymes (CAZymes). To explore novel and efficient CAZymes, we analyzed the 4,142 metagenome-assembled genomes (MAGs) of the horse gut microbiota and found the MAG117.bin13 genome (Bacteroides fragilis) contains the highest number of polysaccharide utilisation loci sites (PULs), indicating its high capability for carbohydrate degradation. Bioinformatics analysis indicate that the PULs region of the MAG117.bin13 genome encodes many hypothetical proteins, which are important sources for exploring novel CAZymes. Interestingly, we discovered a hypothetical protein (595 amino acids). This protein exhibits potential CAZymes activity and has a lower similarity to CAZymes, we named it BfLac2275. We purified the protein using prokaryotic expression technology and studied its enzymatic function. The hydrolysis experiment of the polysaccharide substrate showed that the BfLac2275 protein has the ability to degrade α-lactose (156.94 U/mg), maltose (92.59 U/mg), raffinose (86.81 U/mg), and hyaluronic acid (5.71 U/mg). The enzyme activity is optimal at pH 5.0 and 30 ℃, indicating that the hypothetical protein BfLac2275 is a novel and multifunctional CAZymes in the glycoside hydrolases (GHs). These properties indicate that BfLac2275 has broad application prospects in many fields such as plant polysaccharide decomposition, food industry, animal feed additives and enzyme preparations. This study not only serves as a reference for exploring novel CAZymes encoded by gut microbiota but also provides an example for further studying the functional annotation of hypothetical genes in metagenomic assembly genomes.

MAO inhibiting phytochemicals from the roots of Glycyrrhiza glabra L

Glycyrrhizin, a natural compound that is substantially present in Glycyrrhiza glabra L. (Gg) root. Monoamine oxidase B (MAOB) inhibitor is used for the treatment of several important neuropsychological diseases like Parkinson's disease. Gg is known to possess psychoactive properties which relates to its MAO inhibitory potential. This study sought to determine the MAO inhibition property of glycyrrhizin from Gg root extract. The Aqueous extract containing glycyrrhizin was isolated from the root of Gg and characterized using TLC, HPLC, and LC-MS techniques. In silico docking was conducted using Extra precision Glide 2018, Schrödinger docking suite. In addition, the pharmacokinetic properties of the compounds were predicted using SwissADME. The binding energies of the glycyrrhizin correlated well with their in vitro MAO inhibitory potential. Glycyrrhizin exhibited potent inhibitory activity towards MAOB whereas, an aqueous extract of Gg root inhibits both A and B forms of MAO enzyme. Further, molecular docking and molecular dynamics simulation showed that liquiritigenin and methoxyglabridin showed higher stability than other inhibitor compounds from the Gg root extract. These observations suggest that the phytochemicals from the Gg root extract have potent MAO inhibition properties, which can be exploited for the treatment of neurodegenerative disorders.Communicated by Ramaswamy H. Sarma.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Proteomic analysis of boar seminal plasma: Putative markers for fertility based on capacity of semen preservation at 17°C

Boar seminal plasma (SP) proteins were associated with differences on sperm resistance to cooling at 17°C. However, information about seminal plasma proteins in boars classified by capacity of semen preservation and in vivo fertility remains lacking. Thus, the objective was to evaluate the SP proteome in boars classified by capacity of semen preservation and putative biomarkers for fertility. The ejaculates from high-preservation (HP) showed higher progressive motility during all 5 days than the low-preservation (LP) boars. There was no difference for farrowing rate between ejaculates from LP (89.7%) and HP boars (88.4%). The LP boars presented lower total piglets born (14.0 ± 0.2) than HP (14.8 ± 0.2; p < 0.01). A total of 257 proteins were identified, where 184 were present in both classes of boar, and 41 and 32 were identified only in LP and HP boars, respectively. Nine proteins were differently expressed: five were more abundant in HP (SPMI, ZPBP1, FN1, HPX, and C3) and four in LP boars (B2M, COL1A1, NKX3-2, and MPZL1). The HP boars had an increased abundance of SP proteins related to sperm resistance and fecundation process which explains the better TPB. LP boars had a higher abundance of SP proteins associated with impaired spermatogenesis.

Alkalihalobacillus clausii PA21 transcriptome profiling and functional analysis revealed the metabolic pathway involved in glycoalkaloids degradation

Glycoalkaloids (GAs), including α-solanine and α-chaconine, are secondary metabolites found in potato, which are toxic to higher animals. In a previous study, Alkalihalobacillus clausii PA21 showed the capacity to degrade GAs. Herein, the transcriptome response of PA21 to α-solanine or α-chaconine was evaluated. In total, 3170 and 2783 differential expressed genes (DEGs) were found in α-solanine- and α-chaconine-treated groups, respectively, with most DEGs upregulated. Moreover, GAs activated transmembrane transport, carbohydrate metabolism, transcription, quorum sensing, and bacterial chemotaxis in PA21 to withstand GA-induced stress and promote GAs degradation. Furthermore, qRT-PCR analysis confirmed the upregulation of degrading enzymes and components involved in GA degradation in PA21. In addition, the GAs-degrading enzymes were heterologous expressed, purified, and incubated with GAs to analyze the degradation products. The results showed that α-solanine was degraded to β1-solanine, β2-solanine, γ-solanine, and solanidine by β-glucosidase, α-rhamnosidase, and β-galactosidase. Meanwhile, α-chaconine was degraded to β1-chaconine, β2-chaconine, γ-chaconine, and solanidine by β-glucosidase and α-rhamnosidase. Overall, the molecular mechanism underlying GAs degradation by PA21 was revealed by RNAseq combined with protein expression and function studies, thus providing the basis for the development of engineered bacteria that can efficiently degrade GAs to promote their application in the control of GAs in potatoes.

Insights into the Bioinformatics and Transcriptional Analysis of the Elongator Complexes (ELPs) Gene Family of Wheat: TaELPs Contribute to Wheat Abiotic Stress Tolerance and Leaf Senescence

Elongator complexes (ELPs) are the protein complexes that promote transcription through histone acetylation in eukaryotic cells and interact with elongating RNA polymerase II (RNAPII). ELPs' role in plant growth and development, signal transduction, and response to biotic and abiotic stresses have been confirmed in model plants. However, the functions of the wheat ELP genes are not well documented. The present study identified 18 members of the ELPs from the wheat genome with a homology search. Further, bioinformatics and transcription patterns in response to different stress conditions were analyzed to dissect their potential regulatory mechanisms in wheat. Gene duplication analysis showed that 18 pairs of ELP paralogous genes were derived from segmental duplication, which was divided into six clades by protein phylogenetic and cluster analysis. The orthologous analysis of wheat TaELP genes showed that TaELP genes may have evolved from orthologous genes of other plant species or closely related plants. Moreover, a variety of cis-acting regulatory elements (CAREs) related to growth and development, hormone response, and biotic and abiotic stresses were identified in the TaELPs' promoter regions. The qRT-PCR analysis showed that the transcription of TaELPs was induced under hormone, salt, and drought stress and during leaf senescence. The TaELP2 gene was silenced with BSMV-VIGS, and TaELP2 was preliminarily verified to be involved in the regulation of wheat leaf senescence. Overall, TaELP genes might be regulated by hormone signaling pathways and response to abiotic stress and leaf senescence, which could be investigated further as potential candidate genes for wheat abiotic stress tolerance and yield improvement.

Identification of common genes in obesity and cancer through network interaction and targeting those genes by virtual screening approach

Obesity may have an effect on cancer outcomes, resulting in global inequalities in cancer survival and death. Microarray data analysis was done to identify differentially expressed genes (DEGs) in obese and cancer patients. Total 1977 differentially expressed genes among obesity and gastric cancer, breast cancer, pancreatic cancer, and colorectal cancer were used to build a gene interaction network, which was then analyzed by using Cytoscape software. It has been identified that JUN, CXCL12, and LEP genes show a higher degree and stress, and play an important role in obesity and cancer progression. Further, CXCL12 and LEP were taken for virtual screening study with coumarin and its derivatives to develop a drug against obesity and cancer. The interactions of CXCL12 and LEP with coumarins were studied by molecular docking and it shows good interaction as well as docking score as compared to the standard one. The ADME properties were predicted to check the drug-likeness activity of coumarins and the most of the drug-likeness activities are in admire range. The Binding free energy of the docked complex was calculated by performing MM-GBSA. The molecular docking, ADME properties prediction, and MM-GBSA was performed on Maestro 12.6. The top docked score compounds were further subjected to molecular dynamic simulation to check the stability by using GROMACS. The MM-PBSA study was performed to calculate the binding energy components as well as the energy contributions of specific amino acids. The resultant compounds could be a potent anti-obesity and anti-cancer drug.Communicated by Ramaswamy H. Sarma.

Lactobacillus for ribosome peptide editing cancer

This study reviews newly discovered insect peptide point mutations as new possible cancer research targets. To interpret newly discovered peptide point mutations in insects as new possible cancer research targets, we focused on the numerous peptide changes found in the 'CSP' family on the sex pheromone gland of the female silkworm moth Bombyx mori. We predict that the Bombyx peptide modifications will have a significant effect on cancer CUP (cancers of unknown primary) therapy and that bacterial peptide editing techniques, specifically Lactobacillus combined to CRISPR, will be used to regulate ribosomes and treat cancer in humans.

Birth-and-death evolution of ribonuclease 9 genes in Cetartiodactyla

RNase9 plays a reproductive function and has been recognized as an important member of the ribonuclease (RNase) A superfamily, a gene family that is widely used as a model for molecular evolutionary studies. Here, we identified 178 RNase9 genes from 95 Cetartiodactyla species that represent all four lineages and 21 families of this clade. Unexpectedly, RNase9 experienced an evolutionary scenario of "birth and death" in Ruminantia, and expression analyses showed that duplicated RNase9A and RNase9B genes are expressed in reproductive tissues (epididymis, vas deferens or prostate). This expression pattern combined with the estimate that these genes duplicated during the middle Eocene, a time when Ruminantia become a successful lineage, suggests that the RNase9 gene duplication might have been advantageous for promoting sperm motility and male fertility as an adaptation to climate seasonality changes of this period. In contrast, all RNase9 genes were lost in the Cetacean lineage, which might be associated with their high levels of prostatic lesions and lower reproductive rates as adaptations to a fully aquatic environment and a balance to the demands of ocean resources. This study reveals a complex and intriguing evolutionary history and functional divergence for RNase9 in Cetartiodactyla, providing new insights into the evolution of the RNaseA superfamily and molecular mechanisms for organismal adaptations to the environment.

The nepenthesin insert in the Plasmodium falciparum aspartic protease plasmepsin V is necessary for enzyme function

Plasmepsin V (PM V) is a pepsin-like aspartic protease essential for growth of the malarial parasite Plasmodium falciparum. Previous work has shown PM V to be an endoplasmic reticulum-resident protease that processes parasite proteins destined for export into the host cell. Depletion or inhibition of the enzyme is lethal during asexual replication within red blood cells as well as during the formation of sexual stage gametocytes. The structure of the Plasmodium vivax PM V has been characterized by X-ray crystallography, revealing a canonical pepsin fold punctuated by structural features uncommon to secretory aspartic proteases; however, the function of this unique structure is unclear. Here, we used parasite genetics to probe these structural features by attempting to rescue lethal PM V depletion with various mutant enzymes. We found an unusual nepenthesin 1-type insert in the PM V gene to be essential for parasite growth and PM V activity. Mutagenesis of the nepenthesin insert suggests that both its amino acid sequence and one of the two disulfide bonds that undergird its structure are required for the insert's role in PM V function. Furthermore, molecular dynamics simulations paired with Markov state modeling suggest that mutations to the nepenthesin insert may allosterically affect PM V catalysis through multiple mechanisms. Taken together, these data provide further insights into the structure of the P. falciparum PM V protease.

Enhancing regioselectivity of sucrose phosphorylase by loop engineering for glycosylation of L-ascorbic acid

Sucrose phosphorylase (SPase) has a remarkable capacity to synthesize numerous glucosides from abundantly available sucrose under mild conditions but suffers from specificity and regioselectivity issues. In this study, a loop engineering strategy was introduced to enhance the regioselectivity and substrate specificity of SPase for the efficient synthesis of 2-O-α-D-glucopyranosyl-L-ascorbic acid (AA-2G) via L-ascorbic acid (L-AA). P134, L341, and L343 were identified as "hotspots" for modulating the flexibility of loops, which significantly influenced the H-bonding network of L-AA in the active site, as well as the entrance of the substrate channel, thereby altering the regioselectivity and substrate specificity. Finally, the mutant L341V/L343F, with near-perfect control of the selectivity synthesis of the 2-OH group of L-AA (> 99%), was obtained. The AA-2G production by the mutant reached 244 g L^-1 in a whole-cell biotransformation system, and the conversion rate of L-AA reached 64%, which is the highest level reported to date. Our work also provides a successful loop engineering case for modulating the regioselectivity and specificity of sucrose phosphorylase. KEY POINTS: • "Hotspots" were identified in the flexible loops of sucrose phosphorylase. • Mutants exhibited improved regioselectivity and specificity against L-ascorbic acid. • Synthesized AA-2G with high yield and regioselectivity by whole-cell of mutant.

A Novel Chondroitin AC Lyase With Broad Substrate Specificity From Pedobacter rhizosphaerae: Cloning, Expression, and Characterization

Chondroitin AC lyase (ChSaseAC) is one of the essential polysaccharides lyases in low molecular chondroitin sulfate production. In this work, a novel PrChSaseAC from Pedobacter rhizosphaerae was successfully cloned, expressed in Escherichia coli. After optimizing the induction, the recombinant PrChSaseAC could be expressed efficiently at 0.1 mM IPTG, 25°C, and 12 h induction. Then, it was purified with Ni-NTA affinity chromatography. The characterization of the purified PrChSaseAC showed that it had high specific activity and good storage stability, which would favor the production of low molecular weight chondroitin sulfate. It also displayed activity toward chondroitin sulfate C and hyaluronic acid. PrChSaseAC had the highest activity at pH 7.5, 37°C, 10 mM Ca²⁺, and 5 mg/ml of chondroitin sulfate A. Molecular docking of substrate and enzyme showed the interactions between the enzyme and substrate; it revealed that the enzyme showed high activity to CS-A and hyaluronic acid, but lower activity to CS-C attributed to the structure of the binding pocket. The high stability and specific activity of the enzyme will benefit the industrial production or clinical treatment.

In Silico and Transcription Analysis of Trehalose-6-phosphate Phosphatase Gene Family of Wheat: Trehalose Synthesis Genes Contribute to Salinity, Drought Stress and Leaf Senescence

Trehalose-6-phosphate phosphatase (TPP) genes take part in trehalose metabolism and also in stress tolerance, which has been well documented in many species but poorly understood in wheat. The present research has identified a family of 31 TPP genes in Triticum aestivum L. through homology searches and classified them into five clades by phylogenetic tree analysis, providing evidence of an evolutionary status with Hordeum vulgare, Brachypodium distachyon and Oryza sativa. The exon-intron distribution revealed a discrete evolutionary history and projected possible gene duplication occurrences. Furthermore, different computational approaches were used to analyze the physical and chemical properties, conserved domains and motifs, subcellular and chromosomal localization, and three-dimensional (3-D) protein structures. Cis-regulatory elements (CREs) analysis predicted that TaTPP promoters consist of CREs related to plant growth and development, hormones, and stress. Transcriptional analysis revealed that the transcription levels of TaTPPs were variable in different developmental stages and organs. In addition, qRT-PCR analysis showed that different TaTPPs were induced under salt and drought stresses and during leaf senescence. Therefore, the findings of the present study give fundamental genomic information and possible biological functions of the TaTPP gene family in wheat and will provide the path for a better understanding of TaTPPs involvement in wheat developmental processes, stress tolerance, and leaf senescence.

Structure prediction of honey bee vitellogenin: a multi-domain protein important for insect immunity

Vitellogenin (Vg) has been implicated as a central protein in the immunity of egg‐laying animals. Studies on a diverse set of species suggest that Vg supports health and longevity through binding to pathogens. Specific studies of honey bees (Apis mellifera) further indicate that the vitellogenin (vg) gene undergoes selection driven by local pathogen pressures. Determining the complete 3D structure of full‐length Vg (flVg) protein will provide insights regarding the structure–function relationships underlying allelic variation. Honey bee Vg has been described in terms of function, and two subdomains have been structurally described, while information about the other domains is lacking. Here, we present a structure prediction, restrained by experimental data, of flVg from honey bees. To achieve this, we performed homology modeling and used AlphaFold before using a negative‐stain electron microscopy map to restrict, orient, and validate our 3D model. Our approach identified a highly conserved Ca²⁺‐ion‐binding site in a von Willebrand factor domain that might be central to Vg function. Thereafter, we used rigid‐body fitting to predict the relative position of high‐resolution domains in a flVg model. This mapping represents the first experimentally validated full‐length protein model of a Vg protein and is thus relevant for understanding Vg in numerous species. Our results are also specifically relevant to honey bee health, which is a topic of global concern due to rapidly declining pollinator numbers.

Predicted Structure of Fully Activated Tas1R3/1R3' Homodimer Bound to G Protein and Natural Sugars: Structural Insights into G Protein Activation by a Class C Sweet Taste Homodimer with Natural Sugars

The Tas1R3 G protein-coupled receptor constitutes the main component of sweet taste sensory response in humans via forming a heterodimer with Tas1R2 or a homodimer with Tas1R3. The Tas1R3/1R3' homodimer serves as a low-affinity sweet taste receptor, stimulating gustducin G protein (GGust) signaling in the presence of a high concentration of natural sugars. This provides an additional means to detect the taste of natural sugars, thereby differentiating the flavors between natural sugars and artificial sweeteners. We report here the predicted 3D structure of active state Tas1R3/1R3' homodimer complexed with heterotrimeric GGust and sucrose. We discovered that the GGust makes ionic anchors to intracellular loops 1 and 2 of Tas1R3 while the Gα-α5 helix engages the cytoplasmic region extensively through salt bridge and hydrophobic interactions. We show that in the activation of this complex the Venus flytrap domains of the homodimer undergo a remarkable twist up to ∼100° rotation around the vertical axis to adopt a closed-closed conformation while the intracellular region relaxes to an open-open conformation. We find that binding of sucrose to the homodimer stabilizes a preactivated conformation with a largely open intracellular region that recruits and activates the GGust. Upon activation, the Gα subunit spontaneously opens up the nucleotide-binding site, making nucleotide exchange facile for signaling. This activation of GGust promotes the interdomain twist of the Venus flytrap domains. These structures and transformations could potentially be a basis for the design of new sweeteners with higher activity and less unpleasant flavors.

Mucoadhesive Low Molecular Chitosan Complexes to Protect rHuKGF from Proteolysis: In-vitro Characterization and FHs 74 Int Cell Proliferation Studies

BACKGROUND

Recombinant Keratinocyte Growth Factor (rHuKGF) is a therapeutic protein used widely in oral mucositis after chemotherapy in various cancers, stimulating lung morphogenesis and gastrointestinal tract cell proliferation. In this research study, chitosan-rHuKGF polymeric complex was implemented to improve the stability of rHuKGF and used as rejuvenation therapy for the treatment of oral mucositis in cancer patients.

OBJECTIVE

Complexation of rHuKGF with mucoadhesive low molecular weight chitosan to protect rHuKGF from proteolysis and investigate the effect of chitosan-rHuKGF complex on the proliferation rate of FHs 74 Int cells.

METHODS

The interaction between chitosan and rHuKGF was studied by molecular docking. Malvern ZetaSizer Nano Zs and Fourier-Transform Infrared spectroscopy (FTIR) tests were carried out to characterize the chitosan-rHuKGF complex. In addition, SDS-PAGE was performed to investigate the interaction between chitosan-rHuKGF complex and pepsin. The effect of chitosan-rHuKGF complex on the proliferation rate of FHs 74 Int cells was studied by MTT assay.

RESULTS

Chitosan-rHuKGF complex was formed through the hydrogen bonding proven by the docking studies. A stable chitosan-rHuKGF complex was formed at pH 4.5 and was protected from proteolysis and assessed by SDS PAGE. According to the MTT assay results, chitosan-rHuKGF complex increased the cell proliferation rate of FHs 74 Int cells.

CONCLUSION

The developed complex improved the stability and the biological function of rHuKGF.

Cold Denaturation of Proteins: Where Bioinformatics Meets Thermodynamics to Offer a Mechanistic Understanding: Pea Protein As a Case Study

Protein structure can be altered with heat, but models which predict denaturation show that globular proteins also spontaneously unfold at low temperatures through cold denaturation. By an analysis of the primary structure of pea protein using bioinformatic modeling, a mechanism of pea protein cold denaturation is proposed. Pea protein is then fractionated into partially purified legumin and vicilin components, suspended in ethanol, and subjected to low temperatures (-10 to -20 °C). The structural characterizations of the purified fractions are conducted through FTIR, ζ potential, dynamic light scattering, and oil binding, and these are compared to the results of commercial protein isolates. The observed structural changes suggest that pea protein undergoes changes in structure as the result of low-temperature treatments, which could lead to innovative industrial processing techniques for functionalization by low-temperature processing.

Bioinformatic Analysis of Structure and Function of LIM Domains of Human Zyxin Family Proteins

Members of the human Zyxin family are LIM domain-containing proteins that perform critical cellular functions and are indispensable for cellular integrity. Despite their importance, not much is known about their structure, functions, interactions and dynamics. To provide insights into these, we used a set of in-silico tools and databases and analyzed their amino acid sequence, phylogeny, post-translational modifications, structure-dynamics, molecular interactions, and functions. Our analysis revealed that zyxin members are ohnologs. Presence of a conserved nuclear export signal composed of LxxLxL/LxxxLxL consensus sequence, as well as a possible nuclear localization signal, suggesting that Zyxin family members may have nuclear and cytoplasmic roles. The molecular modeling and structural analysis indicated that Zyxin family LIM domains share similarities with transcriptional regulators and have positively charged electrostatic patches, which may indicate that they have previously unanticipated nucleic acid binding properties. Intrinsic dynamics analysis of Lim domains suggest that only Lim1 has similar internal dynamics properties, unlike Lim2/3. Furthermore, we analyzed protein expression and mutational frequency in various malignancies, as well as mapped protein-protein interaction networks they are involved in. Overall, our comprehensive bioinformatic analysis suggests that these proteins may play important roles in mediating protein-protein and protein-nucleic acid interactions.

A computer-aided approach to identify novel Leishmania major protein disulfide isomerase inhibitors for treatment of leishmaniasis

Leishmaniasis is an infectious disease caused by parasites of the genus Leishmania and transmitted by the bite of a sand fly. To date, most available drugs for treatment are toxic and beyond the economic means of those affected by the disease. Protein disulfide isomerase (PDI) is a chaperone protein that plays a major role in the folding of newly synthesized proteins, specifically assisting in disulfide bond formation, breakage, or rearrangement in all non-native proteins. In previous work, we demonstrated that Leishmania major PDI (LmPDI) has an essential role in pathogen virulence. Furthermore, inhibition of LmPDI further blocked parasite infection in macrophages. In this study, we utilized a computer-aided approach to design a series of LmPDI inhibitors. Fragment-based virtual screening allowed for the understanding of the inhibitors' modes of action on LmPDI active sites. The generated compounds obtained after multiple rounds of virtual screening were synthesized and significantly inhibited target LmPDI reductase activity and were shown to decrease in vitro parasite growth in human monocyte-derived macrophages. This novel cheminformatics and synthetic approach led to the identification of a new series of compounds that might be optimized into novel drugs, likely more specific and less toxic for the treatment of leishmaniasis.

Exploiting the role of stereoelectronic effects to design the antagonists of the human complement C3a receptor

Stereoelectronic effects are crucial in governing the conformational behaviour of small molecules bearing heterocyclic rings adjacent to amides.

A highly active heparinase I from Bacteroides cellulosilyticus: Cloning, high level expression, and molecular characterization

As one of the most extensively studied glycosaminoglycan lyases, heparinase I has been used in producing low or ultra-low molecular weight heparin. Its' important applications are to neutralize the heparin in human blood and analyze heparin structure in the clinic. However, the low productivity and activity of the enzyme have greatly hindered its applications. In this study, a novel Hep-I from Bacteroides cellulosilyticus (BcHep-I) was successfully cloned and heterologously expressed in E. coli BL21 (DE3) as a soluble protein. The molecular mass and isoelectric point (pI) of the enzyme are 44.42 kDa and 9.02, respectively. And the characterization of BcHep-I after purified with Ni-NTA affinity chromatography suggested that it is a mesophilic enzyme. BcHep-I can be activated by 1 mM Ca2+, Mg2+, and Mn2+, while severely inhibited by Zn2+, Co2+, and EDTA. The specific activity of the enzyme was 738.3 U·mg-1 which is the highest activity ever reported. The Km and Vmax were calculated as 0.17 mg·mL-1 and 740.58 U·mg-1, respectively. Besides, the half-life of 300 min at 30°C showed BcHep-I has practical applications. Homology modeling and substrate docking revealed that Gln15, Lys74, Arg76, Lys104, Arg149, Gln208, Tyr336, Tyr342, and Lys338 were mainly involved in the substrate binding of Hep-I, and 11 hydrogen bonds were formed between heparin and the enzyme. These results indicated that BcHep-I with high activity has great potential applications in the industrial production of heparin, especially in the clinic to neutralize heparin.

Multimodal small-molecule screening for human prion protein binders

Prion disease is a rapidly progressive neurodegenerative disorder caused by misfolding and aggregation of the prion protein (PrP), and there are currently no therapeutic options. PrP ligands could theoretically antagonize prion formation by protecting the native protein from misfolding or by targeting it for degradation, but no validated small-molecule binders have been discovered to date. We deployed a variety of screening methods in an effort to discover binders of PrP, including 19F-observed and saturation transfer difference (STD) NMR spectroscopy, differential scanning fluorimetry (DSF), DNA-encoded library selection, and in silico screening. A single benzimidazole compound was confirmed in concentration-response, but affinity was very weak (Kd > 1 mm), and it could not be advanced further. The exceptionally low hit rate observed here suggests that PrP is a difficult target for small-molecule binders. Whereas orthogonal binder discovery methods could yield high-affinity compounds, non-small-molecule modalities may offer independent paths forward against prion disease.

Schistosomal Sulfotransferase Interaction with Oxamniquine Involves Hybrid Mechanism of Induced-fit and Conformational Selection

Background:

Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention.

Objective:

The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations.

Methods:

Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA.

Results:

The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol.

Conclusion:

The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.

The crystal structure of protein-transporting chaperone BCP1 from Saccharomyces cerevisiae

BCP1 is a protein enriched in the nucleus that is required for Mss4 nuclear export and identified as the chaperone of ribosomal protein Rpl23 in Saccharomyces cerevisiae. According to sequence homology, BCP1 is related to the mammalian BRCA2-interacting protein BCCIP and belongs to the BCIP protein family (PF13862) in the Pfam database. However, the BCIP family has no discernible similarity to proteins with known structure. Here, we report the crystal structure of BCP1, presenting an α/β fold in which the central antiparallel β-sheet is flanked by helices. Protein structural classification revealed that BCP1 has similarity to the GNAT superfamily but no conserved substrate-binding residues. Further modeling and protein-protein docking work provide a plausible model to explain the interaction between BCP1 and Rpl23. Our structural analysis presents the first structure of BCIP family and provides a foundation for understanding the molecular basis of BCP1 as a chaperone of Rpl23 for ribosome biosynthesis.

Genetic polymorphisms and transcription profiles associated with intracranial aneurysm: a key role for NOTCH3

Intracranial aneurysm (IA) incidence is about 1~2%. However, the specific mechanisms of IA onset and development need further study. Our objective was to discover novel IA-related genes to determine possible etiologies further. We performed next-generation sequencing on nineteen Chinese patients with familial IA and one patient with sporadic IA. We obtained mRNA expression data of 129 samples from Gene Expression Omnibus (GEO) and made statistical computing to discover differentially expressed genes (DEGs). The screened IA-related gene NOTCH3 was determined by bioinformatic data mining. We verified the IA-related indicators of NOTCH3. Association was found between IA and the NOTCH3 SNPs rs779314594, rs200504060 and rs2285981. Levels of NOTCH3 mRNA were lower in IA tissue than in control tissue, but higher in peripheral blood neutrophils from IA patients than in neutrophils from controls. Levels of NOTCH3 protein were lower in IA tissue than in cerebral artery tissue. NOTCH3 also decreased the expression of angiogenesis factors in human umbilical vein endothelial cells. Variation in NOTCH3 and alteration of its expression in cerebral artery or neutrophils may contribute to IA. Our findings also describe a bioinformatic-experimental approach that may prove useful for probing the pathophysiology of other complex diseases.

Subnanometer structures of HIV-1 envelope trimers on aldrithiol-2-inactivated virus particles

The HIV-1 envelope glycoprotein (Env) trimer, composed of gp120 and gp41 subunits, mediates viral entry into cells. Recombinant Env trimers have been studied structurally, but characterization of Env embedded in intact virus membranes has been limited to low resolution. Here, we deploy cryo-electron tomography and subtomogram averaging to determine the structures of Env trimers on aldrithiol-2 (AT-2)-inactivated virions in ligand-free, antibody-bound and CD4-bound forms at subnanometer resolution. Tomographic reconstructions document molecular features consistent with high-resolution structures of engineered soluble and detergent-solubilized Env trimers. One of three conformational states previously predicted by smFRET was not observed by cryo-ET, potentially owing to AT-2 inactivation. We did observe Env trimers to open in situ in response to CD4 binding, with an outward movement of gp120-variable loops and an extension of a critical gp41 helix. Overall features of Env trimer embedded in AT-2-treated virions appear well-represented by current engineered trimers.

An insertion unique to SARS-CoV-2 exhibits superantigenic character strengthened by recent mutations

Multisystem Inflammatory Syndrome in Children (MIS-C) associated with Coronavirus Disease 2019 (COVID-19) is a newly recognized condition in which children with recent SARS-CoV-2 infection present with a constellation of symptoms including hypotension, multiorgan involvement, and elevated inflammatory markers. These symptoms and the associated laboratory values strongly resemble toxic shock syndrome, an escalation of the cytotoxic adaptive immune response triggered upon the binding of pathogenic superantigens to MHCII molecules and T cell receptors (TCRs). Here, we used structure-based computational models to demonstrate that the SARS-CoV-2 spike (S) exhibits a high-affinity motif for binding TCR, interacting closely with both the α- and β-chains variable domains' complementarity-determining regions. The binding epitope on S harbors a sequence motif unique to SARS-CoV-2 (not present in any other SARS coronavirus), which is highly similar in both sequence and structure to bacterial superantigens. Further examination revealed that this interaction between the virus and human T cells is strengthened in the context of a recently reported rare mutation (D839Y/N/E) from a European strain of SARS-CoV-2. Furthermore, the interfacial region includes selected residues from a motif shared between the SARS viruses from the 2003 and 2019 pandemics, which has intracellular adhesion molecule (ICAM)-like character. These data suggest that the SARS-CoV-2 S may act as a superantigen to drive the development of MIS-C as well as cytokine storm in adult COVID-19 patients, with important implications for the development of therapeutic approaches.

Multifunctional 6-fluoro-3-[3-(pyrrolidin-1-yl)propyl]-1,2-benzoxazoles targeting behavioral and psychological symptoms of dementia (BPSD)

Patients suffering from dementia experience cognitive deficits and 90% of them show non-cognitive behavioral and psychological symptoms of dementia (BPSD). The spectrum of BPSD includes agitation, depression, anxiety and psychosis. Antipsychotics, e.g. quetiapine, have been commonly used off-label to control the burdensome symptoms, though they cause serious side effects and further cognitive impairment. Therefore, the development of targeted therapy for BPSD, suitable for elderly patients, remains relevant. A multitarget-directed ligand, acting on serotonin 5-HT_2A and dopamine D₂ receptors (R) and thus exerting anti-aggressive and antipsychotic activity, as well as on 5-HT₆Rs and 5-HT₇Rs (potential pro-cognitive, antidepressant and anxiolytic activity), poses a promising strategy for the treatment of BPSD. Antitargeting muscarinic M₃R and hERG channel is expected to reduce the risk of side effects. We obtained a series of stereoisomeric compounds by combining 6-fluoro-1,2-benzoxazole moiety and arylsulfonamide fragment through pyrrolidin-1-yl-propyl linker. N-[(3R)-1-[3-(6-fluoro-1,2-benzoxazol-3-yl)propyl]pyrrolidin-3-yl]-1-benzothiophene-2-sulfonamide showed a substantial affinity for the targets of interest (pK_i = 8.32-9.35) and no significant interaction with the antitargets. Functional studies revealed its antagonist efficacy (pK_B = 7.41-9.03). The lead compound showed a promising profile of antipsychotic-like activity in amphetamine- and MK-801-induced hyperlocomotion (MED = 2.5 mg/kg), antidepressant-like, as well as anxiolytic-like activity in mice (MED = 0.312 and 1.25 mg/kg in the forced swim and four-plate tests, respectively). Notably, the novel compound didn't affect spontaneous locomotor activity, nor induced catalepsy or memory deficits (step-through passive avoidance test) in therapeutically relevant doses, which proved its benign safety profile. The overall pharmacological characteristics of the lead compound outperformed the reference drug quetiapine, making it a promising option for evaluation in the treatment of BPSD.

Identification of a Novel Inhibitor of Catabolite Control Protein A from Staphylococcus aureus

Catabolite control protein A is a highly conserved transcriptional regulator in Gram-positive bacteria. Herein, we report a specific small-molecule inhibitor of Staphylococcus aureus catabolite control protein A (SaCcpA). The compound abrogates the regulatory function of SaCcpA, resulting in decreased expression of an S. aureus major cytotoxin, α-hemolysin. The observed synergism between the compound and antibiotics against S. aureus suggests its potential application in a combination therapy to combat antimicrobial resistance.

Homology Modeling of P2X Receptors

Since the X-ray structure of the zebra fish P2X4 receptor in the closed state was published in 2009 homology modeling has been used to generate structural models for P2X receptors. In this chapter, we outline how to use the MODELLER software to generate such structural models for P2X receptors whose structures have not been solved yet.

Conformational States of Exchange Protein Directly Activated by cAMP (EPAC1) Revealed by Ensemble Modeling and Integrative Structural Biology

Exchange proteins directly activated by cAMP (EPAC1 and EPAC2) are important allosteric regulators of cAMP-mediated signal transduction pathways. To understand the molecular mechanism of EPAC activation, we performed detailed Small-Angle X-ray Scattering (SAXS) analysis of EPAC1 in its apo (inactive), cAMP-bound, and effector (Rap1b)-bound states. Our study demonstrates that we can model the solution structures of EPAC1 in each state using ensemble analysis and homology models derived from the crystal structures of EPAC2. The N-terminal domain of EPAC1, which is not conserved between EPAC1 and EPAC2, appears folded and interacts specifically with another component of EPAC1 in each state. The apo-EPAC1 state is a dynamic mixture of a compact (Rg = 32.9 Å, 86%) and a more extended (Rg = 38.5 Å, 13%) conformation. The cAMP-bound form of EPAC1 in the absence of Rap1 forms a dimer in solution; but its molecular structure is still compatible with the active EPAC1 conformation of the ternary complex model with cAMP and Rap1. Herein, we show that SAXS can elucidate the conformational states of EPAC1 activation as it proceeds from the compact, inactive apo conformation through a previously unknown intermediate-state, to the extended cAMP-bound form, and then binds to its effector (Rap1b) in a ternary complex.

HER3-targeted protein chimera forms endosomolytic capsomeres and self-assembles into stealth nucleocapsids for systemic tumor homing of RNA interference in vivo

RNA interference represents a potent intervention for cancer treatment but requires a robust delivery agent for transporting gene-modulating molecules, such as small interfering RNAs (siRNAs). Although numerous molecular approaches for siRNA delivery are adequate in vitro, delivery to therapeutic targets in vivo is limited by payload integrity, cell targeting, efficient cell uptake, and membrane penetration. We constructed nonviral biomaterials to transport small nucleic acids to cell targets, including tumor cells, on the basis of the self-assembling and cell-penetrating activities of the adenovirus capsid penton base. Our recombinant penton base chimera contains polypeptide domains designed for noncovalent assembly with anionic molecules and tumor homing. Here, structural modeling, molecular dynamics simulations, and functional assays suggest that it forms pentameric units resembling viral capsomeres that assemble into larger capsid-like structures when combined with siRNA cargo. Pentamerization forms a barrel lined with charged residues mediating pH-responsive dissociation and exposing masked domains, providing insight on the endosomolytic mechanism. The therapeutic impact was examined on tumors expressing high levels of HER3/ErbB3 that are resistant to clinical inhibitors. Our findings suggest that our construct may utilize ligand mimicry to avoid host attack and target the siRNA to HER3+ tumors by forming multivalent capsid-like structures.

Modeling the Architecture of Depolymerase-Containing Receptor Binding Proteins in Klebsiella Phages

Klebsiella pneumoniae carries a thick polysaccharide capsule. This highly variable chemical structure plays an important role in its virulence. Many Klebsiella bacteriophages recognize this capsule with a receptor binding protein (RBP) that contains a depolymerase domain. This domain degrades the capsule to initiate phage infection. RBPs are highly specific and thus largely determine the host spectrum of the phage. A majority of known Klebsiella phages have only one or two RBPs, but phages with up to 11 RBPs with depolymerase activity and a broad host spectrum have been identified. A detailed bioinformatic analysis shows that similar RBP domains repeatedly occur in K. pneumoniae phages with structural RBP domains for attachment of an RBP to the phage tail (anchor domain) or for branching of RBPs (T4gp10-like domain). Structural domains determining the RBP architecture are located at the N-terminus, while the depolymerase is located in the center of protein. Occasionally, the RBP is complemented with an autocleavable chaperone domain at the distal end serving for folding and multimerization. The enzymatic domain is subjected to an intense horizontal transfer to rapidly shift the phage host spectrum without affecting the RBP architecture. These analyses allowed to model a set of conserved RBP architectures, indicating evolutionary linkages.

The Contribution of the Ankyrin Repeat Domain of TRPV1 as a Thermal Module

The TRPV1 cation nonselective ion channel plays an essential role in thermosensation and perception of other noxious stimuli. TRPV1 can be activated by low extracellular pH, high temperature, or naturally occurring pungent molecules such as allicin, capsaicin, or resiniferatoxin. Its noxious thermal sensitivity makes it an important participant as a thermal sensor in mammals. However, details of the mechanism of channel activation by increases in temperature remain unclear. Here, we used a combination of approaches to try to understand the role of the ankyrin repeat domain (ARD) in channel behavior. First, a computational modeling approach by coarse-grained molecular dynamics simulation of the whole TRPV1 embedded in a phosphatidylcholine and phosphatidylethanolamine membrane provides insight into the dynamics of this channel domain. Global analysis of the structural ensemble shows that the ARD is a region that sustains high fluctuations during dynamics at different temperatures. We then performed biochemical and thermal stability studies of the purified ARD by the means of circular dichroism and tryptophan fluorescence and demonstrate that this region undergoes structural changes at similar temperatures that lead to TRPV1 activation. Our data suggest that the ARD is a dynamic module and that it may participate in controlling the temperature sensitivity of TRPV1.

Stereospecificity control in aminoacyl-tRNA-synthetases: new evidence of d-amino acids activation and editing

The homochirality of amino acids is vital for the functioning of the translation apparatus. l-Amino acids predominate in proteins and d-amino acids usually represent diverse regulatory functional physiological roles in both pro- and eukaryotes. Aminoacyl-tRNA-synthetases (aaRSs) ensure activation of proteinogenic or nonproteinogenic amino acids and attach them to cognate or noncognate tRNAs. Although many editing mechanisms by aaRSs have been described, data about the protective role of aaRSs in d-amino acids incorporation remained unknown. Tyrosyl- and alanyl-tRNA-synthetases were represented as distinct members of this enzyme family. To study the potential to bind and edit noncognate substrates, Thermus thermophilus alanyl-tRNA-synthetase (AlaRS) and tyrosyl-tRNA-synthetase were investigated in the context of d-amino acids recognition. Here, we showed that d-alanine was effectively activated by AlaRS and d-Ala-tRNAAla, formed during the erroneous aminoacylation, was edited by AlaRS. On the other hand, it turned out that d-aminoacyl-tRNA-deacylase (DTD), which usually hydrolyzes d-aminoacyl-tRNAs, was inactive against d-Ala-tRNAAla. To support the finding about DTD, computational docking and molecular dynamics simulations were run. Overall, our work illustrates the novel function of the AlaRS editing domain in stereospecificity control during translation together with trans-editing factor DTD. Thus, we propose different evolutionary strategies for the maintenance of chiral selectivity during translation.

SWISS-MODEL: homology modelling of protein structures and complexes

Homology modelling has matured into an important technique in structural biology, significantly contributing to narrowing the gap between known protein sequences and experimentally determined structures. Fully automated workflows and servers simplify and streamline the homology modelling process, also allowing users without a specific computational expertise to generate reliable protein models and have easy access to modelling results, their visualization and interpretation. Here, we present an update to the SWISS-MODEL server, which pioneered the field of automated modelling 25 years ago and been continuously further developed. Recently, its functionality has been extended to the modelling of homo- and heteromeric complexes. Starting from the amino acid sequences of the interacting proteins, both the stoichiometry and the overall structure of the complex are inferred by homology modelling. Other major improvements include the implementation of a new modelling engine, ProMod3 and the introduction a new local model quality estimation method, QMEANDisCo. SWISS-MODEL is freely available at https://swissmodel.expasy.org.

Substrate-assisted mechanism of catalytic hydrolysis of misaminoacylated tRNA required for protein synthesis fidelity

d-aminoacyl-tRNA-deacylase (DTD) prevents the incorporation of d-amino acids into proteins during translation by hydrolyzing the ester bond between mistakenly attached amino acids and tRNAs. Despite extensive study of this proofreading enzyme, the precise catalytic mechanism remains unknown. Here, a combination of biochemical and computational investigations has enabled the discovery of a new substrate-assisted mechanism of d-Tyr-tRNA^Tyr hydrolysis by Thermus thermophilus DTD. Several functional elements of the substrate, misacylated tRNA, participate in the catalysis. During the hydrolytic reaction, the 2'-OH group of the А76 residue of d-Tyr-tRNA^Tyr forms a hydrogen bond with a carbonyl group of the tyrosine residue, stabilizing the transition-state intermediate. Two water molecules participate in this reaction, attacking and assisting ones, resulting in a significant decrease in the activation energy of the rate-limiting step. The amino group of the d-Tyr aminoacyl moiety is unprotonated and serves as a general base, abstracting the proton from the assisting water molecule and forming a more nucleophilic ester-attacking species. Quantum chemical methodology was used to investigate the mechanism of hydrolysis. The DFT-calculated deacylation reaction is in full agreement with the experimental data. The Gibbs activation energies for the first and second steps were 10.52 and 1.05 kcal/mol, respectively, highlighting that the first step of the hydrolysis process is the rate-limiting step. Several amino acid residues of the enzyme participate in the coordination of the substrate and water molecules. Thus, the present work provides new insights into the proofreading details of misacylated tRNAs and can be extended to other systems important for translation fidelity.

Identification and functional characterization, including cytokine production modulation, of the novel chicken Interleukin-11

Interleukin (IL)-11 plays an important role in the immune system. However, IL-11 has not yet been characterized in avian species, including chickens. This study is the first to clone and functionally characterize chicken IL-11 (chIL-11). Multiple alignments and phylogenetic tree comparisons of chIL-11 with IL-11 proteins from other species revealed high levels of conservation and a close relationship between chicken and Japanese quail IL-11. Our results demonstrate that chIL-11 was a functional ligand of IL-11RA and IL-6ST in chicken HD11 and OU2 cell lines, as well as activated and regulated JAK-STAT, NF-κB, PI3K/AKT, and MAPK signaling pathways in chicken cell lines. In addition, chIL-11 inhibited nitric oxide production, affected proliferation of both tested cell lines, inhibited Type 1 and 17 T helper (Th) cytokine and IL-26, IL-12, and IL-17A-induced interferon-γ production, and enhanced Th2 cytokine (IL-4 and IL-10) production. Taken together, functional analysis of chIL-11 revealed it bound to IL-11RA and IL-6ST and activated the JAK-STAT, NF-κB, and MAPK signaling pathways, which resulted in modulation of Th1/Th17 and Th2 cytokine production in chicken HD11 and OU2 cell lines. Overall, this indicates chIL-11 has a role in both the innate and adaptive immune system.

SWISS-MODEL: homology modelling of protein structures and complexes

Homology modelling has matured into an important technique in structural biology, significantly contributing to narrowing the gap between known protein sequences and experimentally determined structures. Fully automated workflows and servers simplify and streamline the homology modelling process, also allowing users without a specific computational expertise to generate reliable protein models and have easy access to modelling results, their visualization and interpretation. Here, we present an update to the SWISS-MODEL server, which pioneered the field of automated modelling 25 years ago and been continuously further developed. Recently, its functionality has been extended to the modelling of homo- and heteromeric complexes. Starting from the amino acid sequences of the interacting proteins, both the stoichiometry and the overall structure of the complex are inferred by homology modelling. Other major improvements include the implementation of a new modelling engine, ProMod3 and the introduction a new local model quality estimation method, QMEANDisCo. SWISS-MODEL is freely available at https://swissmodel.expasy.org.

Homology threading to generate RNA polymerase structures

Homology threading is a powerful technology for generating structural models based on homologous structures. Here we use threading to generate four complex RNA polymerase models. The models appear to be as useful as x-ray crystal structures or cryo-electron microscopy structures to support research projects.

Anticipating innovations in structural biology

In this review, we describe how the interplay among science, technology and community interests contributed to the evolution of four structural biology data resources. We present the method by which data deposited by scientists are prepared for worldwide distribution, and argue that data archiving in a trusted repository must be an integral part of any scientific investigation.

Sequence-dependent dynamical instability of the human prion protein: a comparative simulation study

The present study aimed to explore the most probable regions of the human prion protein backbone for which the initial steps of conformational transitions as a result of intrinsic and extrinsic perturbing factors on the protein structure can be assigned. A total of 0.3-μs molecular dynamics simulations on several analog structures of the protein have been performed. To mimic the impact of the extrinsic and intrinsic destructive parameters on the dynamical characteristics of the protein, mild acidic conditions and R208H mutation have been simulated. The findings indicated that distribution of conformational flexibilities along the protein chain was almost independent of the induced perturbing factors, and was mostly centralized on certain distinct parts of the structure comprising residues 132-145 and 187-203. Analyses also revealed that the segment comprising residues 187-203 may be considered as a peptide sequence, possessing high potential to start the initial steps of conformational rearrangements due to the induced physicochemical alterations. Sequence alignment and molecular dynamics data also revealed that segment 178-203 prefers to accommodate in extended structures rather than α-helices. Region 178-203 may be considered as a peptide switch capable of initiating the conformational transitions due to the introduced modifications and perturbing parameters.

Prediction on the risk population of idiosyncratic adverse reactions based on molecular docking with mutant proteins

Idiosyncratic adverse drug reactions are drug reactions that occur rarely and unpredictably among the population. These reactions often occur after a drug is marketed, which means that they are strongly related to the genotype of the population. The prediction of such adverse reactions is a major challenge because of the lack of appropriate test models during the drug development process. In this study, we chose withdrawn drugs because the reasons why they were withdrawn and from which countries or regions is easily obtained. We selected Dilevalol and its chiral drug (Labetalol) as the investigatory drugs, as they have been withdrawn from a European market (Britain) because of serious hepatotoxicity. First, we searched for and obtained the Dilevalol-induced- liver-injury related protein, multidrug resistance protein 1 (MDR1), from the Comparative Toxicogenomics Database (CTD). Then, we searched and extracted 477 non-synonymous single nucleotide polymorphisms (nsSNP) on MDR1 in the dbSNP database. Second, we used the VarMod tool to predict the functional changes of MDR1 induced by these nsSNPs, from which we extracted the nsSNPs that significantly change the functions of this protein. Third, we built the three-dimensional structures of those variant proteins and used AutoDock to perform a docking study, choosing the best model to determine the sites of nsSNPs. Finally, we used the data from the 1000 Genomes Project to verify the dominant population distribution of the risk SNP. We applied the same strategy to the post-marketing drug-induced liver injury drugs to further test the feasibility of our method.

Discovery of a proteolytic flagellin family in diverse bacterial phyla that assembles enzymatically active flagella

Bacterial flagella are cell locomotion and occasional adhesion organelles composed primarily of the polymeric protein flagellin, but to date have not been associated with any enzymatic function. Here, we report the bioinformatics-driven discovery of a class of enzymatic flagellins that assemble to form proteolytically active flagella. Originating by a metallopeptidase insertion into the central flagellin hypervariable region, this flagellin family has expanded to at least 74 bacterial species. In the pathogen, Clostridium haemolyticum, metallopeptidase-containing flagellin (which we termed flagellinolysin) is the second most abundant protein in the flagella and is localized to the extracellular flagellar surface. Purified flagellar filaments and recombinant flagellin exhibit proteolytic activity, cleaving nearly 1000 different peptides. With ~ 20,000 flagellin copies per  ~ 10-μm flagella this assembles the largest proteolytic complex known. Flagellum-mediated extracellular proteolysis expands our understanding of the functional plasticity of bacterial flagella, revealing this family as enzymatic biopolymers that mediate interactions with diverse peptide substrates.

So far no enzymatic activity has been attributed to flagellin, the major component of bacterial flagella. Here the authors use bioinformatic analysis and identify a metallopeptidase insertion in flagellins from 74 bacterial species and show that recombinant flagellin and flagellar filaments have proteolytic activity.