Electrostatics of nanosystems: application to microtubules and the ribosome

A lineage-specific protein network at the trypanosome nuclear envelope

The nuclear envelope (NE) separates translation and transcription and is the location of multiple functions, including chromatin organization and nucleocytoplasmic transport. The molecular basis for many of these functions have diverged between eukaryotic lineages. Trypanosoma brucei, a member of the early branching eukaryotic lineage Discoba, highlights many of these, including a distinct lamina and kinetochore composition. Here, we describe a cohort of proteins interacting with both the lamina and NPC, which we term lamina-associated proteins (LAPs). LAPs represent a diverse group of proteins, including two candidate NPC-anchoring pore membrane proteins (POMs) with architecture conserved with S. cerevisiae and H. sapiens, and additional peripheral components of the NPC. While many of the LAPs are Kinetoplastid specific, we also identified broadly conserved proteins, indicating an amalgam of divergence and conservation within the trypanosome NE proteome, highlighting the diversity of nuclear biology across the eukaryotes, increasing our understanding of eukaryotic and NPC evolution.

WUREN: Whole-modal union representation for epitope prediction

B-cell epitope identification plays a vital role in the development of vaccines, therapies, and diagnostic tools. Currently, molecular docking tools in B-cell epitope prediction are heavily influenced by empirical parameters and require significant computational resources, rendering a great challenge to meet large-scale prediction demands. When predicting epitopes from antigen-antibody complex, current artificial intelligence algorithms cannot accurately implement the prediction due to insufficient protein feature representations, indicating novel algorithm is desperately needed for efficient protein information extraction. In this paper, we introduce a multimodal model called WUREN (Whole-modal Union Representation for Epitope predictioN), which effectively combines sequence, graph, and structural features. It achieved AUC-PR scores of 0.213 and 0.193 on the solved structures and AlphaFold-generated structures, respectively, for the independent test proteins selected from DiscoTope3 benchmark. Our findings indicate that WUREN is an efficient feature extraction model for protein complexes, with the generalizable application potential in the development of protein-based drugs. Moreover, the streamlined framework of WUREN could be readily extended to model similar biomolecules, such as nucleic acids, carbohydrates, and lipids.

A map of glycation and glycoxidation sites in collagen I of human cortical bone: Effects of sex and type 2 diabetes

Complications of diabetes is a major health problem affecting multiple organs including bone, where the chronic disease increases the risk of fragility fractures. One hypothesis suggests a pathogenic role for hyperglycemia-induced modification of proteins, a.k.a. advanced glycation end products (AGEs), resulting in structural and functional damage to bone extracellular matrix (ECM). Evidence supporting this hypothesis has been limited by the lack of comprehensive information about the location of AGEs that accumulate in vivo at specific sites within the proteins of bone ECM. Analyzing extracts from cortical bone of cadaveric femurs by liquid chromatography tandem mass spectrometry, we generated a quantitative AGE map of human collagen I for male and female adult donors with and without diabetes. The map describes the chemical nature, sequence position, and levels of four major physiological AGEs, e.g. carboxymethyllysine, and an AGE precursor fructosyllysine within the collagen I triple-helical region. The important features of the map are: 1) high map reproducibility in the individual bone extracts, i.e. 20 male and 20 female donors; 2) localization of modifications to distinct clusters: 10 clusters containing 34 AGE sites in male donors and 9 clusters containing 28 sites in female donors; 3) significant increases in modification levels in diabetes at multiple sites: 26 out of 34 sites in males and in 17 out of 28 sites in females; and 4) generally higher modification levels in male vs. female donors. Moreover, the AGE levels at multiple individual sites correlated with total bone pentosidine levels in male but not in female donors. Molecular dynamics simulations and molecular modeling predicted significant impact of modifications on solvent exposure, charge distribution, and hydrophobicity of the triple helix as well as disruptions to the structure of collagen I fibril. In summary, the AGE map of collagen I revealed diabetes-induced, sex-specific non-enzymatic modifications at distinct triple helical sites that can disrupt collagen structure, thus proposing a specific mechanism of AGE contribution to diabetic complications in human bone.

Structure-based virtual screening of natural compounds against wild and mutant (R1155K, A1156T and D1168A) NS3-4A protease of Hepatitis C virus

NS3-4A, a serine protease, is a primary target for drug development against Hepatitis C Virus (HCV). However, the effectiveness of potent next-generation protease inhibitors is limited by the emergence of mutations and resulting drug resistance. To address this, in this study a structure-based drug design approach is employed to screen a large library of 7320 natural compounds against both wild-type and mutant variants of NS3-4A protease. Telaprevir, a widely used protease inhibitor, was recruited as the control drug. The top 10 compounds with favorable binding affinities underwent drug-likeness evaluation. Based on ADMET studies, complexes of NP_024762 and NP_006776 were selected for molecular dynamic simulations. Principal component analysis (PCA) was employed to explore the conformational space and protein dynamics of the protein-ligand complex using a Free Energy Landscape (FEL) approach. The cosine values obtained from FEL analysis ranged from 0 to 1, and eigenvectors with cosine values below 0.2 were chosen for further analysis. To forecast binding free energies and evaluate energy contributions per residue, the MM-PBSA method was employed. The results highlighted the crucial role of amino acids in the catalytic domain for the binding of the protease with phytochemicals. Stable associations between the top compounds and the target protease were confirmed by the formation of hydrogen bonds in the binding pocket involving residues: His1057, Gly1137, Ser1139, and Ala1157. These findings suggest the potential of these compounds for further validation through biological evaluation.Communicated by Ramaswamy H. Sarma.

Polyphasic approaches to identify and understand α-glucosidase inhibitory potential of secondary metabolites of Withania coagulans fruit

Withania coagulans (WC) is used in traditional and Ayurveda medicine to treat various ailments, including diabetes. Our investigation found that WC fruit hexane extract effectively suppresses α-glucosidase activity (IC50 = 0.013 mg/ml, Ki = 0.012 mg/ml). The purified molecule has an IC50 of 0.004 mg/ml and Ki of 0.0037 mg/ml. FTIR examination indicates distinctive peaks at 3500, 2900, 1770, and 1500 cm-1 corresponding to functional groups OH bending, CH stretching, CO stretching, and CO stretching. GCMS analysis reveals plant secondary metabolites (PSM) such as n-hexadecenoic acid and methyl 9,10-octadecadienoate. NMR confirms the existence of olefinic fatty acids. The bioactive fraction recorded a non-competitive mode of inhibition of α-glucosidase activity. The cytotoxicity exhibited against HELA cell was IC50 0.4 mg/ml and found positive in inhibiting the growth of Bacillus cereus, Escherichia coli, Klebsiella pneumoniae, Staphylococcus aureus, and Pseudomonas aeruginosa. Additionally, ensemble docking and molecular dynamic simulation showed that, out of the four PSMs examined, methyl 12,13-tetradecadienoate interacted with the α-glucosidase enzyme's allosteric site (BE -128.78 kJ/mol) and changed the configurations of the catalytic sites, as demonstrated by the enzyme's decreased affinity for isomaltose. The study found that PSMs from WC fruit may inhibit α-glucosidase, making them viable candidates for antidiabetic medication development.

Exploring high hydrostatic pressure effects on anthocyanin binding to serum albumin and food-derived transferrins

This study investigated the effects of high hydrostatic pressure (HHP) on the binding interactions of cyanindin-3-O-glucoside (C3G) to bovine serum albumin, human serum albumin (HSA), bovine lactoferrin, and ovotransferrin. Fluorescence quenching revealed that HHP reduced C3G-binding affinity to HSA, while having a largely unaffected role for the other proteins. Notably, pretreating HSA at 500 MPa significantly increased its dissociation constant with C3G from 24.7 to 34.3 μM. Spectroscopic techniques suggested that HSA underwent relatively pronounced tertiary structural alterations after HHP treatments. The C3G-HSA binding mechanisms under pressure were further analyzed through molecular dynamics simulation. The localized structural changes in HSA under pressure might weaken its interaction with C3G, particularly polar interactions such as hydrogen bonds and electrostatic forces, consequently leading to a decreased binding affinity. Overall, the importance of pressure-induced structural alterations in proteins influencing their binding with anthocyanins was highlighted, contributing to optimizing HHP processing for anthocyanin-based products.

Biochemical and structural insights of a recombinant AA16 LPMO from the marine and sponge-symbiont Peniophora sp

Lytic polysaccharide monooxygenases (LPMOs) are copper-dependent enzymes that oxidize polysaccharides, leading to their cleavage. LPMOs are classified into eight CAZy families (AA9-11, AA13-17), with the functionality of AA16 being poorly characterized. This study presents biochemical and structural data for an AA16 LPMO (PnAA16) from the marine sponge symbiont Peniophora sp. Phylogenetic analysis revealed that PnAA16 clusters separately from previously characterized AA16s. However, the structural modelling of PnAA16 showed the characteristic immunoglobulin-like fold of LPMOs, with a conserved his-brace motif coordinating a copper ion. The copper-bound PnAA16 showed greater thermal stability than its apo-form, highlighting copper's role in enzyme stability. Functionally, PnAA16 demonstrated oxidase activity, producing 5 μM H₂O₂ after 30 min, but showed 20 times lower peroxidase activity (0.27 U/g) compared to a fungal AA9. Specific activity assays indicated that PnAA16 acts only on cellohexaose, generating native celloligosaccharides (C3 to C5) and oxidized products with regioselective oxidation at C1 and C4 positions. Finally, PnAA16 boosted the activity of a cellulolytic cocktail for cellulose saccharification in the presence of ascorbic acid, hydrogen peroxide, or both. In conclusion, the present work provides insights into the AA16 family, expanding the understanding of their structural and functional relationships and biotechnological potential.

Design and synthesis of novel main protease inhibitors of COVID-19: quinoxalino[2,1-b]quinazolin-12-ones

The COVID-19 pandemic represents a substantial global challenge, being a significant cause of mortality in numerous countries. Thus, it is imperative to conduct research to develop effective therapies to combat COVID-19. The primary aim of this study is to employ a two-step tandem reaction involving 2,3-dichloroquinoxaline and 2-amino-N-substituted benzamides in alkaline media/DMF at an elevated temperature to design and synthesize a series of polycyclic derivatives endowed with quinoxalino[2,1-b]quinazolin-12-one framework. Following synthesis, the newly synthesized heterocycles were evaluated for their potential as inhibitors of the main protease of SARS-CoV-2 by means of molecular docking and dynamic simulation techniques. The in silico investigation demonstrated that all tested compounds effectively establish stable binding interactions, primarily through multiple hydrogen bonding and hydrophobic interactions, at the active site of the enzyme. These findings offer crucial structural insights that can be employed in future endeavors toward designing potent inhibitors targeting the main protease (Mpro). Among the investigated compounds, the p-tolylamino-substituted quinoxalino[2,1-b]quinazolinone derivative exhibited the most promise as an inhibitor of the main protease in COVID-19. Consequently, it warrants further investigation both in vitro and in vivo to identify it as a prospective candidate for anti-SARS-CoV-2 drug development.

Molecular Dynamics Reveals Altered Interactions between Belzutifan and HIF-2 with Natural Variant G323E or Proximal Phosphorylation at T324

In patients with von-Hippel Lindau (VHL) disease, hypoxia-independent accumulation of HIF-2α leads to increased transcriptional activity of HIF-2α:ARNT that drives cancers such as renal cell carcinoma. Belzutifan, a recently FDA-approved drug, is designed to prevent the transcriptional activity of HIF-2α:ARNT, thereby overcoming the consequences of its unnatural accumulation in VHL-dependent cancers. Emerging evidence suggests that the naturally occurring variant G323E located in the HIF-2α drug binding pocket prevents inhibitory activity of belzutifan analogs, though the mechanism of inhibition remains unclear. Interestingly, proximal phosphorylation at neighboring T324, previously shown to regulate HIF-2 protein interactions, has also been proposed to affect HIF-2 drug binding. Here, we used molecular dynamics (MD) simulations to understand and compare the molecular-level effects of G323E and phospho-T324 (pT324) on the belzutifan bound-HIF-2α:ARNT complex. We find that both G323E and pT324 increase structural flexibility within the drug binding site and reduce the apparent binding affinity for belzutifan. Whereas the effects of G323E are concentrated in the binding pocket Fα helix within the HIF-2α PAS-B domain, pT324 decreased the belzutifan binding affinity and stabilized the HIF-2 heterodimer through an alternate mechanism involving polar interactions between the HIF-2α PAS-B and PAS-A domains. Further analysis via ensemble machine learning uncovered important and distinct interchain residue interactions modified by G323E and pT324. These findings reveal a molecular mechanism of G323E-induced drug resistance and suggest that pT324 may also affect the efficacy of HIF-2 drug binding interactions via allosteric effects.

PLK1-mediated phosphorylation cascade activates Mis18 complex to ensure centromere inheritance

Accurate chromosome segregation requires the attachment of microtubules to centromeres, epigenetically defined by the enrichment of CENP-A nucleosomes. During DNA replication, CENP-A nucleosomes undergo dilution. To preserve centromere identity, correct amounts of CENP-A must be restored in a cell cycle-controlled manner orchestrated by the Mis18 complex (Mis18α-Mis18β-Mis18BP1). We demonstrate here that PLK1 interacts with the Mis18 complex by recognizing self-primed phosphorylations of Mis18α (Ser54) and Mis18BP1 (Thr78 and Ser93) through its Polo-box domain. Disrupting these phosphorylations perturbed both centromere recruitment of the CENP-A chaperone HJURP and new CENP-A loading. Biochemical and functional analyses showed that phosphorylation of Mis18α and PLK1 binding were required to activate Mis18α-Mis18β and promote Mis18 complex-HJURP interaction. Thus, our study reveals key molecular events underpinning the licensing role of PLK1 in ensuring accurate centromere inheritance.

Identification of novel broad-spectrum antiviral drugs targeting the N-terminal domain of the FIPV nucleocapsid protein

Coronaviruses pose serious threats to human and animal health worldwide, of which their structural nucleocapsid (N) proteins play multiple key roles in viral replication. However, the structures of animal coronavirus N proteins are poorly understood, posing challenges for research on their functions and pathogenic mechanisms as well as the development of N protein-based antiviral drugs. Therefore, N proteins must be further explored as potential antiviral targets. We determined the structure of the NNTD of feline infectious peritonitis virus (FIPV) and identified 3,6-dihydroxyflavone (3,6- DHF) as an effective N protein inhibitor. 3,6-DHF successfully inhibited FIPV replication in CRFK cells, showing broad-spectrum activity and effectiveness against drugresistant strains. Our study provides important insights for developing novel broadspectrum anti-coronavirus drugs and treating infections caused by drug-resistant mutant strains.

Poisson-Boltzmann-based machine learning model for electrostatic analysis

Electrostatics is of paramount importance to chemistry, physics, biology, and medicine. The Poisson-Boltzmann (PB) theory is a primary model for electrostatic analysis. However, it is highly challenging to compute accurate PB electrostatic solvation free energies for macromolecules due to the nonlinearity, dielectric jumps, charge singularity, and geometric complexity associated with the PB equation. The present work introduces a PB-based machine learning (PBML) model for biomolecular electrostatic analysis. Trained with the second-order accurate MIBPB solver, the proposed PBML model is found to be more accurate and faster than several eminent PB solvers in electrostatic analysis. The proposed PBML model can provide highly accurate PB electrostatic solvation free energy of new biomolecules or new conformations generated by molecular dynamics with much reduced computational cost.

Structure-based drug designing against Leishmania donovani using docking and molecular dynamics simulation studies: exploring glutathione synthetase as a drug target

In the pursuit of developing novel anti-leishmanial agents, we conducted an extensive computational study to screen inhibitors from the FDA-approved ZINC database against Leishmania donovani glutathione synthetase. The three-dimensional structure of Leishmania donovani glutathione synthetase was constructed by homology modeling, using the crystallographic structure of Trypanosoma brucei glutathione synthetase as a template. Subsequently, molecular docking studies were carried out for a large number of compounds using AutoDock Vina. Among the screened compounds, we selected the top five with strong binding affinity to Leishmania donovani glutathione synthetase but having a very low affinity to its human homolog. Further investigations on protein-ligand complexes were done by conducting molecular dynamics (MD) simulation and MM/PBSA analysis. The results revealed that Olysio (Simeprevir) exhibited the lowest binding energy (-89.21 kcal/mol), followed by Telithromycin (-45.34 kcal/mol). These findings showed that these compounds have the potential to act as inhibitors of glutathione synthetase. Hence, our study provides valuable insights for the development of a novel therapeutic strategy against Leishmania donovani by targeting the glutathione synthetase enzyme.Communicated by Ramaswamy H. Sarma.

Identification of candidate target genes of oral squamous cell carcinoma using high-throughput RNA-Seq data and in silico studies of their interaction with naturally occurring bioactive compounds

Oral Squamous Cell Carcinoma (OSCC) accounts for more than 90% of all kinds of oral neoplasms that develop in the oral cavity. It is a type of malignancy that shows high morbidity and recurrence rate, but data on the disease's target genes and biomarkers is still insufficient. In this study, in silico studies have been performed to find out the novel target genes and their potential therapeutic inhibitors for the effective and efficient treatment of OSCC. The DESeq2 package of RStudio was used in the current investigation to screen and identify differentially expressed genes for OSCC. As a result of gene expression analysis, the top 10 novel genes were identified using the Cytohubba plugin of Cytoscape, and among them, the ubiquitin-conjugating enzyme (UBE2D1) was found to be upregulated and playing a significant role in the progression of human oral cancers. Following this, naturally occurring compounds were virtually evaluated and simulated against the discovered novel target as prospective drugs utilizing the Maestro, Schrodinger, and Gromacs software. In a simulated screening of naturally occurring potential inhibitors against the novel target UBE2D1, Epigallocatechin 3-gallate, Quercetin, Luteoline, Curcumin, and Baicalein were identified as potent inhibitors. Novel identified gene UBE2D1 has a significant role in the proliferation of human cancers through suppression of 'guardian of genome' p53 via ubiquitination dependent pathway. Therefore, the treatment of OSCC may benefit significantly from targeting this gene and its discovered naturally occurring inhibitors.Communicated by Ramaswamy H. Sarma.

Unraveling effective extracellular signal-regulated kinase 2 inhibitors: a de novo drug design strategy enhanced by in-depth in silico analyses

Extracellular signal-regulated kinase 2 (ERK-2) is a serine/threonine protein kinase in eukaryotic cells and belongs to the mitogen-activated protein kinase (MAPK) family. An activated form of ERK-2 phosphorylates substrates in the nucleus or cytoplasm and causes specific proteins to be expressed or activated, regulating cell proliferation, differentiation and other functions. Caffeic acid (3,4 - dihydroxy cinnamic acid), as previously reported, directly interacts with ERK-2 and reduces its effects in vitro. It is also reported to have a variety of pharmacological effects, including anti-inflammatory, immunomodulatory, antioxidant and anticancer activities. In the current study, a deep-learning protocol was employed to develop effective 100 compounds by modifying the chemical structure of DHC to improve its inhibitory performance against ERK-2. Calculations of physicochemical properties for those compounds revealed that 20 compounds had drug scores better than DHC (≥ 80%). Following that, molecular docking calculations were performed on the selected compounds and DHC. The obtained data revealed that five compounds had docking scores better than DHC (≥ -5.9 kcal/mol). Moreover, data from molecular mechanics and the Poisson - Boltzmann surface area (MM/PBSA) binding energy over 200 ns MD simulation confirmed that Cmd-1 and Cmd-4 exhibited higher stability with ΔGbinding of -40.8 and -49.1 kcal/mol, respectively, which is better than DHC (-35.1 kcal/mol). Finally, various energetic and structural studies showed the high stability of the two generated compounds within the active site of ERK-2. This study highlights the potential use of Cmd-1 and Cmd-4 as promising anti-ERK-2 drug candidates.Communicated by Ramaswamy H. Sarma.

Tunable crystalline assemblies using surface-engineered protein cages

Assembly of nanoparticles into superlattices yields nanomaterials with novel properties. We have recently shown that engineered protein cages are excellent building blocks for the assembly of inorganic nanoparticles into highly structured hybrid materials, with unprecedented precision. In this study, we show that the protein matrix, composed of surface-charged protein cages, can be readily tuned to achieve a number of different crystalline assemblies. Simply by altering the assembly conditions, different types of crystalline structures were produced, without the need to further modify the cages. Future work can utilize these new protein scaffolds to create nanoparticle superlattices with various assembly geometries and thus tune the functionality of these hybrid materials.

Identification and mechanism of wheat protein disulfide isomerase-promoted gluten network formation

Formation of the gluten network depends on glutenin crosslinking via disulfide bonds, and wheat protein disulfide isomerase (wPDI) plays an important role in this process. Here, we identify a substrate gluten protein of wPDI and the mechanism underlying wPDI-promoted glutenin crosslinking. Farinographic, rheologic, and alveographic analysis unambiguously proves that wPDI improves gluten network formation, which is directly observed by 3D reconstruction of the gluten network. Protein analysis and LC-MS/MS reveal that glutenin subunit 1Dx5 is primarily recruited by wPDI to participate in gluten network formation, and its cysteine-containing N-terminal domain (1Dx5-NTD), which harbors three cysteine residues for crosslinking, is purified. 1Dx5-NTD interacts with wPDI in both redox states, possibly folded by reduced wPDI and then catalyzed by oxidized wPDI, as further evidenced by wPDI-promoted self-crosslinking. Consistent with macroscopic observations, our results suggest that wPDI folds 1Dx5-NTD into β-strand structure that favors disulfide bond formation.

Structure and rational engineering of the PglX methyltransferase and specificity factor for BREX phage defence

Bacteria have evolved a broad range of systems that provide defence against their viral predators, bacteriophages. Bacteriophage Exclusion (BREX) systems recognise and methylate 6 bp non-palindromic motifs within the host genome, and prevent replication of non-methylated phage DNA that encodes these same motifs. How BREX recognises cognate motifs has not been fully understood. In this study we characterise BREX from pathogenic Salmonella and present X-ray crystallographic structures of the conserved BREX protein, PglX. The PglX N-terminal domain encodes the methyltransferase, whereas the C-terminal domain is for motif recognition. We also present the structure of PglX bound to the phage-derived DNA mimic, Ocr, an inhibitor of BREX activity. Our analyses propose modes for DNA-binding by PglX and indicate that both methyltransferase activity and defence require larger BREX complexes. Through rational engineering of PglX we broaden both the range of phages targeted, and the host motif sequences that are methylated by BREX. Our data demonstrate that PglX is used to recognise specific DNA sequences for BREX activity, contributing to motif recognition for both phage defence and host methylation.

Effect of ethanol on the elasticities of double-stranded RNA and DNA revealed by magnetic tweezers and simulations

The elasticities of double-stranded (ds) DNA and RNA, which are critical to their biological functions and applications in materials science, can be significantly modulated by solution conditions such as ions and temperature. However, there is still a lack of a comprehensive understanding of the role of solvents in the elasticities of dsRNA and dsDNA in a comparative way. In this work, we explored the effect of ethanol solvent on the elasticities of dsRNA and dsDNA by magnetic tweezers and all-atom molecular dynamics simulations. We found that the bending persistence lengths and contour lengths of dsRNA and dsDNA decrease monotonically with the increase in ethanol concentration. Furthermore, the addition of ethanol weakens the positive twist-stretch coupling of dsRNA, while promotes the negative twist-stretch coupling of dsDNA. Counter-intuitively, the lower dielectric environment of ethanol causes a significant re-distribution of counterions and enhanced ion neutralization, which overwhelms the enhanced repulsion along dsRNA/dsDNA, ultimately leading to the softening in bending for dsRNA and dsDNA. Moreover, for dsRNA, ethanol causes slight ion-clamping across the major groove, which weakens the major groove-mediated twist-stretch coupling, while for dsDNA, ethanol promotes the stretch-radius correlation due to enhanced ion binding and consequently enhances the helical radius-mediated twist-stretch coupling.

Scorpion α-toxin LqhαIT specifically interacts with a glycan at the pore domain of voltage-gated sodium channels

Voltage-gated sodium (Nav) channels sense membrane potential and drive cellular electrical activity. The deathstalker scorpion α-toxin LqhαIT exerts a strong action potential prolonging effect on Nav channels. To elucidate the mechanism of action of LqhαIT, we determined a 3.9 Å cryoelectron microscopy (cryo-EM) structure of LqhαIT in complex with the Nav channel from Periplaneta americana (NavPas). We found that LqhαIT binds to voltage sensor domain 4 and traps it in an "S4 down" conformation. The functionally essential C-terminal epitope of LqhαIT forms an extensive interface with the glycan scaffold linked to Asn330 of NavPas that augments a small protein-protein interface between NavPas and LqhαIT. A combination of molecular dynamics simulations, structural comparisons, and prior mutagenesis experiments demonstrates the functional importance of this toxin-glycan interaction. These findings establish a structural basis for the specificity achieved by scorpion α-toxins and reveal the conserved glycan as an essential component of the toxin-binding epitope.

Structural study of a light chain mispaired bispecific predicts mechanism of downstream separation

Bispecific antibodies expressed and assembled from a single upstream culture require the correct balance and pairing of four different heavy and light chains (HC and LC). The increased potential for chain-mispaired species challenges the downstream purification of this new format. While clearance of HC-mispaired species, including homodimers and half-antibodies, has been assessed, removal of LC mispairs requires a more stringent approach. Here, we report two case studies in which separation is achieved, as well as the structural basis of these separations: (A) In the first case, a main species with a positively charged patch in the correctly formed variable fragment (Fv) is disrupted when paired with the wrong LC. This LC-mispaired variant binds more weakly to a cation exchange resin and can be washed off in a chromatography step. (B) A second molecule whose LC mispair introduces a negative-charge patch and hydrophobic patch in close proximity, presenting increased binding to a multimodal anion exchange resin. This LC-mispaired variant can be retained on the column under conditions in which the bispecific is recovered. In both case studies, the molecular structural analysis by protein surface properties models correlated well with the chromatography experiments. The comprehensive interpretation of experimental and computational results has provided a better understanding of strategies and potential applications for predicting the downstream purification of complex molecules.

Chemical Synthesis and Structure-Activity Relationship Studies of the Coagulation Factor Xa Inhibitor Tick Anticoagulant Peptide from the Hematophagous Parasite Ornithodoros moubata

Tick Anticoagulant Peptide (TAP), a 60-amino acid protein from the soft tick Ornithodoros moubata, inhibits activated coagulation factor X (fXa) with almost absolute specificity. Despite TAP and Bovine Pancreatic Trypsin Inhibitor (BPTI) (i.e., the prototype of the Kunitz-type protease inhibitors) sharing a similar 3D fold and disulphide bond topology, they have remarkably different amino acid sequence (only ~24% sequence identity), thermal stability, folding pathways, protease specificity, and even mechanism of protease inhibition. Here, fully active and correctly folded TAP was produced in reasonably high yields (~20%) by solid-phase peptide chemical synthesis and thoroughly characterised with respect to its chemical identity, disulphide pairing, folding kinetics, conformational dynamics, and fXa inhibition. The versatility of the chemical synthesis was exploited to perform structure-activity relationship studies on TAP by incorporating non-coded amino acids at positions 1 and 3 of the inhibitor. Using Hydrogen-Deuterium Exchange Mass Spectrometry, we found that TAP has a remarkably higher conformational flexibility compared to BPTI, and propose that these different dynamics could impact the different folding pathway and inhibition mechanisms of TAP and BPTI. Hence, the TAP/BPTI pair represents a nice example of divergent evolution, while the relative facility of TAP synthesis could represent a good starting point to design novel synthetic analogues with improved pharmacological profiles.

Structural modifications and kinetic effects of KRAS interactions with HRAS and NRAS: an in silico comparative analysis of KRAS mutants

The RAS genes which code for KRAS, HRAS, and NRAS are three of the most frequently mutated oncogenes responsible for cancer deaths. Tumorigenesis is one of the most significant outcomes of deregulation of RAS GTPases. Although the structures have been extensively studied, there is still more to be discovered about the actual binding conformations of the three isoforms, especially when mutated, to design an inhibitory drug. Recent studies have identified important interactions between the three isoforms that affect the oncogenic strength of the others when they are mutated. In this study, we utilize molecular dynamics simulations to examine the modifications of the structural property, mechanism, and kinetic energy of KRAS when interacting individually and with HRAS and NRAS. Notably, we found that WT-KRAS' orientation when bound to WT-HRAS vs. WT-NRAS is rotated 180°, with mutants demonstrating a similar binding pattern. The binding sites of the isoforms with KRAS share similarities with those involved in the GDP/GTP active site and site of KRAS dimerization. Thus, the isoform interaction can serve as an inhibitory method of KRAS actions. This study advances the understanding of inhibiting RAS-driven cancers through a novel isoform interaction approach only recently discovered, which has been proven to be an effective alternate therapeutic approach. We developed a blueprint of the interaction which would be beneficial in the development of KRAS mutant-specific and pan-KRAS mutant inhibitory drugs that mimic the isoform interactions. Our results support the direct interaction inhibition mechanism of mutant KRAS when bound to WT-HRAS and WT-NRAS by the isoforms' hypervariable region binding to the G-domain of KRAS. Furthermore, our results support the approach of reducing the effects of oncogenic KRAS by altering the concentration of the isoforms or a drug alternative based on the overall structural and kinetic stability, as well as the binding strength of the mutant-isoform complexes.

The molecular basis of Human FN3K mediated phosphorylation of glycated substrate

Glycation, a non-enzymatic post-translational modification occurring on proteins, can be actively reversed via site-specific phosphorylation of the fructose-lysine moiety by FN3K kinase, to impact the cellular function of target protein. A regulatory axis between FN3K and glycated protein targets has been associated with conditions like diabetes and cancer. However the molecular basis of this relationship has not been explored so far. Here, we determined a series of crystal structures of HsFN3K in apo-state, and in complex with different nucleotide analogs together with a sugar substrate mimic to reveal the features important for its kinase activity and substrate recognition. Additionally, the dynamics in sugar substrate binding during the kinase catalytic cycle provide important mechanistic insights into HsFN3K function. Our structural work provides the molecular basis for rationale small molecule design targeting FN3K.

Competitive Interaction of the SGFRKMAF Peptide with 3CLpro Dimerization Intermediates: A Brownian Dynamics Investigation

The SGFRKMAF peptide is known to inhibit the dimerization of 3CLpro monomers, which is essential for SARS-CoV-2 replication. The mechanism behind this, however, is largely unknown. In this work, we used Brownian dynamics simulations to compare and contrast 3CLpro monomer-monomer interactions and 3CLpro monomer-SGFRKMAF peptide interactions. We found that formation of the 3CLpro wild-type dimer could potentially involve formation of three intermediates that are primarily stabilized by G11-G124, S1-S301, and T118-G278 interactions. Analysis of 3CLpro monomer interaction with the SGFRKMAF peptide, however, revealed the presence of eight basins of interactions where the peptide assumes the highest local densities at the 3CLPro monomer surface. The second highest-density basin was found to coincide with the interface region of the wild-type 3CLpro dimer, thereby directly blocking the 3CLpro dimer-dimer interactions. The other basins, however, were found to lie far from the interface region. Notably, we found that only 6% of the BD trajectories end up directly into the basin at the interface region and ∼39% of the trajectories end up into those basins lying away from the interface region, indicating a greater role for peptide binding at sites away from the dimer interface region. Importantly, the locations of the basins lying away from the interface were found to coincide with the 3CLpro residues involved in stabilization of the 3CLpro monomer-monomer intermediates. Given that the rate constant of the peptide reaching the monomer surface was found to be almost an order of magnitude higher than the rate constant of monomer-monomer association, the SGFRKMAF peptide has the potential to inhibit dimerization of 3CLpro monomers not only through blocking the interface region but also through blocking the formation of the intermediates involved in the dimerization process. This could potentially open new avenues for 3CLpro dimerization inhibitors that transcend traditional X-ray-based discovery approaches.

Bioinformatics-based identification of GH12 endoxyloglucanases in citrus-pathogenic Penicillium spp

Millions of tons of citrus peel waste are produced every year as a byproduct of the juice industry. Citrus peel is rich in pectin and xyloglucan, but while the pectin is extracted for use in the food industry, the xyloglucan is currently not valorized. To target hydrolytic degradation of citrus peel xyloglucan into oligosaccharides, we have used bioinformatics to identify three glycoside hydrolase 12 (GH12) endoxyloglucanases (EC 3.2.1.151) from the citrus fruit pathogens Penicillium italicum GL-Gan1 and Penicillium digitatum Pd1 and characterized them on xyloglucan obtained by alkaline extraction from citrus peel. The enzymes displayed pH-temperature optima of pH 4.6-5.3 and 35-37°C. PdGH12 from P. digitatum and PiGH12A from P. italicum share 84% sequence identity and displayed similar kinetics, although kcat was highest for PdGH12. In contrast, PiGH12B from P. italicum, which has the otherwise conserved Trp in subsite -4 replaced with a Tyr, displayed a 3 times higher KM and a 4 times lower kcat/KM than PiGH12A, but was the most thermostable enzyme of the three Penicillium-derived endoxyloglucanases. The benchmark enzyme AnGH12 from Aspergillus nidulans was more thermally stable and had a higher pH-temperature optimum than the enzymes from Penicillum spp. The difference in structure of the xyloglucan oligosaccharides extracted from citrus peel xyloglucan and tamarind xyloglucan by the new endoxyloglucanases was determined by LC-MS. The inclusion of citrus peel xyloglucan demonstrated that the endoxyloglucanases liberated fucosylated xyloglucan oligomers, implying that these enzymes have the potential to upgrade citrus peel residues to produce oligomers useful as intermediates or bioactive compounds.

Effect of E134K pathogenic mutation of SMN protein on SMN-SmD1 interaction, with implication in spinal muscular atrophy: A molecular dynamics study

Spinal muscular atrophy (SMA) is a disease that results from mutations in the Survival of Motor Neuron (SMN) gene 1, leading to muscle atrophy due to motor neurons degeneration. SMN plays a crucial role in the assembly of spliceosomal small nuclear ribonucleoprotein complexes via binding to the arginine-glycine rich C-terminal tails of Sm proteins recognized by SMN Tudor domain. E134K Tudor mutation, cause of the more severe type I SMA, compromises the SMN-Sm interaction without a perturbation of the domain fold. By molecular dynamics simulations, we investigated the mechanism of Tudor-SmD1 interaction, and the effects on it of E134K mutation. It was observed that E134 is crucial to catch the positive dimethylated arginines (DMRs) of the SmD1 tail that, wrapping around the acidic Tudor surface, enters a central DMR into an aromatic cage. The flexible cage residue Y130 must be blocked from the wrapped tail to assure a stable binding. The charge inversion in E134K mutation causes the loss of a critical anchor point, disfavoring the tail wrapping and leaving Y130 free to swing, leading to DMR detachments and exposition of the C-terminal region of the tail. This could suggest new hypotheses regarding a possible autoimmune response by anti-Sm autoantibodies.

Structural basis for Mis18 complex assembly and its implications for centromere maintenance

The centromere, defined by the enrichment of CENP-A (a Histone H3 variant) containing nucleosomes, is a specialised chromosomal locus that acts as a microtubule attachment site. To preserve centromere identity, CENP-A levels must be maintained through active CENP-A loading during the cell cycle. A central player mediating this process is the Mis18 complex (Mis18α, Mis18β and Mis18BP1), which recruits the CENP-A-specific chaperone HJURP to centromeres for CENP-A deposition. Here, using a multi-pronged approach, we characterise the structure of the Mis18 complex and show that multiple hetero- and homo-oligomeric interfaces facilitate the hetero-octameric Mis18 complex assembly composed of 4 Mis18α, 2 Mis18β and 2 Mis18BP1. Evaluation of structure-guided/separation-of-function mutants reveals structural determinants essential for cell cycle controlled Mis18 complex assembly and centromere maintenance. Our results provide new mechanistic insights on centromere maintenance, highlighting that while Mis18α can associate with centromeres and deposit CENP-A independently of Mis18β, the latter is indispensable for the optimal level of CENP-A loading required for preserving the centromere identity.

SARS-CoV-2: Prediction of critical ionic amino acid mutations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), that caused coronavirus disease 2019 (COVID-19), has been studied thoroughly, and several variants are revealed across the world with their corresponding mutations. Studies and vaccines development focus on the genetic mutations of the S protein due to its vital role in allowing the virus attach and fuse with the membrane of a host cell. In this perspective, we study the effects of all ionic amino acid mutations of the SARS-CoV-2 viral spike protein S1 when bound to Antibody CC12.1 within the SARS-CoV-2:CC12.1 complex model. Binding free energy calculations between SARS-CoV-2 and antibody CC12.1 are based on the Analysis of Electrostatic Similarities of Proteins (AESOP) framework, where the electrostatic potentials are calculated using Adaptive Poisson-Boltzmann Solver (APBS). The atomic radii and charges that feed into the APBS calculations are calculated using the PDB2PQR software. Our results are the first to propose in silico potential life-threatening mutations of SARS-CoV-2 beyond the present mutations found in the five common variants worldwide. We find each of the following mutations: K378A, R408A, K424A, R454A, R457A, K458A, and K462A, to play significant roles in the binding to Antibody CC12.1, since they are turned into strong inhibitors on both chains of the S1 protein, whereas the mutations D405A, D420A, and D427A, show to play important roles in this binding, as they are turned into mild inhibitors on both chains of the S1 protein.

A prion-like domain is required for phase separation and chloroplast RNA processing during cold acclimation in Arabidopsis

Arabidopsis (Arabidopsis thaliana) plants can produce photosynthetic tissue with active chloroplasts at temperatures as low as 4°C, and this process depends on the presence of the nuclear-encoded, chloroplast-localized RNA-binding protein CP29A. In this study, we demonstrate that CP29A undergoes phase separation in vitro and in vivo in a temperature-dependent manner, which is mediated by a prion-like domain (PLD) located between the two RNA recognition motif domains of CP29A. The resulting droplets display liquid-like properties and are found near chloroplast nucleoids. The PLD is required to support chloroplast RNA splicing and translation in cold-treated tissue. Together, our findings suggest that plant chloroplast gene expression is compartmentalized by inducible condensation of CP29A at low temperatures, a mechanism that could play a crucial role in plant cold resistance.

Gadolinium-Based NMR Spin Relaxation Measurements of Near-Surface Electrostatic Potentials of Biomolecules

NMR spectroscopy is an important tool for the measurement of the electrostatic properties of biomolecules. To this point, paramagnetic relaxation enhancements (PREs) of 1H nuclei arising from nitroxide cosolutes in biomolecular solutions have been used to measure effective near-surface electrostatic potentials (ϕENS) of proteins and nucleic acids. Here, we present a gadolinium (Gd)-based NMR method, exploiting Gd chelates with different net charges, for measuring ϕENS values and demonstrate its utility through applications to a number of biomolecular systems. The use of Gd-based cosolutes offers several advantages over nitroxides for ϕENS measurements. First, unlike nitroxide compounds, Gd chelates enable electrostatic potential measurements on oxidation-sensitive proteins that require reducing agents. Second, the large electron spin quantum number of Gd (7/2) results in notably larger PREs for Gd chelates when used at the same concentrations as nitroxide radicals. Thus, it is possible to measure ϕENS values exclusively from + and - charged compounds even for highly charged biomolecules, avoiding the use of neutral cosolutes that, as we further establish here, limits the accuracy of the measured electrostatic potentials. In addition, the smaller concentrations of cosolutes required minimize potential binding to sites on macromolecules. Fourth, the closer proximity of the paramagnetic center and charged groups within Gd chelates, in comparison to the corresponding nitroxide compounds, enables more accurate predictions of ϕENS potentials for cross-validation of the experimental results. The Gd-based method described here, thus, broadens the applicability of studies of biomolecular electrostatics using solution NMR spectroscopy.

Molecular mechanisms of proteoglycan-mediated semaphorin signaling in axon guidance

The precise assembly of a functional nervous system relies on axon guidance cues. Beyond engaging their cognate receptors and initiating signaling cascades that modulate cytoskeletal dynamics, guidance cues also bind components of the extracellular matrix, notably proteoglycans, yet the role and mechanisms of these interactions remain poorly understood. We found that Drosophila secreted semaphorins bind specifically to glycosaminoglycan (GAG) chains of proteoglycans, showing a preference based on the degree of sulfation. Structural analysis of Sema2b unveiled multiple GAG-binding sites positioned outside canonical plexin-binding site, with the highest affinity binding site located at the C-terminal tail, characterized by a lysine-rich helical arrangement that appears to be conserved across secreted semaphorins. In vivo studies revealed a crucial role of the Sema2b C-terminal tail in specifying the trajectory of olfactory receptor neurons. We propose that secreted semaphorins tether to the cell surface through interactions with GAG chains of proteoglycans, facilitating their presentation to cognate receptors on passing axons.

ROS1 kinase inhibition reimagined: identifying repurposed drug via virtual screening and molecular dynamics simulations for cancer therapeutics

Precision medicine has revolutionized modern cancer therapeutic management by targeting specific molecular aberrations responsible for the onset and progression of tumorigenesis. ROS proto-oncogene 1 (ROS1) is a receptor tyrosine kinase (RTK) that can induce tumorigenesis through various signaling pathways, such as cell proliferation, survival, migration, and metastasis. It has emerged as a promising therapeutic target in various cancer types. However, there is very limited availability of specific ROS1 inhibitors for therapeutic purposes. Exploring repurposed drugs for rapid and effective treatment is a useful approach. In this study, we utilized an integrated approach of virtual screening and molecular dynamics (MD) simulations of repurposing existing drugs for ROS1 kinase inhibition. Using a curated library of 3648 FDA-approved drugs, virtual screening identified drugs capable of binding to ROS1 kinase domain. The results unveil two hits, Midostaurin and Alectinib with favorable binding profiles and stable interactions with the active site residues of ROS1. These hits were subjected to stability assessment through all-atom MD simulations for 200 ns. MD results showed that Midostaurin and Alectinib were stable with ROS1. Taken together, the study showed a rational framework for the selection of repurposed Midostaurin and Alectinib with ROS1 inhibitory potential for therapeutic management after further validation.

Ribosome External Electric Field Regulates Metabolic Enzyme Activity: The RAMBO Effect

Ribosomes bind to many metabolic enzymes and change their activity. A general mechanism for ribosome-mediated amplification of metabolic enzyme activity, RAMBO, was formulated and elucidated for the glycolytic enzyme triosephosphate isomerase, TPI. The RAMBO effect results from a ribosome-dependent electric field-substrate dipole interaction energy that can increase or decrease the ground state of the reactant and product to regulate catalytic rates. NMR spectroscopy was used to determine the interaction surface of TPI binding to ribosomes and to measure the corresponding kinetic rates in the absence and presence of intact ribosome particles. Chemical cross-linking and mass spectrometry revealed potential ribosomal protein binding partners of TPI. Structural results and related changes in TPI energetics and activity show that the interaction between TPI and ribosomal protein L11 mediate the RAMBO effect.

Influence of the Water Model on the Structure and Interactions of the GPR40 Protein with the Lipid Membrane and the Solvent: Rigid versus Flexible Water Models

G protein-coupled receptors (GPCR) are responsible for modulating various physiological functions and are thus related to the pathophysiology of different diseases. Being potential therapeutic targets, multiple computational methodologies have been developed to analyze their behavior and interactions with other species. The solvent, on the other hand, has received much less attention. In this work, we analyzed the effect of four explicit water models on the structure and interactions of the GPR40 receptor in its apo form. We employed the rigid SPC/E and TIP4P models, and their flexible versions, the FBA/ϵ and TIP4P/ϵflex. We explored the structural changes and their correlation with some bulk dynamic properties of water. Our results showed an adverse effect on the conservation of the secondary structure of the receptor with all the models due to the breaking of the intramolecular hydrogen bond network, being more evident for the TIP4P models. Notably, all four models brought the receptor to states similar to the active one, modifying the intracellular part of the TM5 and TM6 domains in a "hinge" type movement, allowing the opening of the structure. Regarding the dynamic properties, the rigid models showed results comparable to those obtained in other studies on membrane systems. However, flexible models exhibit disparities in the molecular representation of systems. Surprisingly, the FBA/ϵ model improves the molecular picture of several properties, even though their agreement with bulk diffusion is poorer. These findings reinforce our idea that exploring other water models or improving the current ones, to better represent the membrane interface, can lead to a positive impact on the description of the signal transduction mechanisms and the search of new drugs by targeting these receptors.

An Unfitted Finite Element Poisson-Boltzmann Solver with Automatic Resolving of Curved Molecular Surface

So far, the existing Poisson-Boltzmann (PB) solvers that accurately take into account the interface jump conditions need a pregenerated body-fitted mesh (molecular surface mesh). However, qualified biomolecular surface meshing and its implementation into numerical methods remains a challenging and laborious issue, which practically hinders the progress of further developments and applications of a bunch of numerical methods in this field. In addition, even with a molecular surface mesh, it is only a low-order approximation of the original curved surface. In this article, an interface-penalty finite element method (IPFEM), which is a typical unfitted finite element method, is proposed to solve the Poisson-Boltzmann equation (PBE) without requiring the user to generate a molecular surface mesh. The Gaussian molecular surface is used to represent the molecular surface and can be automatically resolved with a high-order approximation within our method. Theoretical convergence rates of the IPFEM for the linear PB equation have been provided and are well validated on a benchmark problem with an analytical solution (we also noticed from numerical examples that the IPFEM has similar convergence rates for the nonlinear PBE). Numerical results on a set of different-sized biomolecules demonstrate that the IPFEM is numerically stable and accurate in the calculation of biomolecular electrostatic solvation energy.

Conformational Properties of Poly(A)-Binding Protein Complexed with Poly(A) RNA

Microscopic understanding of protein-RNA interactions is important for different biological activities, such as RNA transport, translation, splicing, silencing, etc. Polyadenine (Poly(A)) binding proteins (PABPs) make up a class of regulatory proteins that play critical roles in protecting the poly(A) tails of cellular mRNAs from nuclease degradation. In this work, we performed molecular dynamics simulations to investigate the conformational modifications of human PABP protein and poly(A) RNA that occur during complexation. It is demonstrated that the intermediate linker domain of the protein transforms from a disordered coil-like structure to a helical form during the recognition process, leading to the formation of the complex. On the other hand, disordered collapsed coil-like RNA on complexation has been found to transform into a rigid extended conformation. Importantly, the binding free energy calculation showed that the thermodynamic stability of the complex is primarily guided by favorable hydrophobic interactions between the protein and the RNA.

Anti-Staphy Peptides Rationally Designed from Cry10Aa Bacterial Protein

Bacterial infections pose a significant threat to human health, constituting a major challenge for healthcare systems. Antibiotic resistance is particularly concerning in the context of treating staphylococcal infections. In addressing this challenge, antimicrobial peptides (AMPs), characterized by their hydrophobic and cationic properties, unique mechanism of action, and remarkable bactericidal and immunomodulatory capabilities, emerge as promising alternatives to conventional antibiotics for tackling bacterial multidrug resistance. This study focuses on the Cry10Aa protein as a template for generating AMPs due to its membrane-penetrating ability. Leveraging the Joker algorithm, six peptide variants were derived from α-helix 3 of Cry10Aa, known for its interaction with lipid bilayers. In vitro, antimicrobial assays determined the minimum inhibitory concentration (MIC) and minimum bactericidal concentration (MBC) required for inhibiting the growth of Staphylococcus aureus, Escherichia coli, Acinetobacter baummanii, Enterobacter cloacae, Enterococcus facallis, Klebsiella pneumonia, and Pseudomonas aeruginosa. Time-kill kinetics were performed using the parental peptide AMPCry10Aa, as well as AMPCry10Aa_1 and AMPCry10Aa_5, against E. coli ATCC, S. aureus 111 and S. aureus ATCC strains showing that AMPCry10Aa_1 and AMPCry10Aa_5 peptides can completely reduce the initial bacterial load with less than 2 h of incubation. AMPCry10Aa_1 and AMPCry 10Aa_5 present stability in human serum and activity maintenance up to 37 °C. Cytotoxicity assays, conducted using the MTT method, revealed that all of the tested peptides exhibited cell viability >50% (IC50). The study also encompassed evaluations of the structure and physical-chemical properties. The three-dimensional structures of AMPCry10Aa and AMPCry10Aa_5 were determined through nuclear magnetic resonance (NMR) spectroscopy, indicating the adoption of α-helical segments. Electron paramagnetic resonance (EPR) spectroscopy elucidated the mechanism of action, demonstrating that AMPCry10Aa_5 enters the outer membranes of E. coli and S. aureus, causing substantial increases in lipid fluidity, while AMPCry10Aa slightly increases lipid fluidity in E. coli. In conclusion, the results obtained underscore the potential of Cry10Aa as a source for developing antimicrobial peptides as alternatives to conventional antibiotics, offering a promising avenue in the battle against antibiotic resistance.

Mechanism of Abnormal Activation of MEK1 Induced by Dehydroalanine Modification

Mitogen-activated protein kinase kinase 1 (MAPK kinase 1, MEK1) is a key kinase in the mitogen-activated protein kinase (MAPK) signaling pathway. MEK1 mutations have been reported to lead to abnormal activation that is closely related to the malignant growth and spread of various tumors, making it an important target for cancer treatment. Targeting MEK1, four small-molecular drugs have been approved by the FDA, including Trametinib, Cobimetinib, Binimetinib, and Selumetinib. Recently, a study showed that modification with dehydroalanine (Dha) can also lead to abnormal activation of MEK1, which has the potential to promote tumor development. In this study, we used molecular dynamics simulations and metadynamics to explore the mechanism of abnormal activation of MEK1 caused by the Dha modification and predicted the inhibitory effects of four FDA-approved MEK1 inhibitors on the Dha-modified MEK1. The results showed that the mechanism of abnormal activation of MEK1 caused by the Dha modification is due to the movement of the active segment, which opens the active pocket and exposes the catalytic site, leading to sustained abnormal activation of MEK1. Among four FDA-approved inhibitors, only Selumetinib clearly blocks the active site by changing the secondary structure of the active segment from α-helix to disordered loop. Our study will help to explain the mechanism of abnormal activation of MEK1 caused by the Dha modification and provide clues for the development of corresponding inhibitors.

Nicotine Enantioselectively Targets Myeloid Differentiation Protein 2 and Inhibits the Toll-like Receptor 4 Signaling

Psychoactive substances, including morphine and methamphetamine, have been shown to interact with the classic innate immune receptor Toll-like receptor 4 (TLR4) and its partner protein myeloid differentiation protein 2 (MD2) in a nonenantioselective manner. (-)-Nicotine, the primary alkaloid in tobacco and a key component of highly addictive cigarettes, targets the TLR4/MD2, influencing TLR4 signaling pathways. Existing as two enantiomers, the stereoselective recognition of nicotine by TLR4/MD2 in the context of the innate immune response remains unclear. In this study, we synthesized (+)-nicotine and investigated its effects alongside (-)-nicotine on lipopolysaccharide (LPS)-induced TLR4 signaling. (-)-Nicotine dose-dependently inhibited proinflammatory factors such as tumor necrosis factor α (TNF-α), interleukin 6 (IL-6), and cyclooxygenase-2 (COX-2). In contrast, (+)-nicotine showed no such inhibitory effects. Molecular dynamics simulations revealed that (-)-nicotine exhibited a stronger affinity with the TLR4 coreceptor MD2 than (+)-nicotine. Additionally, in silico simulations revealed that both nicotine enantiomers initially attach to the entrance of the MD2 cavity, creating a metastable state before they fully enter the cavity. In the metastable state, (-)-nicotine established more stable interactions with the surrounding residues at the entrance of the MD2 cavity compared to those of (+)-nicotine. This highlights the crucial role of the MD2 cavity entrance in the chiral recognition of nicotine. These findings provide valuable insights into the distinct interactions between nicotine enantiomers and the TLR4 coreceptor MD2, underscoring the enantioselective effect of nicotine on modulating TLR4 signaling.

A Membrane-Assisted Mechanism for the Release of Ceramide from the CERT START Domain

Ceramide transfer protein CERT is the mediator of nonvesicular transfer of ceramide from the ER to Golgi. In CERT, START is the domain responsible for the binding and transport of ceramide. A wealth of structural data has revealed a helix-grip fold surrounding a large hydrophobic cavity holding the ceramide. Yet, little is known about the mechanisms by which START releases the ceramide through the polar region and into the packed environment of cellular membranes. As such events do not lend themselves easily to experimental investigations, we used multiple unbiased microsecond-long molecular simulations. We propose a membrane-assisted mechanism in which the membrane acts as an allosteric effector initiating the release of ceramide and where the passage of the ceramide acyl chains is facilitated by the intercalation of a single phosphatidylcholine lipid in the cavity, practically greasing the ceramide way out. We verify using free energy calculation and experimental lipidomics data that CERT forms stable complexes with phosphatidylcholine lipids, in addition to ceramide, thus providing validation for the proposed mechanism.

Structural study for substrate recognition of human N-terminal glutamine amidohydrolase 1 in the arginine N-degron pathway

The N-degron pathway determines the half-life of proteins by selectively destabilizing the proteins bearing N-degrons. N-terminal glutamine amidohydrolase 1 (NTAQ1) plays an essential role in the arginine N-degron (Arg/N-degron) pathway as an initializing enzyme via the deamidation of the N-terminal (Nt) glutamine (Gln). However, the Nt-serine-bound conformation of hNTAQ1 according to the previously identified crystal structure suggests the possibility of other factors influencing the recognition of Nt residues by hNTAQ1. Hence, in the current study, we aimed to further elucidate the substrate recognition of hNTAQ1; specifically, we explored 12 different substrate-binding conformations of hNTAQ1 depending on the subsequent residue of Nt-Gln. Results revealed that hNTAQ1 primarily interacts with the protein Nt backbone, instead of the side chain, for substrate recognition. Here, we report that the Nt backbone of proteins appears to be a key component of hNTAQ1 function and is the main determinant of substrate recognition. Moreover, not all second residues from Nt-Gln, but rather distinctive and charged residues, appeared to aid in detecting substrate recognition. These new findings define the substrate-recognition process of hNTAQ1 and emphasize the importance of the subsequent Gln residue in the Nt-Gln degradation system. Our extensive structural and biochemical analyses provide insights into the substrate specificity of the N-degron pathway and shed light on the mechanism underlying hNTAQ1 substrate recognition. An improved understanding of the protein degradation machinery could aid in developing therapies to promote overall health through enhanced protein regulation, such as targeted protein therapies.

Characterization of the Charge Heterogeneity of a Monoclonal Antibody That Binds to Both Cation Exchange and Anion Exchange Columns under the Same Binding Conditions

Therapeutic antibodies play an important role in the public healthcare system to treat patients with a variety of diseases. Protein characterization using an array of analytical tools provides in-depth information for drug quality, safety, efficacy, and the further understanding of the molecule. A therapeutic antibody candidate MAB1 exhibits unique binding properties to both cation and anion exchange columns at neutral pH. This uniqueness disrupts standard purification processes and necessitates adjustments in manufacturing. This study identifies that the charge heterogeneity of MAB1 is primarily due to the N-terminal cyclization of glutamine to pyroglutamine and, to a lesser extent, succinimide intermediate, deamidation, and C-terminal lysine. Using three approaches, i.e., deferential chemical labeling, H/D exchange, and molecular modeling, the binding to anion exchange resins is attributed to negatively charged patches on the antibody's surface, involving specific carboxylic acid residues. The methodologies shown here can be extended to study protein binding orientation in column chromatography.

Inclusion of Water Multipoles into the Implicit Solvation Framework Leads to Accuracy Gains

The current practical "workhorses" of the atomistic implicit solvation─the Poisson-Boltzmann (PB) and generalized Born (GB) models─face fundamental accuracy limitations. Here, we propose a computationally efficient implicit solvation framework, the Implicit Water Multipole GB (IWM-GB) model, that systematically incorporates the effects of multipole moments of water molecules in the first hydration shell of a solute, beyond the dipole water polarization already present at the PB/GB level. The framework explicitly accounts for coupling between polar and nonpolar contributions to the total solvation energy, which is missing from many implicit solvation models. An implementation of the framework, utilizing the GAFF force field and AM1-BCC atomic partial charges model, is parametrized and tested against the experimental hydration free energies of small molecules from the FreeSolv database. The resulting accuracy on the test set (RMSE ∼ 0.9 kcal/mol) is 12% better than that of the explicit solvation (TIP3P) treatment, which is orders of magnitude slower. We also find that the coupling between polar and nonpolar parts of the solvation free energy is essential to ensuring that several features of the IWM-GB model are physically meaningful, including the sign of the nonpolar contributions.

Long-range charge transfer mechanism of the III2IV2 mycobacterial supercomplex

Aerobic life is powered by membrane-bound redox enzymes that shuttle electrons to oxygen and transfer protons across a biological membrane. Structural studies suggest that these energy-transducing enzymes operate as higher-order supercomplexes, but their functional role remains poorly understood and highly debated. Here we resolve the functional dynamics of the 0.7 MDa III2IV2 obligate supercomplex from Mycobacterium smegmatis, a close relative of M. tuberculosis, the causative agent of tuberculosis. By combining computational, biochemical, and high-resolution (2.3 Å) cryo-electron microscopy experiments, we show how the mycobacterial supercomplex catalyses long-range charge transport from its menaquinol oxidation site to the binuclear active site for oxygen reduction. Our data reveal proton and electron pathways responsible for the charge transfer reactions, mechanistic principles of the quinone catalysis, and how unique molecular adaptations, water molecules, and lipid interactions enable the proton-coupled electron transfer (PCET) reactions. Our combined findings provide a mechanistic blueprint of mycobacterial supercomplexes and a basis for developing drugs against pathogenic bacteria.

Differences in the Membrane-Binding Properties of Flaviviral Nonstructural 1 (NS1) Protein: Comparative Simulations of Zika and Dengue Virus NS1 Proteins in Explicit Bilayers

NS1 in flaviviruses is the only nonstructural protein that is secretory and interacts with different cellular components of the host cell membrane. NS1 is localized in the ER as a dimer to facilitate viral replication. Crystal structures of NS1 homologues from zika (ZIKV) and dengue (DENV) viruses have revealed the organization of different domains in NS1 dimers. The β-roll and the connector and intertwined loop regions of wing domains of NS1 have been shown to interact with the membranes. In this study, we have performed multiple molecular dynamics (MD) simulations of ZIKV and DENV NS1 systems in apo and in POPE bilayers with different cholesterol concentrations (0, 20 and 40%). The NS1 protein was placed just above the membrane surface, and for each NS1-membrane system two to three independent simulations with 600 ns production run were performed. At the end of the production runs, ZIKV NS1 inserts deeper inside the membrane compared to the DENV counterpart. Unlike ZIKV NS1, the orientation of DENV NS1 is asymmetric in which one of the chains in the dimer interacts with the membrane while the other is more exposed to the solvent. The β-roll region in ZIKV NS1 penetrates beyond the headgroup region and interacts with the lipid acyl chains while the C-terminal region barely interacts with the headgroup. Specific residues in the intertwined region deeply penetrate inside the membrane. The role of charged and aromatic residues of ZIKV NS1 in strongly interacting with the membrane components is revealed. The presence of cholesterol affects the extent of insertion in the membrane and interaction of individual residues. Overall, membrane-binding properties of ZIKV NS1 significantly differ from its counterpart in DENV. The differences found in the binding and insertion of NS1 can be used to design drugs and novel antibodies that can be flavivirus specific.

Ca2+ permeation through C-terminal cleaved, but not full-length human Pannexin1 hemichannels, mediates cell death

Pannexin1 hemichannels (Panx1 HCs) are found in the membrane of most mammalian cells and communicate the intracellular and extracellular spaces, enabling the passive transfer of ions and small molecules. They are involved in physiological and pathophysiological conditions. During apoptosis, the C-terminal tail of Panx1 is proteolytically cleaved, but the permeability features of hemichannels and their role in cell death remain elusive. To address these topics, HeLa cells transfected with full-length human Panx1 (fl-hPanx1) or C-terminal truncated hPanx1 (Δ371hPanx1) were exposed to alkaline extracellular saline solution, increasing the activity of Panx1 HCs. The Δ371hPanx1 HC was permeable to DAPI and Etd+, but not to propidium iodide, whereas fl-hPanx1 HC was only permeable to DAPI. Furthermore, the cytoplasmic Ca2+ signal increased only in Δ371hPanx1 cells, which was supported by bioinformatics approaches. The influx of Ca2+ through Δ371hPanx1 HCs was necessary to promote cell death up to about 95% of cells, whereas the exposure to alkaline saline solution without Ca2+ failed to induce cell death, and the Ca2+ ionophore A23187 promoted more than 80% cell death even in fl-hPanx1 transfectants. Moreover, cell death was prevented with carbenoxolone or 10Panx1 in Δ371hPanx1 cells, whereas it was undetectable in HeLa Panx1-/- cells. Pretreatment with Ferrostatin-1 and necrostatin-1 did not prevent cell death, suggesting that ferroptosis or necroptosis was not involved. In comparison, zVAD-FMK, a pancaspase inhibitor, reduced death by ~60%, suggesting the involvement of apoptosis. Therefore, alkaline pH increases the activity of Δ371hPanx1HCs, leading to a critical intracellular free-Ca2+ overload that promotes cell death.

Rational Approach toward COVID-19's Main Protease Inhibitors: A Hierarchical Biochemoinformatics Analysis

This study investigated the potential of selected compounds as inhibitors of SARS-CoV-2 Mpro through pharmacokinetic and toxicological analyses, molecular docking, and molecular dynamics simulations. In silico molecular docking simulations revealed promising ligands with favorable binding affinities for Mpro, ranging from -6.2 to -9.5 kcal/mol. Moreover, molecular dynamics simulations demonstrated the stability of protein-ligand complexes over 200 ns, maintaining protein secondary structures. MM-PBSA analysis revealed favorable interactions between ligands and Mpro, with negative binding energy values. Hydrogen bond formation capacity during molecular dynamics was confirmed, indicating consistent interactions with Mpro catalytic residues. Based on these findings, selected ligands show promise for future studies in developing COVID-19 treatments.

Antibiotic discovery with artificial intelligence for the treatment of Acinetobacter baumannii infections

Global challenges presented by multidrug-resistant Acinetobacter baumannii infections have stimulated the development of new treatment strategies. We reported that outer membrane protein W (OmpW) is a potential therapeutic target in A. baumannii. Here, a library of 11,648 natural compounds was subjected to a primary screening using quantitative structure-activity relationship (QSAR) models generated from a ChEMBL data set with >7,000 compounds with their reported minimal inhibitory concentration (MIC) values against A. baumannii followed by a structure-based virtual screening against OmpW. In silico pharmacokinetic evaluation was conducted to assess the drug-likeness of these compounds. The ten highest-ranking compounds were found to bind with an energy score ranging from -7.8 to -7.0 kcal/mol where most of them belonged to curcuminoids. To validate these findings, one lead compound exhibiting promising binding stability as well as favorable pharmacokinetics properties, namely demethoxycurcumin, was tested against a panel of A. baumannii strains to determine its antibacterial activity using microdilution and time-kill curve assays. To validate whether the compound binds to the selected target, an OmpW-deficient mutant was studied and compared with the wild type. Our results demonstrate that demethoxycurcumin in monotherapy and in combination with colistin is active against all A. baumannii strains. Finally, the compound was found to significantly reduce the A. baumannii interaction with host cells, suggesting its anti-virulence properties. Collectively, this study demonstrates machine learning as a promising strategy for the discovery of curcuminoids as antimicrobial agents for combating A. baumannii infections.

IMPORTANCE

Acinetobacter baumannii presents a severe global health threat, with alarming levels of antimicrobial resistance rates resulting in significant morbidity and mortality in the USA, ranging from 26% to 68%, as reported by the Centers for Disease Control and Prevention (CDC). To address this threat, novel strategies beyond traditional antibiotics are imperative. Computational approaches, such as QSAR models leverage molecular structures to predict biological effects, expediting drug discovery. We identified OmpW as a potential therapeutic target in A. baumannii and screened 11,648 natural compounds. We employed QSAR models from a ChEMBL bioactivity data set and conducted structure-based virtual screening against OmpW. Demethoxycurcumin, a lead compound, exhibited promising antibacterial activity against A. baumannii, including multidrug-resistant strains. Additionally, demethoxycurcumin demonstrated anti-virulence properties by reducing A. baumannii interaction with host cells. The findings highlight the potential of artificial intelligence in discovering curcuminoids as effective antimicrobial agents against A. baumannii infections, offering a promising strategy to address antibiotic resistance.

In vitro and in silico investigation of effects of antimicrobial peptides from Solanaceae plants against rice sheath blight pathogen Rhizoctinia solani

Rhizoctonia solani, the causative agent of sheath blight disease in rice, poses a significant threat to agricultural productivity. Traditional management approaches involving chemical fungicides have been effective but come with detrimental consequences for the ecosystem. This study aimed to investigate sustainable alternatives in the form of antifungal peptides derived from Solanaceous plant species as potential agents against R. solani. Peptide extracts were obtained using an optimized antimicrobial peptide (AMP) extraction method and desalted using the solid-phase extraction technique. The antifungal potential of peptide-rich extracts from Solanum tuberosum and Capsicum annum was assessed through in vitro tests employing the agar well diffusion method. Furthermore, peptide-protein docking analysis was performed on HPEPDOCK and HDOCK server; and molecular dynamics simulations (MDS) of 100 ns period were performed using the Gromacs 2020.4. The results demonstrated significant inhibition zones for both extracts at concentrations of 100 mg/mL. Additionally, the extracts of Solanum tuberosum and Capsicum annum had minimum inhibitory concentrations of 50 mg/mL and 25 mg/mL, respectively with minimum fungicidal concentrations of 25 mg/mL. Insights into the potential mechanisms of key peptides inhibiting R. solani targets were gleaned from in-silico studies. Notably, certain AMPs exhibited favorable free energy of binding against pathogenicity-related targets, including histone demethylase, sortin nexin, and squalene synthase, in protein-peptide docking simulations. Extended molecular dynamics simulations lasting 100 ns and MM-PBSA calculations were performed on select protein-peptide complexes. AMP10 displayed the most favorable binding free energy against all target proteins, with AMP3, AMP12b, AMP6, and AMP15 also exhibiting promising results against specific targets of R. solani. These findings underscore the potential of peptide extracts from S. tuberosum and C. annum as effective antifungal agents against rice sheath blight caused by R. solani.