GATSol, an enhanced predictor of protein solubility through the synergy of 3D structure graph and large language modeling.

BACKGROUND: Protein solubility is a critically important physicochemical property closely related to protein expression. For example, it is one of the main factors to be considered in the design and production of antibody drugs and a prerequisite for realizing various protein functions. Although several solubility prediction models have emerged in recent years, many of these models are limited to capturing information embedded in one-dimensional amino acid sequences, resulting in unsatisfactory predictive performance. RESULTS: In this study, we introduce a novel Graph Attention network-based protein Solubility model, GATSol, which represents the 3D structure of proteins as a protein graph. In addition to the node features of amino acids extracted by the state-of-the-art protein large language model, GATSol utilizes amino acid distance maps generated using the latest AlphaFold technology. Rigorous testing on independent eSOL and the Saccharomyces cerevisiae test datasets has shown that GATSol outperforms most recently introduced models, especially with respect to the coefficient of determination R2, which reaches 0.517 and 0.424, respectively. It outperforms the current state-of-the-art GraphSol by 18.4% on the S. cerevisiae_test set. CONCLUSIONS: GATSol captures 3D dimensional features of proteins by building protein graphs, which significantly improves the accuracy of protein solubility prediction. Recent advances in protein structure modeling allow our method to incorporate spatial structure features extracted from predicted structures into the model by relying only on the input of protein sequences, which simplifies the entire graph neural network prediction process, making it more user-friendly and efficient. As a result, GATSol may help prioritize highly soluble proteins, ultimately reducing the cost and effort of experimental work. The source code and data of the GATSol model are freely available at https://github.com/binbinbinv/GATSol .

Explore new quinoxaline pharmacophore tethered sulfonamide fragments as in vitro α-glucosidase, α-amylase, and acetylcholinesterase inhibitors with ADMET and molecular modeling simulation.

A new series of quinoxaline-sulfonamide derivatives 3-12 were synthesized using fragment-based drug design by reaction of quinoxaline sulfonyl chloride (QSC) with different amines and hydrazines. The quinoxaline-sulfonamide derivatives were evaluated for antidiabetic and anti-Alzheimer's potential against α-glucosidase, α-amylase, and acetylcholinesterase enzymes. These derivatives showed good to moderate potency against α-amylase and α-glucosidase with inhibitory percentages between 24.34 ± 0.01%-63.09 ± 0.02% and 28.95 ± 0.04%-75.36 ± 0.01%, respectively. Surprisingly, bis-sulfonamide quinoxaline derivative 4 revealed the most potent activity with inhibitory percentages of 75.36 ± 0.01% and 63.09 ± 0.02% against α-glucosidase and α-amylase compared to acarbose (IP = 57.79 ± 0.01% and 67.33 ± 0.01%), respectively. Moreover, the quinoxaline derivative 3 exhibited potency as α-glucosidase and α-amylase inhibitory with a minute decline from compound 4 and acarbose with inhibitory percentages of 44.93 ± 0.01% and 38.95 ± 0.01%. Additionally, in vitro acetylcholinesterase inhibitory activity for designed derivatives exhibited weak to moderate activity. Still, sulfonamide-quinoxaline derivative 3 emerged as the most active member with inhibitory percentage of 41.92 ± 0.02% compared with donepezil (IP = 67.27 ± 0.60%). The DFT calculations, docking simulation, target prediction, and ADMET analysis were performed and discussed in detail.

Structural basis for human mitochondrial tRNA maturation.

The human mitochondrial genome is transcribed into two RNAs, containing mRNAs, rRNAs and tRNAs, all dedicated to produce essential proteins of the respiratory chain. The precise excision of tRNAs by the mitochondrial endoribonucleases (mt-RNase), P and Z, releases all RNA species from the two RNA transcripts. The tRNAs then undergo 3'-CCA addition. In metazoan mitochondria, RNase P is a multi-enzyme assembly that comprises the endoribonuclease PRORP and a tRNA methyltransferase subcomplex. The requirement for this tRNA methyltransferase subcomplex for mt-RNase P cleavage activity, as well as the mechanisms of pre-tRNA 3'-cleavage and 3'-CCA addition, are still poorly understood. Here, we report cryo-EM structures that visualise four steps of mitochondrial tRNA maturation: 5' and 3' tRNA-end processing, methylation and 3'-CCA addition, and explain the defined sequential order of the tRNA processing steps. The methyltransferase subcomplex recognises the pre-tRNA in a distinct mode that can support tRNA-end processing and 3'-CCA addition, likely resulting from an evolutionary adaptation of mitochondrial tRNA maturation complexes to the structurally-fragile mitochondrial tRNAs. This subcomplex can also ensure a tRNA-folding quality-control checkpoint before the sequential docking of the maturation enzymes. Altogether, our study provides detailed molecular insight into RNA-transcript processing and tRNA maturation in human mitochondria.

Reactivities of acrylamide warheads toward cysteine targets: a QM/ML approach to covalent inhibitor design.

Covalent inhibition offers many advantages over non-covalent inhibition, but covalent warhead reactivity must be carefully balanced to maintain potency while avoiding unwanted side effects. While warhead reactivities are commonly measured with assays, a computational model to predict warhead reactivities could be useful for several aspects of the covalent inhibitor design process. Studies have shown correlations between covalent warhead reactivities and quantum mechanic (QM) properties that describe important aspects of the covalent reaction mechanism. However, the models from these studies are often linear regression equations and can have limitations associated with their usage. Applications of machine learning (ML) models to predict covalent warhead reactivities with QM descriptors are not extensively seen in the literature. This study uses QM descriptors, calculated at different levels of theory, to train ML models to predict reactivities of covalent acrylamide warheads. The QM/ML models are compared with linear regression models built upon the same QM descriptors and with ML models trained on structure-based features like Morgan fingerprints and RDKit descriptors. Experiments show that the QM/ML models outperform the linear regression models and the structure-based ML models, and literature test sets demonstrate the power of the QM/ML models to predict reactivities of unseen acrylamide warhead scaffolds. Ultimately, these QM/ML models are effective, computationally feasible tools that can expedite the design of new covalent inhibitors.

The step-by-step assembly mechanism of secreted flavivirus NS1 tetramer and hexamer captured at atomic resolution.

Flaviviruses encode a conserved, membrane-associated nonstructural protein 1 (NS1) with replication and immune evasion functions. The current knowledge of secreted NS1 (sNS1) oligomers is based on several low-resolution structures, thus hindering the development of drugs and vaccines against flaviviruses. Here, we revealed that recombinant sNS1 from flaviviruses exists in a dynamic equilibrium of dimer-tetramer-hexamer states. Two DENV4 hexameric NS1 structures and several tetrameric NS1 structures from multiple flaviviruses were solved at atomic resolution by cryo-EM. The stacking of the tetrameric NS1 and hexameric NS1 is facilitated by the hydrophobic ß-roll and connector domains. Additionally, a triacylglycerol molecule located within the central cavity may play a role in stabilizing the hexamer. Based on differentiated interactions between the dimeric NS1, two distinct hexamer models (head-to-head and side-to-side hexamer) and the step-by-step assembly mechanisms of NS1 dimer into hexamer were proposed. We believe that our study sheds light on the understanding of the NS1 oligomerization and contributes to NS1-based therapies.

TRPV3 activation by different agonists accompanied by lipid dissociation from the vanilloid site.

TRPV3 represents both temperature- and ligand-activated transient receptor potential (TRP) channel. Physiologically relevant opening of TRPV3 channels by heat has been captured structurally, while opening by agonists has only been observed in structures of mutant channels. Here, we present cryo-EM structures that illuminate opening and inactivation of wild-type human TRPV3 in response to binding of two types of agonists: either the natural cannabinoid tetrahydrocannabivarin (THCV) or synthetic agonist 2-aminoethoxydiphenylborane (2-APB). We found that THCV binds to the vanilloid site, while 2-APB binds to the S1-S4 base and ARD-TMD linker sites. Despite binding to distally located sites, both agonists induce similar pore opening and cause dissociation of a lipid that occupies the vanilloid site in their absence. Our results uncover different but converging allosteric pathways through which small-molecule agonists activate TRPV3 and provide a framework for drug design and understanding the role of lipids in ion channel function.

Calcium-induced structural transitions are central to the folding, function, and processing of serratiopeptidase zymogen into mature form.

Serratia marcescens is an emerging health-threatening, gram-negative opportunistic pathogen associated with a wide variety of localized and life-threatening systemic infections. One of the most crucial virulence factors produced by S. marcescens is serratiopeptidase, a 50.2-kDa repeats-in-toxin (RTX) family broad-specificity zinc metalloprotease. RTX family proteins are functionally diverse exoproteins of gram-negative bacteria that exhibit calcium-dependent structural dynamicity and are secreted through a common type-1 secretion system (T1SS) machinery. To evaluate the impact of various divalent ligands on the folding and maturation of serratiopeptidase zymogen, the protein was purified and a series of structural and functional investigations were undertaken. The results indicate that calcium binding to the C-terminal RTX domain acts as a folding switch, triggering a disordered-to-ordered transition in the enzyme's conformation. Further, the auto-processing of the 16-amino acid N-terminal pro-peptide results in the maturation of the enzyme. The binding of calcium ions to serratiopeptidase causes a highly cooperative conformational transition in its structure, which is essential for the enzyme's activation and maturation. This conformational change is accompanied by an increase in solubility and enzymatic activity. For efficient secretion and to minimize intracellular toxicity, the enzyme needs to be in an unfolded extended form. The calcium-rich extracellular environment favors the folding and processing of zymogen into mature serratiopeptidase, i.e., the holo-form required by S. marcescens to establish infections and survive in different environmental niches.

Identifying and avoiding radiation damage in macromolecular crystallography.

Radiation damage remains one of the major impediments to accurate structure solution in macromolecular crystallography. The artefacts of radiation damage can manifest as structural changes that result in incorrect biological interpretations being drawn from a model, they can reduce the resolution to which data can be collected and they can even prevent structure solution entirely. In this article, we discuss how to identify and mitigate against the effects of radiation damage at each stage in the macromolecular crystal structure-solution pipeline.

Structural mechanisms of Tad pilus assembly and its interaction with an RNA virus.

Caulobacter crescentus Tad (tight adherence) pili, part of the type IV pili family, are crucial for mechanosensing, surface adherence, bacteriophage (phage) adsorption, and cell-cycle regulation. Unlike other type IV pilins, Tad pilins lack the typical globular ß sheet domain responsible for pilus assembly and phage binding. The mechanisms of Tad pilus assembly and its interaction with phage ΦCb5 have been elusive. Using cryo-electron microscopy, we unveiled the Tad pilus assembly mechanism, featuring a unique network of hydrogen bonds at its core. We then identified the Tad pilus binding to the ΦCb5 maturation protein (Mat) through its ß region. Notably, the amino terminus of ΦCb5 Mat is exposed outside the capsid and phage/pilus interface, enabling the attachment of fluorescent and affinity tags. These engineered ΦCb5 virions can be efficiently assembled and purified in Escherichia coli, maintaining infectivity against C. crescentus, which presents promising applications, including RNA delivery and phage display.

New-generation biofilm effective antimicrobial peptides and a real-time anti-biofilm activity assay: CoMIC.

Nowadays, it is very important to produce new-generation drugs with antimicrobial properties that will target biofilm-induced infections. The first target for combating these microorganisms, which are the source itself. Antimicrobial peptides, which are more effective than antibiotics due to their ability to kill microorganisms and use a different metabolic pathway, are among the new options today. The aim of this study is to develop new-generation antibiotics that inhibit both biofilm-producing bacteria and the biofilm itself. For this purpose, we designed four different peptides by combining two amino acid forms (D- and L-) with the same sequence having alpha helix structures. It was found that the combined use of these two forms can increase antimicrobial efficacy more than 30-fold. These results are supported by molecular modeling and scanning electron microscopy (SEM), at the same time cytotoxicity (IC50) and hemotoxicity (HC50) values remained within the safe range. Furthermore, antibiofilm activities of these peptides were investigated. Since the existing biofilm inhibition methods in the literature do not technically simulate the exact situation, in this study, we have developed a real-time observable biofilm model and a new detection method based on it, which we call the CoMIC method. Findings have shown that the NET1 peptide with D-leucine amino acid in its structure and the NET3 peptide with D-arginine amino acid in its structure are effective in inhibiting biofilm. As a conclusion, our peptides can be considered as potential next-generation broad-spectrum antibiotic molecule/drug candidates that might be used in biofilm and clinical important bacteria. KEY POINTS: • Antimicrobial peptides were developed to inhibit both biofilms producing bacteria and the biofilm itself. • CoMIC will fill a very crucial gap in understanding biofilms and conducting the necessary quantitative studies. • Molecular modelling studies, NET1 peptide molecules tends to move towards and adhere to the membrane within nanoseconds.

Structural basis and synergism of ATP and Na+ activation in bacterial K+ uptake system KtrAB.

The K+ uptake system KtrAB is essential for bacterial survival in low K+ environments. The activity of KtrAB is regulated by nucleotides and Na+. Previous studies proposed a putative gating mechanism of KtrB regulated by KtrA upon binding to ATP or ADP. However, how Na+ activates KtrAB and the Na+ binding site remain unknown. Here we present the cryo-EM structures of ATP- and ADP-bound KtrAB from Bacillus subtilis (BsKtrAB) both solved at 2.8 Å. A cryo-EM density at the intra-dimer interface of ATP-KtrA was identified as Na+, as supported by X-ray crystallography and ICP-MS. Thermostability assays and functional studies demonstrated that Na+ binding stabilizes the ATP-bound BsKtrAB complex and enhances its K+ flux activity. Comparing ATP- and ADP-BsKtrAB structures suggests that BsKtrB Arg417 and Phe91 serve as a channel gate. The synergism of ATP and Na+ in activating BsKtrAB is likely applicable to Na+-activated K+ channels in central nervous system.

Identification of domains in Plasmodium falciparum proteins of unknown function using DALI search on AlphaFold predictions.

Plasmodium falciparum, the causative agent of malaria, poses a significant global health challenge, yet much of its biology remains elusive. A third of the genes in the P. falciparum genome lack annotations regarding their function, impeding our understanding of the parasite's biology. In this study, we employ structure predictions and the DALI search algorithm to analyse proteins encoded by uncharacterized genes in the reference strain 3D7 of P. falciparum. By comparing AlphaFold predictions to experimentally determined protein structures in the Protein Data Bank, we found similarities to known domains in 353 proteins of unknown function, shedding light on their potential functions. The lowest-scoring 5% of similarities were additionally validated using the size-independent TM-align algorithm, confirming the detected similarities in 88% of the cases. Notably, in over 70 P. falciparum proteins the presence of domains resembling heptatricopeptide repeats, which are typically involvement in RNA binding and processing, was detected. This suggests this family, which is important in transcription in mitochondria and apicoplasts, is much larger in Plasmodium parasites than previously thought. The results of this domain search provide a resource to the malaria research community that is expected to inform and enable experimental studies.

Structural basis for the intracellular regulation of ferritin degradation.

The interaction between nuclear receptor coactivator 4 (NCOA4) and the iron storage protein ferritin is a crucial component of cellular iron homeostasis. The binding of NCOA4 to the FTH1 subunits of ferritin initiates ferritinophagy-a ferritin-specific autophagic pathway leading to the release of the iron stored inside ferritin. The dysregulation of NCOA4 is associated with several diseases, including neurodegenerative disorders and cancer, highlighting the NCOA4-ferritin interface as a prime target for drug development. Here, we present the cryo-EM structure of the NCOA4-FTH1 interface, resolving 16 amino acids of NCOA4 that are crucial for the interaction. The characterization of mutants, designed to modulate the NCOA4-FTH1 interaction, is used to validate the significance of the different features of the binding site. Our results explain the role of the large solvent-exposed hydrophobic patch found on the surface of FTH1 and pave the way for the rational development of ferritinophagy modulators.

HighFold: accurately predicting structures of cyclic peptides and complexes with head-to-tail and disulfide bridge constraints.

In recent years, cyclic peptides have emerged as a promising therapeutic modality due to their diverse biological activities. Understanding the structures of these cyclic peptides and their complexes is crucial for unlocking invaluable insights about protein target-cyclic peptide interaction, which can facilitate the development of novel-related drugs. However, conducting experimental observations is time-consuming and expensive. Computer-aided drug design methods are not practical enough in real-world applications. To tackles this challenge, we introduce HighFold, an AlphaFold-derived model in this study. By integrating specific details about the head-to-tail circle and disulfide bridge structures, the HighFold model can accurately predict the structures of cyclic peptides and their complexes. Our model demonstrates superior predictive performance compared to other existing approaches, representing a significant advancement in structure-activity research. The HighFold model is openly accessible at https://github.com/hongliangduan/HighFold.

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

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

Leptospiral LipL45 lipoprotein undergoes processing and shares structural similarities with bacterial sigma regulators.

Leptospirosis is a widespread zoonotic infectious disease of human and veterinary concern caused by pathogenic spirochetes of the genus Leptospira. To date, little progress towards understanding leptospiral pathogenesis and identification of virulence factors has been made, which is the main bottleneck for developing effective measures against the disease. Some leptospiral proteins, including LipL32, Lig proteins, LipL45, and LipL21, are being considered as potential virulence factors or vaccine candidates. However, their function remains to be established. LipL45 is the most expressed membrane lipoprotein in leptospires, upregulated when the bacteria are transferred to temperatures resembling the host, expressed during infection, suppressed after culture attenuation, and known to suffer processing in vivo and in vitro, generating fragments. Based on body of evidence, we hypothesized that the LipL45 processing might occur by an auto-cleavage event, deriving two fragments. The results presented here, based on bioinformatics, structure modeling analysis, and experimental data, corroborate that LipL45 processing probably includes a self-catalyzed non-proteolytic event and suggest the participation of LipL45 in cell-surface signaling pathways, as the protein shares structural similarities with bacterial sigma regulators. Our data indicate that LipL45 might play an important role in response to environmental conditions, with possible function in the adaptation to the host.

Novel structure of secreted small molecular weight antigen Mtb12 from Mycobacterium tuberculosis.

Mtb12, a small protein secreted by Mycobacterium tuberculosis, is known to elicit immune responses in individuals infected with the pathogen. It serves as an antigen recognized by the host's immune system. Due to its immunogenic properties and pivotal role in tuberculosis (TB) pathogenesis, Mtb12 is considered a promising candidate for TB diagnosis and vaccine development. However, the structural and functional properties of Mtb12 are largely unexplored, representing a significant gap in our understanding of M. tuberculosis biology. In this study, we present the first structure of Mtb12, which features a unique tertiary configuration consisting of four beta strands and four alpha helices. Structural analysis reveals that Mtb12 has a surface adorned with a negatively charged pocket adjacent to a central cavity. The features of these structural elements and their potential effects on the function of Mtb12 warrant further exploration. These findings offer valuable insights for vaccine design and the development of diagnostic tools.

Modulating anti-inflammatory and anticancer properties by designing a family of metal-complexes based on 5-nitropicolinic acid.

A new family of six complexes based on 5-nitropicolinic acid (5-npic) and transition metals has been obtained: [M(5-npic)2]n (MII = Mn (1) and Cd (2)), [Cu(5-npic)2]n (3), and [M(5-npic)2(H2O)2] (MII = Co (4), Ni (5), and Zn (6)), which display 1D, 2D, and mononuclear structures, respectively, thanks to different coordination modes of 5-npic. After their physicochemical characterization by single-crystal X-ray diffraction (SCXRD), elemental analyses (EA), and spectroscopic techniques, quantum chemical calculations using Time-Dependent Density Functional Theory (TD-DFT) were performed to further study the luminescence properties of compounds 2 and 6. The potential anticancer activity of all complexes was tested against three tumor cell lines, B16-F10, HT29, and HepG2, which are models widely used for studying melanoma, colon cancer, and liver cancer, respectively. The best results were found for compounds 2 and 4 against B16-F10 (IC50 = 26.94 and 45.10 µg mL-1, respectively). In addition, anti-inflammatory studies using RAW 264.7 cells exhibited promising activity for 2, 3, and 6 (IC50 NO = 5.38, 24.10, and 17.63 µg mL-1, respectively). This multidisciplinary study points to complex 2, based on CdII, as a promising anticancer and anti-inflammatory material.

Transformation of a Viral Capsid from Nanocages to Nanotubes and Then to Hydrogels: Redirected Self-Assembly and Effects on Immunogenicity.

The ability to manipulate the self-assembly of proteins is essential to understanding the mechanisms of life and beneficial to fabricating advanced nanomaterials. Here, we report the transformation of the MS2 phage capsid from nanocages to nanotubes and then to nanotube hydrogels through simple point mutations guided by interfacial interaction redesign. We demonstrate that site 70, which lies in the flexible FG loop of the capsid protein (CP), is a "magic" site that can largely dictate the final morphology of assemblies. By varying the amino acid at site 70, with the aid of a cysteine-to-alanine mutation at site 46, we achieved the assembly of double-helical or single-helical nanotubes in addition to nanocages. Furthermore, an additional cysteine substitution on the surface of nanotubes mediated their cross-linking to form hydrogels with reducing agent responsiveness. The hierarchical self-assembly system allowed for the investigation of morphology-related immunogenicity of MS2 CPs, which revealed dramatic differences among nanocages, nanotubes, and nanotube hydrogels in terms of immune response types, antibody levels and T cell functions. This study provides insights into the assembly manipulation of protein nanomaterials and the customized design of nanovaccines and drug delivery systems.

Molecular modeling of some commercially available antiviral drugs and their derivatives against SARS-CoV-2 infection.

Numerous prior studies have identified therapeutic targets that could effectively combat severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, including the angiotensin-converting enzyme 2 (ACE2) receptor, RNA-dependent RNA polymerase (RdRp), and Main protease (Mpro). In parallel, antiviral compounds like abacavir, acyclovir, adefovir, amantadine, amprenavir, darunavir, didanosine, oseltamivir, penciclovir, and tenofovir are under investigation for their potential in drug repurposing to address this infection. The aim of the study was to determine the effect of modifying the functional groups of the aforementioned antivirals in silico. Using the genetic optimization for ligand docking algorithm on software Maestro (version 11.1), the modified antivirals were docked onto ACE2 receptor, RdRp, and Mpro. Using QuickProp (Maestro v11.1), PASS (prediction of activity spectra for the substances), and altogether with SwissADME, the ADMET (absorption, distribution, metabolism, excretion, and toxicity) of the modified antivirals, as well as their bioavailability and the predicted activity spectra, were determined. Discovery studio software was used to undertake post-docking analysis. Among the 10 antivirals, N(CH3)2 derivative of darunavir, N(CH3)2 derivative of amprenavir and NCH3 derivative of darunavir exhibited best binding affinities with ACE2 receptor (docking scores: -10.333, -9.527 and -9.695 kJ/mol, respectively). Moreover, NCH3 derivative of abacavir (-6.506 kJ/mol), NO2 derivative of didanosine (-6.877 kJ/mol), NCH3 derivative of darunavir (-7.618 kJ/mol) exerted promising affinity to Mpro. In conclusion, the results of the in silico screenings can serve as a useful information for future experimental works.