Stoichiometry and architecture of the human pyruvate dehydrogenase complex

The pyruvate dehydrogenase complex (PDHc) is a key megaenzyme linking glycolysis with the citric acid cycle. In mammalian PDHc, dihydrolipoamide acetyltransferase (E2) and the dihydrolipoamide dehydrogenase-binding protein (E3BP) form a 60-subunit core that associates with the peripheral subunits pyruvate dehydrogenase (E1) and dihydrolipoamide dehydrogenase (E3). The structure and stoichiometry of the fully assembled, mammalian PDHc or its core remained elusive. Here, we demonstrate that the human PDHc core is formed by 48 E2 copies that bind 48 E1 heterotetramers and 12 E3BP copies that bind 12 E3 homodimers. Cryo-electron microscopy, together with native and cross-linking mass spectrometry, confirmed a core model in which 8 E2 homotrimers and 12 E2-E2-E3BP heterotrimers assemble into a pseudoicosahedral particle such that the 12 E3BP molecules form six E3BP-E3BP intertrimer interfaces distributed tetrahedrally within the 60-subunit core. The even distribution of E3 subunits in the peripheral shell of PDHc guarantees maximum enzymatic activity of the megaenzyme.

Study of individual domains contributing to MALT1 dimerization in BCL10-independent and dependent assembly

The CARMA-BCL10-MALT1 (CBM) signalosome functions as a pivotal supramolecular module, integrating diverse receptor-induced signaling pathways to regulate BCL10-dependent NF-kB activation in innate and adaptive immunity. Conversely, the API2-MALT1 fusion protein in t(11; 18)(q21; q21) MALT lymphoma constitutively induces BCL10-independent NF-kB activation. MALT1 dimer formation is indispensable for the requisite proteolytic activity and is critical for NF-kB activation regulation in both scenarios. However, the molecular assembly of MALT1 individual domains in CBM activation remains elusive. Here we report the crystal structure of the MALT1 death domain (DD) at a resolution of 2.1 Å, incorporating reconstructed residues in previously disordered loops 1 and 2. Additionally, we observe a conformational regulation element (CRE) regulating stem-helix formation in NLRPs pyrin (PYD) within the MALT1 DD structure. The structure reveals a stem-helix-mediated dimer further corroborated in solution. To elucidate how the BCL10 filament facilitates MALT1 dimerization, we reconstitute a BCL10-CARD-MALT1-DD-IG1-IG2 complex model. We propose a N+7 rule for BCL10-dependent MALT1 dimerization via the IG1-IG2 domain and for MALT1-dependent cleavage in trans. Biochemical data further indicates concentration-dependent dimerization of the MALT1 IG1-IG2 domain, facilitating MALT1 dimerization in BCL10-independent manner. Our findings provide a structural and biochemical foundation for understanding MALT1 dimeric mechanisms, shedding light on potential BCL10-independent MALT1 dimer formation and high-order BCL10-MALT1 assembly.

Bioengineering a plant NLR immune receptor with a robust binding interface toward a conserved fungal pathogen effector

Bioengineering of plant immune receptors has emerged as a key strategy for generating novel disease resistance traits to counteract the expanding threat of plant pathogens to global food security. However, current approaches are limited by rapid evolution of plant pathogens in the field and may lack durability when deployed. Here, we show that the rice nucleotide-binding, leucine-rich repeat (NLR) immune receptor Pik-1 can be engineered to respond to a conserved family of effectors from the multihost blast fungus pathogen Magnaporthe oryzae. We switched the effector binding and response profile of the Pik NLR from its cognate rice blast effector AVR-Pik to the host-determining factor pathogenicity toward weeping lovegrass 2 (Pwl2) by installing a putative host target, OsHIPP43, in place of the native integrated heavy metal-associated domain (generating Pikm-1OsHIPP43). This chimeric receptor also responded to other PWL alleles from diverse blast isolates. The crystal structure of the Pwl2/OsHIPP43 complex revealed a multifaceted, robust interface that cannot be easily disrupted by mutagenesis, and may therefore provide durable, broad resistance to blast isolates carrying PWL effectors in the field. Our findings highlight how the host targets of pathogen effectors can be used to bioengineer recognition specificities that have more robust properties compared to naturally evolved disease resistance genes.

Enhanced stabilisation and reduced fibril forming potential of an amyloidogenic light chain using a variable heavy domain to mimic the homodimer complex

Light chain amyloidosis (AL), is classified as a plasma cell dyscrasia, whereby a mutant plasma cell multiplies uncontrollably and secretes enormous amounts of immunoglobulin-free light chain (FLC) fragments. These FLCs undergo a process of misfolding and aggregation into amyloid fibrils, that can cause irreversible system-wide damage. Current treatments that focus on depleting the underlying plasma cell clone are often poorly tolerated, particularly in patients with severe cardiac involvement, meaning patient prognosis is poor. An alternative treatment approach currently being explored is the inhibition of FLC aggregation by stabilisation of the native conformer. Here, we aimed to identify and characterise antibody fragments that target FLC domains and promote their stabilisation. Using phage-display screening methods, we identified a variable heavy (VH) domain, termed VH1, targeted towards the FLC. Using differential scanning fluorimetry and surface plasmon resonance, VH1 was characterised to bind and kinetically stabilise an amyloidogenic FLC, whereby a > 5.5 °C increase in thermal stability was noted. This improved stability corresponded to the inhibition of fibril formation, where 10 : 1 LC : VH1 concentration reduced aggregation to baseline levels. X-ray crystallographic structures of the LC : VH1 complex at atomic resolution revealed binding in a 1 : 1 ratio, mimicking the dimeric antigen binding sites of the native immunoglobulin molecule and the native LC homodimer.

The oligomeric states of dye-decolorizing peroxidases from Streptomyces lividans and their implications for mechanism of substrate oxidation

A common evolutionary mechanism in biology to drive function is protein oligomerization. In prokaryotes, the symmetrical assembly of repeating protein units to form homomers is widespread, yet consideration in vitro of whether such assemblies have functional or mechanistic consequences is often overlooked. Dye-decolorizing peroxidases (DyPs) are one such example, where their dimeric α + β barrel units can form various oligomeric states, but the oligomer influence, if any, on mechanism and function has received little attention. In this work, we have explored the oligomeric state of three DyPs found in Streptomyces lividans, each with very different mechanistic behaviors in their reactions with hydrogen peroxide and organic substrates. Using analytical ultracentrifugation, we reveal that except for one of the A-type DyPs where only a single sedimenting species is detected, oligomer states ranging from homodimers to dodecamers are prevalent in solution. Using cryo-EM on preparations of the B-type DyP, we determined a 3.02 Å resolution structure of a hexamer assembly that corresponds to the dominant oligomeric state in solution as determined by analytical ultracentrifugation. Furthermore, cryo-EM data detected sub-populations of higher-order oligomers, with one of these formed by an arrangement of two B-type DyP hexamers to give a dodecamer assembly. Our solution and structural insights of these oligomer states provide a new framework to consider previous mechanistic studies of these DyP members and are discussed in terms of long-range electron transfer for substrate oxidation and in the "storage" of oxidizable equivalents on the heme until a two-electron donor is available.

In silico prediction of CD8+ and CD4+ T cell epitopes in Leishmania major proteome: Using immunoinformatics

The leishmaniases are NDTs (neglected tropical diseases) that affect people all over the world. They are brought on by protozoans from the genus Leishmania and disseminated by phlebotomine flies that are afflicted with the disease. The best option to manage and lower the incidence of these diseases has been thought by the creation of a safe and effective vaccination. This research used an in silico based mining approach to look for high potential epitopes that might bind to MHC Class I and MHC Class II molecules (mainly; HLA-A*02:01 & HLA-DRB1*03:01) from human population in order to promote vaccine development. Based on the presence of signal peptides, GPI anchors, antigenicity predictions, and a subtractive proteomic technique, we have screened 17 putative antigenic proteins from the 8083 total proteins of L. major. After that thorough immunogenic epitope prediction were done using IEDB-AR tools. We isolated five immunogenic epitopes (three 9-mer & two 15-mer) from five antigenic proteins through docking and MD simulation analysis. Finally, these five anticipated epitopes, viz., TLPEIPVNV, ELMAPVFGL, TLAAAVALL, NSINIRLDGVTSAGF and NVPLVVDASSLFRVA have considerably stronger binding potential with their respective alleles and may trigger immunological responses. The goal of this work was to identify MHC restricted epitopes for CD8+ and CD4+ T cells activation using immunoinformatics in order to identify potential vaccine candidates against L. major parasites.

Structure of a novel form of phosphopantetheine adenylyltransferase from Klebsiella pneumoniae at 2.59 Å resolution

Phosphopantetheine adenylyltransferase (EC. 2.7.7.3, PPAT) catalyzes the penultimate step of the multistep reaction in the coenzyme A (CoA) biosynthesis pathway. In this step, an adenylyl group from adenosine triphosphate (ATP) is transferred to 4'-phosphopantetheine (PNS) yielding 3'-dephospho-coenzyme A (dpCoA) and pyrophosphate (PPi). PPAT from strain C3 of Klebsiella pneumoniae (KpPPAT) was cloned, expressed and purified. It was crystallized using 0.1 M HEPES buffer and PEG10000 at pH 7.5. The crystals belonged to tetragonal space group P41212 with cell dimensions of a = b = 72.82 Å and c = 200.37 Å. The structure was determined using the molecular replacement method and refined to values of 0.208 and 0.255 for Rcryst and Rfree factors, respectively. The structure determination showed the presence of three crystallographically independent molecules A, B and C in the asymmetric unit. The molecules A and B are observed in the form of a dimer in the asymmetric unit while molecule C belongs to the second dimer whose partner is related by crystallographic twofold symmetry. The polypeptide chain of KpPPAT folds into a β/α structure. The conformations of the side chains of several residues in the substrate binding site in KpPPAT are significantly different from those reported in other PPATs. As a result, the modes of binding of substrates, phosphopantetheine (PNS) and adenosine triphosphate (ATP) differ considerably. The binding studies using fluorescence spectroscopy indicated a KD value of 3.45 × 10-4 M for ATP which is significantly lower than the corresponding values reported for PPAT from other species.

Activation of Thoeris antiviral system via SIR2 effector filament assembly

To survive bacteriophage (phage) infections, bacteria developed numerous anti-phage defence systems1-7. Some of them (for example, type III CRISPR-Cas, CBASS, Pycsar and Thoeris) consist of two modules: a sensor responsible for infection recognition and an effector that stops viral replication by destroying key cellular components8-12. In the Thoeris system, a Toll/interleukin-1 receptor (TIR)-domain protein, ThsB, acts as a sensor that synthesizes an isomer of cyclic ADP ribose, 1''-3' glycocyclic ADP ribose (gcADPR), which is bound in the Smf/DprA-LOG (SLOG) domain of the ThsA effector and activates the silent information regulator 2 (SIR2)-domain-mediated hydrolysis of a key cell metabolite, NAD+ (refs. 12-14). Although the structure of ThsA has been solved15, the ThsA activation mechanism remained incompletely understood. Here we show that 1''-3' gcADPR, synthesized in vitro by the dimeric ThsB' protein, binds to the ThsA SLOG domain, thereby activating ThsA by triggering helical filament assembly of ThsA tetramers. The cryogenic electron microscopy (cryo-EM) structure of activated ThsA revealed that filament assembly stabilizes the active conformation of the ThsA SIR2 domain, enabling rapid NAD+ depletion. Furthermore, we demonstrate that filament formation enables a switch-like response of ThsA to the 1''-3' gcADPR signal.

In silico and structural investigation of sulfonamides targeting VraSR two component system in methicillin-resistant Staphylococcus aureus

Drug-resistant Staphylococcus aureus strains are global health concerns. Several studies have shown that these strains can develop defences against cell wall antibiotics such as β-lactams, glycopeptides and daptomycin which target cell wall biosynthesis. The coordination of these responses have been associated with two component system (TCS) regulated by histidine kinase protein (VraS) and its cognate regulator VraR which influences the target DNA upon signal recognition. Computer-based screening methods, predictions and simulations have emerged as more efficient and quick ways to identify promising new compound leads from large databases against emerging drug targets thus allowing prediction of small select set of molecules for further validations. These combined approaches conserve valuable time and resources. Due to methicillin resistance, sulfonamide-derivative medications have been found to be effective treatment strategy to treat S. aureus infections. The current study used ligand-based virtual screening (LBVS) to identify powerful sulfonamide derivative inhibitors from an antibacterial compound library against VraSR signaling components, VraS and VraR. We identified promising sulfonamide derivative [compound 5: (4-[(1-{[(3,5-Dimethoxyphenyl)Carbamoyl]Methyl}-2,4-Dioxo-1,2,3,4-Tetrahydroquinazolin-3-Yl)Methyl]-N-[(Furan-2-Yl)Methyl]Benzamide)] with reasonable binding parameters of -31.38 kJ/mol and ΔGbind score of -294.32 kJ/mol against ATP binding domain of sensor kinase VraS. We further identified four compounds N1 (PCID83276726), N3 (PCID83276757), N9 (PCID3672584), and N10 (PCID20900589) against VraR DNA binding domain (VraRC) with ΔGbind energies of -190.27, -237.54, -165.21, and -190.88 kJ/mol, respectively. Structural and simulation analyses further suggest their stable interactions with DNA interacting residues and potential to disrupt DNA binding domain dimerization; therefore, it is prudent to further investigate and characterize them as VraR dimer disruptors and inhibit other promoter binding site. Interestingly, the discovery of drugs that target VraS and VraR may open new therapeutic avenues for drug-resistant S. aureus. These predictions based on screening, simulations and binding affinities against VraSR components hold promise for opening novel therapeutic avenues against drug-resistant S. aureus and present opportunities for repositioning efforts. These efforts aim to create analogs with enhanced potency and selectivity against two-component signaling systems that significantly contribute to virulence in MRSA or VRSA. These analyses contribute valuable insights into potential avenues for combating antibiotic-resistant S. aureus through computationally driven drug discovery strategies.Communicated by Ramaswamy H. Sarma.

Reconciling ASPP-p53 binding mode discrepancies through an ensemble binding framework that bridges crystallography and NMR data

ASPP2 and iASPP bind to p53 through their conserved ANK-SH3 domains to respectively promote and inhibit p53-dependent cell apoptosis. While crystallography has indicated that these two proteins employ distinct surfaces of their ANK-SH3 domains to bind to p53, solution NMR data has suggested similar surfaces. In this study, we employed multi-scale molecular dynamics (MD) simulations combined with free energy calculations to reconcile the discrepancy in the binding modes. We demonstrated that the binding mode based solely on a single crystal structure does not enable iASPP's RT loop to engage with p53's C-terminal linker-a verified interaction. Instead, an ensemble of simulated iASPP-p53 complexes facilitates this interaction. We showed that the ensemble-average inter-protein contacting residues and NMR-detected interfacial residues qualitatively overlap on ASPP proteins, and the ensemble-average binding free energies better match experimental KD values compared to single crystallgarphy-determined binding mode. For iASPP, the sampled ensemble complexes can be grouped into two classes, resembling the binding modes determined by crystallography and solution NMR. We thus propose that crystal packing shifts the equilibrium of binding modes towards the crystallography-determined one. Lastly, we showed that the ensemble binding complexes are sensitive to p53's intrinsically disordered regions (IDRs), attesting to experimental observations that these IDRs contribute to biological functions. Our results provide a dynamic and ensemble perspective for scrutinizing these important cancer-related protein-protein interactions (PPIs).

Architecture and potential roles of a delta-class glutathione S-transferase in protecting honey bee from agrochemicals

The European honey bee, Apis mellifera, serves as the principle managed pollinator species globally. In recent decades, honey bee populations have been facing serious health threats from combined biotic and abiotic stressors, including diseases, limited nutrition, and agrochemical exposure. Understanding the molecular mechanisms underlying xenobiotic adaptation of A. mellifera is critical, considering its extensive exposure to phytochemicals and agrochemicals present in the environment. In this study, we conducted a comprehensive structural and functional characterization of AmGSTD1, a delta class glutathione S-transferase (GST), to unravel its roles in agrochemical detoxification and antioxidative stress responses. We determined the 3-dimensional (3D) structure of a honey bee GST using protein crystallography for the first time, providing new insights into its molecular structure. Our investigations revealed that AmGSTD1 metabolizes model substrates, including 1-chloro-2,4-dinitrobenzene (CDNB), p-nitrophenyl acetate (PNA), phenylethyl isothiocyanate (PEITC), propyl isothiocyanate (PITC), and the oxidation byproduct 4-hydroxynonenal (HNE). Moreover, we discovered that AmGSTD1 exhibits binding affinity with the fluorophore 8-Anilinonaphthalene-1-sulfonic acid (ANS), which can be inhibited with various herbicides, fungicides, insecticides, and their metabolites. These findings highlight the potential contribution of AmGSTD1 in safeguarding honey bee health against various agrochemicals, while also mitigating oxidative stress resulting from exposure to these substances.

Structure-based insights into the mechanism of [4Fe-4S]-dependent sulfur insertase LarE

Several essential cellular metabolites, such as enzyme cofactors, contain sulfur atoms and their biosynthesis requires specific thiolation enzymes. LarE is an ATP-dependent sulfur insertase, which catalyzes the sequential conversion of the two carboxylate groups of the precursor of the lactate racemase cofactor into thiocarboxylates. Two types of LarE enzymes are known, one that uses a catalytic cysteine as a sacrificial sulfur donor, and the other one that uses a [4Fe-4S] cluster as a cofactor. Only the crystal structure of LarE from Lactobacillus plantarum (LpLarE) from the first class has been solved. We report here the crystal structure of LarE from Methanococcus maripaludis (MmLarE), belonging to the second class, in the cluster-free (apo-) and cluster-bound (holo-) forms. The structure of holo-MmLarE shows that the [4Fe-4S] cluster is chelated by three cysteines only, leaving an open coordination site on one Fe atom. Moreover, the fourth nonprotein-bonded iron atom was able to bind an anionic ligand such as a phosphate group or a chloride ion. Together with the spectroscopic analysis of holo-MmLarE and the previously reported biochemical investigations of holo-LarE from Thermotoga maritima, these crystal structures support the hypothesis of a reaction mechanism, in which the [4Fe-4S] cluster binds a hydrogenosulfide ligand in place of the chloride anion, thus generating a [4Fe-5S] intermediate, and transfers it to the substrate, as in the case of [4Fe-4S]-dependent tRNA thiolation enzymes.

Assembly Requirements for the Construction of Large-Scale Binary Protein Structures

The precise assembly of multiple biomacromolecules into well-defined structures and materials is of great importance for various biomedical and nanobiotechnological applications. In this study, we investigate the assembly requirements for two-component materials using charged protein nanocages as building blocks. To achieve this, we designed several variants of ferritin nanocages to determine the surface characteristics necessary for the formation of large-scale binary three-dimensional (3D) assemblies. These nanocage variants were employed in protein crystallization experiments and macromolecular crystallography analyses, complemented by computational methods. Through the screening of nanocage variant combinations at various ionic strengths, we identified three essential features for successful assembly: (1) the presence of a favored crystal contact region, (2) the presence of a charged patch not involved in crystal contacts, and (3) sufficient distinctiveness between the nanocages. Surprisingly, the absence of noncrystal contact mediating patches had a detrimental effect on the assemblies, highlighting their unexpected importance. Intriguingly, we observed the formation of not only binary structures but also both negatively and positively charged unitary structures under previously exclusively binary conditions. Overall, our findings will inform future design strategies by providing some design rules, showcasing the utility of supercharging symmetric building blocks in facilitating the assembly of biomacromolecules into large-scale binary 3D assemblies.

Structural analysis and functional evaluation of the disordered ß-hexosyltransferase region from Hamamotoa (Sporobolomyces) singularis

Hamamotoa (Sporobolomyces) singularis codes for an industrially important membrane bound ß-hexosyltransferase (BHT), (BglA, UniprotKB: Q564N5) that has applications in the production of natural fibers such as galacto-oligosaccharides (GOS) and natural sugars found in human milk. When heterologously expressed by Komagataella phaffii GS115, BHT is found both membrane bound and soluble secreted into the culture medium. In silico structural predictions and crystal structures support a glycosylated homodimeric enzyme and the presence of an intrinsically disordered region (IDR) with membrane binding potential within its novel N-terminal region (1-110 amino acids). Additional in silico analysis showed that the IDR may not be essential for stable homodimerization. Thus, we performed progressive deletion analyses targeting segments within the suspected disordered region, to determine the N-terminal disorder region's impact on the ratio of membrane-bound to secreted soluble enzyme and its contribution to enzyme activity. The ratio of the soluble secreted to membrane-bound enzyme shifted from 40% to 53% after the disordered N-terminal region was completely removed, while the specific activity was unaffected. Furthermore, functional analysis of each glycosylation site found within the C-terminal domain revealed reduced total secreted protein activity by 58%-97% in both the presence and absence of the IDR, indicating that glycosylation at all four locations is required by the host for the secretion of active enzyme and independent of the removed disordered N-terminal region. Overall, the data provides evidence that the disordered region only partially influences the secretion and membrane localization of BHT.

Advanced interactional characterization of the inhibitory effect of anthocyanin extract from Hibiscus sabdariffa L. on apple polyphenol oxidase

In this study, a comprehensive approach to advance the inhibitory effect of Hibiscus sabdariffa extract on apple polyphenol oxidase (PPO) was performed. PPO was extracted, purified, and characterized for optimal activity, whereas response surface methodology generated a quadratic polynomial model to fit the experimental results for hibiscus extraction. The optimum conditions allowed to predict a maximum recovery of anthocyanins (256.11 mg delphinidin-3-O-glucoside/g), with a validated value of 272.87 mg delphinidin-3-O-glucoside/g dry weight (DW). The chromatographic methods highlighted the presence of gallic acid (36,812.90 µg/g DW extract), myricetin (141,933.84 µg/g DW extract), caffeic acid (101,394.07 µg/g DW extract), sinapic acid (1157.46 µg/g DW extract), kaempferol (2136.76 µg/g DW extract), and delphinidin 3-O-β-d-glucoside (226,367.08 µg/g DW extract). The inactivation of PPO followed a first-order kinetic model. A temperature-mediated flexible fit between PPO and anthocyanins was suggested, whereas the molecular docking tests indicated that PPO is a good receptor for cafestol, gallic acid, and catechin, involving hydrophobic and hydrogen bond interactions. PRACTICAL APPLICATION: It is well known that enzymatic browning is one of the most important challenges in the industrial minimal processing of selected fruit and vegetable products. Novel inhibitors for polyphenol oxidase are proposed in this study by using an anthocyanin-enriched extract from Hibiscus sabdariffa L. Based on our results, combining the chemical effect of phytochemicals from hibiscus extract with different functional groups with minimal heating could be an interesting approach for the development of a new strategy to inhibit apple polyphenol oxidase.

Peptide aptamer-based time-resolved fluoroimmunoassay for CHIKV diagnosis

BACKGROUND

Chikungunya virus (CHIKV) and Dengue virus (DENV) have similar clinical symptoms, which often induce misdiagnoses. Therefore, an antigen detection diagnostic system that can clearly identify these two viruses is desirable.

METHODS

In this study, we developed a novel peptide with high affinity and specificity to CHIKV, and further constructed peptide aptamer-based TRFIA assay to efficiently detect CHIKV. Peptide aptamer B2 (ITPQSSTTEAEL) and B3 (DTQGSNWI) were obtained through computer-aided design and selected as CHIKV-specific peptide aptamers based on their high binding affinity, strong hydrogen bonding, and RMSD of molecular docking. Then, a sandwich-Time-Resolved Fluoroimmunoassay (TRFIA) was successfully constructed for the detection of the interaction between peptide aptamers and viruses.

RESULTS

When using B2 as the detection element, highly specific detection of CHIKV E2 was achieved with detection limits of 8.5 ng/ml in PBS solution. Variation coefficient between inter-assay showed the disturbances received from the detection of clinical fluid specimens (including serum and urine), were also within acceptable limits. The detection limits for 10-fold dilution serum and urine were 57.8 ng/mL and 147.3 ng/mL, respectively. The fluorescent signal intensity exhibited a good linear correlation with E2 protein concentration in the range of 0-1000 ng/mL, indicating the potential for quantitative detection of E2 protein.

CONCLUSIONS

These results demonstrate that the construction of peptide aptamers with high affinity and specificity provides an excellent method for rapid diagnostic element screening, and the developed peptide aptamer B2 contributed to better detection of CHIKV viral particles compared to traditional antibodies.

A critical evaluation of protein kinase regulation by activation loop autophosphorylation

Phosphorylation of proteins is a ubiquitous mechanism of regulating their function, localization, or activity. Protein kinases, enzymes that use ATP to phosphorylate protein substrates are, therefore, powerful signal transducers in eukaryotic cells. The mechanism of phosphoryl-transfer is universally conserved among protein kinases, which necessitates the tight regulation of kinase activity for the orchestration of cellular processes with high spatial and temporal fidelity. In response to a stimulus, many kinases enhance their own activity by autophosphorylating a conserved amino acid in their activation loop, but precisely how this reaction is performed is controversial. Classically, kinases that autophosphorylate their activation loop are thought to perform the reaction in trans, mediated by transient dimerization of their kinase domains. However, motivated by the recently discovered regulation mechanism of activation loop cis-autophosphorylation by a kinase that is autoinhibited in trans, we here review the various mechanisms of autoregulation that have been proposed. We provide a framework for critically evaluating biochemical, kinetic, and structural evidence for protein kinase dimerization and autophosphorylation, and share some thoughts on the implications of these mechanisms within physiological signaling networks.

The Modulation of Sucrose Nonfermenting 1-Related Protein Kinase 2.6 State by Persulfidation and Phosphorylation: Insights from Molecular Dynamics Simulations

SnRK2.6 (SUCROSE NONFERMENTING 1-RELATED PROTEIN KINASE2.6) has been characterized as a molecular switch for the intracellular abscisic acid (ABA) signal-transduction pathway. Normally, SnRK2.6 is kept in an "off" state, forming a binary complex with protein phosphatase type 2Cs (PP2Cs). Upon stressful conditions, SnRK2.6 turns into an "on" state by its release from PP2Cs and then phosphorylation at Ser175. However, how the "on" and "off" states for SnRK2.6 are fine-tuned, thereby controlling the initiation and braking processes of ABA signaling, is still largely unclear. SnRK2.6 activity was tightly regulated through protein post-translational modifications (PTM), such as persulfidation and phosphorylation. Taking advantage of molecular dynamics simulations, our results showed that Cys131/137 persulfidation on SnRK2.6 induces destabilized binding and weakened interactions between SnRK2.6 and HAB1 (HYPERSENSITIVE TO ABA1), an important PP2C family protein. This unfavorable effect on the association of the SnRK2.6-HAB1 complex suggests that persulfidation functions are a positive regulator of ABA signaling initiation. In addition, Ser267 phosphorylation in persulfidated SnRK2.6 renders a stable physical association between SnRK2.6 and HAB1, a key characterization for SnRK2.6 inhibition. Rather than Ser175, HAB1 cannot dephosphorylate Ser267 in SnRK2.6, which implies that the retained phosphorylation status of Ser267 could ensure that the activated SnRK2.6 reforms the binary complex to cease ABA signaling. Taken together, our findings expand current knowledge concerning the regulation of persulfidation and phosphorylation on the state transition of SnRK2.6 and provide insights into the fine-tuned mechanism of ABA signaling.

Toll-Like Receptor 3 Mediates Aortic Stenosis Through a Conserved Mechanism of Calcification

BACKGROUND

Calcific aortic valve disease (CAVD) is characterized by a phenotypic switch of valvular interstitial cells to bone-forming cells. Toll-like receptors (TLRs) are evolutionarily conserved pattern recognition receptors at the interface between innate immunity and tissue repair. Type I interferons (IFNs) are not only crucial for an adequate antiviral response but also implicated in bone formation. We hypothesized that the accumulation of endogenous TLR3 ligands in the valvular leaflets may promote the generation of osteoblast-like cells through enhanced type I IFN signaling.

METHODS

Human valvular interstitial cells isolated from aortic valves were challenged with mechanical strain or synthetic TLR3 agonists and analyzed for bone formation, gene expression profiles, and IFN signaling pathways. Different inhibitors were used to delineate the engaged signaling pathways. Moreover, we screened a variety of potential lipids and proteoglycans known to accumulate in CAVD lesions as potential TLR3 ligands. Ligand-receptor interactions were characterized by in silico modeling and verified through immunoprecipitation experiments. Biglycan (Bgn), Tlr3, and IFN-α/β receptor alpha chain (Ifnar1)-deficient mice and a specific zebrafish model were used to study the implication of the biglycan (BGN)-TLR3-IFN axis in both CAVD and bone formation in vivo. Two large-scale cohorts (GERA [Genetic Epidemiology Research on Adult Health and Aging], n=55 192 with 3469 aortic stenosis cases; UK Biobank, n=257 231 with 2213 aortic stenosis cases) were examined for genetic variation at genes implicated in BGN-TLR3-IFN signaling associating with CAVD in humans.

RESULTS

Here, we identify TLR3 as a central molecular regulator of calcification in valvular interstitial cells and unravel BGN as a new endogenous agonist of TLR3. Posttranslational BGN maturation by xylosyltransferase 1 (XYLT1) is required for TLR3 activation. Moreover, BGN induces the transdifferentiation of valvular interstitial cells into bone-forming osteoblasts through the TLR3-dependent induction of type I IFNs. It is intriguing that Bgn-/-, Tlr3-/-, and Ifnar1-/- mice are protected against CAVD and display impaired bone formation. Meta-analysis of 2 large-scale cohorts with >300 000 individuals reveals that genetic variation at loci relevant to the XYLT1-BGN-TLR3-interferon-α/β receptor alpha chain (IFNAR) 1 pathway is associated with CAVD in humans.

CONCLUSIONS

This study identifies the BGN-TLR3-IFNAR1 axis as an evolutionarily conserved pathway governing calcification of the aortic valve and reveals a potential therapeutic target to prevent CAVD.

Evolutionary relevance of single nucleotide variants within the forebrain exclusive human accelerated enhancer regions

Background

Human accelerated regions (HARs) are short conserved genomic sequences that have acquired significantly more nucleotide substitutions than expected in the human lineage after divergence from chimpanzees. The fast evolution of HARs may reflect their roles in the origin of human-specific traits. A recent study has reported positively-selected single nucleotide variants (SNVs) within brain-exclusive human accelerated enhancers (BE-HAEs) hs1210 (forebrain), hs563 (hindbrain) and hs304 (midbrain/forebrain). By including data from archaic hominins, these SNVs were shown to be Homo sapiens-specific, residing within transcriptional factors binding sites (TFBSs) for SOX2 (hs1210), RUNX1/3 (hs563), and FOS/JUND (hs304). Although these findings suggest that the predicted modifications in TFBSs may have some role in present-day brain structure, work is required to verify the extent to which these changes translate into functional variation.

Results

To start to fill this gap, we investigate the SOX2 SNV, with both forebrain expression and strong signal of positive selection in humans. We demonstrate that the HMG box of SOX2 binds in vitro with Homo sapiens-specific derived A-allele and ancestral T-allele carrying DNA sites in BE-HAE hs1210. Molecular docking and simulation analysis indicated highly favourable binding of HMG box with derived A-allele containing DNA site when compared to site carrying ancestral T-allele.

Conclusion

These results suggest that adoptive changes in TF affinity within BE-HAE hs1210 and other HAR enhancers in the evolutionary history of Homo sapiens might.

have brought about changes in gene expression patterns and have functional consequences on forebrain formation and evolution.

Methods

The present study employ electrophoretic mobility shift assays (EMSA) and molecular docking and molecular dynamics simulations approaches.

Is BF.7 more infectious than other Omicron subtypes: Insights from structural and simulation studies of BF.7 spike RBD variant

Fear of a fresh infection wave and a global health issue in the ongoing COVID-19 pandemic have been rekindled by the appearance of two new novel variants BF.7 and BA.4/5 of Omicron lineages. Predictions of increased antibody evasion capabilities and transmissibility have been recognised in addition to the existing lineages (BA.1.1, BA.2, BA.2.12.1 and BA.3) as cause for worry. In comparison to Omicron, BA.4 and BF.7 share nine mutations in the spike protein, Leu371Phe, Thr376Ala, Asp405Asn, Arg408Ser, Ser446Gly, Leu452Arg, Phe486Val, Arg493Gln, Ser496Gly, whereas BF.7 contains an additional mutation, Arg346Thr, in the receptor binding domain (RBD) region. Due to the critical need for analysis and data on the BA.4 and BF.7 variants, we have computationally analyzed the interaction pattern between the Omicron, BA.4 and BF.7 RBD and angiotensin-converting enzyme 2 (ACE2) to determine the influence of these unique mutations on the structures, functions, and binding affinity of RBD towards ACE2. These analyses also allow to compare molecular models to previously reported data to evaluate the robustness of our methods for quick prediction of emerging future variants. The docking results reveal that BA.4 and BF.7 have particularly strong interactions with ACE2 when compared to Omicron, as shown by several parameters such as salt bridge, hydrogen bond, and non-bonded interactions. In addition, the estimations of binding free energy corroborated the findings further. BA.4 and BF.7 were found to bind to ACE2 with similar affinities (-72.14 and - 71.54 kcal/mol, respectively) and slightly stronger than Omicron (-70.04 kcal/mol). The differences in the binding pattern between the Omicron, BA.4 and BF.7 variant complexes indicated that the BA.4 and BF.7 RBD substitutions Asp405Asn, Ser446Gly, Leu452Arg, Phe486Val and Arg493Gln caused additional interactions with ACE2. In addition, normal mode analyses also indicate more stable conformations of BA.4 and BF.7 RBDs against human ACE2. Based on these structural and simulation analyses, we hypothesized that these changes may affect the binding affinity of BA.4 and BF.7 with ACE2.

Unveiling mutation effects on the structural dynamics of the main protease from SARS-CoV-2 with hybrid simulation methods

The main protease of SARS-CoV-2 (called M^pro or 3CL^pro) is essential for processing polyproteins encoded by viral RNA. Several M^pro mutations were found in SARS-CoV-2 variants, which are related to higher transmissibility, pathogenicity, and resistance to neutralization antibodies. Macromolecules adopt several favored conformations in solution depending on their structure and shape, determining their dynamics and function. In this study, we used a hybrid simulation method to generate intermediate structures along the six lowest frequency normal modes and sample the conformational space and characterize the structural dynamics and global motions of WT SARS-CoV-2 M^pro and 48 mutations, including mutations found in P.1, B.1.1.7, B.1.351, B.1.525 and B.1.429+B.1.427 variants. We tried to contribute to the elucidation of the effects of mutation in the structural dynamics of SARS-CoV-2 M^pro. A machine learning analysis was performed following the investigation regarding the influence of the K90R, P99L, P108S, and N151D mutations on the dimeric interface assembling of the SARS-CoV-2 M^pro. The parameters allowed the selection of potential structurally stable dimers, which demonstrated that some single surface aa substitutions not located at the dimeric interface (K90R, P99L, P108S, and N151D) are able to induce significant quaternary changes. Furthermore, our results demonstrated, by a Quantum Mechanics method, the influence of SARS-CoV-2 M^pro mutations on the catalytic mechanism, confirming that only one of the chains of the WT and mutant SARS-CoV-2 M^pros are prone to cleave substrates. Finally, it was also possible to identify the aa residue F140 as an important factor related to the increasing enzymatic reactivity of a significant number of SARS-CoV-2 M^pro conformations generated by the normal modes-based simulations.

Mutation pattern of cancer suppressor p53 gene in factory workers exposed to polypropylene and impact on their reproductive hormones

BACKGROUND

Polypropylene is a thermoplastic polymer playing the role of an endocrine disruptor that interferes with the union, emission, transport or elimination of normal hormones. Epidemiological information indicated the relation of endocrine-disturbing chemicals with prostate cancer, testis tumor and diminished fertility. p53 is a key tumor silencer gene. The present study aimed to evaluate luteinizing hormone (LH) and follicle-stimulating hormone (FSH) levels and the risk of p53 mutations as a result of exposure to polypropylene in non-tumorous adult male factory workers.

METHODS

In total, 150 (controls = 35, workers = 115) subjects were recruited. Groups were maintained according to the tenure of exposure G1 (1-5 years), G2 (6-10 years), G3 (11-15 years) and G4 (16-20 years). Concentrations of LH and FSH were determined through an enzyme-linked immunosorbent assay. Genotyping analysis was performed by polymerase chain reaction based gel electrophoresis followed by DNA sequencing. The structural and functional impact of the mutation on the p53 structure was evaluated using 50-ns molecular dynamics (MD) simulations and protein-DNA docking.

RESULTS

Mean plasma LH levels were significantly decreased in G1 (p > 0.05) as well as the G2, G3 and G4 (p > 0.001) groups. Similarly, FSH levels were significant decrease in G1 (p > 0.05), G2 (p > 0.01), G3 (p > 0.001) and G4 (p > 0.001) compared to the control group. Sequencing results found three variants i.e. g.13450 T>G, g.13430C>T and g.13737G>A. One of them was predicted to be disease-causing others are polymorphisms. MD simulation of missense mutation R273H showed no structural impact on the protein structure in MD simulation, but it resulted in weaker binding of p53 with the DNA that might lower the gene expression of cell cycle regulatory proteins.

CONCLUSIONS

These findings predict decreased fertility and risk of malignancies in the future. The spectrum of p53 mutations as a result of polypropylene exposure in the Pakistani population has not been investigated before. Further studies and meta-analyses are required to elucidate the role of different plasticizers in reproduction and cancer-causing risk factors in a larger population.

Vaccine target and carrier molecule nontypeable Haemophilus influenzae protein D dimerizes like the close Escherichia coli GlpQ homolog but unlike other known homolog dimers

We have determined the 1.8 Å X-ray crystal structure of nonlipidated (i.e., N-terminally truncated) nontypeable Haemophilus influenzae (NTHi; H. influenzae) protein D. Protein D exists on outer membranes of H. influenzae strains and acts as a virulence factor that helps invade human cells. Protein D is a proven successful antigen in animal models to treat obstructive pulmonary disease (COPD) and otitis media (OM), and when conjugated to polysaccharides also has been used as a carrier molecule for human vaccines, for example in GlaxoSmithKline Synflorix™. NTHi protein D shares high sequence and structural identify to the Escherichia coli (E. coli) glpQ gene product (GlpQ). E. coli GlpQ is a glycerophosphodiester phosphodiesterase (GDPD) with a known dimeric structure in the Protein Structural Database, albeit without an associated publication. We show here that both structures exhibit similar homodimer organization despite slightly different crystal lattices. Additionally, we have observed both the presence of weak dimerization and the lack of dimerization in solution during size exclusion chromatography (SEC) experiments yet have distinctly observed dimerization in native mass spectrometry analyses. Comparison of NTHi protein D and E. coli GlpQ with other homologous homodimers and monomers shows that the E. coli and NTHi homodimer interfaces are distinct. Despite this distinction, NTHi protein D and E. coli GlpQ possess a triose-phosphate isomerase (TIM) barrel domain seen in many of the other homologs. The active site of NTHi protein D is located near the center of this TIM barrel. A putative glycerol moiety was modeled in two different conformations (occupancies) in the active site of our NTHi protein D structure and we compared this to ligands modeled in homologous structures. Our structural analysis should aid in future efforts to determine structures of protein D bound to substrates, analog intermediates, and products, to fully appreciate this reaction scheme and aiding in future inhibitor design.

A novel consensus-based computational pipeline for screening of antibody therapeutics for efficacy against SARS-CoV-2 variants of concern including Omicron variant

Multiple severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continue to evolve carrying flexible amino acid substitutions in the spike protein's receptor binding domain (RBD). These substitutions modify the binding of the SARS-CoV-2 to human angiotensin-converting enzyme 2 (hACE2) receptor and have been implicated in altered host fitness, transmissibility, and efficacy against antibody therapeutics and vaccines. Reliably predicting the binding strength of SARS-CoV-2 variants RBD to hACE2 receptor and neutralizing antibodies (NAbs) can help assessing their fitness, and rapid deployment of effective antibody therapeutics, respectively. Here, we introduced a two-step computational framework with 3-fold validation that first identified dissociation constant as a reliable predictor of binding affinity in hetero- dimeric and trimeric protein complexes. The second step implements dissociation constant as descriptor of the binding strengths of SARS-CoV-2 variants RBD to hACE2 and NAbs. Then, we examined several variants of concerns (VOCs) such as Alpha, Beta, Gamma, Delta, and Omicron and demonstrated that these VOCs RBD bind to the hACE2 with enhanced affinity. Furthermore, the binding affinity of Omicron variant's RBD was reduced with majority of the RBD-directed NAbs, which is highly consistent with the experimental neutralization data. By studying the atomic contacts between RBD and NAbs, we revealed the molecular footprints of four NAbs (GH-12, P2B-1A1, Asarnow_3D11, and C118)-that may likely neutralize the recently emerged Omicron variant-facilitating enhanced binding affinity. Finally, our findings suggest a computational pathway that could aid researchers identify a range of current NAbs that may be effective against emerging SARS-CoV-2 variants.

Structure-guided mutagenesis of the capsid protein indicates that a nanovirus requires assembled viral particles for systemic infection

Nanoviruses are plant multipartite viruses with a genome composed of six to eight circular single-stranded DNA segments. The distinct genome segments are encapsidated individually in icosahedral particles that measure ≈18 nm in diameter. Recent studies on the model species Faba bean necrotic stunt virus (FBNSV) revealed that complete sets of genomic segments rarely occur in infected plant cells and that the function encoded by a given viral segment can complement the others across neighbouring cells, presumably by translocation of the gene products through unknown molecular processes. This allows the viral genome to replicate, assemble into viral particles and infect anew, even with the distinct genome segments scattered in different cells. Here, we question the form under which the FBNSV genetic material propagates long distance within the vasculature of host plants and, in particular, whether viral particle assembly is required. Using structure-guided mutagenesis based on a 3.2 Å resolution cryogenic-electron-microscopy reconstruction of the FBNSV particles, we demonstrate that specific site-directed mutations preventing capsid formation systematically suppress FBNSV long-distance movement, and thus systemic infection of host plants, despite positive detection of the mutated coat protein when the corresponding segment is agroinfiltrated into plant leaves. These results strongly suggest that the viral genome does not propagate within the plant vascular system under the form of uncoated DNA molecules or DNA:coat-protein complexes, but rather moves long distance as assembled viral particles.

Integration of immunoinformatics and cheminformatics to design and evaluate a multitope vaccine against Klebsiella pneumoniae and Pseudomonas aeruginosa coinfection

Introduction: Klebsiella pneumoniae (K. pneumoniae) and Pseudomonas aeruginosa (P. aeruginosa) are the most common Gram-negative bacteria associated with pneumonia and coinfecting the same patient. Despite their high virulence, there is no effective vaccine against them. Methods: In the current study, the screening of several proteins from both pathogens highlighted FepA and OmpK35 for K. pneumonia in addition to HasR and OprF from P. aeruginosa as promising candidates for epitope mapping. Those four proteins were linked to form a multitope vaccine, that was formulated with a suitable adjuvant, and PADRE peptides to finalize the multitope vaccine construct. The final vaccine's physicochemical features, antigenicity, toxicity, allergenicity, and solubility were evaluated for use in humans. Results: The output of the computational analysis revealed that the designed multitope construct has passed these assessments with satisfactory scores where, as the last stage, we performed a molecular docking study between the potential vaccine construct and K. pneumonia associated immune receptors, TLR4 and TLR2, showing affinitive to both targets with preferentiality for the TLR4 receptor protein. Validation of the docking studies has proceeded through molecular dynamics simulation, which estimated a strong binding and supported the nomination of the designed vaccine as a putative solution for K. pneumoniae and P. aeruginosa coinfection. Here, we describe the approach for the design and assessment of our potential vaccine.

Soluble domains of cytochrome c-556 and Rieske iron-sulfur protein from Chlorobaculum tepidum: Crystal structures and interaction analysis

In photosynthetic green sulfur bacteria, the electron transfer reaction from menaquinol:cytochrome c oxidoreductase to the P840 reaction center (RC) complex occurs directly without any involvement of soluble electron carrier protein(s). X-ray crystallography has determined the three-dimensional structures of the soluble domains of the CT0073 gene product and Rieske iron-sulfur protein (ISP). The former is a mono-heme cytochrome c with an α-absorption peak at 556 nm. The overall fold of the soluble domain of cytochrome c-556 (designated as cyt c-556_sol) consists of four α-helices and is very similar to that of water-soluble cyt c-554 that independently functions as an electron donor to the P840 RC complex. However, the latter's remarkably long and flexible loop between the α3 and α4 helices seems to make it impossible to be a substitute for the former. The structure of the soluble domain of the Rieske ISP (Rieske_sol protein) shows a typical β-sheets-dominated fold with a small cluster-binding and a large subdomain. The architecture of the Rieske_sol protein is bilobal and belongs to those of b₆f-type Rieske ISPs. Nuclear magnetic resonance (NMR) measurements revealed weak non-polar but specific interaction sites on Rieske_sol protein when mixed with cyt c-556_sol. Therefore, menaquinol:cytochrome c oxidoreductase in green sulfur bacteria features a Rieske/cytb complex tightly associated with membrane-anchored cyt c-556.

HBD-2 variants and SARS-CoV-2: New insights into inter-individual susceptibility

Background

A deep understanding of the causes of liability to SARS-CoV-2 is essential to develop new diagnostic tests and therapeutics against this serious virus in order to overcome this pandemic completely. In the light of the discovered role of antimicrobial peptides [such as human b-defensin-2 (hBD-2) and cathelicidin LL-37] in the defense against SARS-CoV-2, it became important to identify the damaging missense mutations in the genes of these molecules and study their role in the pathogenesis of COVID-19.

Methods

We conducted a comprehensive analysis with multiple in silico approaches to identify the damaging missense SNPs for hBD-2 and LL-37; moreover, we applied docking methods and molecular dynamics analysis to study the impact of the filtered mutations.

Results

The comprehensive analysis reveals the presence of three damaging SNPs in hBD-2; these SNPs were predicted to decrease the stability of hBD-2 with a damaging impact on hBD-2 structure as well. G51D and C53G mutations were located in highly conserved positions and were associated with differences in the secondary structures of hBD-2. Docking-coupled molecular dynamics simulation analysis revealed compromised binding affinity for hBD-2 SNPs towards the SARS-CoV-2 spike domain. Different protein-protein binding profiles for hBD-2 SNPs, in relation to their native form, were guided through residue-wise levels and differential adopted conformation/orientation.

Conclusions

The presented model paves the way for identifying patients prone to COVID-19 in a way that would guide the personalization of both the diagnostic and management protocols for this serious disease.

Mapping of Methyl Epitopes of a Peptide-Drug with Its Receptor by 2D STDD-Methyl TROSY NMR Spectroscopy

The current trend in the biopharmaceutical market has boosted the development and production of biological drugs with high efficacy and fidelity for receptor binding. While high-resolution structural insights into binding epitopes of the receptor are indispensable for better therapeutic design, it is tedious and costly. In this work, we develop a protocol by integrating two well-known NMR-based solution-state methods. Saturation transfer double-difference with methyl-TROSY (STDD-Methyl TROSY NMR) was used to probe methyl binding epitopes of the ligand in a label-free environment. This study was carried out with Human insulin as a model peptide drug, with the insulin growth factor receptor (IGFR), which is an off-target receptor for insulin. Methyl epitopes identified from STDD-Methyl TROSY NMR spectroscopy were validated through the HADDOCK platform to generate a drug-receptor model. Since this method can be applied at natural abundance, it has the potential to screen a large set of peptide-drug interactions for optimum receptor binding. Thus, we propose STDD-Methyl TROSY NMR spectroscopy as a technique for rapid screening of biologics for the development of optimized biopharmaceutics.

Densest subgraph-based methods for protein-protein interaction hot spot prediction

Background

Hot spots play an important role in protein binding analysis. The residue interaction network is a key point in hot spot prediction, and several graph theory-based methods have been proposed to detect hot spots. Although the existing methods can yield some interesting residues by network analysis, low recall has limited their abilities in finding more potential hot spots.

Result

In this study, we develop three graph theory-based methods to predict hot spots from only a single residue interaction network. We detect the important residues by finding subgraphs with high densities, i.e., high average degrees. Generally, a high degree implies a high binding possibility between protein chains, and thus a subgraph with high density usually relates to binding sites that have a high rate of hot spots. By evaluating the results on 67 complexes from the SKEMPI database, our methods clearly outperform existing graph theory-based methods on recall and F-score. In particular, our main method, Min-SDS, has an average recall of over 0.665 and an f2-score of over 0.364, while the recall and f2-score of the existing methods are less than 0.400 and 0.224, respectively.

Conclusion

The Min-SDS method performs best among all tested methods on the hot spot prediction problem, and all three of our methods provide useful approaches for analyzing bionetworks. In addition, the densest subgraph-based methods predict hot spots with only one residue interaction network, which is constructed from spatial atomic coordinate data to mitigate the shortage of data from wet-lab experiments.

Glycosylation of a key cubilin Asn residue results in reduced binding to albumin

Kidney disease often manifests with an increase in proteinuria, which can result from both glomerular and/or proximal tubule injury. The proximal tubules are the major site of protein and peptide endocytosis of the glomerular filtrate, and cubilin is the proximal tubule brush border membrane glycoprotein receptor that binds filtered albumin and initiates its processing in proximal tubules. Albumin also undergoes multiple modifications depending upon the physiologic state. We previously documented that carbamylated albumin had reduced cubilin binding, but the effects of cubilin modifications on binding albumin remain unclear. Here, we investigate the cubilin-albumin binding interaction to define the impact of cubilin glycosylation and map the key glycosylation sites while also targeting specific changes in a rat model of proteinuria. We identified a key Asn residue, N1285, that when glycosylated reduced albumin binding. In addition, we found a pH-induced conformation change may contribute to ligand release. To further define the albumin-cubilin binding site, we determined the solution structure of cubilin's albumin-binding domain, CUB7,8, using small-angle X-ray scattering and molecular modeling. We combined this information with mass spectrometry crosslinking experiments of CUB7,8 and albumin that provides a model of the key amino acids required for cubilin-albumin binding. Together, our data supports an important role for glycosylation in regulating the cubilin interaction with albumin, which is altered in proteinuria and provides new insight into the binding interface necessary for the cubilin-albumin interaction.

Design of immunogens for eliciting antibody responses that may protect against SARS-CoV-2 variants

The rise of SARS-CoV-2 variants and the history of outbreaks caused by zoonotic coronaviruses point to the need for next-generation vaccines that confer protection against variant strains. Here, we combined analyses of diverse sequences and structures of coronavirus spikes with data from deep mutational scanning to design SARS-CoV-2 variant antigens containing the most significant mutations that may emerge. We trained a neural network to predict RBD expression and ACE2 binding from sequence, which allowed us to determine that these antigens are stable and bind to ACE2. Thus, they represent viable variants. We then used a computational model of affinity maturation (AM) to study the antibody response to immunization with different combinations of the designed antigens. The results suggest that immunization with a cocktail of the antigens is likely to promote evolution of higher titers of antibodies that target SARS-CoV-2 variants than immunization or infection with the wildtype virus alone. Finally, our analysis of 12 coronaviruses from different genera identified the S2’ cleavage site and fusion peptide as potential pan-coronavirus vaccine targets.

SARS-CoV-2 variants have already emerged and future variants may pose greater threats to the efficacy of current vaccines. Rather than using a reactive approach to vaccine development that would lag behind the evolution of the virus, such as updating the sequence in the vaccine with a current variant, we sought to use a proactive approach that predicts some of the mutations that could arise that could evade current immune responses. Then, by including these mutations in a new vaccine antigen, we might be able to protect against those potential variants before they appear. Toward this end, we used various computational methods including sequence analysis and machine learning to design such antigens. We then used simulations of antibody development, and the results suggest that immunization with our designed antigens is likely to result in an antibody response that is better able to target SARS-CoV-2 variants than current vaccines. We also leveraged our sequence analysis to suggest that a particular site on the spike protein could serve as a useful target for a pan-coronavirus vaccine.

Prediction of the structural interface between fibroblast growth factor23 and Burosumab using alanine scanning and molecular docking

Burosumab, an FGF23 targeting monoclonal antibody, was approved by the FDA in 2018 for use in children and adults with X-linked hypophosphatemia (or XLH). While several clinical studies have demonstrated the long-term safety and efficacy of Burosumab, the molecular basis of FGF23-Burosumab interaction which underpins its mechanism of action remains unknown. In this study, we employed molecular docking combined with alanine scanning of epitope and paratope to predict a model of FGF23-Burosumab interaction. Then, we used the model to understand the species-species cross-reactivity of Burosumab and to reverse engineer mouse FGF23 with 'back to human' mutations to bind Burosumab. Finally, we redesigned the CDRs with two mutations to engineer an affinity enhanced variant of the antibody. Our study provides insights into the FGF23-Burosumab interaction and demonstrates that alanine-scanning coupled with molecular docking can be used to optimize antibody candidates (e.g., structure-guided affinity maturation) for therapeutic use.

The Hidden Intricacies of Aquaporins: Remarkable Details in a Common Structural Scaffold

Evolution turned aquaporins (AQPs) into the most efficient facilitators of passive water flow through cell membranes at no expense of solute discrimination. In spite of a plethora of solved AQP structures, many structural details remain hidden. Here, by combining extensive sequence- and structural-based analysis of a unique set of 20 non-redundant high-resolution structures and molecular dynamics simulations of four representatives, key aspects of AQP stability, gating, selectivity, pore geometry, and oligomerization, with a potential impact on channel functionality, are identified. The general view of AQPs possessing a continuous open water pore is challenged and it is depicted that AQPs' selectivity is not exclusively shaped by pore-lining residues but also by the relative arrangement of transmembrane helices. Moreover, this analysis reveals that hydrophobic interactions constitute the main determinant of protein thermal stability. Finally, a numbering scheme of the conserved AQP scaffold is established, facilitating direct comparison of, for example, disease-causing mutations and prediction of potential structural consequences. Additionally, the results pave the way for the design of optimized AQP water channels to be utilized in biotechnological applications.

Bivalent BET Bromodomain Inhibitors Confer Increased Potency and Selectivity for BRDT via Protein Conformational Plasticity

Bromodomain and extraterminal domain (BET) proteins are important regulators of gene transcription and chromatin remodeling. BET family members BRD4 and BRDT are validated targets for cancer and male contraceptive drug development, respectively. Due to the high structural similarity of the acetyl-lysine binding sites, most reported inhibitors lack intra-BET selectivity. We surmised that protein-protein interactions induced by bivalent inhibitors may differ between BRD4 and BRDT, conferring an altered selectivity profile. Starting from nonselective monovalent inhibitors, we developed cell-active bivalent BET inhibitors with increased activity and selectivity for BRDT. X-ray crystallographic and solution studies revealed unique structural states of BRDT and BRD4 upon interaction with bivalent inhibitors. Varying spacer lengths and symmetric vs unsymmetric connections resulted in the same dimeric states, whereas different chemotypes induced different dimers. The findings indicate that the increased intra-BET selectivity of bivalent inhibitors is due to the differential plasticity of BET bromodomains upon inhibitor-induced dimerization.

HUG Domain Is Responsible for Active Dimer Stabilization in an NrdJd Ribonucleotide Reductase

Ribonucleotide reductases (RNRs) catalyze the reduction of ribonucleotides to the corresponding deoxyribonucleotides. The catalytic activity of most RNRs depends on the formation of a dimer of the catalytic subunits. The active site is located at the interface, and part of the substrate binding site and regulatory mechanisms work across the subunit in the dimer. In this study, we describe and characterize a novel domain responsible for forming the catalytic dimer in several class II RNRs. The 3D structure of the class II RNR from Rhodobacter sphaeroides reveals a so far undescribed α-helical domain in the dimer interface, which is embracing the other subunit. Genetic removal of this HUG domain leads to a severe reduction of activity paired with reduced dimerization capability. In comparison with other described RNRs, the enzyme with this domain is less dependent on the presence of nucleotides to act as allosteric effectors in the formation of dimers. The HUG domain appears to serve as an interlock to keep the dimer intact and functional even at low enzyme and/or effector concentrations.

Brachypodium Antifreeze Protein Gene Products Inhibit Ice Recrystallisation, Attenuate Ice Nucleation, and Reduce Immune Response

Antifreeze proteins (AFPs) from the model crop, Brachypodium distachyon, allow freeze survival and attenuate pathogen-mediated ice nucleation. Intriguingly, Brachypodium AFP genes encode two proteins, an autonomous AFP and a leucine-rich repeat (LRR). We present structural models which indicate that ice-binding motifs on the ~13 kDa AFPs can “spoil” nucleating arrays on the ~120 kDa bacterial ice nucleating proteins used to form ice at high sub-zero temperatures. These models are consistent with the experimentally demonstrated decreases in ice nucleating activity by lysates from wildtype compared to transgenic Brachypodium lines. Additionally, the expression of Brachypodium LRRs in transgenic Arabidopsis inhibited an immune response to pathogen flagella peptides (flg22). Structural models suggested that this was due to the affinity of the LRR domains to flg22. Overall, it is remarkable that the Brachypodium genes play multiple distinctive roles in connecting freeze survival and anti-pathogenic systems via their encoded proteins’ ability to adsorb to ice as well as to attenuate bacterial ice nucleation and the host immune response.

Structural Characterization of Two Short Unspecific Peroxygenases: Two Different Dimeric Arrangements

Unspecific peroxygenases (UPOs) are extracellular fungal enzymes of biotechnological interest as self-sufficient (and more stable) counterparts of cytochrome P450 monooxygenases, the latter being present in most living cells. Expression hosts and structural information are crucial for exploiting UPO diversity (over eight thousand UPO-type genes were identified in sequenced genomes) in target reactions of industrial interest. However, while many thousands of entries in the Protein Data Bank include molecular coordinates of P450 enzymes, only 19 entries correspond to UPO enzymes, and UPO structures from only two species (Agrocybe aegerita and Hypoxylon sp.) have been published to date. In the present study, two UPOs from the basidiomycete Marasmius rotula (rMroUPO) and the ascomycete Collariella virescens (rCviUPO) were crystallized after sequence optimization and Escherichia coli expression as active soluble enzymes. Crystals of rMroUPO and rCviUPO were obtained at sufficiently high resolution (1.45 and 1.95 Å, respectively) and the corresponding structures were solved by molecular replacement. The crystal structures of the two enzymes (and two mutated variants) showed dimeric proteins. Complementary biophysical and molecular biology studies unveiled the diverse structural bases of the dimeric nature of the two enzymes. Intermolecular disulfide bridge and parallel association between two α-helices, among other interactions, were identified at the dimer interfaces. Interestingly, one of the rCviUPO variants incorporated the ability to produce fatty acid diepoxides—reactive compounds with valuable cross-linking capabilities—due to removal of the enzyme C-terminal tail located near the entrance of the heme access channel. In conclusion, different dimeric arrangements could be described in (short) UPO crystal structures.

Over the rainbow: structural characterization of the chromoproteins gfasPurple, amilCP, spisPink and eforRed

Anthozoan chromoproteins are highly pigmented, diversely coloured and readily produced in recombinant expression systems. While they are a versatile and powerful building block in synthetic biology for applications such as biosensor development, they are not widely used in comparison to the related fluorescent proteins, partly due to a lack of structural characterization to aid protein engineering. Here, high-resolution X-ray crystal structures of four open-source chromoproteins, gfasPurple, amilCP, spisPink and eforRed, are presented. These proteins are dimers in solution, and mutation at the conserved dimer interface leads to loss of visible colour development in gfasPurple. The chromophores are trans and noncoplanar in gfasPurple, amilCP and spisPink, while that in eforRed is cis and noncoplanar, and also emits fluorescence. Like other characterized chromoproteins, gfasPurple, amilCP and eforRed contain an sp 2-hybridized N-acylimine in the peptide bond preceding the chromophore, while spisPink is unusual and demonstrates a true sp 3-hybridized trans-peptide bond at this position. It was found that point mutations at the chromophore-binding site in gfasPurple that substitute similar amino acids to those in amilCP and spisPink generate similar colours. These features and observations have implications for the utility of these chromoproteins in protein engineering and synthetic biology applications.

High-valency Anti-CD99 Antibodies Toward the Treatment of T Cell Acute Lymphoblastic Leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive form of leukemia that currently requires intensive chemotherapy. While childhood T-ALL is associated with high cure rates, adult T-ALL is not, and both are associated with significant short- and long-term morbidities. Thus, less toxic and effective strategies to treat T-ALL are needed. CD99 is overexpressed on T-ALL blasts at diagnosis and at relapse. Although targeting CD99 with cytotoxic antibodies has been proposed, the molecular features required for their activity are undefined. We identified human antibodies that selectively bound to the extracellular domain of human CD99, and the most potent clone, 10A1, shared an epitope with a previously described cytotoxic IgM antibody. We engineered clone 10A1 in bivalent, trivalent, tetravalent, and dodecavalent formats. Increasing the antibody valency beyond two had no effects on binding to T-ALL cells. In contrast, a valency of ≥3 was required for cytotoxicity, suggesting a mechanism of action in which an antibody clusters ≥3 CD99 molecules to induce cytotoxicity. We developed a human IgG-based tetravalent version of 10A1 that exhibited cytotoxic activity to T-ALL cells but not to healthy peripheral blood cells. The crystal structure of the 10A1 Fab in complex with a CD99 fragment revealed that the antibody primarily recognizes a proline-rich motif (PRM) of CD99 in a manner reminiscent of SH3-PRM interactions. This work further validates CD99 as a promising therapeutic target in T-ALL and defines a pathway toward the development of a selective therapy against T-ALL.

Mechanisms of Cre recombinase synaptic complex assembly and activation illuminated by Cryo-EM

Cre recombinase selectively recognizes DNA and prevents non-specific DNA cleavage through an orchestrated series of assembly intermediates. Cre recombines two loxP DNA sequences featuring a pair of palindromic recombinase binding elements and an asymmetric spacer region, by assembly of a tetrameric synaptic complex, cleavage of an opposing pair of strands, and formation of a Holliday junction intermediate. We used Cre and loxP variants to isolate the monomeric Cre-loxP (54 kDa), dimeric Cre₂-loxP (110 kDa), and tetrameric Cre₄-loxP₂ assembly intermediates, and determined their structures using cryo-EM to resolutions of 3.9, 4.5 and 3.2 Å, respectively. Progressive and asymmetric bending of the spacer region along the assembly pathway enables formation of increasingly intimate interfaces between Cre protomers and illuminates the structural bases of biased loxP strand cleavage order and half-the-sites activity. Application of 3D variability analysis to the tetramer data reveals constrained conformational sampling along the pathway between protomer activation and Holliday junction isomerization. These findings underscore the importance of protein and DNA flexibility in Cre-mediated site selection, controlled activation of alternating protomers, the basis for biased strand cleavage order, and recombination efficiency. Such considerations may advance development of site-specific recombinases for use in gene editing applications.

Structural basis of cyclic oligoadenylate binding to the transcription factor Csa3 outlines cross talk between type III and type I CRISPR systems

RNA interference by type III CRISPR systems results in the synthesis of cyclic oligoadenylate (cOA) second messengers, which are known to bind and regulate various CARF domain-containing nuclease receptors. The CARF domain-containing Csa3 family of transcriptional factors associated with the DNA-targeting type I CRISPR systems regulate expression of various CRISPR and DNA repair genes in many prokaryotes. In this study, we extend the known receptor repertoire of cOA messengers to include transcriptional factors by demonstrating specific binding of cyclic tetra-adenylate (cA4) to Saccharolobus solfataricus Csa3 (Csa3Sso). Our 2.0-Å resolution X-ray crystal structure of cA4-bound full-length Csa3Sso reveals the binding of its CARF domain to an elongated conformation of cA4. Using cA4 binding affinity analyses of Csa3Sso mutants targeting the observed Csa3Sso•cA4 structural interface, we identified a Csa3-specific cA4 binding motif distinct from a more widely conserved cOA-binding CARF motif. Using a rational surface engineering approach, we increased the cA4 binding affinity of Csa3Sso up to ∼145-fold over the wildtype, which has potential applications for future second messenger-driven CRISPR gene expression and editing systems. Our in-solution Csa3Sso structural analysis identified cA4-induced allosteric and asymmetric conformational rearrangement of its C-terminal winged helix-turn-helix effector domains, which could potentially be incompatible to DNA binding. However, specific in vitro binding of the purified Csa3Sso to its putative promoter (PCas4a) was found to be cA4 independent, suggesting a complex mode of Csa3Sso regulation. Overall, our results support cA4-and Csa3-mediated cross talk between type III and type I CRISPR systems.

Thermal Stability and Inhibitory Action of Red Grape Skin Phytochemicals against Enzymes Associated with Metabolic Syndrome

The present study focuses on heat-induced structural changes and the degradation kinetics of phytochemicals and antioxidant activity of red grape skin extract. The thermal degradation of anthocyanins, flavonoids, polyphenols, and antioxidant activity followed a first-order kinetic model, increasing with temperature due to the intensification of the degradation process. The activation energy (Ea) highlighted this phenomenon. Likewise, the kinetic and thermodynamic parameters certified the irreversible degradation of the bioactive compounds from the skin of the Băbească neagră grape variety. Both temperature and duration of heating had a significant impact on the content of bioactive compounds. In addition, the red grape skin extract inhibited certain enzymes such as α-amylase, α-glucosidase, lipase, and lipoxygenase, which are associated with metabolic syndrome and inflammation. Further knowledge on the possible inhibition mechanisms exerted by the major anthocyanins found in red grape skin extract on the metabolic syndrome-associated enzymes was gathered upon running molecular docking tests. Detailed analysis of the resulting molecular models revealed that malvidin 3-O-glucoside binds in the vicinity of the catalytic site of α-amylase and lipase, whereas no direct contact with catalytic amino acids was identified in the case of α-glucosidase and lipoxygenase.

In silico study on the effects of disulfide bonds in ORF8 of SARS-CoV-2

The COVID-19 epidemic, caused by virus SARS-CoV-2, has been a pandemic and threatening everyone's health in the past two years. In SARS-CoV-2, ORF8 is one of the most important accessory...

Protein-Protein Docking: Past, Present, and Future

The biological significance of proteins attracted the scientific community in exploring their characteristics. The studies shed light on the interaction patterns and functions of proteins in a living body. Due to their practical difficulties, reliable experimental techniques pave the way for introducing computational methods in the interaction prediction. Automated methods reduced the difficulties but could not yet replace experimental studies as the field is still evolving. Interaction prediction problem being critical needs highly accurate results, but none of the existing methods could offer reliable performance that can parallel with experimental results yet. This article aims to assess the existing computational docking algorithms, their challenges, and future scope. Blind docking techniques are quite helpful when no information other than the individual structures are available. As more and more complex structures are being added to different databases, information-driven approaches can be a good alternative. Artificial intelligence, ruling over the major fields, is expected to take over this domain very shortly.

Naturally occurring cancer-associated mutations disrupt oligomerization and activity of protein arginine methyltransferase 1 (PRMT1)

Protein arginine methylation is a posttranslational modification catalyzed by the protein arginine methyltransferase (PRMT) enzyme family. Dysregulated protein arginine methylation is linked to cancer and a variety of other human diseases. PRMT1 is the predominant PRMT isoform in mammalian cells and acts in pathways regulating transcription, DNA repair, apoptosis, and cell proliferation. PRMT1 dimer formation, which is required for methyltransferase activity, is mediated by interactions between a structure called the dimerization arm on one monomer and a surface of the Rossman Fold of the other monomer. Given the link between PRMT1 dysregulation and disease and the link between PRMT1 dimerization and activity, we searched the Catalogue of Somatic Mutations in Cancer (COSMIC) database to identify potential inactivating mutations occurring in the PRMT1 dimerization arm. We identified three mutations that correspond to W215L, Y220N, and M224V substitutions in human PRMT1V2 (isoform 1) (W197L, Y202N, M206V in rat PRMT1V1). Using a combination of site-directed mutagenesis, analytical ultracentrifugation, native PAGE, and activity assays, we found that these conservative substitutions surprisingly disrupt oligomer formation and substantially impair both S-adenosyl-L-methionine (AdoMet) binding and methyltransferase activity. Molecular dynamics simulations suggest that these substitutions introduce novel interactions within the dimerization arm that lock it in a conformation not conducive to dimer formation. These findings provide a clear, if putative, rationale for the contribution of these mutations to impaired arginine methylation in cells and corresponding health consequences.

DYW domain structures imply an unusual regulation principle in plant organellar RNA editing catalysis

RNA editosomes selectively deaminate cytidines to uridines in plant organellar transcripts–mostly to restore protein functionality and consequently facilitate mitochondrial and chloroplast function. The RNA editosomal pentatricopeptide repeat proteins serve target RNA recognition, whereas the intensively studied DYW domain elicits catalysis. Here we present structures and functional data of a DYW domain in an inactive ground state and activated. DYW domains harbour a cytidine deaminase fold and a C-terminal DYW motif, with catalytic and structural zinc atoms, respectively. A conserved gating domain within the deaminase fold regulates the active site sterically and mechanistically in a process that we termed gated zinc shutter. Based on the structures, an autoinhibited ground state and its activation are cross-validated by RNA editing assays and differential scanning fluorimetry. We anticipate that, in vivo, the framework of an active plant RNA editosome triggers the release of DYW autoinhibition to ensure a controlled and coordinated cytidine deamination playing a key role in mitochondrial and chloroplast homeostasis.

Whey Protein Isolate-Xylose Maillard-Based Conjugates with Tailored Microencapsulation Capacity of Flavonoids from Yellow Onions Skins

The objective of this study is to encapsulate flavonoids from yellow onion skins in whey protein isolates (WPI) and xylose (X), by Maillard-based conjugates, as an approach to improve the ability to entrap flavonoids and to develop powders with enhanced antioxidant activity. WPI (0.6%, w/v) was conjugated to X (0.3%, w/v) through the Maillard reaction at 90 °C for 120 min, in the presence of a flavonoid-enriched extract. Two variants of powders were obtained by freeze-drying. The glycation of WPI allowed a better encapsulation efficiency, up to 90.53 ± 0.29%, corresponding to a grafting degree of 30.38 ± 1.55%. The molecular modelling approach was used to assess the impact of X interactions with α-lactalbumin and β-lactoglobulin on the ability of these proteins to bind the main flavonoids from the yellow onion skins. The results showed that X might compete with quercetin glucosides to bind with α-lactalbumin. No interference was found in the case of β-lactoglobulin. The microstructural appearance of the powders revealed finer spherosomes in powder with WPI-X conjugates via the Maillard reaction. The powders were added to nachos, followed by a phytochemical characterization, in order to test their potential added value. An increase in antioxidant activity was observed, with no significant changes during storage.