Statistical potential for assessment and prediction of protein structures

Clinical features, functional consequences, and rescue pharmacology of missense GRID1 and GRID2 human variants

GRID1 and GRID2 encode the enigmatic GluD1 and GluD2 proteins, which form tetrameric receptors that play important roles in synapse organization and development of the central nervous system. Variation in these genes has been implicated in neurodevelopmental phenotypes. We evaluated GRID1 and GRID2 human variants from the literature, ClinVar, and clinical laboratories and found that many of these variants reside in intolerant domains, including the amino terminal domain of both GRID1 and GRID2. Other conserved regions, such as the M3 transmembrane domain, show different intolerance between GRID1 and GRID2. We introduced these variants into GluD1 and GluD2 cDNA and performed electrophysiological and biochemical assays to investigate the mechanisms of dysfunction of GRID1/2 variants. One variant in the GRID1 distal amino terminal domain resides at a position predicted to interact with Cbln2/Cbln4, and the variant disrupts complex formation between GluD1 and Cbln2, which could perturb its role in synapse organization. We also discovered that, like the lurcher mutation (GluD2-A654T), other rare variants in the GRID2 M3 domain create constitutively active receptors that share similar pathogenic phenotypes. We also found that the SCHEMA schizophrenia M3 variant GluD1-A650T produced constitutively active receptors. We tested a variety of compounds for their ability to inhibit constitutive currents of GluD receptor variants and found that pentamidine potently inhibited GluD2-T649A constitutive channels (IC50 50 nM). These results identify regions of intolerance to variation in the GRID genes, illustrate the functional consequences of GRID1 and GRID2 variants, and suggest how these receptors function normally and in disease.

Viper Venom Phospholipase A2 Database: The Structural and Functional Anatomy of a Primary Toxin in Envenomation

Viper venom phospholipase A2 enzymes (vvPLA2s) and phospholipase A2-like (PLA2-like) proteins are two of the principal toxins in viper venom that are responsible for the severe myotoxic and neurotoxic effects caused by snakebite envenoming, among other pathologies. As snakebite envenoming is the deadliest neglected tropical disease, a complete understanding of these proteins' properties and their mechanisms of action is urgently needed. Therefore, we created a database comprising information on the holo-form, cofactor-bound 3D structure of 217 vvPLA2 and PLA2-like proteins in their physiologic environment, as well as 79 membrane-bound viper species from 24 genera, which we have made available to the scientific community to accelerate the development of new anti-snakebite drugs. In addition, the analysis of the sequenced, 3D structure of the database proteins reveals essential aspects of the anatomy of the proteins, their toxicity mechanisms, and the conserved binding site areas that may anchor universal interspecific inhibitors. Moreover, it pinpoints hypotheses for the molecular origin of the myotoxicity of the PLA2-like proteins. Altogether, this study provides an understanding of the diversity of these toxins and how they are conserved, and it indicates how to develop broad, interspecies, efficient small-molecule inhibitors to target the toxin's many mechanisms of action.

Genomic surveillance of omicron B.1.1.529 SARS-CoV-2 and its variants between December 2021 and March 2023 in Tamil Nadu, India-A state-wide prospective longitudinal study

A state-wide prospective longitudinal investigation of the genomic surveillance of the omicron B.1.1.529 SARS-CoV-2 variant and its sublineages in Tamil Nadu, India, was conducted between December 2021 and March 2023. The study aimed to elucidate their mutational patterns and their genetic interrelationship in the Indian population. The study identified several unique mutations at different time-points, which likely could attribute to the changing disease characteristics, transmission, and pathogenicity attributes of omicron variants. The study found that the omicron variant is highly competent in its mutating potentials, and that it continues to evolve in the general population, likely escaping from natural as well as vaccine-induced immune responses. Our findings suggest that continuous surveillance of viral variants at the global scenario is warranted to undertake intervention measures against potentially precarious SARS-CoV-2 variants and their evolution.

Computational Design of α-Conotoxins to Target Specific Nicotinic Acetylcholine Receptor Subtypes

Nicotinic acetylcholine receptors (nAChRs) are drug targets for neurological diseases and disorders, but selective targeting of the large number of nAChR subtypes is challenging. Marine cone snail α-conotoxins are potent blockers of nAChRs and some have been engineered to achieve subtype selectivity. This engineering effort would benefit from rapid computational methods able to predict mutational energies, but current approaches typically require high-resolution experimental structures, which are not widely available for α-conotoxin complexes. Herein, five mutational energy prediction methods were benchmarked using crystallographic and mutational data on two acetylcholine binding protein/α-conotoxin systems. Molecular models were developed for six nAChR subtypes in complex with five α-conotoxins that were studied through 150 substitutions. The best method was a combination of FoldX and molecular dynamics simulations, resulting in a predictive Matthews Correlation Coefficient (MCC) of 0.68 (85 % accuracy). Novel α-conotoxin mutants designed using this method were successfully validated by experimental assay with improved pharmaceutical properties. This work paves the way for the rapid design of subtype-specific nAChR ligands and potentially accelerated drug development.

Making Use of Averaging Methods in MODELLER for Protein Structure Prediction

Recent advances in protein structure prediction, driven by AlphaFold 2 and machine learning, demonstrate proficiency in static structures but encounter challenges in capturing essential dynamic features crucial for understanding biological function. In this context, homology-based modeling emerges as a cost-effective and computationally efficient alternative. The MODELLER (version 10.5, accessed on 30 November 2023) algorithm can be harnessed for this purpose since it computes intermediate models during simulated annealing, enabling the exploration of attainable configurational states and energies while minimizing its objective function. There have been a few attempts to date to improve the models generated by its algorithm, and in particular, there is no literature regarding the implementation of an averaging procedure involving the intermediate models in the MODELLER algorithm. In this study, we examined MODELLER's output using 225 target-template pairs, extracting the best representatives of intermediate models. Applying an averaging procedure to the selected intermediate structures based on statistical potentials, we aimed to determine: (1) whether averaging improves the quality of structural models during the building phase; (2) if ranking by statistical potentials reliably selects the best models, leading to improved final model quality; (3) whether using a single template versus multiple templates affects the averaging approach; (4) whether the "ensemble" nature of the MODELLER building phase can be harnessed to capture low-energy conformations in holo structures modeling. Our findings indicate that while improvements typically fall short of a few decimal points in the model evaluation metric, a notable fraction of configurations exhibit slightly higher similarity to the native structure than MODELLER's proposed final model. The averaging-building procedure proves particularly beneficial in (1) regions of low sequence identity between the target and template(s), the most challenging aspect of homology modeling; (2) holo protein conformations generation, an area in which MODELLER and related tools usually fall short of the expected performance.

In silico screening of phytoconstituents as potential anti-inflammatory agents targeting NF-κB p65: an approach to promote burn wound healing

Chronic burn wounds are frequently characterised by a prolonged and dysregulated inflammatory phase that is mediated by over-activation of NF-κB p65. Synthetic wound healing drugs used for treatment of inflammation are primarily associated with several shortcomings which reduce their therapeutic index. In this scenario, phytoconstituents that exhibit multifaceted biological activities including anti-inflammatory effects have emerged as a promising therapeutic alternative. However, identification and isolation of phytoconstituents from medicinal herbs is a cumbersome method that is linked to profound uncertainty. Hence, present study aimed to identify prospective phytoconstituents as inhibitors of RHD of NF-κB p65 by utilizing in silico approach. Virtual screening of 2821 phytoconstituents was performed against protein model. Out of 2821 phytoconstituents, 162 phytoconstituents displayed a higher binding affinity (≤ -8.0 kcal/mol). These 162 phytoconstituents were subjected to ADMET predictions, and 15 of them were found to satisfy Lipinski's rule of five and showed favorable pharmacokinetic properties. Among these 15 phytoconstituents, 5 phytoconstituents with high docking scores i.e. silibinin, bismurrayaquinone A, withafastuosin B, yuccagenin, (+)-catechin 3-gallate were selected for molecular dynamics (MD) simulation analysis. Results of MD simulation indicated that withafastuosin B, (+)-catechin 3-gallate and yuccagenin produced a compact and stable complex with protein without significant variations in conformation. Relative binding energy analysis of best hit molecules indicate that withafastuosin B, and (+)-catechin 3-gallate exhibit high binding affinity with target protein among other lead molecules. Findings of study suggest that these phytoconstituents could serve as promising anti-inflammatory agents for treatment of burn wounds by inhibiting the RHD of NF-κB p65.Communicated by Ramaswamy H. Sarma.

Delivery to, and Reactivation of, the p53 Pathway in Cancer Cells Using a Grafted Cyclotide Conjugated with a Cell-Penetrating Peptide

Peptides are promising drug modalities that can modulate protein-protein interactions, but their application is hampered by their limited ability to reach intracellular targets. Here, we improved the cytosolic delivery of a peptide blocking p53:MDM2/X interactions using a cyclotide as a stabilizing scaffold. We applied several design strategies to improve intracellular delivery and found that the conjugation of the lead cyclotide to the cyclic cell-penetrating peptide cR10 was the most effective. Conjugation allowed cell internalization at micromolar concentration and led to elevated intracellular p53 levels in A549, MCF7, and MCF10A cells, as well as inducing apoptosis in A549 cells without causing membrane disruption. The lead peptide had >35-fold improvement in inhibitory activity and increased cellular uptake compared to a previously reported cyclotide p53 activator. In summary, we demonstrated the delivery of a large polar cyclic peptide in the cytosol and confirmed its ability to modulate intracellular protein-protein interactions involved in cancer.

RNAdvisor: a comprehensive benchmarking tool for the measure and prediction of RNA structural model quality

RNA is a complex macromolecule that plays central roles in the cell. While it is well known that its structure is directly related to its functions, understanding and predicting RNA structures is challenging. Assessing the real or predictive quality of a structure is also at stake with the complex 3D possible conformations of RNAs. Metrics have been developed to measure model quality while scoring functions aim at assigning quality to guide the discrimination of structures without a known and solved reference. Throughout the years, many metrics and scoring functions have been developed, and no unique assessment is used nowadays. Each developed assessment method has its specificity and might be complementary to understanding structure quality. Therefore, to evaluate RNA 3D structure predictions, it would be important to calculate different metrics and/or scoring functions. For this purpose, we developed RNAdvisor, a comprehensive automated software that integrates and enhances the accessibility of existing metrics and scoring functions. In this paper, we present our RNAdvisor tool, as well as state-of-the-art existing metrics, scoring functions and a set of benchmarks we conducted for evaluating them. Source code is freely available on the EvryRNA platform: https://evryrna.ibisc.univ-evry.fr.

Ligand coupling mechanism of the human serotonin transporter differentiates substrates from inhibitors

The presynaptic serotonin transporter (SERT) clears extracellular serotonin following vesicular release to ensure temporal and spatial regulation of serotonergic signalling and neurotransmitter homeostasis. Prescription drugs used to treat neurobehavioral disorders, including depression, anxiety, and obsessive-compulsive disorder, trap SERT by blocking the transport cycle. In contrast, illicit drugs of abuse like amphetamines reverse SERT directionality, causing serotonin efflux. Both processes result in increased extracellular serotonin levels. By combining molecular dynamics simulations with biochemical experiments and using a homologous series of serotonin analogues, we uncovered the coupling mechanism between the substrate and the transporter, which triggers the uptake of serotonin. Free energy analysis showed that only scaffold-bound substrates could initiate SERT occlusion through attractive long-range electrostatic interactions acting on the bundle domain. The associated spatial requirements define substrate and inhibitor properties, enabling additional possibilities for rational drug design approaches.

HumanLectome, an update of UniLectin for the annotation and prediction of human lectins

The UniLectin portal (https://unilectin.unige.ch/) was designed in 2019 with the goal of centralising curated and predicted data on carbohydrate-binding proteins known as lectins. UniLectin is also intended as a support for the study of lectomes (full lectin set) of organisms or tissues. The present update describes the inclusion of several new modules and details the latest (https://unilectin.unige.ch/humanLectome/), covering our knowledge of the human lectome and comprising 215 unevenly characterised lectins, particularly in terms of structural information. Each HumanLectome entry is protein-centric and compiles evidence of carbohydrate recognition domain(s), specificity, 3D-structure, tissue-based expression and related genomic data. Other recent improvements regarding interoperability and accessibility are outlined.

The Application of MD Simulation to Lead Identification, Vaccine Design, and Structural Studies in Combat against Leishmaniasis - A Review

Drug discovery, vaccine design, and protein interaction studies are rapidly moving toward the routine use of molecular dynamics simulations (MDS) and related methods. As a result of MDS, it is possible to gain insights into the dynamics and function of identified drug targets, antibody-antigen interactions, potential vaccine candidates, intrinsically disordered proteins, and essential proteins. The MDS appears to be used in all possible ways in combating diseases such as cancer, however, it has not been well documented as to how effectively it is applied to infectious diseases such as Leishmaniasis. As a result, this review aims to survey the application of MDS in combating leishmaniasis. We have systematically collected articles that illustrate the implementation of MDS in drug discovery, vaccine development, and structural studies related to Leishmaniasis. Of all the articles reviewed, we identified that only a limited number of studies focused on the development of vaccines against Leishmaniasis through MDS. Also, the PCA and FEL studies were not carried out in most of the studies. These two were globally accepted utilities to understand the conformational changes and hence it is recommended that this analysis should be taken up in similar approaches in the future.

Assessment of the bound conformation of bombesin to the BB1 and BB2 receptors

Bombesin is an endogenous peptide involved in a wide spectrum of physiological activities ranging from satiety, control of circadian rhythm and thermoregulation in the central nervous system, to stimulation of gastrointestinal hormone release, activation of macrophages and effects on development in peripheral tissues. Actions of the peptide are mediated through the two high affinity G-protein coupled receptors BB1R and BB2R. Under pathophysiological conditions, these receptors are overexpressed in many different types of tumors, such as prostate cancer, breast cancer, small and non-small cell lung cancer and pancreatic cancer. This observation has been used for designing cell markers, but it has not been yet exploited for therapeutical purposes. Despite the enormous biological interest of the peptide, little is known about the stereochemical features that contribute to their activity. On the one hand, mutagenesis studies identified a few receptor residues important for high bombesin affinity and on the other, a few studies focused on the relevance of diverse residues of the peptide for receptor activation. Models of the peptide bound to BB1R and BB2R can be helpful to improve our understanding of the stereochemical features granting bombesin activity. Accordingly, the present study describes the computational process followed to construct such models by means of Steered Molecular Dynamics, using models of the peptide and its receptors. Present results provide new insights into the structure-activity relationships of bombesin and its receptors, as well as render an explanation for the differential binding affinity observed towards BB1R and BB2R. Finally, these models can be further exploited to help for designing novel small molecule peptidomimetics with improved pharmacokinetics profile.

Interaction between Cfd1 and Nbp35 proteins involved in cytosolic FeS cluster assembly machinery deciphers a stable complexation in Leishmania donovani

Leishmania donovani is the causative unicellular parasite for visceral leishmaniasis (VL); and FeS proteins are likely to be very essential for their survival and viability. Cytosolic FeS cluster assembly (CIA) machinery is one of the four systems for the biosynthesis and transfer of FeS clusters among eukaryotes; Cfd1 and Nbp35 are the scaffold components for cytosolic FeS cluster biogenesis. We investigated the role of CIA machinery components and purified Cfd1 and Nbp35 proteins of L. donovani. We also investigated the interactive nature between LdCfd1 and LdNbp35 proteins by in silico analysis, in vitro co-purification, pull down assays along with in vivo immuno-precipitation; which inferred that both LdCfd1 and LdNbp35 proteins are interacting with each other. Thus, our collective data revealed the interaction between these two proteins which forms a stable complex that can be attributed to the cellular process of FeS clusters biogenesis, and transfer to target apo-proteins of L. donovani. The expression of Cfd1 and Nbp35 proteins in Amp B resistant parasites is up-regulated leading to increased amount of FeS proteins. Hence, it favors increased tolerance towards ROS level, which helps parasites survival under drug pressure contributing in Amphotericin B resistance.

Atomistic ensemble of active SHP2 phosphatase

SHP2 phosphatase plays an important role in regulating several intracellular signaling pathways. Pathogenic mutations of SHP2 cause developmental disorders and are linked to hematological malignancies and cancer. SHP2 comprises two tandemly-arranged SH2 domains, a catalytic PTP domain, and a disordered C-terminal tail. Under physiological, non-stimulating conditions, the catalytic site of PTP is occluded by the N-SH2 domain, so that the basal activity of SHP2 is low. Whereas the autoinhibited structure of SHP2 has been known for two decades, its active, open structure still represents a conundrum. Since the oncogenic mutant SHP2E76K almost completely populates the active, open state, this mutant has been extensively studied as a model for activated SHP2. By molecular dynamics simulations and accurate explicit-solvent SAXS curve predictions, we present the heterogeneous atomistic ensemble of constitutively active SHP2E76K in solution, encompassing a set of conformational arrangements and radii of gyration in agreement with experimental SAXS data.

The cytokinin receptor OHK4/OsHK4 regulates inflorescence architecture in rice via an IDEAL PLANT ARCHITECTURE1/WEALTHY FARMER'S PANICLE-mediated positive feedback circuit

Inflorescence architecture is important for rice (Oryza sativa) grain yield. The phytohormone cytokinin (CK) has been shown to regulate rice inflorescence development; however, the underlying mechanism mediated by CK perception is still unclear. Employing a forward genetic approach, we isolated an inactive variant of the CK receptor OHK4/OsHK4 gene named panicle length1, which shows decreased panicle size due to reduced inflorescence meristem (IM) activity. A 2-amino acid deletion in the long α-helix stalk of the sensory module of OHK4 impairs the homodimerization and ligand-binding capacity of the receptor, even though the residues do not touch the ligand-binding domain or the dimerization interface. This deletion impairs CK signaling that occurs through the type-B response regulator OsRR21, which acts downstream of OHK4 in controlling inflorescence size. Meanwhile, we found that IDEAL PLANT ARCHITECTURE1(IPA1)/WEALTHY FARMER'S PANICLE (WFP), encoding a positive regulator of IM development, acts downstream of CK signaling and is directly activated by OsRR21. Additionally, we revealed that IPA1/WFP directly binds to the OHK4 promoter and upregulates its expression through interactions with 2 TCP transcription factors, forming a positive feedback circuit. Altogether, we identified the OHK4-OsRR21-IPA1 regulatory module, providing important insights into the role of CK signaling in regulating rice inflorescence architecture.

Exploring the conformational changes of the Munc18-1/syntaxin 1a complex

Neurotransmitters are released from synaptic vesicles, the membrane of which fuses with the plasma membrane upon calcium influx. This membrane fusion reaction is driven by the formation of a tight complex comprising the plasma membrane N-ethylmaleimide-sensitive factor attachment receptor (SNARE) proteins syntaxin-1a and SNAP-25 with the vesicle SNARE protein synaptobrevin. The neuronal protein Munc18-1 forms a stable complex with syntaxin-1a. Biochemically, syntaxin-1a cannot escape the tight grip of Munc18-1, so formation of the SNARE complex is inhibited. However, Munc18-1 is essential for the release of neurotransmitters in vivo. It has therefore been assumed that Munc18-1 makes the bound syntaxin-1a available for SNARE complex formation. Exactly how this occurs is still unclear, but it is assumed that structural rearrangements occur. Here, we used a series of mutations to specifically weaken the complex at different positions in order to induce these rearrangements biochemically. Our approach was guided through sequence and structural analysis and supported by molecular dynamics simulations. Subsequently, we created a homology model showing the complex in an altered conformation. This conformation presumably represents a more open arrangement of syntaxin-1a that permits the formation of a SNARE complex to be initiated while still bound to Munc18-1. In the future, research should investigate how this central reaction for neuronal communication is controlled by other proteins.

Resolving binding pathways and solvation thermodynamics of plant hormone receptors

Plant hormones are small molecules that regulate plant growth, development, and responses to biotic and abiotic stresses. They are specifically recognized by the binding site of their receptors. In this work, we resolved the binding pathways for eight classes of phytohormones (auxin, jasmonate, gibberellin, strigolactone, brassinosteroid, cytokinin, salicylic acid, and abscisic acid) to their canonical receptors using extensive molecular dynamics simulations. Furthermore, we investigated the role of water displacement and reorganization at the binding site of the plant receptors through inhomogeneous solvation theory. Our findings predict that displacement of water molecules by phytohormones contributes to free energy of binding via entropy gain and is associated with significant free energy barriers for most systems analyzed. Also, our results indicate that displacement of unfavorable water molecules in the binding site can be exploited in rational agrochemical design. Overall, this study uncovers the mechanism of ligand binding and the role of water molecules in plant hormone perception, which creates new avenues for agrochemical design to target plant growth and development.

Pharmacoinformatic approach to identify potential phytochemicals against SARS-CoV-2 spike receptor-binding domain in native and variants of concern

Severe acute respiratory syndrome-related coronavirus-2 (SARS-CoV-2) pathogenesis is initiated by the binding of SARS-CoV-2 spike (S) protein with the angiotensin-converting enzyme 2 receptor (ACE2R) on the host cell surface. The receptor-binding domain (RBD) of the S protein mediates the binding and is more prone to mutations resulting in the generation of different variants. Recently, molecules with the potential to inhibit the interaction of S protein with ACE2R have been of interest due to their therapeutic value. In this context, the present work was performed to identify potential RBD binders from the Indian medicinal plant's phytochemical database through virtual screening, molecular docking, and molecular dynamic simulation. Briefly, 1578 compounds filtered from 9596 phytochemicals were chosen for screening against the RBD of the native SARS-CoV-2 S protein. Based on the binding energy, the top 30 compounds were selected and re-docked individually against the native and five variants of concern (VOCs: alpha, beta, gamma, delta, and omicron) of SARS-CoV-2. Four phytochemicals, namely withanolide F, serotobenine, orobanchol, and gibberellin A51, were found to be potential RBD binders in native and all SARS-CoV-2 VOCs. Among the four, withanolide F exhibited lower binding energy (- 10.84 to - 8.56 kcal/mol) and better ligand efficiency (- 0.3 to - 0.25) against all forms of RBD and hence was subjected to a 100 ns MD simulation which confirmed its stringent binding to the RBDs in native and VOCs. The study prioritizes withanolide F as a prospective COVID-19 (Coronavirus disease) therapeutic agent based on the observations. It warrants deeper investigations into the four promising leads for understanding their precise therapeutic value.

Regulation of STING activity in DNA sensing by ISG15 modification

Sensing of human immunodeficiency virus type 1 (HIV-1) DNA is mediated by the cyclic GMP-AMP synthase-stimulator of interferon genes (cGAS-STING) signaling axis. Signal transduction and regulation of this cascade is achieved by post-translational modifications. Here we show that cGAS-STING-dependent HIV-1 sensing requires interferon-stimulated gene 15 (ISG15). ISG15 deficiency inhibits STING-dependent sensing of HIV-1 and STING agonist-induced antiviral response. Upon external stimuli, STING undergoes ISGylation at residues K224, K236, K289, K347, K338, and K370. Inhibition of STING ISGylation at K289 suppresses STING-mediated type Ⅰ interferon induction by inhibiting its oligomerization. Of note, removal of STING ISGylation alleviates gain-of-function phenotype in STING-associated vasculopathy with onset in infancy (SAVI). Molecular modeling suggests that ISGylation of K289 is an important regulator of oligomerization. Taken together, our data demonstrate that ISGylation at K289 is crucial for STING activation and represents an important regulatory step in DNA sensing of viruses and autoimmune responses.

Computational Modeling Study of the Binding of Aging and Non-Aging Inhibitors with Neuropathy Target Esterase

Neuropathy target esterase (NTE) is a serine hydrolase with phospholipase B activity, which is involved in maintaining the homeostasis of phospholipids. It can be inhibited by aging inhibitors such as some organophosphorus (OP) compounds, which leads to delayed neurotoxicity with distal degeneration of axons. However, the detailed binding conformation of aging and non-aging inhibitors with NTE is not known. In this study, new computational models were constructed by using MODELLER 10.3 and AlphaFold2 to further investigate the inhibition mechanism of aging and non-aging compounds using molecular docking. The results show that the non-aging compounds bind the hydrophobic pocket much deeper than aging compounds and form the hydrophobic interaction with Phe1066. Therefore, the unique binding conformation of non-aging compounds may prevent the aging reaction. These important differences of the binding conformations of aging and non-aging inhibitors with NTE may help explain their different inhibition mechanism and the protection of non-aging NTE inhibitors against delayed neuropathy.

Improving antibody optimization ability of generative adversarial network through large language model

Generative adversarial networks (GANs) have successfully generated functional protein sequences. However, traditional GANs often suffer from inherent randomness, resulting in a lower probability of obtaining desirable sequences. Due to the high cost of wet-lab experiments, the main goal of computer-aided antibody optimization is to identify high-quality candidate antibodies from a large range of possibilities, yet improving the ability of GANs to generate these desired antibodies is a challenge. In this study, we propose and evaluate a new GAN called the Language Model Guided Antibody Generative Adversarial Network (AbGAN-LMG). This GAN uses a language model as an input, harnessing such models' powerful representational capabilities to improve the GAN's generation of high-quality antibodies. We conducted a comprehensive evaluation of the antibody libraries and sequences generated by AbGAN-LMG for COVID-19 (SARS-CoV-2) and Middle East Respiratory Syndrome (MERS-CoV). Results indicate that AbGAN-LMG has learned the fundamental characteristics of antibodies and that it improved the diversity of the generated libraries. Additionally, when generating sequences using AZD-8895 as the target antibody for optimization, over 50% of the generated sequences exhibited better developability than AZD-8895 itself. Through molecular docking, we identified 70 antibodies that demonstrated higher affinity for the wild-type receptor-binding domain (RBD) of SARS-CoV-2 compared to AZD-8895. In conclusion, AbGAN-LMG demonstrates that language models used in conjunction with GANs can enable the generation of higher-quality libraries and candidate sequences, thereby improving the efficiency of antibody optimization. AbGAN-LMG is available at http://39.102.71.224:88/.

In silico identification of five binding sites on the SARS-CoV-2 spike protein and selection of seven ligands for such sites

To contribute to the development of products capable of complexing with the SARS-CoV-2 spike protein, and thus preventing the virus from entering the host cell, this work aimed at discovering binding sites in the whole protein structure, as well as selecting substances capable of binding efficiently to such sites. Initially, the three-dimensional structure of the protein, with all receptor binding domains in the closed state, underwent blind docking with 38 substances potentially capable of binding to this protein according to the literature. This allowed the identification of five binding sites. Then, those substances with more affinities for these sites underwent pharmacophoric search in the ZINC15 database. The 14,329 substances selected from ZINC15 were subjected to docking to the five selected sites of the spike protein. The ligands with more affinities for the protein sites, as well as the selected sites themselves, were used in the de novo design of new ligands that were also docked to the binding sites of the protein. The best ligands, regardless of their origins, were used to form complexes with the spike protein, which were subsequently used in molecular dynamics simulations and calculations of ligands affinities to the protein through the molecular mechanics/Poisson-Boltzmann surface area method (MMPBSA). Seven substances with good affinities to the spike protein (-12.9 to -20.6 kcal/mol), satisfactory druggability (Bioavailability score: 0.17 to 0.55), and low acute toxicity to mice (LD50: 751 to 1421 mg/kg) were selected as potentially useful for the future development of new products to manage COVID-19 infections.Communicated by Ramaswamy H. Sarma.

EasyModel: a user-friendly web-based interface based on MODELLER

Three-dimensional protein structures are invaluable sources of information for the functional annotation of protein molecules. Describing the function of a protein sequence is one of the most common problems in biology. Generally, this problem can be facilitated by studying the tertiary structure of proteins. In the lack of protein structures, comparative modeling often provides a useful three-dimensional model of the protein associated with at least one known protein structure. Comparative modeling predicts the tertiary structure of a certain protein sequence (target) mainly based on its homological sequence to the sequence of one or more proteins with known structures (templates). MODELLER is one of the most widely used tools for homology or comparative modeling of three-dimensional protein structures. However, most users find it challenging to start with MODELLER as it is a command line based and requires knowledge of basic Python scripting to use it efficiently. In this study, a web-based interface has been designed to predict the tertiary structure of proteins based on Modeller, which does the comparative modeling automatically, and uses PHP and Python programming languages. This tool is called "EasyModel" and is available at http://bioinf.modares.ac.ir/software/easymodel/ . EasyModel provides a straightforward graphical interface for Modeller that can be used in only one browser.

Identification of PPA1 inhibitor candidates for potential repurposing in cancer medicine

Inorganic pyrophosphatase 1 (PPA1) is pivotal to cellular metabolism as it facilitates the hydrolysis of PPi-a by-product of various metabolic processes that influence cell growth and differentiation. Overexpression of PPA1 enzyme has been linked to diminished patient survival and was shown to influence tumor cell dynamics, thereby positioning it as a potential therapy target for a variety of cancers including colorectal cancer, diffuse large B-cell lymphoma, and lung adenocarcinoma. Despite this therapeutic promise, there are no known inhibitors of PPA1 as of today. In this study, we searched for potential PPA1 inhibitors using a molecular docking screen of 30 470 compounds with a history of clinical trials and/or US Food and Drug Administration approval. We specifically targeted the active pocket that coincides with the established catalytic domain. Our screen identified promising hits, which we further subjected to ADMET (absorption, distribution, metabolism, excretion, and toxicity) filtering. Subsequent molecular dynamics (MD) analyses were conducted on devazepide, quinotolast, and tarazepide-the three substances that successfully navigated all filters. MD analyses reinforced the stability of the protein-ligand complexes and confirmed ligand binding, as substantiated by our root mean square deviation, radius of gyration and secondary structures of proteins analyses. Furthermore, Molecular Mechanics Poisson-Boltzmann Surface Area calculations post-MD identified devazepide and quinotolast as showing higher binding affinities; being supported by principal component analysis, free energy landscape, and dynamic cross-correlation matrix results. Overall, our study reveals devazepide and quinotolast as potential candidates for PPA1 inhibition which could be considered for repurposing studies that need further experimental validation. These results not only reveal a potential for clinical repurposing for PPA1 inhibition but they also offer valuable insights into the development of future compounds for targeting the crucial PPA1 enzyme.

Targeting Toll-like receptor-driven systemic inflammation by engineering an innate structural fold into drugs

There is a clinical need for conceptually new treatments that target the excessive activation of inflammatory pathways during systemic infection. Thrombin-derived C-terminal peptides (TCPs) are endogenous anti-infective immunomodulators interfering with CD14-mediated TLR-dependent immune responses. Here we describe the development of a peptide-based compound for systemic use, sHVF18, expressing the evolutionarily conserved innate structural fold of natural TCPs. Using a combination of structure- and in silico-based design, nuclear magnetic resonance spectroscopy, biophysics, mass spectrometry, cellular, and in vivo studies, we here elucidate the structure, CD14 interactions, protease stability, transcriptome profiling, and therapeutic efficacy of sHVF18. The designed peptide displays a conformationally stabilized, protease resistant active innate fold and targets the LPS-binding groove of CD14. In vivo, it shows therapeutic efficacy in experimental models of endotoxin shock in mice and pigs and increases survival in mouse models of systemic polymicrobial infection. The results provide a drug class based on Nature´s own anti-infective principles.

A computationally designed antigen eliciting broad humoral responses against SARS-CoV-2 and related sarbecoviruses

The threat of spillovers of coronaviruses associated with the severe acute respiratory syndrome (SARS) from animals to humans necessitates vaccines that offer broader protection from sarbecoviruses. By leveraging a viral-genome-informed computational method for selecting immune-optimized and structurally engineered antigens, here we show that a single antigen based on the receptor binding domain of the spike protein of sarbecoviruses elicits broad humoral responses against SARS-CoV-1, SARS-CoV-2, WIV16 and RaTG13 in mice, rabbits and guinea pigs. When administered as a DNA immunogen or by a vector based on a modified vaccinia virus Ankara, the optimized antigen induced vaccine protection from the Delta variant of SARS-CoV-2 in mice genetically engineered to express angiotensin-converting enzyme 2 and primed by a viral-vector vaccine (AZD1222) against SARS-CoV-2. A vaccine formulation incorporating mRNA coding for the optimized antigen further validated its broad immunogenicity. Vaccines that elicit broad immune responses across subgroups of coronaviruses may counteract the threat of zoonotic spillovers of betacoronaviruses.

Stilbenoid compounds inhibit NF-κB-mediated inflammatory responses in the Drosophila intestine

Introduction

Stilbenoid compounds have been described to have anti-inflammatory properties in animal models in vivo, and have been shown to inhibit Ca2+-influx through the transient receptor potential ankyrin 1 (TrpA1).

Methods

To study how stilbenoid compounds affect inflammatory signaling in vivo, we have utilized the fruit fly, Drosophila melanogaster, as a model system. To induce intestinal inflammation in the fly, we have fed flies with the intestinal irritant dextran sodium sulphate (DSS).

Results

We found that DSS induces severe changes in the bacteriome of the Drosophila intestine, and that this dysbiosis causes activation of the NF-κB transcription factor Relish. We have taken advantage of the DSS-model to study the anti-inflammatory properties of the stilbenoid compounds pinosylvin (PS) and pinosylvin monomethyl ether (PSMME). With the help of in vivo approaches, we have identified PS and PSMME to be transient receptor ankyrin 1 (TrpA1)-dependent antagonists of NF-κB-mediated intestinal immune responses in Drosophila. We have also computationally predicted the putative antagonist binding sites of these compounds at Drosophila TrpA1.

Discussion

Taken together, we show that the stilbenoids PS and PSMME have anti-inflammatory properties in vivo in the intestine and can be used to alleviate chemically induced intestinal inflammation in Drosophila.

Modeling the Characteristic Residues of Chlorophyll f Synthase (ChlF) from Halomicronema hongdechloris to Determine Its Reaction Mechanism

Photosystem II (PSII) is a quinone-utilizing photosynthetic system that converts light energy into chemical energy and catalyzes water splitting. PsbA (D1) and PsbD (D2) are the core subunits of the reaction center that provide most of the ligands to redox-active cofactors and exhibit photooxidoreductase activities that convert quinone and water into quinol and dioxygen. The performed analysis explored the putative uncoupled electron transfer pathways surrounding P680+ induced by far-red light (FRL) based on photosystem II (PSII) complexes containing substituted D1 subunits in Halomicronema hongdechloris. Chlorophyll f-synthase (ChlF) is a D1 protein paralog. Modeling PSII-ChlF complexes determined several key protein motifs of ChlF. The PSII complexes included a dysfunctional Mn4CaO5 cluster where ChlF replaced the D1 protein. We propose the mechanism of chlorophyll f synthesis from chlorophyll a via free radical chemistry in an oxygenated environment created by over-excited pheophytin a and an inactive water splitting reaction owing to an uncoupled Mn4CaO5 cluster in PSII-ChlF complexes. The role of ChlF in the formation of an inactive PSII reaction center is under debate, and putative mechanisms of chlorophyll f biosynthesis are discussed.

Construction of Histone-Protein Complex Structures by Peptide Growing

The structures of histone complexes are master keys to epigenetics. Linear histone peptide tails often bind to shallow pockets of reader proteins via weak interactions, rendering their structure determination challenging. In the present study, a new protocol, PepGrow, is introduced. PepGrow uses docked histone fragments as seeds and grows the full peptide tails in the reader-binding pocket, producing atomic-resolution structures of histone-reader complexes. PepGrow is able to handle the flexibility of histone peptides, and it is demonstrated to be more efficient than linking pre-docked peptide fragments. The new protocol combines the advantages of popular program packages and allows fast generation of solution structures. AutoDock, a force-field-based program, is used to supply the docked peptide fragments used as structural seeds, and the building algorithm of Modeller is adopted and tested as a peptide growing engine. The performance of PepGrow is compared to ten other docking methods, and it is concluded that in situ growing of a ligand from a seed is a viable strategy for the production of complex structures of histone peptides at atomic resolution.

Protons taken hostage: Dynamic H-bond networks of the pH-sensing GPR68

Proton-sensing G Protein Coupled Receptors (GPCRs) sense changes in the extracellular pH to effect cell signaling for cellular homeostasis. They tend to be overexpressed in solid tumors associated with acidic extracellular pH, and are of direct interest as drug targets. How proton-sensing GPCRs sense extracellular acidification and activate upon protonation change is important to understand, because it may guide the design of therapeutics. Lack of publicly available experimental structures make it challenging to discriminate between conflicting mechanisms proposed for proton-binding, as main roles have been assigned to either an extracellular histidine cluster or to an internal carboxylic triad. Here we present a protocol to derive and evaluate structural models of the proton-sensing GPR68. This approach integrates state-of-the-art homology modeling with microsecond-timescale atomistic simulations, and with a detailed assessment of the compatibility of the structural models with known structural features of class A GPCRs. To decipher structural elements of potential interest for protonation-coupled conformational changes of GPR68, we used the best-compatible model as a starting point for independent atomistic simulations of GPR68 with different protonation states, and graph computations to characterize the response of GPR68 to changes in protonation. We found that GPR68 hosts an extended hydrogen-bond network that inter-connects the extracellular histidine cluster to the internal carboxylic triad, and which can even reach groups at the cytoplasmic G-protein binding site. Taken together, results suggest that GPR68 relies on dynamic, hydrogen-bond networks to inter-connect extracellular and internal proton-binding sites, and to elicit conformational changes at the cytoplasmic G-protein binding site.

Identification of novel inhibitors of Neisseria gonorrhoeae MurI using homology modeling, structure-based pharmacophore, molecular docking, and molecular dynamics simulation-based approach

MurI is one of the most significant role players in the biosynthesis of the peptidoglycan layer in Neisseria gonorrhoeae (Ng). We attempted to highlight the structural and functional relationship between Ng-MurI and D-glutamate to design novel molecules targeting this interaction. The three-dimensional (3D) model of the protein was constructed by homology modeling and the quality and consistency of generated model were assessed. The binding site of the protein was identified by molecular docking studies and a pharmacophore was identified using the interactions of the control ligand. The structure-based pharmacophore model was validated and employed for high-throughput virtual screening and molecular docking to identify novel Ng-MurI inhibitors. Finally, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with the substrate glutamate and novel molecules facilitated us to confirm the stability of the protein-ligand docked complexes. The 100 ns MD simulations of the potential lead compounds with protein confirmed that the modeled complexes were stable. This study identifies novel potential compounds with good fitness and docking scores, which made the interactions of biological significance within the protein active site. Hence, the identified compounds may act as new leads to design and develop Ng-MurI inhibitors.Communicated by Ramaswamy H. Sarma.

M Protein from Dengue virus oligomerizes to pentameric channel protein: in silico analysis study

The Dengue virus M protein is a 75 amino acid polypeptide with two helical transmembranes (TM). The TM domain oligomerizes to form an ion channel, facilitating viral release from the host cells. The M protein has a critical role in the virus entry and life cycle, making it a potent drug target. The oligomerization of the monomeric protein was studied using ab initio modeling and molecular dynamics (MD) simulation in an implicit membrane environment. The representative structures obtained showed pentamer as the most stable oligomeric state, resembling an ion channel. Glutamic acid, threonine, serine, tryptophan, alanine, isoleucine form the pore-lining residues of the pentameric channel, conferring an overall negative charge to the channel with approximate length of 51.9 Å. Residue interaction analysis (RIN) for M protein shows that Ala94, Leu95, Ser112, Glu124, and Phe155 are the central hub residues representing the physicochemical interactions between domains. The virtual screening with 165 different ion channel inhibitors from the ion channel library shows monovalent ion channel blockers, namely lumacaftor, glipizide, gliquidone, glisoxepide, and azelnidipine to be the inhibitors with high docking scores. Understanding the three-dimensional structure of M protein will help design therapeutics and vaccines for Dengue infection.

Improving Geometric Validation Metrics and Ensuring Consistency with Experimental Data through TrioSA: An NMR Refinement Protocol

Protein model refinement a the crucial step in improving the quality of a predicted protein model. This study presents an NMR refinement protocol called TrioSA (torsion-angle and implicit-solvation-optimized simulated annealing) that improves the accuracy of backbone/side-chain conformations and the overall structural quality of proteins. TrioSA was applied to a subset of 3752 solution NMR protein structures accompanied by experimental NMR data: distance and dihedral angle restraints. We compared the initial NMR structures with the TrioSA-refined structures and found significant improvements in structural quality. In particular, we observed a reduction in both the maximum and number of NOE (nuclear Overhauser effect) violations, indicating better agreement with experimental NMR data. TrioSA improved geometric validation metrics of NMR protein structure, including backbone accuracy and the secondary structure ratio. We evaluated the contribution of each refinement element and found that the torsional angle potential played a significant role in improving the geometric validation metrics. In addition, we investigated protein-ligand docking to determine if TrioSA can improve biological outcomes. TrioSA structures exhibited better binding prediction compared to the initial NMR structures. This study suggests that further development and research in computational refinement methods could improve biomolecular NMR structural determination.

Paralog-dependent isogenic cell assay cascade generates highly selective SLC16A3 inhibitors

Despite being considered druggable and attractive therapeutic targets, most of the solute carrier (SLC) membrane transporters remain pharmacologically underexploited. One of the reasons for this is a lack of reliable chemical screening assays, made difficult by functional redundancies among SLCs. In this study we leveraged synthetic lethality between the lactate transporters SLC16A1 and SLC16A3 in a screening strategy that we call paralog-dependent isogenic cell assay (PARADISO). The system involves five isogenic cell lines, each dependent on various paralog genes for survival/fitness, arranged in a screening cascade tuned for the identification of SLC16A3 inhibitors. We screened a diversity-oriented library of ∼90,000 compounds and further developed our hits into slCeMM1, a paralog-selective and potent SLC16A3 inhibitor. By implementing chemoproteomics, we showed that slCeMM1 is selective also at the proteome-wide level, thus fulfilling an important criterion for chemical probes. This study represents a framework for the development of specific cell-based drug discovery assays.

Impaired tissue homing by the Ikzf3N159S variant is mediated by interfering with Ikaros function

AIOLOS, encoded by IKZF3, is a member of the IKZF family of proteins that plays an important role in regulating late B-cell differentiation. Human individuals heterozygous for the AIOLOS p.N160S variant displayed impaired humoral immune responses as well as impaired B and T cell development. We have previously reported that a mouse strain harboring an Ikzf3N159S allele that corresponds to human IKZF3N160S recapitulated immune-deficient phenotypes, such as impaired B cell development and loss of CD23 expression. In this study, we investigated the effect of the Ikzf3N159S variant and found that B1a cell development was impaired in Ikzf3N159S/N159S mice. In addition, CD62L expression was severely decreased in both B and T lymphocytes by the Ikzf3N159S mutation, in a dose-dependent manner. Mixed bone marrow chimera experiments have revealed that most immunodeficient phenotypes, including low CD62L expression, occur in intrinsic cells. Interestingly, while Ikzf3N159S/N159S lymphocytes were still present in the spleen, they were completely outcompeted by control cells in the lymph nodes, suggesting that the capacity for homing or retention in the lymph nodes was lost due to the Ikzf3N159S mutation. The homing assay confirmed severely decreased homing abilities to lymph nodes of Ikzf3N159S/N159S B and T lymphocytes but selective enrichment of CD62L expressing Ikzf3N159S/N159S lymphocytes in lymph nodes. This finding suggests that impaired CD62L expression is the major reason for the impaired homing capacity caused by the Ikzf3N159S mutation. Interestingly, an excess amount of Ikaros, but not Aiolos, restored CD62L expression in Ikzf3N159S/N159S B cells. Together with the loss of CD62L expression due to Ikaros deficiency, the AiolosN159S mutant protein likely interferes with Ikaros function through heterodimerization, at least in activating the Sell gene encoding CD62L expression. Thus, our results revealed that AiolosN159S causes some immunodeficient phenotypes via the pathogenesis referred to as the heterodimeric interference as observed for AiolosG158R variant.

Structural basis of allosteric modulation of metabotropic glutamate receptor activation and desensitization

The metabotropic glutamate receptors (mGluRs) are neuromodulatory family C G protein coupled receptors which assemble as dimers and allosterically couple extracellular ligand binding domains (LBDs) to transmembrane domains (TMDs) to drive intracellular signaling. Pharmacologically, mGluRs can be targeted either at the LBDs by glutamate and synthetic orthosteric compounds or at the TMDs by allosteric modulators. Despite the potential of allosteric TMD-targeting compounds as therapeutics, an understanding of the functional and structural basis of their effects on mGluRs is limited. Here we use a battery of approaches to dissect the distinct functional and structural effects of orthosteric versus allosteric ligands. We find using electrophysiological and live cell imaging assays that both agonists and positive allosteric modulators (PAMs) can drive activation and desensitization of mGluRs. The effects of PAMs are pleiotropic, including both the ability to boost the maximal response to orthosteric agonists and to serve independently as desensitization-biased agonists across mGluR subtypes. Conformational sensors reveal PAM-driven inter-subunit re-arrangements at both the LBD and TMD. Motivated by this, we determine cryo-electron microscopy structures of mGluR3 in the presence of either an agonist or antagonist alone or in combination with a PAM. These structures reveal PAM-driven re-shaping of intra- and inter-subunit conformations and provide evidence for a rolling TMD dimer interface activation pathway that controls G protein and beta-arrestin coupling.

Highlights

-Agonists and PAMs drive mGluR activation, desensitization, and endocytosis-PAMs are desensitization-biased and synergistic with agonists-Four combinatorial ligand conditions reveal an ensemble of full-length mGluR structures with novel interfaces-Activation and desensitization involve rolling TMD interfaces which are re-shaped by PAM.

Fluorescent protein lifetimes report densities and phases of nuclear condensates during embryonic stem-cell differentiation

Fluorescent proteins (FP) are frequently used for studying proteins inside cells. In advanced fluorescence microscopy, FPs can report on additional intracellular variables. One variable is the local density near FPs, which can be useful in studying densities within cellular bio-condensates. Here, we show that a reduction in fluorescence lifetimes of common monomeric FPs reports increased levels of local densities. We demonstrate the use of this fluorescence-based variable to report the distribution of local densities within heterochromatin protein 1α (HP1α) in mouse embryonic stem cells (ESCs), before and after early differentiation. We find that local densities within HP1α condensates in pluripotent ESCs are heterogeneous and cannot be explained by a single liquid phase. Early differentiation, however, induces a change towards a more homogeneous distribution of local densities, which can be explained as a liquid-like phase. In conclusion, we provide a fluorescence-based method to report increased local densities and apply it to distinguish between homogeneous and heterogeneous local densities within bio-condensates.

Acetylcholinesterase Inhibition Activity of Hippeastrum papilio (Ravenna) Van Scheepen (Amaryllidaceae) Using Zebrafish Brain Homogenates

The Amaryllidaceae family constitutes an interesting source of exclusive alkaloids with a broad spectrum of biological activity. Galanthamine, the most relevant one, has been commercialized for the palliative treatment of Alzheimer's disease symptoms since 2001 due to its potential as an acetylcholinesterase (AChE) inhibitor. In vitro screenings against AChE by applying different Amaryllidaceae species and alkaloids have been reported in the literature; however, they are usually carried out using purified market enzymes. The main goal of this work is to evaluate the AChE inhibitory potential of Hippeastrum papilio (Amaryllidaceae) extracts using zebrafish brain homogenates. The biological assays show that the H. papilio bulb extracts present an interesting AChE inhibitory activity in comparison with the positive reference control galanthamine (IC50 values of 1.20 ± 0.10 and 0.79 ± 0.15 μg/mL, respectively). The chemical profile of H. papilio shows that this species has a high amount of galanthamine, which may contribute to the inhibitory effect on AChE activity of zebrafish brains. Computational experiments were used to build the model for zebrafish AChE and to evaluate the interactions between galanthamine and the enzymic active site. This work suggests that zebrafish could represent an important model in the search for bioactive molecules from the Amaryllidaceae family for the treatment of Alzheimer's disease.

A dynamic biomimetic model of the membrane-bound CD4-CD3-TCR complex during pMHC disengagement

The coordinated (dis)engagement of the membrane-bound T cell receptor (TCR)-CD3-CD4 complex from the peptide-major histocompatibility complex (pMHC) is fundamental to TCR signal transduction and T cell effector function. As such, an atomic-scale understanding would not only enhance our basic understanding of the adaptive immune response but would also accelerate the rational design of TCRs for immunotherapy. In this study, we explore the impact of the CD4 coreceptor on the TCR-pMHC (dis)engagement by constructing a molecular-level biomimetic model of the CD3-TCR-pMHC and CD4-CD3-TCR-pMHC complexes within a lipid bilayer. After allowing the system complexes to equilibrate (engage), we use steered molecular dynamics to dissociate (disengage) the pMHC. We find that 1) the CD4 confines the pMHC closer to the T cell by 1.8 nm at equilibrium; 2) CD4 confinement shifts the TCR along the MHC binding groove engaging a different set of amino acids and enhancing the TCR-pMHC bond lifetime; 3) the CD4 translocates under load increasing the interaction strength between the CD4-pMHC, CD4-TCR, and CD4-CD3; and 4) upon dissociation, the CD3-TCR complex undergoes structural oscillation and increased energetic fluctuation between the CD3-TCR and CD3-lipids. These atomic-level simulations provide mechanistic insight on how the CD4 coreceptor impacts TCR-pMHC (dis)engagement. More specifically, our results provide further support (enhanced bond lifetime) for a force-dependent kinetic proofreading model and identify an alternate set of amino acids in the TCR that dominate the TCR-pMHC interaction and could thus impact the design of TCRs for immunotherapy.

Sequestration of histidine kinases by non-cognate response regulators establishes a threshold level of stimulation for bacterial two-component signaling

Bacterial two-component systems (TCSs) consist of a sensor histidine kinase (HK) that perceives a specific signal, and a cognate response regulator (RR) that modulates the expression of target genes. Positive autoregulation improves TCS sensitivity to stimuli, but may trigger disproportionately large responses to weak signals, compromising bacterial fitness. Here, we combine experiments and mathematical modelling to reveal a general design that prevents such disproportionate responses: phosphorylated HKs (HK~Ps) can be sequestered by non-cognate RRs. We study five TCSs of Mycobacterium tuberculosis and find, for all of them, non-cognate RRs that show higher affinity than cognate RRs for HK~Ps. Indeed, in vitro assays show that HK~Ps preferentially bind higher affinity non-cognate RRs and get sequestered. Mathematical modelling indicates that this sequestration would introduce a 'threshold' stimulus strength for eliciting responses, thereby preventing responses to weak signals. Finally, we construct tunable expression systems in Mycobacterium bovis BCG to show that higher affinity non-cognate RRs suppress responses in vivo.

Two point mutations in protocadherin-1 disrupt hantavirus recognition and afford protection against lethal infection

Andes virus (ANDV) and Sin Nombre virus (SNV) are the etiologic agents of severe hantavirus cardiopulmonary syndrome (HCPS) in the Americas for which no FDA-approved countermeasures are available. Protocadherin-1 (PCDH1), a cadherin-superfamily protein recently identified as a critical host factor for ANDV and SNV, represents a new antiviral target; however, its precise role remains to be elucidated. Here, we use computational and experimental approaches to delineate the binding surface of the hantavirus glycoprotein complex on PCDH1's first extracellular cadherin repeat domain. Strikingly, a single amino acid residue in this PCDH1 surface influences the host species-specificity of SNV glycoprotein-PCDH1 interaction and cell entry. Mutation of this and a neighboring residue substantially protects Syrian hamsters from pulmonary disease and death caused by ANDV. We conclude that PCDH1 is a bona fide entry receptor for ANDV and SNV whose direct interaction with hantavirus glycoproteins could be targeted to develop new interventions against HCPS.

Unveiling the interdomain dynamics of type II DNA topoisomerase through all-atom simulations: Implications for understanding its catalytic cycle

Type IIA DNA topoisomerases are complex molecular nanomachines that manage topological states of the DNA molecule in the cell and play a crucial role in cellular processes such as cell division and transcription. They are also established targets of cancer chemotherapy. Starting from the available crystal structure of a fully catalytic topoisomerase IIA homodimer from Saccharomyces cerevisiae, we constructed three states of this molecular motor primarily changing the configurations of the DNA segment bound in the DNA gate and performed μs-long all-atom molecular simulations. A comprehensive analysis revealed a sliding motion within the DNA gate and a teamwork between the N-gate and DNA gate that may be associated with the necessary molecular events that allow passage of the T-segment of DNA. The observed movement of the ATPase dimer relative to the DNA domain was reflected in different interaction patterns between the K-loops of the transducer domain and the B-A-B form of the bound DNA. Based on the obtained results, we mapped simulated configurations to the structures in the proposed catalytic cycle through which type IIA topoisomerases exert their function and discussed the possible transition events. The results extend our understanding of the mechanism of action of type IIA topoisomerases and provide an atomistic interpretation of some of the observed features of these molecular motors.

Evolutional insights into the interaction between Rab7 and RILP in lysosome motility

Lysosome motility is critical for the cellular function. However, Rab7-related transport elements showed genetic differences between vertebrates and invertebrates, making the mechanism of lysosomal motility mysterious. We suggested that Rab7 interacted with RILP as a feature of highly evolved organisms since they could interact with each other in Spodoptera frugiperda but not in Drosophila melanogaster. The N-terminus of Sf-RILP was identified to be necessary for their interaction, and Glu61 was supposed to be the key point for the stability of the interaction. A GC-rich domain on the C-terminal parts of Sf-RILP hampered the expression of Sf-RILP and its interaction with Sf-Rab7. Although the corresponding vital amino acids in the mammalian model at the C-terminus of Sf-RILP turned to be neutral, the C-terminus would also help with the homologous interactions between RILP fragments in insects. The significantly different interactions in invertebrates shed light on the biodiversity and complexity of lysosomal motility.

RNA 3D Structure Prediction: Progress and Perspective

Ribonucleic acid (RNA) molecules play vital roles in numerous important biological functions such as catalysis and gene regulation. The functions of RNAs are strongly coupled to their structures or proper structure changes, and RNA structure prediction has been paid much attention in the last two decades. Some computational models have been developed to predict RNA three-dimensional (3D) structures in silico, and these models are generally composed of predicting RNA 3D structure ensemble, evaluating near-native RNAs from the structure ensemble, and refining the identified RNAs. In this review, we will make a comprehensive overview of the recent advances in RNA 3D structure modeling, including structure ensemble prediction, evaluation, and refinement. Finally, we will emphasize some insights and perspectives in modeling RNA 3D structures.

Canonical and Alternative Auxin Signaling Systems in Mono-, Di-, and Tetraploid Potatoes

It has long been known that the phytohormone auxin plays a promoting role in tuber formation and stress tolerance in potatoes. Our study aimed to identify and characterize the complete sets of auxin-related genes that presumably constitute the entire auxin signaling system in potato (Solanum tuberosum L.). The corresponding genes were retrieved from sequenced genomes of the doubled monoploid S. tuberosum DM1-3-516-R44 (DM) of the Phureja group, the heterozygous diploid line RH89-039-16 (RH), and the autotetraploid cultivar Otava. Both canonical and noncanonical auxin signaling pathways were considered. Phylogenetic and domain analyses of deduced proteins were supplemented by expression profiling and 3D molecular modeling. The canonical and ABP1-mediated pathways of auxin signaling appeared to be well conserved. The total number of potato genes/proteins presumably involved in canonical auxin signaling is 46 and 108 in monoploid DM and tetraploid Otava, respectively. Among the studied potatoes, spectra of expressed genes obviously associated with auxin signaling were partly cultivar-specific and quite different from analogous spectrum in Arabidopsis. Most of the noncanonical pathways found in Arabidopsis appeared to have low probability in potato. This was equally true for all cultivars used irrespective of their ploidy. Thus, some important features of the (noncanonical) auxin signaling pathways may be variable and species-specific.

Individual bat virome analysis reveals co-infection and spillover among bats and virus zoonotic potential

Bats are reservoir hosts for many zoonotic viruses. Despite this, relatively little is known about the diversity and abundance of viruses within individual bats, and hence the frequency of virus co-infection and spillover among them. We characterize the mammal-associated viruses in 149 individual bats sampled from Yunnan province, China, using an unbiased meta-transcriptomics approach. This reveals a high frequency of virus co-infection (simultaneous infection of bat individuals by multiple viral species) and spillover among the animals studied, which may in turn facilitate virus recombination and reassortment. Of note, we identify five viral species that are likely to be pathogenic to humans or livestock, based on phylogenetic relatedness to known pathogens or in vitro receptor binding assays. This includes a novel recombinant SARS-like coronavirus that is closely related to both SARS-CoV and SARS-CoV-2. In vitro assays indicate that this recombinant virus can utilize the human ACE2 receptor such that it is likely to be of increased emergence risk. Our study highlights the common occurrence of co-infection and spillover of bat viruses and their implications for virus emergence.

Effects of Bacterial Metabolites on the Wnt4 Protein in Dental-Pulp-Stem-Cells-Based Endodontic Pulpitis Treatment

Porphyromonas gingivalis is associated with endodontic pulpitis, causing damage to the dental pulp, leading to severe pain and a decline in quality of life. Regenerative pulp treatments using dental pulp stem cells (DPSCs) can be hindered by interactions between DPSCs and the infecting bacteria. The protein WNT family member 4 (Wnt4) plays a critical role in the differentiation of DPSCs and the regeneration of odontogenic tissue. However, the specific influence of P. gingivalis on Wnt4 remains unclear. In this study, we employed a computational approach to investigate the underlying mechanisms through which P. gingivalis-produced metabolites inhibit the Wnt4 protein, thereby diminishing the regenerative potential and therapeutic efficacy of odontogenic tissue. Among the metabolites examined, C₂₉H₄₆N₇O₁₈P₃S^-4 exhibited the strongest inhibitory effect on the Wnt4 protein, as evidenced by the lowest binding energy score of -6782 kcal/mol. Molecular dynamic simulation trajectories revealed that the binding of C₂₉H₄₆N₇O₁₈P₃S^-4 significantly altered the structural dynamics and stability of the Wnt4 protein. These alterations in protein trajectories may have implications for the molecular function of Wnt4 and its associated pathways. Overall, our findings shed light on the inhibitory impact of P. gingivalis-produced metabolites on the Wnt4 protein. Further in vitro, in vivo, and clinical studies are necessary to validate and expand upon our findings.

Conservation and divergence of the G-interfaces of Drosophila melanogaster septins

Septins possess a conserved guanine nucleotide-binding (G) domain that participates in the stabilization of organized hetero-oligomeric complexes which assemble into filaments, rings and network-like structures. The fruit fly, Drosophila melanogaster, has five such septin genes encoding Sep1, Sep2, Sep4, Sep5 and Pnut. Here, we report the crystal structure of the heterodimer formed between the G-domains of Sep1 and Sep2, the first from an insect to be described to date. A G-interface stabilizes the dimer (in agreement with the expected arrangement for the Drosophila hexameric particle) and this bears significant resemblance to its human counterparts, even down to the level of individual amino acid interactions. On the other hand, a model for the G-interface formed between the two copies of Pnut which occupy the centre of the hexamer, shows important structural differences, including the loss of a highly favourable bifurcated salt-bridge network. Whereas wild-type Pnut purifies as a monomer, the reintroduction of the salt-bridge network results in stabilizing the dimeric interface in solution as shown by size exclusion chromatography and thermal stability measurements. Adaptive steered molecular dynamics reveals an unzipping mechanism for dimer dissociation which initiates at a point of electrostatic repulsion within the switch II region. Overall, the data contribute to a better understanding of the molecular interactions involved in septin assembly/disassembly.

Computational design of the temperature optimum of an enzyme reaction

Cold-adapted enzymes are characterized both by a higher catalytic activity at low temperatures and by having their temperature optimum down-shifted, compared to mesophilic orthologs. In several cases, the optimum does not coincide with the onset of protein melting but reflects some other type of inactivation. In the psychrophilic α-amylase from an Antarctic bacterium, the inactivation is thought to originate from a specific enzyme-substrate interaction that breaks around room temperature. Here, we report a computational redesign of this enzyme aimed at shifting its temperature optimum upward. A set of mutations designed to stabilize the enzyme-substrate interaction were predicted by computer simulations of the catalytic reaction at different temperatures. The predictions were verified by kinetic experiments and crystal structures of the redesigned α-amylase, showing that the temperature optimum is indeed markedly shifted upward and that the critical surface loop controlling the temperature dependence approaches the target conformation observed in a mesophilic ortholog.