Structural basis for antiepileptic drugs and botulinum neurotoxin recognition of SV2A.

More than one percent of people have epilepsy worldwide. Levetiracetam (LEV) is a successful new-generation antiepileptic drug (AED), and its derivative, brivaracetam (BRV), shows improved efficacy. Synaptic vesicle glycoprotein 2a (SV2A), a putative membrane transporter in the synaptic vesicles (SVs), has been identified as a target of LEV and BRV. SV2A also serves as a receptor for botulinum neurotoxin (BoNT), which is the most toxic protein and has paradoxically emerged as a potent reagent for therapeutic and cosmetic applications. Nevertheless, no structural analysis on AEDs and BoNT recognition by full-length SV2A has been available. Here we describe the cryo-electron microscopy structures of the full-length SV2A in complex with the BoNT receptor-binding domain, BoNT/A2 HC, and either LEV or BRV. The large fourth luminal domain of SV2A binds to BoNT/A2 HC through protein-protein and protein-glycan interactions. LEV and BRV occupy the putative substrate-binding site in an outward-open conformation. A propyl group in BRV creates additional contacts with SV2A, explaining its higher binding affinity than that of LEV, which was further supported by label-free spectral shift assay. Numerous LEV derivatives have been developed as AEDs and positron emission tomography (PET) tracers for neuroimaging. Our work provides a structural framework for AEDs and BoNT recognition of SV2A and a blueprint for the rational design of additional AEDs and PET tracers.

QM/MM Study of the Catalytic Mechanism and Substrate Specificity of the Aromatic Substrate C-Methyltransferase Fur6.

In the field of medical chemistry and other organic chemistry, introducing a methyl group into a designed position has been difficult to achieve. However, owing to the vigorous developments in the field of enzymology, methyltransferases are considered potential tools for addressing this problem. Within the methyltransferase family, Fur6 catalyzes the methylation of C3 of 1,2,4,5,7-pentahydroxynaphthalene (PHN) using S-adenosyl-l-methionine (SAM) as the methyl donor. Here, we report the catalytic mechanism and substrate specificity of Fur6 based on computational studies. Our molecular dynamics (MD) simulation studies reveal the reactive form of PHN and its interactions with the enzyme. Our hybrid quantum mechanics/molecular mechanics (QM/MM) calculations suggest the reaction pathway of the methyl transfer step in which the energy barrier is 8.6 kcal mol-1. Our free-energy calculations with a polarizable continuum model (PCM) indicate that the final deprotonation step of the methylated intermediate occurs after it is ejected into the water solvent from the active center pocket of Fur6. Additionally, our studies on the protonation states, the highest occupied molecular orbital (HOMOs), and the energy barriers of the methylation reaction for the analogs of PHN demonstrate the mechanism of the specificity to PHN. Our study provides valuable insights into Fur6 chemistry, contributing to a deeper understanding of molecular mechanisms and offering an opportunity to engineer the enzyme to achieve high yields of the desired product(s).

Enhancement of SARS-CoV-2 Infection via Crosslinking of Adjacent Spike Proteins by N-Terminal Domain-Targeting Antibodies.

The entry of SARS-CoV-2 into host cells is mediated by the interaction between the spike receptor-binding domain (RBD) and host angiotensin-converting enzyme 2 (ACE2). Certain human antibodies, which target the spike N-terminal domain (NTD) at a distant epitope from the host cell binding surface, have been found to augment ACE2 binding and enhance SARS-CoV-2 infection. Notably, these antibodies exert their effect independently of the antibody fragment crystallizable (Fc) region, distinguishing their mode of action from previously described antibody-dependent infection-enhancing (ADE) mechanisms. Building upon previous hypotheses and experimental evidence, we propose that these NTD-targeting infection-enhancing antibodies (NIEAs) achieve their effect through the crosslinking of neighboring spike proteins. In this study, we present refined structural models of NIEA fragment antigen-binding region (Fab)-NTD complexes, supported by molecular dynamics simulations and hydrogen-deuterium exchange mass spectrometry (HDX-MS). Furthermore, we provide direct evidence confirming the crosslinking of spike NTDs by NIEAs. Collectively, our findings advance our understanding of the molecular mechanisms underlying NIEAs and their impact on SARS-CoV-2 infection.

Comprehensive computational analysis of the SRK-SP11 molecular interaction underlying self-incompatibility in Brassicaceae using improved structure prediction for cysteine-rich proteins.

Plants employ self-incompatibility (SI) to promote cross-fertilization. In Brassicaceae, this process is regulated by the formation of a complex between the pistil determinant S receptor kinase (SRK) and the pollen determinant S-locus protein 11 (SP11, also known as S-locus cysteine-rich protein, SCR). In our previous study, we used the crystal structures of two eSRK-SP11 complexes in Brassica rapa S8 and S9 haplotypes and nine computationally predicted complex models to demonstrate that only the SRK ectodomain (eSRK) and SP11 pairs derived from the same S haplotype exhibit high binding free energy. However, predicting the eSRK-SP11 complex structures for the other 100 + S haplotypes and genera remains difficult because of SP11 polymorphism in sequence and structure. Although protein structure prediction using AlphaFold2 exhibits considerably high accuracy for most protein monomers and complexes, 46% of the predicted SP11 structures that we tested showed < 75 mean per-residue confidence score (pLDDT). Here, we demonstrate that the use of curated multiple sequence alignment (MSA) for cysteine-rich proteins significantly improved model accuracy for SP11 and eSRK-SP11 complexes. Additionally, we calculated the binding free energies of the predicted eSRK-SP11 complexes using molecular dynamics (MD) simulations and observed that some Arabidopsis haplotypes formed a binding mode that was critically different from that of B. rapa S8 and S9. Thus, our computational results provide insights into the haplotype-specific eSRK-SP11 binding modes in Brassicaceae at the residue level. The predicted models are freely available at Zenodo, https://doi.org/10.5281/zenodo.8047768.

Site-specific photo-crosslinking/cleavage for protein-protein interface identification reveals oligomeric assembly of lysosomal-associated membrane protein type 2A in mammalian cells.

Genetic code expansion enables site-specific photo-crosslinking by introducing photo-reactive non-canonical amino acids into proteins at defined positions during translation. This technology is widely used for analyzing protein-protein interactions and is applicable in mammalian cells. However, the identification of the crosslinked region still remains challenging. Here, we developed a new method to identify the crosslinked region by pre-installing a site-specific cleavage site, an α-hydroxy acid (Nε -allyloxycarbonyl-α-hydroxyl-l-lysine acid, AllocLys-OH), into the target protein. Alkaline treatment cleaves the crosslinked complex at the position of the α-hydroxy acid residue and thus helps to identify which side of the cleavage site, either closer to the N-terminus or C-terminus, the crosslinked site is located within the target protein. A series of AllocLys-OH introductions narrows down the crosslinked region. By applying this method, we identified the crosslinked regions in lysosomal-associated membrane protein type 2A (LAMP2A), a receptor of chaperone-mediated autophagy, in mammalian cells. The results suggested that at least two interfaces are involved in the homophilic interaction, which requires a trimeric or higher oligomeric assembly of adjacent LAMP2A molecules. Thus, the combination of site-specific crosslinking and site-specific cleavage promises to be useful for revealing binding interfaces and protein complex geometries.

Chemigenetic indicators based on synthetic chelators and green fluorescent protein.

Molecular fluorescent indicators are versatile tools for dynamic imaging of biological systems. We now report a class of indicators that are based on the chemigenetic combination of a synthetic ion-recognition motif and a protein-based fluorophore. Specifically, we have developed a calcium ion (Ca2+) indicator that is based on genetic insertion of circularly permuted green fluorescent protein into HaloTag protein self-labeled with a ligand containing the Ca2+ chelator 1,2-bis(2-aminophenoxy)ethane-N,N,N',N'-tetraacetic acid. We have demonstrated the versatility of this design by also developing a sodium ion (Na+) indicator using a crown-ether-containing ligand. This approach affords bright and sensitive ion indicators that can be applicable to cell imaging. This design can enable the development of chemigenetic indicators with ion or molecular specificities that have not been realized with fully protein-based indicators.

Calculations of the binding free energies of the Comprehensive in vitro Proarrhythmia Assay (CiPA) reference drugs to cardiac ion channels.

The evaluation of the inhibitory activities of drugs on multiple cardiac ion channels is required for the accurate assessment of proarrhythmic risks. Moreover, the in silico prediction of such inhibitory activities of drugs on cardiac channels can improve the efficiency of the drug-development process. Here, we performed molecular docking simulations to predict the complex structures of 25 reference drugs that were proposed by the Comprehensive in vitro Proarrhythmia Assay consortium using two cardiac ion channels, the human ether-a-go-go-related gene (hERG) potassium channel and human NaV1.5 (hNaV1.5) sodium channel, with experimentally available structures. The absolute binding free energy (ΔGbind) values of the predicted structures were calculated by a molecular dynamics-based method and compared with the experimental half-maximal inhibitory concentration (IC50) data. Furthermore, the regression analysis between the calculated values and negative of the common logarithm of the experimental IC50 values (pIC50) revealed that the calculated values of four and ten drugs deviated significantly from the regression lines of the hERG and hNaV1.5 channels, respectively. We reconsidered the docking poses and protonation states of the drugs based on the experimental data and recalculated their ΔGbind values. Finally, the calculated ΔGbind values of 24 and 19 drugs correlated with their experimental pIC50 values (coefficients of determination=0.791 and 0.613 for the hERG and hNaV1.5 channels, respectively). Thus, the regression analysis between the calculated ΔGbind and experimental IC50 data ensured the realization of an increased number of reliable complex structures.

Feedback-AVPGAN: Feedback-guided generative adversarial network for generating antiviral peptides.

In this study, we propose Feedback-AVPGAN, a system that aims to computationally generate novel antiviral peptides (AVPs). This system relies on the key premise of the Generative Adversarial Network (GAN) model and the Feedback method. GAN, a generative modeling approach that uses deep learning methods, comprises a generator and a discriminator. The generator is used to generate peptides; the generated proteins are fed to the discriminator to distinguish between the AVPs and non-AVPs. The original GAN design uses actual data to train the discriminator. However, not many AVPs have been experimentally obtained. To solve this problem, we used the Feedback method to allow the discriminator to learn from the existing as well as generated synthetic data. We implemented this method using a classifier module that classifies each peptide sequence generated by the GAN generator as AVP or non-AVP. The classifier uses the transformer network and achieves high classification accuracy. This mechanism enables the efficient generation of peptides with a high probability of exhibiting antiviral activity. Using the Feedback method, we evaluated various algorithms and their performance. Moreover, we modeled the structure of the generated peptides using AlphaFold2 and determined the peptides having similar physicochemical properties and structures to those of known AVPs, although with different sequences.

Brassinosteroid-induced gene repression requires specific and tight promoter binding of BIL1/BZR1 via DNA shape readout.

BRZ-INSENSITIVE-LONG 1 (BIL1)/BRASSINAZOLE-RESISTANT 1 (BZR1) and its homologues are plant-specific transcription factors that convert the signalling of the phytohormones brassinosteroids (BRs) to transcriptional responses, thus controlling various physiological processes in plants. Although BIL1/BZR1 upregulates some BR-responsive genes and downregulates others, the molecular mechanism underlying the dual roles of BIL1/BZR1 is still poorly understood. Here we show that BR-responsive transcriptional repression by BIL1/BZR1 requires the tight binding of BIL1/BZR1 alone to the 10 bp elements of DNA fragments containing the known 6 bp core-binding motifs at the centre. Furthermore, biochemical and structural evidence demonstrates that the selectivity for two nucleobases flanking the core motifs is realized by the DNA shape readout of BIL1/BZR1 without direct recognition of the nucleobases. These results elucidate the molecular and structural basis of transcriptional repression by BIL1/BZR1 and contribute to further understanding of the dual roles of BIL1/BZR1 in BR-responsive gene regulation.

Uptake mechanism of iron-phytosiderophore from the soil based on the structure of yellow stripe transporter.

Calcareous soils cover one-third of all land and cause severe growth defects in plants due to the poor water solubility of iron at high pH. Poaceae species use a unique chelation strategy, whereby plants secrete a high-affinity metal chelator, known as phytosiderophores (mugineic acids), and reabsorb the iron-phytosiderophore complex by the yellow stripe 1/yellow stripe 1-like (YS1/YSL) transporter for efficient uptake of iron from the soil. Here, we present three cryo-electron microscopy structures of barley YS1 (HvYS1) in the apo state, in complex with an iron-phytosiderophore complex, Fe(III)-deoxymugineic acid (Fe(III)-DMA), and in complex with the iron-bound synthetic DMA analog (Fe(III)-PDMA). The structures reveal a homodimeric assembly mediated through an anti-parallel ß-sheet interaction with cholesterol hemisuccinate. Each protomer adopts an outward open conformation, and Fe(III)-DMA is bound near the extracellular space in the central cavity. Fe(III)-PDMA occupies the same binding site as Fe(III)-DMA, demonstrating that PDMA can function as a potent fertilizer in an essentially identical manner to DMA. Our results provide a structural framework for iron-phytosiderophore recognition and transport by YS1/YSL transporters, which will enable the rational design of new, high-potency fertilizers.

Prediction of protein mononucleotide binding sites using AlphaFold2 and machine learning.

In this study, we developed a system that predicts the binding sites of proteins for five mononucleotides (AMP, ADP, ATP, GDP, and GTP). The system comprises two machine learning (ML)-based predictors using a convolutional neural network and a gradient boosting machine, two template-based predictors based on sequence and structure alignment, and a predictor that performs ensemble learning of these four predictors. In this study, data augmentation of ligand binding sites with similar ligand structures was performed. For example, in the prediction of ADP-binding sites using ML methods, the binding sites of AMP and ATP, which have similar structures, are considered. In addition, we constructed the structure models using AlphaFold2, a highly accurate protein prediction method. The secondary structure and dihedral angle information obtained using the model structures were used as ML predictor features. Additionally, in the template-based predictor, the structures of the binding sites were used as templates to be explored based on structure alignment to identify the binding site of the target. Consequently, the template-based predictor based on structure alignment showed the best performance among the four individual predictors, and the ensemble predictor achieved the best performance, with an area under the curve of 0.958 for all mononucleotides.

The α- and ß-Subunit Boundary at the Stem of the Mushroom-Like α3ß3-Type Oxygenase Component of Rieske Non-Heme Iron Oxygenases Is the Rieske-Type Ferredoxin-Binding Site.

Cumene dioxygenase (CumDO) is an initial enzyme in the cumene degradation pathway of Pseudomonas fluorescens IP01 and is a Rieske non-heme iron oxygenase (RO) that comprises two electron transfer components (reductase [CumDO-R] and Rieske-type ferredoxin [CumDO-F]) and one catalytic component (α3ß3-type oxygenase [CumDO-O]). Catalysis is triggered by electrons that are transferred from NAD(P)H to CumDO-O by CumDO-R and CumDO-F. To investigate the binding mode between CumDO-F and CumDO-O and to identify the key CumDO-O amino acid residues for binding, we simulated docking between the CumDO-O crystal structure and predicted model of CumDO-F and identified two potential binding sites: one is at the side-wise site and the other is at the top-wise site in mushroom-like CumDO-O. Then, we performed alanine mutagenesis of 16 surface amino acid residues at two potential binding sites. The results of reduction efficiency analyses using the purified components indicated that CumDO-F bound at the side-wise site of CumDO-O, and K117 of the α-subunit and R65 of the ß-subunit were critical for the interaction. Moreover, these two positively charged residues are well conserved in α3ß3-type oxygenase components of ROs whose electron donors are Rieske-type ferredoxins. Given that these residues were not conserved if the electron donors were different types of ferredoxins or reductases, the side-wise site of the mushroom-like structure is thought to be the common binding site between Rieske-type ferredoxin and α3ß3-type oxygenase components in ROs. IMPORTANCE We clarified the critical amino acid residues of the oxygenase component (Oxy) of Rieske non-heme iron oxygenase (RO) for binding with Rieske-type ferredoxin (Fd). Our results showed that Rieske-type Fd-binding site is commonly located at the stem (side-wise site) of the mushroom-like α3ß3 quaternary structure in many ROs. The resultant binding site was totally different from those reported at the top-wise site of the doughnut-like α3-type Oxy, although α3-type Oxys correspond to the cap (α3 subunit part) of the mushroom-like α3ß3-type Oxys. Critical amino acid residues detected in this study were not conserved if the electron donors of Oxys were different types of Fds or reductases. Altogether, we can suggest that unique binding modes between Oxys and electron donors have evolved, depending on the nature of the electron donors, despite Oxy molecules having shared α3ß3 quaternary structures.

A unique peptide-based pharmacophore identifies an inhibitory compound against the A-subunit of Shiga toxin.

Shiga toxin (Stx), a major virulence factor of enterohemorrhagic Escherichia coli (EHEC), can cause fatal systemic complications. Recently, we identified a potent inhibitory peptide that binds to the catalytic A-subunit of Stx. Here, using biochemical structural analysis and X-ray crystallography, we determined a minimal essential peptide motif that occupies the catalytic cavity and is required for binding to the A-subunit of Stx2a, a highly virulent Stx subtype. Molecular dynamics simulations also identified the same motif and allowed determination of a unique pharmacophore for A-subunit binding. Notably, a series of synthetic peptides containing the motif efficiently inhibit Stx2a. In addition, pharmacophore screening and subsequent docking simulations ultimately identified nine Stx2a-interacting molecules out of a chemical compound database consisting of over 7,400,000 molecules. Critically, one of these molecules markedly inhibits Stx2a both in vitro and in vivo, clearly demonstrating the significance of the pharmacophore for identifying therapeutic agents against EHEC infection.

Mutagenesis analysis of GMN motif in Arabidopsis thaliana Mg2+ transporter MRS2-1.

Magnesium is an important nutrient for plants, but much is still unknown about plant Mg2+ transporters. Combining with the structural prediction of AlphaFold2, we used mutagenesis and 28Mg uptake assay to study the highly conserved "GMN" motif of Arabidopsis thaliana MRS2-1 (AtMRS2-1) transporter. We demonstrated that the glycine and methionine in GMN motif are essential for AtMRS2-1 to transport Mg2+.

Structural basis for the assembly and quinone transport mechanisms of the dimeric photosynthetic RC-LH1 supercomplex.

The reaction center (RC) and light-harvesting complex 1 (LH1) form a RC-LH1 core supercomplex that is vital for the primary reactions of photosynthesis in purple phototrophic bacteria. Some species possess the dimeric RC-LH1 complex with a transmembrane polypeptide PufX, representing the largest photosynthetic complex in anoxygenic phototrophs. However, the details of the architecture and assembly mechanism of the RC-LH1 dimer are unclear. Here we report seven cryo-electron microscopy (cryo-EM) structures of RC-LH1 supercomplexes from Rhodobacter sphaeroides. Our structures reveal that two PufX polypeptides are positioned in the center of the S-shaped RC-LH1 dimer, interlocking association between the components and mediating RC-LH1 dimerization. Moreover, we identify another transmembrane peptide, designated PufY, which is located between the RC and LH1 subunits near the LH1 opening. PufY binds a quinone molecule and prevents LH1 subunits from completely encircling the RC, creating a channel for quinone/quinol exchange. Genetic mutagenesis, cryo-EM structures, and computational simulations provide a mechanistic understanding of the assembly and electron transport pathways of the RC-LH1 dimer and elucidate the roles of individual components in ensuring the structural and functional integrity of the photosynthetic supercomplex.

Cryptic Oxidative Transamination of Hydroxynaphthoquinone in Natural Product Biosynthesis.

Pyridoxal 5'-phosphate (PLP)-dependent enzymes are a group of versatile enzymes that catalyze various reactions, but only a small number of them react with O2. Here, we report an unprecedented PLP-dependent enzyme, NphE, that catalyzes both transamination and two-electron oxidation using O2 as an oxidant. Our intensive analysis reveals that NphE transfers the l-glutamate-derived amine to 1,3,6,8-tetrahydroxynaphthalene-derived mompain to form 8-amino-flaviolin (8-AF) via a highly conjugated quinonoid intermediate that is reactive with O2. During the NphE reaction, O2 is reduced to yield H2O2. An integrated technique involving NphE structure prediction by AlphaFold v2.0 and molecular dynamics simulation suggested the O2-accessible cavity. Our in vivo results demonstrated that 8-AF is a genuine biosynthetic intermediate for the 1,3,6,8-tetrahydroxynaphthalene-derived meroterpenoid naphterpin without an amino group, which was supported by site-directed mutagenesis. This study clearly establishes the NphE reaction product 8-AF as a common intermediate with a cryptic amino group for the biosynthesis of terpenoid-polyketide hybrid natural products.

Structurally-discovered KLF4 variants accelerate and stabilize reprogramming to pluripotency.

Non-genetically modified somatic cells can only be inefficiently and stochastically reprogrammed to pluripotency by exogenous expression of reprogramming factors. Low competence of natural reprogramming factors may prevent the majority of cells to successfully and synchronously reprogram. Here we screened DNA-interacting amino acid residues in the zinc-finger domain of KLF4 for enhanced reprogramming efficiency using alanine-substitution scanning methods. Identified KLF4 L507A mutant accelerated and stabilized reprogramming to pluripotency in both mouse and human somatic cells. By testing all the variants of L507 position, variants with smaller amino acid residues in the KLF4 L507 position showed higher reprogramming efficiency. L507A bound more to promoters or enhancers of pluripotency genes, such as KLF5, and drove gene expression of these genes during reprogramming. Molecular dynamics simulations predicted that L507A formed additional interactions with DNA. Our study demonstrates how modifications in amino acid residues of DNA-binding domains enable next-generation reprogramming technology with engineered reprogramming factors.

Structures of human pannexin-1 in nanodiscs reveal gating mediated by dynamic movement of the N terminus and phospholipids.

Pannexin (PANX) family proteins form large-pore channels that mediate purinergic signaling. We analyzed the cryo-EM structures of human PANX1 in lipid nanodiscs to elucidate the gating mechanism and its regulation by the amino terminus in phospholipids. The wild-type channel has an amino-terminal funnel in the pore, but in the presence of the inhibitor probenecid, a cytoplasmically oriented amino terminus and phospholipids obstruct the pore. Functional analysis using whole-cell patch-clamp and oocyte voltage clamp showed that PANX1 lacking the amino terminus did not open and had a dominant negative effect on channel activity, thus confirming that the amino-terminal domain played an essential role in channel opening. These observations suggest that dynamic conformational changes in the amino terminus of human PANX1 are associated with lipid movement in and out of the pore. Moreover, the data provide insight into the gating mechanism of PANX1 and, more broadly, other large-pore channels.

Structural and Molecular Basis of the Catalytic Mechanism of Geranyl Pyrophosphate C6-Methyltransferase: Creation of an Unprecedented Farnesyl Pyrophosphate C6-Methyltransferase.

Prenyl pyrophosphate methyltransferases enhance the structural diversity of terpenoids. However, the molecular basis of their catalytic mechanisms is poorly understood. In this study, using multiple strategies, we characterized a geranyl pyrophosphate (GPP) C6-methyltransferase, BezA. Biochemical analysis revealed that BezA requires Mg2+ and solely methylates GPP. The crystal structures of BezA and its complex with S-adenosyl homocysteine were solved at 2.10 and 2.56 Å, respectively. Further analyses using site-directed mutagenesis, molecular docking, molecular dynamics simulations, and quantum mechanics/molecular mechanics calculations revealed the molecular basis of the methylation reaction. Importantly, the function of E170 as a catalytic base to complete the methylation reaction was established. We also succeeded in switching the substrate specificity by introducing a W210A substitution, resulting in an unprecedented farnesyl pyrophosphate C6-methyltransferase.