Structural insights into phosphatidylethanolamine N-methyltransferase PmtA mediating bacterial phosphatidylcholine synthesis.

Phosphatidylethanolamine N-methyltransferase (PmtA) catalyzes the biosynthesis of phosphatidylcholine (PC) from phosphatidylethanolamine (PE). Although PC is one of the major phospholipids constituting bilayer membranes in eukaryotes, certain bacterial species encode PmtA, a membrane-associated methyltransferase, to produce PC, which is correlated with cellular stress responses, adaptability to environmental changes, and symbiosis or virulence with eukaryotic hosts. Depending on the organism, multiple PmtAs may be required for producing monomethyl- and dimethyl-PE derivatives along with PC, whereas in organisms such as Rubellimicrobium thermophilum, a single enzyme is sufficient to direct all three methylation steps. In this study, we present the x-ray crystal structures of PmtA from R. thermophilum in complex with dimethyl-PE and S-adenosyl-l-homocysteine, as well as in its lipid-free form. Moreover, we demonstrate that the enzyme associates with the cellular membrane via electrostatic interactions facilitated by a group of critical basic residues and can successively methylate PE and its methylated derivatives, culminating in the production of PC.

Structural study of the intrinsically disordered tardigrade damage suppressor protein (Dsup) and its complex with DNA.

Studies of proteins, found in one of the most stress-resistant animals tardigrade Ramazzottius varieornatus, aim to reveal molecular principles of extreme tolerance to various types of stress and developing applications based on them for medicine, biotechnology, pharmacy, and space research. Tardigrade DNA/RNA-binding damage suppressor protein (Dsup) reduces DNA damage caused by reactive oxygen spices (ROS) produced upon irradiation and oxidative stresses in Dsup-expressing transgenic organisms. This work is focused on the determination of structural features of Dsup protein and Dsup-DNA complex, which refines details of protective mechanism. For the first time, intrinsically disordered nature of Dsup protein with highly flexible structure was experimentally proven and characterized by the combination of small angle X-ray scattering (SAXS) technique, circular dichroism spectroscopy, and computational methods. Low resolution models of Dsup protein and an ensemble of conformations were presented. In addition, we have shown that Dsup forms fuzzy complex with DNA.

Structural studies of ß-glucosidase from the thermophilic bacterium Caldicellulosiruptor saccharolyticus.

ß-Glucosidase from the thermophilic bacterium Caldicellulosiruptor saccharolyticus (Bgl1) has been denoted as having an attractive catalytic profile for various industrial applications. Bgl1 catalyses the final step of in the decomposition of cellulose, an unbranched glucose polymer that has attracted the attention of researchers in recent years as it is the most abundant renewable source of reduced carbon in the biosphere. With the aim of enhancing the thermostability of Bgl1 for a broad spectrum of biotechnological processes, it has been subjected to structural studies. Crystal structures of Bgl1 and its complex with glucose were determined at 1.47 and 1.95 Šresolution, respectively. Bgl1 is a member of glycosyl hydrolase family 1 (GH1 superfamily, EC 3.2.1.21) and the results showed that the 3D structure of Bgl1 follows the overall architecture of the GH1 family, with a classical (ß/α)8 TIM-barrel fold. Comparisons of Bgl1 with sequence or structural homologues of ß-glucosidase reveal quite similar structures but also unique structural features in Bgl1 with plausible functional roles.

The cAMP receptor protein from Gardnerella vaginalis is not regulated by ligands.

Overgrowth of Gardnerella vaginalis causes an imbalance in vaginal microecology. The pathogenicity of G. vaginalis is directly regulated by the cAMP receptor protein (CRP). In this study, we resolve the crystal structure of CRPGv at a resolution of 2.22 Å and find some significant differences from homologous proteins. The first 23 amino acids of CRPGv are inserted into the ligand binding pocket, creating a strong steric barrier to ligand entry that has not been seen previously in its homologues. In the absence of ligands, the two α helices used by CRPGv to bind oligonucleotide chains are exposed and can specifically bind TGTGA-N6-TCACA sequences. cAMP and other ligands of CRP homologs are not cofactors of CRPGv. There is no coding gene of the adenylate cyclase, and cAMP could not be identified in G. vaginalis by liquid chromatography tandem mass spectrometry. We speculate that CRPGv may achieve fine regulation through a conformational transformation different from that of its homologous proteins, and this conformational transformation is no longer dependent on small molecules, but may be aided by accessory proteins. CRPGv is the first discovered CRP that is not ligand-regulated, and its active conformation provides a structural basis for drug screening.

Structural characterization of green fluorescent protein in the I-state.

Green fluorescent protein (GFP) is widely utilized as a fluorescent tag in biochemical fields. Whereas the intermediate (I) state has been proposed in the photoreaction cycle in addition to the A and B states, until now the structure of I has only been estimated by computational studies. In this paper, we report the crystal structures of the I stabilizing variants of GFP at high resolutions where respective atoms can be observed separately. Comparison with the structures in the other states highlights the structural feature of the I state. The side chain of one of the substituted residues, Val203, adopts the gauche- conformation observed for Thr203 in the A state, which is different from the B state. On the other hand, His148 interacts with the chromophore by ordinary hydrogen bonding with a distance of 2.85 Å, while the weaker interaction by longer distances is observed in the A state. Therefore, it was indicated that it is possible to distinguish three states A, B and I by the two hydrogen bond distances Oγ-Thr203···Oη-chromophore and Nδ1-His148···Oη-chromophore. We discuss the characteristics of the I intermediate of wild-type GFP on the bases of the structure estimated from the variant structures by quantum chemical calculations.

Structure and dynamics of cholesterol-mediated aquaporin-0 arrays and implications for lipid rafts.

Aquaporin-0 (AQP0) tetramers form square arrays in lens membranes through a yet unknown mechanism, but lens membranes are enriched in sphingomyelin and cholesterol. Here, we determined electron crystallographic structures of AQP0 in sphingomyelin/cholesterol membranes and performed molecular dynamics (MD) simulations to establish that the observed cholesterol positions represent those seen around an isolated AQP0 tetramer and that the AQP0 tetramer largely defines the location and orientation of most of its associated cholesterol molecules. At a high concentration, cholesterol increases the hydrophobic thickness of the annular lipid shell around AQP0 tetramers, which may thus cluster to mitigate the resulting hydrophobic mismatch. Moreover, neighboring AQP0 tetramers sandwich a cholesterol deep in the center of the membrane. MD simulations show that the association of two AQP0 tetramers is necessary to maintain the deep cholesterol in its position and that the deep cholesterol increases the force required to laterally detach two AQP0 tetramers, not only due to protein-protein contacts but also due to increased lipid-protein complementarity. Since each tetramer interacts with four such 'glue' cholesterols, avidity effects may stabilize larger arrays. The principles proposed to drive AQP0 array formation could also underlie protein clustering in lipid rafts.

AI-Predicted Protein Deformation Encodes Energy Landscape Perturbation.

AI algorithms have proven to be excellent predictors of protein structure, but whether and how much these algorithms can capture the underlying physics remains an open question. Here, we aim to test this question using the Alphafold2 (AF) algorithm: We use AF to predict the subtle structural deformation induced by single mutations, quantified by strain, and compare with experimental datasets of corresponding perturbations in folding free energy ΔΔG. Unexpectedly, we find that physical strain alone-without any additional data or computation-correlates almost as well with ΔΔG as state-of-the-art energy-based and machine-learning predictors. This indicates that the AF-predicted structures alone encode fine details about the energy landscape. In particular, the structures encode significant information on stability, enough to estimate (de-)stabilizing effects of mutations, thus paving the way for the development of novel, structure-based stability predictors for protein design and evolution.

Phalloidin and DNase I-bound F-actin pointed end structures reveal principles of filament stabilization and disassembly.

Actin filament turnover involves subunits binding to and dissociating from the filament ends, with the pointed end being the primary site of filament disassembly. Several molecules modulate filament turnover, but the underlying mechanisms remain incompletely understood. Here, we present three cryo-EM structures of the F-actin pointed end in the presence and absence of phalloidin or DNase I. The two terminal subunits at the undecorated pointed end adopt a twisted conformation. Phalloidin can still bind and bridge these subunits, inducing a conformational shift to a flattened, F-actin-like state. This explains how phalloidin prevents depolymerization at the pointed end. Interestingly, two DNase I molecules simultaneously bind to the phalloidin-stabilized pointed end. In the absence of phalloidin, DNase I binding would disrupt the terminal actin subunit packing, resulting in filament disassembly. Our findings uncover molecular principles of pointed end regulation and provide structural insights into the kinetic asymmetry between the actin filament ends.

Revealing the pH-dependent conformational changes in sol g 2.1 protein and potential ligands binding.

Sol g 2, a major protein found in the venom of the tropical fire ant (Solenopsis geminata), is well-known for its ability to bind various hydrophobic molecules. In this study, we investigate the binding activity of recombinant Sol g 2.1 protein (rSol g 2.1) with potential molecules, including (E)-ß-Farnesene, α-Caryophyllene, and 1-Octen-3-ol at different pH levels (pH 7.4 and 5.5) using fluorescence competitive binding assays (FCBA). Our results revealed that Sol g 2.1 protein has higher affinity binding with these ligands at neutral pH. Relevance to molecular docking and molecular dynamics simulations were utilized to provide insights into the stability and conformational dynamics of Sol g 2.1 and its ligand complexes. After simulation, we found that Sol g 2.1 protein has higher affinity binding with these ligands as well as high structural stability at pH 7.4 than at an acidic pH level, indicating by RMSD, RMSF, Rg, SASA, and principal component analysis (PCA). Additionally, the Sol g 2.1 protein complexes at pH 7.4 showed significantly lower binding free energy (∆Gbind) and higher total residue contributions, particularly from key non-polar amino acids such as Trp36, Met40, Cys62, and Ile104, compared to the lower pH environment. These explain why they exhibited higher binding affinity than the lower pH. Therefore, we suggested that Sol g 2.1 protein is a pH-responsive carrier protein. These findings also expand our understanding of protein-ligand interactions and offer potential avenues for the development of innovative drug delivery strategies targeting Sol g 2.1 protein.

Modelling protein complexes with crosslinking mass spectrometry and deep learning.

Scarcity of structural and evolutionary information on protein complexes poses a challenge to deep learning-based structure modelling. We integrate experimental distance restraints obtained by crosslinking mass spectrometry (MS) into AlphaFold-Multimer, by extending AlphaLink to protein complexes. Integrating crosslinking MS data substantially improves modelling performance on challenging targets, by helping to identify interfaces, focusing sampling, and improving model selection. This extends to single crosslinks from whole-cell crosslinking MS, opening the possibility of whole-cell structural investigations driven by experimental data. We demonstrate this by revealing the molecular basis of iron homoeostasis in Bacillus subtilis.

Structural insight into the DNMT1 reaction cycle by cryo-electron microscopy.

DNMT1 is an essential DNA methyltransferase that catalyzes the transfer of methyl groups to CpG islands in DNA and generates a prominent epigenetic mark. The catalytic activity of DNMT1 relies on its conformational plasticity and ability to change conformation from an auto-inhibited to an activated state. Here, we present four cryo-EM reconstructions of apo DNMT1 and DNTM1: non-productive DNA, DNTM1: H3Ub2-peptide, DNTM1: productive DNA complexes. Our structures demonstrate the flexibility of DNMT1's N-terminal regulatory domains during the transition from an apo 'auto-inhibited' to a DNA-bound 'non-productive' and finally a DNA-bound 'productive' state of DNMT1. Furthermore, we address the regulation of DNMT1's methyltransferase activity by a DNMT1-selective small-molecule inhibitor and ubiquitinated histone H3. We observe that DNMT1 binds DNA in a 'non-productive' state despite the presence of the inhibitor and present the cryo-EM reconstruction of full-length DNMT1 in complex with a di-ubiquitinated H3 peptide analogue. Taken together, our results provide structural insights into the reaction cycle of DNMT1.

Structural and functional effects of phosphopriming and scaffolding in the kinase GSK-3ß.

Glycogen synthase kinase 3ß (GSK-3ß) targets specific signaling pathways in response to distinct upstream signals. We used structural and functional studies to dissect how an upstream phosphorylation step primes the Wnt signaling component ß-catenin for phosphorylation by GSK-3ß and how scaffolding interactions contribute to this reaction. Our crystal structure of GSK-3ß bound to a phosphoprimed ß-catenin peptide confirmed the expected binding mode of the phosphoprimed residue adjacent to the catalytic site. An aspartate phosphomimic in the priming site of ß-catenin adopted an indistinguishable structure but reacted approximately 1000-fold slower than the native phosphoprimed substrate. This result suggests that substrate positioning alone is not sufficient for catalysis and that native phosphopriming interactions are necessary. We also obtained a structure of GSK-3ß with an extended peptide from the scaffold protein Axin that bound with greater affinity than that of previously crystallized Axin fragments. This structure neither revealed additional contacts that produce the higher affinity nor explained how substrate interactions in the GSK-3ß active site are modulated by remote Axin binding. Together, our findings suggest that phosphopriming and scaffolding produce small conformational changes or allosteric effects, not captured in the crystal structures, that activate GSK-3ß and facilitate ß-catenin phosphorylation. These results highlight limitations in our ability to predict catalytic activity from structure and have potential implications for the role of natural phosphomimic mutations in kinase regulation and phosphosite evolution.

Structures of the Mycobacterium tuberculosis efflux pump EfpA reveal the mechanisms of transport and inhibition.

As the first identified multidrug efflux pump in Mycobacterium tuberculosis (Mtb), EfpA is an essential protein and promising drug target. However, the functional and inhibitory mechanisms of EfpA are poorly understood. Here we report cryo-EM structures of EfpA in outward-open conformation, either bound to three endogenous lipids or the inhibitor BRD-8000.3. Three lipids inside EfpA span from the inner leaflet to the outer leaflet of the membrane. BRD-8000.3 occupies one lipid site at the level of inner membrane leaflet, competitively inhibiting lipid binding. EfpA resembles the related lysophospholipid transporter MFSD2A in both overall structure and lipid binding sites and may function as a lipid flippase. Combining AlphaFold-predicted EfpA structure, which is inward-open, we propose a complete conformational transition cycle for EfpA. Together, our results provide a structural and mechanistic foundation to comprehend EfpA function and develop EfpA-targeting anti-TB drugs.

Teaching old docks new tricks with machine learning enhanced ensemble docking.

We here introduce Ensemble Optimizer (EnOpt), a machine-learning tool to improve the accuracy and interpretability of ensemble virtual screening (VS). Ensemble VS is an established method for predicting protein/small-molecule (ligand) binding. Unlike traditional VS, which focuses on a single protein conformation, ensemble VS better accounts for protein flexibility by predicting binding to multiple protein conformations. Each compound is thus associated with a spectrum of scores (one score per protein conformation) rather than a single score. To effectively rank and prioritize the molecules for further evaluation (including experimental testing), researchers must select which protein conformations to consider and how best to map each compound's spectrum of scores to a single value, decisions that are system-specific. EnOpt uses machine learning to address these challenges. We perform benchmark VS to show that for many systems, EnOpt ranking distinguishes active compounds from inactive or decoy molecules more effectively than traditional ensemble VS methods. To encourage broad adoption, we release EnOpt free of charge under the terms of the MIT license.

Empowering AlphaFold2 for protein conformation selective drug discovery with AlphaFold2-RAVE.

Small-molecule drug design hinges on obtaining co-crystallized ligand-protein structures. Despite AlphaFold2's strides in protein native structure prediction, its focus on apo structures overlooks ligands and associated holo structures. Moreover, designing selective drugs often benefits from the targeting of diverse metastable conformations. Therefore, direct application of AlphaFold2 models in virtual screening and drug discovery remains tentative. Here, we demonstrate an AlphaFold2-based framework combined with all-atom enhanced sampling molecular dynamics and Induced Fit docking, named AF2RAVE-Glide, to conduct computational model-based small-molecule binding of metastable protein kinase conformations, initiated from protein sequences. We demonstrate the AF2RAVE-Glide workflow on three different mammalian protein kinases and their type I and II inhibitors, with special emphasis on binding of known type II kinase inhibitors which target the metastable classical DFG-out state. These states are not easy to sample from AlphaFold2. Here, we demonstrate how with AF2RAVE these metastable conformations can be sampled for different kinases with high enough accuracy to enable subsequent docking of known type II kinase inhibitors with more than 50% success rates across docking calculations. We believe the protocol should be deployable for other kinases and more proteins generally.

​Fusarium Protein Toolkit: a web-based resource for structural and variant analysis of Fusarium species.

BACKGROUND: ​​The genus Fusarium poses significant threats to food security and safety worldwide because numerous species of the fungus cause destructive diseases and/or mycotoxin contamination in crops. The adverse effects of climate change are exacerbating some existing threats and causing new problems. These challenges highlight the need for innovative solutions, including the development of advanced tools to identify targets for control strategies. DESCRIPTION: In response to these challenges, we developed the Fusarium Protein Toolkit (FPT), a web-based tool that allows users to interrogate the structural and variant landscape within the Fusarium pan-genome. The tool displays both AlphaFold and ESMFold-generated protein structure models from six Fusarium species. The structures are accessible through a user-friendly web portal and facilitate comparative analysis, functional annotation inference, and identification of related protein structures. Using a protein language model, FPT predicts the impact of over 270 million coding variants in two of the most agriculturally important species, Fusarium graminearum and F. verticillioides. To facilitate the assessment of naturally occurring genetic variation, FPT provides variant effect scores for proteins in a Fusarium pan-genome based on 22 diverse species. The scores indicate potential functional consequences of amino acid substitutions and are displayed as intuitive heatmaps using the PanEffect framework. CONCLUSION: FPT fills a knowledge gap by providing previously unavailable tools to assess structural and missense variation in proteins produced by Fusarium. FPT has the potential to deepen our understanding of pathogenic mechanisms in Fusarium, and aid the identification of genetic targets for control strategies that reduce crop diseases and mycotoxin contamination. Such targets are vital to solving the agricultural problems incited by Fusarium, particularly evolving threats resulting from climate change. Thus, FPT has the potential to contribute to improving food security and safety worldwide.

Structure and mechanism of a phosphotransferase system glucose transporter.

Glucose is the primary source of energy for many organisms and is efficiently taken up by bacteria through a dedicated transport system that exhibits high specificity. In Escherichia coli, the glucose-specific transporter IICBGlc serves as the major glucose transporter and functions as a component of the phosphoenolpyruvate-dependent phosphotransferase system. Here, we report cryo-electron microscopy (cryo-EM) structures of the glucose-bound IICBGlc protein. The dimeric transporter embedded in lipid nanodiscs was captured in the occluded, inward- and occluded, outward-facing conformations. Together with biochemical and biophysical analyses, and molecular dynamics (MD) simulations, we provide insights into the molecular basis and dynamics for substrate recognition and binding, including the gates regulating the binding sites and their accessibility. By combination of these findings, we present a mechanism for glucose transport across the plasma membrane. Overall, this work provides molecular insights into the structure, dynamics, and mechanism of the IICBGlc transporter in a native-like lipid environment.

ModFOLD9: A Web Server for Independent Estimates of 3D Protein Model Quality.

Accurate models of protein tertiary structures are now available from numerous advanced prediction methods, although the accuracy of each method often varies depending on the specific protein target. Additionally, many models may still contain significant local errors. Therefore, reliable, independent model quality estimates are essential both for identifying errors and selecting the very best models for further biological investigations. ModFOLD9 is a leading independent server for detecting the local errors in models produced by any method, and it can accurately discriminate between high-quality models from multiple alternative approaches. ModFOLD9 incorporates several new scores from deep learning-based approaches, leading to greatly improved prediction accuracy compared with earlier versions of the server. ModFOLD9 is continuously independently benchmarked, and it is shown to be highly competitive with other public servers. ModFOLD9 is freely available at https://www.reading.ac.uk/bioinf/ModFOLD/.

flDPnn2: Accurate and Fast Predictor of Intrinsic Disorder in Proteins.

Prediction of the intrinsic disorder in protein sequences is an active research area, with well over 100 predictors that were released to date. These efforts are motivated by the functional importance and high levels of abundance of intrinsic disorder, combined with relatively low amounts of experimental annotations. The disorder predictors are periodically evaluated by independent assessors in the Critical Assessment of protein Intrinsic Disorder prediction (CAID) experiments. The recently completed CAID2 experiment assessed close to 40 state-of-the-art methods demonstrating that some of them produce accurate results. In particular, flDPnn2 method, which is the successor of flDPnn that performed well in the CAID1 experiment, secured the overall most accurate results on the Disorder-NOX dataset in CAID2. flDPnn2 implements a number of improvements when compared to its predecessor including changes to the inputs, increased size of the deep network model that we retrained on a larger training set, and addition of an alignment module. Using results from CAID2, we show that flDPnn2 produces accurate predictions very quickly, modestly improving over the accuracy of flDPnn and reducing the runtime by half, to about 27 s per protein. flDPnn2 is freely available as a convenient web server at http://biomine.cs.vcu.edu/servers/flDPnn2/.

CAPRI-Q: The CAPRI resource evaluating the quality of predicted structures of protein complexes.

Protein interactions are essential for cellular processes. In recent years there has been significant progress in computational prediction of 3D structures of individual protein chains, with the best-performing algorithms reaching sub-Ångström accuracy. These techniques are now finding their way into the prediction of protein interactions, adding to the existing modeling approaches. The community-wide Critical Assessment of Predicted Interactions (CAPRI) has been a catalyst for the development of procedures for the structural modeling of protein assemblies by organizing blind prediction experiments. The predicted structures are assessed against unpublished experimentally determined structures using a set of metrics with proven robustness that have been established in the CAPRI community. In addition, several advanced benchmarking databases provide targets against which users can test docking and assembly modeling software. These include the Protein-Protein Docking Benchmark, the CAPRI Scoreset, and the Dockground database, all developed by members of the CAPRI community. Here we present CAPRI-Q, a stand-alone model quality assessment tool, which can be freely downloaded or used via a publicly available web server. This tool applies the CAPRI metrics to assess the quality of query structures against given target structures, along with other popular quality metrics such as DockQ, TM-score and l-DDT, and classifies the models according to the CAPRI model quality criteria. The tool can handle a variety of protein complex types including those involving peptides, nucleic acids, and oligosaccharides. The source code is freely available from https://gitlab.in2p3.fr/cmsb-public/CAPRI-Q and its web interface through the Dockground resource at https://dockground.compbio.ku.edu/assessment/.