Comparative protein structure modeling using MODELLER

The Structure and Function of the Bacterial Osmotically Inducible Protein Y

The ability to adapt to osmotically diverse and fluctuating environments is critical to the survival and resilience of bacteria that colonize the human gut and urinary tract. Environmental stress often provides cross-protection against other challenges and increases antibiotic tolerance of bacteria. Thus, it is critical to understand how E. coli and other microbes survive and adapt to stress conditions. The osmotically inducible protein Y (OsmY) is significantly upregulated in response to hypertonicity. Yet its function remains unknown for decades. We determined the solution structure and dynamics of OsmY by nuclear magnetic resonance spectroscopy, which revealed that the two Bacterial OsmY and Nodulation (BON) domains of the protein are flexibly linked under low- and high-salinity conditions. In-cell solid-state NMR further indicates that there are no gross structural changes in OsmY as a function of osmotic stress. Using cryo-electron and super-resolution fluorescence microscopy, we show that OsmY attenuates plasmolysis-induced structural changes in E. coli and improves the time to growth resumption after osmotic upshift. Structure-guided mutational and functional studies demonstrate that exposed hydrophobic residues in the BON1 domain are critical for the function of OsmY. We find no evidence for membrane interaction of the BON domains of OsmY, contrary to current assumptions. Instead, at high ionic strength, we observe an interaction with the water channel, AqpZ. Thus, OsmY does not play a simple structural role in E. coli but may influence a cascade of osmoregulatory functions of the cell.

Nanomechanical collective vibration of SARS-CoV-2 spike proteins

The development of effective therapeutics against COVID-19 requires a thorough understanding of the receptor recognition mechanism of the SARS-CoV-2 spike (S) protein. Here the multidomain collective dynamics on the trimer of the spike protein has been analyzed using normal mode analysis (NMA). A common nanomechanical profile was identified in the spike proteins of SARS-CoV-2 and its variants. The profile involves collective vibrations of the receptor-binding domain (RBD) and the N-terminal domain (NTD), which may mediate the physical interaction process. Quantitative analysis of the collective modes suggests a nanomechanical property involving large-scale conformational changes, which explains the difference in receptor binding affinity among different variants. These results support the use of intrinsic global dynamics as a valuable perspective for studying the allosteric and functional mechanisms of the S protein. This approach also provides a low-cost theoretical toolkit for screening potential pathogenic mutations and drug targets.

Identification of phytochemical inhibitors targeting phosphate acetyltransferase of Mycoplasma genitalium: insights from virtual screening and molecular dynamics studies

Mycoplasma genitalium (M. genitalium) has evolved as a superbug, and the developing antimicrobial resistance with just a few treatment options available is an imminent concern. Due to the emergence of antibiotic resistance, a new antibiotic class or medications are required to combat this pathogen. The phosphate acetyltransferase (PTA) enzyme can be a suitable drug target which is essential for M. genitalium survival and involves in acetate metabolism. To efficiently find potent inhibitors, structure-based drug design approaches targeting the PTA of M. genitalium have been established. In this study, the three most potent phytochemical inhibitors were predicted from virtual screening and these are sitostanyl ferulate, beta-sitosterol-beta-D-glucoside, and brassinolide, with binding energies of - 9.66, - 9.60, and - 9.48 kcal/mol, respectively. The active site residues Thr-125, Arg-300, Ser-299, Tyr-272, and Lys-273 appear to be critical in binding the three predicted potent inhibitors. The results of the molecular dynamics study indicate that the three predicted phytochemical inhibitors have formed stable bonds with PTA. Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) was utilized for the estimation of binding free energy of PTA-phytochemical complexes. Taken together, the findings of our computational work might aid in the development of possible potential drugs to treat and ameliorate the severity of M. genitalium infection.

Firefly toxin lucibufagins evolved after the origin of bioluminescence

Fireflies were believed to originally evolve their novel bioluminescence as warning signals to advertise their toxicity to predators, which was later adopted in adult mating. Although the evolution of bioluminescence has been investigated extensively, the warning signal hypothesis of its origin has not been tested. In this study, we test this hypothesis by systematically determining the presence or absence of firefly toxin lucibufagins (LBGs) across firefly species and inferring the time of origin of LBGs. We confirm the presence of LBGs in the subfamily Lampyrinae, but more importantly, we reveal the absence of LBGs in other lineages, including the subfamilies of Luciolinae, Ototretinae, and Psilocladinae, two incertae sedis lineages, and the Rhagophthalmidae family. Ancestral state reconstructions for LBGs based on firefly phylogeny constructed using genomic data suggest that the presence of LBGs in the common ancestor of the Lampyrinae subfamily is highly supported but unsupported in more ancient nodes, including firefly common ancestors. Our results suggest that firefly LBGs probably evolved much later than the evolution of bioluminescence. We thus conclude that firefly bioluminescence did not originally evolve as direct warning signals for toxic LBGs and advise that future studies should focus on other hypotheses. Moreover, LBG toxins are known to directly target and inhibit the α subunit of Na+, K+-ATPase (ATPα). We further examine the effects of amino acid substitutions in firefly ATPα on its interactions with LBGs. We find that ATPα in LBG-containing fireflies is relatively insensitive to LBGs, which suggests that target-site insensitivity contributes to LBG-containing fireflies' ability to deal with their own toxins.

Targeting the Plasmodium falciparum UCHL3 ubiquitin hydrolase using chemically constrained peptides

The ubiquitin-proteasome system is essential to all eukaryotes and has been shown to be critical to parasite survival as well, including Plasmodium falciparum, the causative agent of the deadliest form of malarial disease. Despite the central role of the ubiquitin-proteasome pathway to parasite viability across its entire life-cycle, specific inhibitors targeting the individual enzymes mediating ubiquitin attachment and removal do not currently exist. The ability to disrupt P. falciparum growth at multiple developmental stages is particularly attractive as this could potentially prevent both disease pathology, caused by asexually dividing parasites, as well as transmission which is mediated by sexually differentiated parasites. The deubiquitinating enzyme PfUCHL3 is an essential protein, transcribed across both human and mosquito developmental stages. PfUCHL3 is considered hard to drug by conventional methods given the high level of homology of its active site to human UCHL3 as well as to other UCH domain enzymes. Here, we apply the RaPID mRNA display technology and identify constrained peptides capable of binding to PfUCHL3 with nanomolar affinities. The two lead peptides were found to selectively inhibit the deubiquitinase activity of PfUCHL3 versus HsUCHL3. NMR spectroscopy revealed that the peptides do not act by binding to the active site but instead block binding of the ubiquitin substrate. We demonstrate that this approach can be used to target essential protein-protein interactions within the Plasmodium ubiquitin pathway, enabling the application of chemically constrained peptides as a novel class of antimalarial therapeutics.

C-type inactivation and proton modulation mechanisms of the TASK3 channel

The TWIK-related acid-sensitive K+ channel 3 (TASK3) belongs to the two-pore domain (K2P) potassium channel family, which regulates cell excitability by mediating a constitutive "leak" potassium efflux in the nervous system. Extracellular acidification inhibits TASK3 channel, but the molecular mechanism by which channel inactivation is coupled to pH decrease remains unclear. Here, we report the cryo-electron microscopy structures of human TASK3 at neutral and acidic pH. Structural comparison revealed selectivity filter (SF) rearrangements upon acidification, characteristic of C-type inactivation, but with a unique structural basis. The extracellular mouth of the SF was prominently dilated and simultaneously blocked by a hydrophobic gate. His98 protonation shifted the conformational equilibrium between the conductive and C-type inactivated SF toward the latter by engaging a cation-π interaction with Trp78, consistent with molecular dynamics simulations and electrophysiological experiments. Our work illustrated how TASK3 is gated in response to extracellular pH change and implies how physiological stimuli might directly modulate the C-type gating of K2P channels.

Interplay between charge distribution and DNA in shaping HP1 paralog phase separation and localization

The heterochromatin protein 1 (HP1) family is a crucial component of heterochromatin with diverse functions in gene regulation, cell cycle control, and cell differentiation. In humans, there are three paralogs, HP1α, HP1β, and HP1γ, which exhibit remarkable similarities in their domain architecture and sequence properties. Nevertheless, these paralogs display distinct behaviors in liquid-liquid phase separation (LLPS), a process linked to heterochromatin formation. Here, we employ a coarse-grained simulation framework to uncover the sequence features responsible for the observed differences in LLPS. We highlight the significance of the net charge and charge patterning along the sequence in governing paralog LLPS propensities. We also show that both highly conserved folded and less-conserved disordered domains contribute to the observed differences. Furthermore, we explore the potential co-localization of different HP1 paralogs in multicomponent assemblies and the impact of DNA on this process. Importantly, our study reveals that DNA can significantly reshape the stability of a minimal condensate formed by HP1 paralogs due to competitive interactions of HP1α with HP1β and HP1γ versus DNA. In conclusion, our work highlights the physicochemical nature of interactions that govern the distinct phase-separation behaviors of HP1 paralogs and provides a molecular framework for understanding their role in chromatin organization.

The Molecular Architecture of the Nuclear Basket

The nuclear pore complex (NPC) is the sole mediator of nucleocytoplasmic transport. Despite great advances in understanding its conserved core architecture, the peripheral regions can exhibit considerable variation within and between species. One such structure is the cage-like nuclear basket. Despite its crucial roles in mRNA surveillance and chromatin organization, an architectural understanding has remained elusive. Using in-cell cryo-electron tomography and subtomogram analysis, we explored the NPC's structural variations and the nuclear basket across fungi (yeast; S. cerevisiae), mammals (mouse; M. musculus), and protozoa (T. gondii). Using integrative structural modeling, we computed a model of the basket in yeast and mammals that revealed how a hub of Nups in the nuclear ring binds to basket-forming Mlp/Tpr proteins: the coiled-coil domains of Mlp/Tpr form the struts of the basket, while their unstructured termini constitute the basket distal densities, which potentially serve as a docking site for mRNA preprocessing before nucleocytoplasmic transport.

Characterization of haloacid dehalogenase superfamily acid phosphatase from Staphylococcus lugdunensis

The haloacid dehalogenase superfamily implicated in bacterial pathogenesis comprises different enzymes having roles in many metabolic pathways. Staphylococcus lugdunensis, a Gram-positive bacterium, is an opportunistic human pathogen causing infections in the central nervous system, urinary tract, bones, peritoneum, systemic conditions and cutaneous infection. The haloacid dehalogenase superfamily proteins play a significant role in the pathogenicity of certain bacteria, facilitating invasion, survival, and proliferation within host cells. The genome of S. lugdunensis encodes more than ten proteins belonging to this superfamily. However, none of them have been characterized. The present work reports the characterization of one of the haloacid dehalogenase superfamily proteins (SLHAD1) from Staphylococcus lugdunensis. The functional analysis revealed that SLHAD1 is a metal-dependent acid phosphatase, which catalyzes the dephosphorylation of phosphorylated metabolites of cellular pathways, including glycolysis, gluconeogenesis, nucleotides, and thiamine metabolism. Based on the substrate specificity and genomic analysis, the physiological function of SLHAD1 in thiamine metabolism has been tentatively assigned. The crystal structure of SLHAD1, lacking 49 residues at the C-terminal, was determined at 1.7 Å resolution with a homodimer in the asymmetric unit. It was observed that SLHAD1 exhibited time-dependent cleavage at a specific point, occurring through a self-initiated process. A combination of bioinformatics, biochemical, biophysical, and structural studies explored unique features of SLHAD1. Overall, the study revealed a detailed characterization of a critical enzyme of the human pathogen Staphylococcus lugdunensis, associated with several life-threatening infections.

Dissecting the Allosteric Fine-Tuning of Enzyme Catalysis

Fully understanding the mechanism of allosteric regulation in biomolecules requires separating and examining all of the involved factors. In enzyme catalysis, allosteric effector binding shifts the structure and dynamics of the active site, leading to modified energetic (e.g., energy barrier) and dynamical (e.g., diffusion coefficient) factors underlying the catalyzed reaction rate. Such modifications can be subtle and dependent on the type of allosteric effector, representing a fine-tuning of protein function. The microscopic description of allosteric regulation at the level of function-dictating factors has prospective applications in fundamental and pharmaceutical sciences, which is, however, largely missing so far. Here, we characterize the allosteric fine-tuning of enzyme catalysis, using human Pin1 as an example, by performing more than half-millisecond all-atom molecular dynamics simulations. Changes of reaction kinetics and the dictating factors, including the free energy surface along the reaction coordinate and the diffusion coefficient of the reaction dynamics, under various enzyme and allosteric effector binding conditions are examined. Our results suggest equal importance of the energetic and dynamical factors, both of which can be modulated allosterically, and the combined effect determines the final allosteric output. We also reveal the potential dynamic basis for allosteric modulation using an advanced statistical technique to detect function-related conformational dynamics. Methods developed in this work can be applied to other allosteric systems.

Role of the leucine-rich repeat protein kinase 2 C-terminal tail in domain cross-talk

Leucine-rich repeat protein kinase 2 (LRRK2) is a multi-domain protein encompassing two of biology's most critical molecular switches, a kinase and a GTPase, and mutations in LRRK2 are key players in the pathogenesis of Parkinson's disease (PD). The availability of multiple structures (full-length and truncated) has opened doors to explore intra-domain cross-talk in LRRK2. A helix extending from the WD40 domain and stably docking onto the kinase domain is common in all available structures. This C-terminal (Ct) helix is a hub of phosphorylation and organelle-localization motifs and thus serves as a multi-functional protein : protein interaction module. To examine its intra-domain interactions, we have recombinantly expressed a stable Ct motif (residues 2480-2527) and used peptide arrays to identify specific binding sites. We have identified a potential interaction site between the Ct helix and a loop in the CORB domain (CORB loop) using a combination of Gaussian accelerated molecular dynamics simulations and peptide arrays. This Ct-Motif contains two auto-phosphorylation sites (T2483 and T2524), and T2524 is a 14-3-3 binding site. The Ct helix, CORB loop, and the CORB-kinase linker together form a part of a dynamic 'CAP' that regulates the N-lobe of the kinase domain. We hypothesize that in inactive, full-length LRRK2, the Ct-helix will also mediate interactions with the N-terminal armadillo, ankyrin, and LRR domains (NTDs) and that binding of Rab substrates, PD mutations, or kinase inhibitors will unleash the NTDs.

Antibodies from dengue patients with prior exposure to Japanese encephalitis virus are broadly neutralizing against Zika virus

Exposure to multiple mosquito-borne flaviviruses within a lifetime is not uncommon; however, how sequential exposures to different flaviviruses shape the cross-reactive humoral response against an antigen from a different serocomplex has yet to be explored. Here, we report that dengue-infected individuals initially primed with the Japanese encephalitis virus (JEV) showed broad, highly neutralizing potencies against Zika virus (ZIKV). We also identified a rare class of ZIKV-cross-reactive human monoclonal antibodies with increased somatic hypermutation and broad neutralization against multiple flaviviruses. One huMAb, K8b, binds quaternary epitopes with heavy and light chains separately interacting with overlapping envelope protein dimer units spanning domains I, II, and III through cryo-electron microscopy and structure-based mutagenesis. JEV virus-like particle immunization in mice further confirmed that such cross-reactive antibodies, mainly IgG3 isotype, can be induced and proliferate through heterologous dengue virus (DENV) serotype 2 virus-like particle stimulation. Our findings highlight the role of prior immunity in JEV and DENV in shaping the breadth of humoral response and provide insights for future vaccination strategies in flavivirus-endemic countries.

Designed allosteric protein logic

The regulation of protein function by external or internal signals is one of the key features of living organisms. The ability to directly control the function of a selected protein would represent a valuable tool for regulating biological processes. Here, we present a generally applicable regulation of proteins called INSRTR, based on inserting a peptide into a loop of a target protein that retains its function. We demonstrate the versatility and robustness of coiled-coil-mediated regulation, which enables designs for either inactivation or activation of selected protein functions, and implementation of two-input logic functions with rapid response in mammalian cells. The selection of insertion positions in tested proteins was facilitated by using a predictive machine learning model. We showcase the robustness of the INSRTR strategy on proteins with diverse folds and biological functions, including enzymes, signaling mediators, DNA binders, transcriptional regulators, reporters, and antibody domains implemented as chimeric antigen receptors in T cells. Our findings highlight the potential of INSRTR as a powerful tool for precise control of protein function, advancing our understanding of biological processes and developing biotechnological and therapeutic interventions.

Effect of single nucleotide polymorphisms on the structure of long noncoding RNAs and their interaction with RNA binding proteins

Long non-coding RNAs (lncRNA) are emerging as a new class of regulatory RNAs with remarkable potential to be utilized as therapeutic targets against many human diseases. Several genome-wide association studies (GWAS) have catalogued Single Nucleotide Polymorphisms (SNPs) present in the noncoding regions of the genome from where lncRNAs originate. In this study, we have selected 67 lncRNAs with GWAS-tagged SNPs and have also investigated their role in affecting the local secondary structures. Majority of the SNPs lead to changes in the secondary structure of lncRNAs to a different extent by altering the base pairing patterns. These structural changes in lncRNA are also manifested in form of alteration in the binding site for RNA binding proteins (RBPs) along with affecting their binding efficacies. Ultimately, these structural modifications may influence the transcriptional and post-transcriptional pathways of these RNAs, leading to the causation of diseases. Hence, it is important to understand the possible underlying mechanism of RBPs in association with GWAS-tagged SNPs in human diseases.

Structural and functional validation of a cloned parasporin from Bacillus thuringiensis isolate KAU 41 native to Western Ghats of India

Parasporins of Bacillus thuringiensis (Bt) exhibit specific toxicity to cancer cells. PCR based mining has identified apoptosis inducing parasporin in KAU41 Bt isolate from the Western Ghats of India. The study aimed to clone and overexpress the parasporin of native KAU41 Bt isolate for determining structural and functional characteristics of the protein. Parasporin gene was cloned in pGEM-T, sequenced, sub-cloned in pET30+ and overexpressed in Escherichia coli. The expressed protein was characterized by SDS-PAGE and in silico methods. Cytotoxicity of cleaved peptide was assessed by MTT assay. SDS-PAGE displayed a 31 kDa protein (rp-KAU41) overexpressed. Upon proteinase K digestion, the protein was cleaved into 29 kDa peptide which was found to be cytotoxic to HeLa cells. The protein has a deduced sequence of 267 amino acids with β-strands folding pattern of crystal protein. Even though rp-KAU41 shared a 99.15% identity to chain-A of non-toxic crystal protein, it only showed a less similarity to the existing parasporins like PS4 (38%) and PS5 (24%) in UPGMA analysis, emphasizing the novelty of rp-KAU41. The protein is predicted to have more similarity to the pore forming toxins of Aerolysin superfamily and an additional loop in rp-KAU41 may be contributing towards its cytotoxicity. The molecular docking with caspase 3 resulted in higher Z dock and Z rank score substantiating its role in the activation of intrinsic pathway of apoptosis. The recombinant parasporin protein, rp-KAU41 is presumed to belong to the Aerolysin superfamily. An interaction with caspase 3 substantiates its role in activating the intrinsic pathway of apoptosis in cancer cells.

Direct interaction between the transmembrane helices stabilize cytochrome P450 2B4 and cytochrome b5 redox complex

The catalytic activity of cytochrome P450 2B4 (CYP2B4) is moderated by its cognate redox partner cytochrome b5 (Cyt-b5). The endoplasmic reticulum (ER) membrane and intermolecular transmembrane (TM) interaction between CYP2B4 and Cyt-b5 regulate the substrate catalysis and the reaction rate. This emphasizes the significance of elucidating the molecular basis of CYP2B4 and Cyt-b5 complexation in a membrane environment to better understand the enzymatic activity of CYP2B4. Our previous solid-state NMR studies revealed the membrane topology of the transmembrane domains of these proteins in the free and complex forms. Here, we show the cross-angle complex formation by the single-pass TM domains of CYP2B4 and Cyt-b5, which is mainly driven by several salt-bridges (E2-R128, R21-D104 and K25-D104), using a multi-microsecond molecular dynamic simulation. Additionally, the leucine-zipper residues (L8, L12, L15, L18 and L19 from CYP2B4) and π-stacking between H23 and F20 residues of CYP2B4 and W110 of Cyt-b5 are identified to stabilize the TM-TM complex in the ER membrane. The simulated tilts of the helices in the free and in the complex are in excellent agreement with solid-state NMR results. The TM-TM packing influences a higher order structural stability when compared to the complex formed by the truncated soluble domains of these two proteins. MM/PBSA based binding free energy estimates nearly 100-fold higher binding affinity (ΔG = -2810.68 ± 696.44 kJ/mol) between the soluble domains of the full-length CYP2B4 and Cyt-b5 when embedded in lipid membrane as compared to the TM-domain-truncated soluble domains (ΔG = -27.406 ± 10.32 kJ/mol). The high-resolution full-length CYP2B4-Cyt-b5 complex structure and its dynamics in a native ER membrane environment reported here could aid in the development of approaches to effectively modulate the drug-metabolism activity of CYP2B4.

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.

Molecular Dynamics-Based Design and Biophysical Evaluation of Thermostable Single-Chain Fv Antibody Mutants Derived from Pharmaceutical Antibodies

Antibody drugs are denatured under physical stress, e.g., friction, heat, and freezing, which triggers formation of aggregates and resultant allergic reactions. Design of a stable antibody is thus critical for the development of antibody drugs. Here, we obtained a thermostable single-chain Fv (scFv) antibody clone by rigidifying the flexible region. We first conducted a short molecular dynamics (MD) simulation (3 runs of 50 ns) to search for weak spots in the scFv antibody, i.e., flexible regions located outside the CDR (complementarity determining region) and the interface between the heavy-chain and light-chain variable regions. We then designed a thermostable mutant and evaluated it by means of a short MD simulation (3 runs of 50 ns) based on reductions in the root-mean-square fluctuation (RMSF) values and formation of new hydrophilic interactions around the weak spot. Finally, we designed the VL-R66G mutant by applying our strategy to scFv derived from trastuzumab. Trastuzumab scFv variants were prepared by using an Escherichia coli expression system, and the melting temperature-measured as a thermostability index-was 5 °C higher than that of the wild-type trastuzumab scFv, while the antigen-binding affinity was unchanged. Our strategy required few computational resources, and would be applicable to antibody drug discovery.

Molecular Mechanism of Nucleosome Recognition by the Pioneer Transcription Factor Sox

Pioneer transcription factors (PTFs) have the remarkable ability to directly bind to chromatin to stimulate vital cellular processes. In this work, we dissect the universal binding mode of Sox PTF by combining extensive molecular simulations and physiochemistry approaches, along with DNA footprinting techniques. As a result, we show that when Sox consensus DNA is located at the solvent-facing DNA strand, Sox binds to the compact nucleosome without imposing any significant conformational changes. We also reveal that the base-specific Sox:DNA interactions (base reading) and Sox-induced DNA changes (shape reading) are concurrently required for sequence-specific nucleosomal DNA recognition. Among three different nucleosome positions located on the positive DNA arm, a sequence-specific reading mechanism is solely satisfied at the superhelical location 2 (SHL2). While SHL2 acts transparently for solvent-facing Sox binding, among the other two positions, SHL4 permits only shape reading. The final position, SHL0 (dyad), on the other hand, allows no reading mechanism. These findings demonstrate that Sox-based nucleosome recognition is essentially guided by intrinsic nucleosome properties, permitting varying degrees of DNA recognition.

An optimized thermodynamics integration protocol for identifying beneficial mutations in antibody design

Accurate identification of beneficial mutations is central to antibody design. Many knowledge-based (KB) computational approaches have been developed to predict beneficial mutations, but their accuracy leaves room for improvement. Thermodynamic integration (TI) is an alchemical free energy algorithm that offers an alternative technique for identifying beneficial mutations, but its performance has not been evaluated. In this study, we developed an efficient TI protocol with high accuracy for predicting binding free energy changes of antibody mutations. The improved TI method outperforms KB methods at identifying both beneficial and deleterious mutations. We observed that KB methods have higher accuracies in predicting deleterious mutations than beneficial mutations. A pipeline using KB methods to efficiently exclude deleterious mutations and TI to accurately identify beneficial mutations was developed for high-throughput mutation scanning. The pipeline was applied to optimize the binding affinity of a broadly sarbecovirus neutralizing antibody 10-40 against the circulating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) omicron variant. Three identified beneficial mutations show strong synergy and improve both binding affinity and neutralization potency of antibody 10-40. Molecular dynamics simulation revealed that the three mutations improve the binding affinity of antibody 10-40 through the stabilization of an altered binding mode with increased polar and hydrophobic interactions. Above all, this study presents an accurate and efficient TI-based approach for optimizing antibodies and other biomolecules.

Unconventional interactions of the TRPV4 ion channel with beta-adrenergic receptor ligands

The transient receptor potential vanilloid 4 (TRPV4) ion channel is present in different tissues including those of the airways. This channel is activated in response to stimuli such as changes in temperature, hypoosmotic conditions, mechanical stress, and chemicals from plants, lipids, and others. TRPV4's overactivity and/or dysfunction has been associated with several diseases, such as skeletal dysplasias, neuromuscular disorders, and lung pathologies such as asthma and cardiogenic lung edema and COVID-19-related respiratory malfunction. TRPV4 antagonists and blockers have been described; nonetheless, the mechanisms involved in achieving inhibition of the channel remain scarce, and the search for safe use of these molecules in humans continues. Here, we show that the widely used bronchodilator salbutamol and other ligands of β-adrenergic receptors inhibit TRPV4's activation. We also demonstrate that inhibition of TRPV4 by salbutamol is achieved through interaction with two residues located in the outer region of the pore and that salbutamol leads to channel closing, consistent with an allosteric mechanism. Our study provides molecular insights into the mechanisms that regulate the activity of this physiopathologically important ion channel.

Biofortification of Crops to Fight Anemia: Role of Vacuolar Iron Transporters

Plant-based foods provide all the crucial nutrients for human health. Among these, iron (Fe) is one of the essential micronutrients for plants and humans. A lack of Fe is a major limiting factor affecting crop quality, production, and human health. There are people who suffer from various health problems due to the low intake of Fe in their plant-based foods. Anemia has become a serious public health issue due to Fe deficiency. Enhancing Fe content in the edible part of food crops is a major thrust area for scientists worldwide. Recent progress in nutrient transporters has provided an opportunity to resolve Fe deficiency or nutritional problems in plants and humans. Understanding the structure, function, and regulation of Fe transporters is essential to address Fe deficiency in plants and to improve Fe content in staple food crops. In this review, we summarized the role of Fe transporter family members in the uptake, cellular and intercellular movement, and long-distance transport of Fe in plants. We draw insights into the role of vacuolar membrane transporters in the crop for Fe biofortification. We also provide structural and functional insights into cereal crops' vacuolar iron transporters (VITs). This review will help highlight the importance of VITs for improving the Fe biofortification of crops and alleviating Fe deficiency in humans.

The Intrinsic Radius as a Key Parameter in the Generalized Born Model to Adjust Protein-Protein Electrostatic Interaction

The generalized Born (GB) model is an extension of the continuum dielectric theory of Born solvation energy and is a powerful method for accelerating the molecular dynamic (MD) simulations of charged biological molecules in water. While the effective dielectric constant of water that varies as a function of the separation distance between solute molecules is incorporated into the GB model, adjustment of the parameters is indispensable for accurate calculation of the Coulomb (electrostatic) energy. One of the key parameters is the lower limit of the spatial integral of the energy density of the electric field around a charged atom, known as the intrinsic radius ρ. Although ad hoc adjustment of ρ has been conducted to improve the Coulombic (ionic) bond stability, the physical mechanism by which ρ affects the Coulomb energy remains unclear. Via energetic analysis of three differently sized systems, here, we clarify that the Coulomb bond stability increases with increasing ρ and that the increased stability is caused by the interaction energy term, not by the self-energy (desolvation energy) term, as was supposed previously. Our results suggest that the use of larger values for the intrinsic radii of hydrogen and oxygen atoms, together with the use of a relatively small value for the spatial integration cutoff in the GB model, can better reproduce the Coulombic attraction between protein molecules.

Molecular Dynamics Refinement of Open State Serotonin 5-HT3A Receptor Structures

Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT₃ receptors are a clinical target for treating nausea and vomiting associated with chemotherapy and radiation therapy (Thompson and Lummis, 2006). Multiple cryo-electron microscopy (cryo-EM) structures of 5-HT₃ receptors have been determined in distinct functional states (e.g., open, closed, etc.) (Basak et al., 2018; Basak et al., 2018; Polovinkin et al., 2018; Zhang et al., 2015). However, recent work has shown that the transmembrane pores of the open 5-HT₃ receptor structures rapidly collapse and become artificially asymmetric in molecular dynamics (MD) simulations. To avoid this hydrophobic collapse, Dämgen and Biggin developed an equilibration protocol that led to a stable open state structure of the glycine receptor in MD simulations (Dämgen and Biggin, 2020). However, the protocol failed to yield open-like structures of the 5-HT₃ receptor in our simulations. Here, we present a refined equilibration protocol that involves the rearrangement of the transmembrane helices to achieve stable open state structures of the 5-HT₃ receptor that allow both water and ion permeation through the channel. Notably, channel gating is mediated through collective movement of the transmembrane helices, involving not only pore lining M2 helices but also their cross-talk with the adjacent M1 and M3 helices. Thus, the successful application of our refined equilibration protocol underscores the importance of the conformational coupling between the transmembrane helices in stabilizing open-like structures of the 5-HT₃ receptor.

Conformational stability of the deamidated and mutated human βB2-crystallin

Previous studies propose that genetic mutations and post-translational modifications in protein crystallins promote protein aggregation and are considered significant risk factors for cataract formation. The βB2-crystallin (HβB2C) forms a high proportion of proteins in the human eye lens. Different congenital mutations and post-translational deamidations in βB2-crystallin have been reported and linked to cataract formation. In this work, we employed extensive all-atom molecular dynamics simulations to evaluate the conformational stability of deamidated and mutated HβB2C. Our results show critical changes in the protein surface and its native contacts due to a modification in the conformational equilibrium of these proteins. The double deamidated (Q70E/Q162E) and single deamidated (Q70E) impact the well compact conformation of the HβB2C. These post-translational modifications allow the exposure of the protein hydrophobic interface, which lead to the exposure of electronegative residues. On the other hand, our mutational studies showed that the S143F mutation modifies the hydrogen-bond network of an antiparallel β-sheet, unfolding the C-terminal domain. Interestingly, the chain termination mutation (Q155X) does not unfold the N-terminal domain. However, the resultant conformation is more compact and avoids the exposure of the hydrophobic interface. Our results provide valuable information about the first steps of HβB2C unfolding in the presence of deamidated amino acids that have been reported to appear during aging. The findings reported in this work are essential for the general knowledge of the initial steps in the cataract formation mechanism, which may be helpful for the further development of molecules with pharmacological potential against cataract disease.

I329L: A Dual Action Viral Antagonist of TLR Activation Encoded by the African Swine Fever Virus (ASFV)

The African Swine Fever Virus (ASFV) is an economically important, large DNA virus which causes a highly contagious and frequently fatal disease in domestic pigs. Due to the acute nature of the infection and the complexity of the protective porcine anti-ASFV response, there is no accepted vaccine in use. As resistance to ASFV is known to correlate with a robust IFN response, the virus is predicted to have evolved strategies to inhibit innate immunity by modulating the IFN response. The deletion of virus host evasion gene(s) inhibiting IFN is a logical solution to develop an attenuated virus vaccine. One such candidate, the ASFV ORF I329L gene, is highly conserved in pathogenic and non-pathogenic virus isolates and in this study we confirm and extend the conclusion that it has evolved for the inhibition of innate immunity initiated through Toll-like receptors (TLRs). Specifically, the ASFV I329L extracellular (ECD) and intracellular (ICD) domains inhibit TLR signalling by two entirely different mechanisms. Bioinformatics modelling suggests that the ECD inhibits several TLR signalling pathways through a short sequence homologous to the conserved TLR dimerization domain, here termed the putative dimerization domain (PDD). Remarkably, both full length and PDD constructs of I329L were demonstrated to inhibit activation, not only of TLR3, but also TLR4, TLR5, TLR8 and TLR9. Additionally, the demonstration of a weak association of I329L with TLR3 is consistent with the formation of a non-signalling I329L-TLR3 heterodimer, perhaps mediated through the PDD of I329L. Finally, the ICD associates with TRIF, thereby impacting on both TLR3 and TLR4 signalling. Thus, I329L offers potential as a general inhibitor of TLR responses and is a rational candidate for construction and testing of an I329L deletion mutant vaccine.

The Biofilm Regulatory Network from Bacillus subtilis: A Structure-Function Analysis

Bacterial biofilms are notorious for their ability to protect bacteria from environmental challenges, most importantly the action of antibiotics. Bacillus subtilis is an extensively studied model organism used to understand the process of biofilm formation. A complex network of principal regulatory proteins including Spo0A, AbrB, AbbA, Abh, SinR, SinI, SlrR, and RemA, work in concert to transition B. subtilis from the free-swimming planktonic state to the biofilm state. In this review, we explore, connect, and summarize decades worth of structural and biochemical studies that have elucidated this protein signaling network. Since structure dictates function, unraveling aspects of protein molecular mechanisms will allow us to devise ways to exploit critical features of the biofilm regulatory pathway, such as possible therapeutic intervention. This review pools our current knowledge base of B. subtilis biofilm regulatory proteins and highlights potential therapeutic intervention points.

Identification of Quorum Quenching N-Acyl Homoserine Lactonases from Priestia aryabhattai J1D and Bacillus cereus G Isolated from the Rhizosphere

Several pathogenic bacteria communicate using N-acyl homoserine lactone (AHL) as a quorum sensing (QS) molecule. The process of interfering with the QS system is known as quorum quenching (QQ), it is an effective tool to control QS-dependent virulence in pathogens. In the present study, rhizosphere bacterial isolates were screened for their ability to produce AHL lactonase enzyme as QQ molecules, which hydrolyses AHL signalling molecules and consequently blocks the QS system. Potent N-hexanoyl-l-homoserine lactone (C6HSL) hydrolytic QQ activity was detected in rhizosphere isolates namely Bacillus cereus G and Priestia aryabhattai J1D. The cell-free supernatant of the bacterial isolates indicated a reduction in biofilm formation in the human pathogens Vibrio cholerae, Pseudomonas aeruginosa, and Staphylococcus aureus without inhibiting cells, signifying their biocontrol property. Furthermore, liquid chromatography high resolution mass spectrometry analysis confirmed C6HSL hydrolytic activity by AHL lactonase produced by these rhizosphere isolates. Also, the aiiA homologous gene from the bacterial isolates showed similarity with the aiiA lactonase gene from Bacillus species, which was further confirmed by homology modelling. In silico structure analysis by comparing with the structure of Bacillus revealed the similarity in the active site, indicating the same degradation pattern. Based on available reported data, the present study indicates the first report of the presence of the aiiA lactonase gene in P. aryabhattai.

Conformational transitions induced by NADH binding promote reduction half-reaction in 2-hydroxybiphenyl-3-monooxygenase catalytic cycle

2-Hydroxybiphenyl-3-monoxygenase from Pseudomonas azelaica is an effective catalyst of the regiospecific conversions of various aromatic compounds. A comprehensive understanding of the complete catalytic cycle, including the as yet unclear details of NADH binding, NADH/FAD interaction as well as related conformational changes could facilitate the rational design of improved enzyme variants for practical applications. Induced fit formation of a specific pocket for the nicotinamide ring at NADH binding has been revealed using advanced molecular simulation methods including metadynamics and QM/MM modeling. The resulting triple stacking interaction of the nicotinamide as well as isoalloxazine rings and evolutionarily correlated amino acid residues of the active site greatly contributes to the stabilization of the charge-transfer complex and determines the Pro-S stereospecificity of the hydride transfer and the low energy barrier 11 kcal/mol. Then the resulting FADH^- anion undergoes the consequent conformational transition of the FAD isoalloxazine ring from the open out to the closed in position which is followed by the binding of an oxygen molecule what is crucial for the next step of substrate oxidation and the completion of the catalytic cycle.

Modulation of the IKS channel by PIP2 requires two binding sites per monomer

The phosphatidyl-inositol-4,5-bisphosphate (PIP₂) lipid has been shown to be crucial for the coupling between the voltage sensor and the pore of the potassium voltage-gated K_V7 channel family, especially the K_V7.1 channel. Expressed in the myocardium membrane, K_V7.1 forms a complex with KCNE1 auxiliary subunits to generate the I_KS current. Here we present molecular models of the transmembrane region of this complex in its three known states, namely the Resting/Closed (RC), the Intermediate/Closed (IC), and the Activated/Open (AO), robustness of which is assessed by agreement with a range of biophysical data. Molecular Dynamics (MD) simulations of these models embedded in a lipid bilayer including phosphatidyl-inositol-4,5-bisphosphate (PIP₂) lipids show that in presence of KCNE1, two PIP₂ lipids are necessary to stabilize each state. The simulations also show that KCNE1 interacts with both PIP₂ binding sites, forming a tourniquet around the pore and preventing its opening. The present investigation provides therefore key molecular elements that govern the role of PIP₂ in KCNE1 modulation of I_KS channels, possibly a common mechanism by which auxiliary KCNE subunits might modulate a variety of other ion channels.

Coiled-Coil Protein Origami: Design, Isolation, and Characterization

Coiled-coil protein origami (CCPO) is a rationally designed de novo protein fold, constructed by concatenating coiled-coil forming segments into a polypeptide chain, that folds into polyhedral nano-cages. To date, nanocages in the shape of a tetrahedron, square pyramid, trigonal prism, and trigonal bipyramid have been successfully designed and extensively characterized following the design principles of CCPO. These designed protein scaffolds and their favorable biophysical properties are suitable for functionalization and other various biotechnological applications. To further facilitate the development, we are presenting a detailed guide to the world of CCPO, starting from design (CoCoPOD, an integrated platform for designing CCPO strictures) and cloning (modified Golden-gate assembly) to fermentation and isolation (NiNTA, Strep-trap, IEX, and SEC) concluding with standard characterization techniques (CD, SEC-MALS, and SAXS).

Molecular interactions underlying the phase separation of HP1α: role of phosphorylation, ligand and nucleic acid binding

Heterochromatin protein 1α (HP1α) is a crucial element of chromatin organization. It has been proposed that HP1α functions through liquid-liquid phase separation (LLPS), which allows it to compact chromatin into transcriptionally repressed heterochromatin regions. In vitro, HP1α can undergo phase separation upon phosphorylation of its N-terminus extension (NTE) and/or through interactions with DNA and chromatin. Here, we combine computational and experimental approaches to elucidate the molecular interactions that drive these processes. In phosphorylation-driven LLPS, HP1α can exchange intradimer hinge-NTE interactions with interdimer contacts, which also leads to a structural change from a compacted to an extended HP1α dimer conformation. This process can be enhanced by the presence of positively charged HP1α peptide ligands and disrupted by the addition of negatively charged or neutral peptides. In DNA-driven LLPS, both positively and negatively charged peptide ligands can perturb phase separation. Our findings demonstrate the importance of electrostatic interactions in HP1α LLPS where binding partners can modulate the overall charge of the droplets and screen or enhance hinge region interactions through specific and non-specific effects. Our study illuminates the complex molecular framework that can fine-tune the properties of HP1α and that can contribute to heterochromatin regulation and function.

Molecular dynamics simulations reveal membrane lipid interactions of the full-length lymphocyte specific kinase (Lck)

The membrane-bound lymphocyte-specific protein-tyrosine kinase (Lck) triggers T cell antigen receptor signalling to initiate adaptive immune responses. Despite many structure-function studies, the mode of action of Lck and the potential role of plasma membrane lipids in regulating Lck's activity remains elusive. Advances in molecular dynamics simulations of membrane proteins in complex lipid bilayers have opened a new perspective in gathering such information. Here, we have modelled the full-length Lck open and closed conformations  using data available from different crystalographic studies and simulated its interaction with the inner leaflet of the T cell plasma membrane. In both conformations, we found that the unstructured unique domain and the structured domains including the kinase interacted with the membrane with a preference for PIP lipids. Interestingly, our simulations suggest that the Lck-SH2 domain interacts with lipids differently in the open and closed Lck conformations, demonstrating that lipid interaction can potentially regulate Lck's conformation and in turn modulate T cell signalling. Additionally, the Lck-SH2 and kinase domain residues that significantly contacted PIP lipids are found to be conserved among the Src family of kinases, thereby potentially representing similar PIP interactions within the family.

Structure of MotA, a flagellar stator protein, from hyperthermophile

Many motile bacteria swim and swarm toward favorable environments using the flagellum, which is rotated by a motor embedded in the inner membrane. The motor is composed of the rotor and the stator, and the motor torque is generated by the change of the interaction between the rotor and the stator induced by the ion flow through the stator. A stator unit consists of two types of membrane proteins termed A and B. Recent cryo-EM studies on the stators from mesophiles revealed that the stator consists of five A and two B subunits, whereas the low-resolution EM analysis showed that purified hyperthermophilic MotA forms a tetramer. To clarify the assembly formation and factors enhancing thermostability of the hyperthermophilic stator, we determined the cryo-EM structure of MotA from Aquifex aeolicus (Aa-MotA), a hyperthermophilic bacterium, at 3.42 Å resolution. Aa-MotA forms a pentamer with pseudo C5 symmetry. A simulated model of the Aa-MotA₅MotB₂ stator complex resembles the structures of mesophilic stator complexes, suggesting that Aa-MotA can assemble into a pentamer equivalent to the stator complex without MotB. The distribution of hydrophobic residues of MotA pentamers suggests that the extremely hydrophobic nature in the subunit boundary and the transmembrane region is a key factor to stabilize hyperthermophilic Aa-MotA.

Computer-Aided Drug Design Boosts RAS Inhibitor Discovery

The Rat Sarcoma (RAS) family (NRAS, HRAS, and KRAS) is endowed with GTPase activity to regulate various signaling pathways in ubiquitous animal cells. As proto-oncogenes, RAS mutations can maintain activation, leading to the growth and proliferation of abnormal cells and the development of a variety of human cancers. For the fight against tumors, the discovery of RAS-targeted drugs is of high significance. On the one hand, the structural properties of the RAS protein make it difficult to find inhibitors specifically targeted to it. On the other hand, targeting other molecules in the RAS signaling pathway often leads to severe tissue toxicities due to the lack of disease specificity. However, computer-aided drug design (CADD) can help solve the above problems. As an interdisciplinary approach that combines computational biology with medicinal chemistry, CADD has brought a variety of advances and numerous benefits to drug design, such as the rapid identification of new targets and discovery of new drugs. Based on an overview of RAS features and the history of inhibitor discovery, this review provides insight into the application of mainstream CADD methods to RAS drug design.

Molecular dynamics of DNA translocation by FtsK

The bacterial FtsK motor harvests energy from ATP to translocate double-stranded DNA during cell division. Here, we probe the molecular mechanisms underlying coordinated DNA translocation in FtsK by performing long timescale simulations of its hexameric assembly and individual subunits. From these simulations we predict signaling pathways that connect the ATPase active site to DNA-gripping residues, which allows the motor to coordinate its translocation activity with its ATPase activity. Additionally, we utilize well-tempered metadynamics simulations to compute free-energy landscapes that elucidate the extended-to-compact transition involved in force generation. We show that nucleotide binding promotes a compact conformation of a motor subunit, whereas the apo subunit is flexible. Together, our results support a mechanism whereby each ATP-bound subunit of the motor conforms to the helical pitch of DNA, and ATP hydrolysis/product release causes a subunit to lose grip of DNA. By ordinally engaging and disengaging with DNA, the FtsK motor unidirectionally translocates DNA.

Genotypic and phenotypic study of antiviral resistance mutations in refractory cytomegalovirus infection

This study describes the genotypic and phenotypic characterisation of novel human cytomegalovirus (HCMV) genetic variants of a cohort of 94 clinically-resistant HCMV patients. Antiviral-resistant mutations were detected in the UL97, UL54 and UL56 target genes of 25/94 (26.6%) patients. The genotype-phenotype correlation study resolved the status of 5 uncharacterised UL54 DNA polymerase (G441S, A543V, F460S, R512C, A928T) and 2 UL56 terminase (F345L, P800L) mutations found in clinical isolates. A928T conferred high triple-resistance to ganciclovir, foscarnet and cidofovir, and A543V had 10-fold reduced susceptibility to cidofovir. Viral growth assays showed G441S, A543V, F345L and P800L impaired viral growth capacities compared with wild-type AD169 HCMV. 3D modelling predicted A543V and A928T phenotypes but not R512C, reinforcing the need for individual characterisation of mutations by recombinant phenotyping. Extending mutation databases is crucial to optimize treatments and to improve the assessment of patients with resistant/refractory HCMV infection.

Cancer-related Mutations with Local or Long-range Effects on an Allosteric Loop of p53

The tumor protein 53 (p53) is involved in transcription-dependent and independent processes. Several p53 variants related to cancer have been found to impact protein stability. Other variants, on the contrary, might have little impact on structural stability and have local or long-range effects on the p53 interactome. Our group previously identified a loop in the DNA binding domain (DBD) of p53 (residues 207-213) which can recruit different interactors. Experimental structures of p53 in complex with other proteins strengthen the importance of this interface for protein-protein interactions. We here characterized with structure-based approaches somatic and germline variants of p53 which could have a marginal effect in terms of stability and act locally or allosterically on the region 207-213 with consequences on the cytosolic functions of this protein. To this goal, we studied 1132 variants in the p53 DBD with structure-based approaches, accounting also for protein dynamics. We focused on variants predicted with marginal effects on structural stability. We then investigated each of these variants for their impact on DNA binding, dimerization of the p53 DBD, and intramolecular contacts with the 207-213 region. Furthermore, we identified variants that could modulate long-range the conformation of the region 207-213 using a coarse-grain model for allostery and all-atom molecular dynamics simulations. Our predictions have been further validated using enhanced sampling methods for 15 variants. The methodologies used in this study could be more broadly applied to other p53 variants or cases where conformational changes of loop regions are essential in the function of disease-related proteins.

In Silico Analysis and Functional Characterization of Antimicrobial and Insecticidal Vicilin from Moth Bean (Vigna aconitifolia (Jacq.) Marechal) Seeds

Vicilin has nutraceutical potential and different noteworthy medicative health-promoting biotic diversions, and it is remarkable against pathogenic microorganisms and insects. In this study, Vigna aconitifolia vicilin (VacV) has been identified and characterized from the seed of Vigna aconitifolia (Jacq.) Marechal (Moth beans). LC-MS/MS analysis of VacV provided seven random fragmented sequences comprising 238 residues, showing significant homology with already reported Vigna radiata vicilin (VraV). VacV was purified using ammonium sulfate precipitation (60%) followed by size exclusion chromatography on Hi-Load 16/60 Superdex 200 pg column and anion-exchange chromatography (Hi trap Q FF column). Purified VacV showed a major ~50 kDa band and multiple lower bands on 12% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) under both reduced and non-reduced conditions. After all, a three-dimensional molecular structure of VacV was predicted, which showed β-sheeted molecular conformation similar to crystallographic structure of VraV. All Vicilins from V. aconitifolia and other plants were divided into six sub-groups by phylogenetic analysis, and VacV shared a high degree of similarity with vicilins of Vigna radiata, Pisum sativum, Lupinus albus, Cicer arietinum and Glycine max. Additionally, VacV (20 μg) has significant growth inhibition against different pathogenic bacteria along strong antifungal activity (50 μg). Likewise, VacV (3.0 mg) produced significant growth reduction in Rice Weevil Sitophilus oryzae larvae after 9 days compared with control. Furthermore, by using MMT assay, the cytotoxicity effect of VacV on the growth of HepG2 liver cancerous cells was tested. VacV showed cytotoxicity against the HepG-2 line and the acquired value was 180 µg after 48 h. Finally, we performed molecular docking against caspase-3 protein (PDB ID: 3DEI) for VacV bioactive receptor interface residues. Hence, our results reveal that VacV, has nutraceutical potential and moth beans can be used as a rich resource of functional foods.

Modulation of Tumor Immune Microenvironment and Prognostic Value of Ferroptosis-Related Genes, and Candidate Target Drugs in Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is the most common type of malignant brain tumor, among which IDH1-wild type GBM has a poor prognosis. Recent studies have shown that ferroptosis-related genes (FRGs) are correlated with the development and progression of cancer. In GBM, the role of FRGs associated with IDH1 status as biological indicators and therapeutic targets remains to be clarified. Ten of FRGs (STEAP3, HSPB1, MAP1LC3A, SOCS1, LOX, CAPG, CP, GDF15, CDKN1A, and CD44) associated with IDH1 status in GBM were identified as key genes through screening by survival analysis and Random Forest using The Cancer Genome Atlas (TCGA) datasets, and the protein expressions of key genes were verified. Transwell and qPCR results showed that ferroptosis promoted the migration of glioblastoma cells and affected the expression of key genes. Our study established the ferroptosis-related prognostic model for GBM patients based on ten key genes by a different modeling method from previous study, the GSVA algorithm. Further, we took the methods of functional enrichment analysis, clinical characteristics, immune cell infiltration, immunomodulator, ESTIMATE and single nucleotide variant (SNV) analysis to study the molecular mechanisms of prognostic model and key genes. The results showed that ten key genes were strongly associated with immune-related factors and were significantly involved in the p53 signaling pathway, senescence and autophagy in cancer, and in the negative regulation of protein kinase activity. Moreover, potential therapeutic drugs were identified by Virtual Screening and Molecular Docking. Our study indicated that the novel ferrotosis-related prognostic model for GBM patients and key genes possessed the prognostic and therapeutic values.

Carbohydrate Metabolism in Bacteria: Alternative Specificities in ADP-Glucose Pyrophosphorylases Open Novel Metabolic Scenarios and Biotechnological Tools

We explored the ability of ADP-glucose pyrophosphorylase (ADP-Glc PPase) from different bacteria to use glucosamine (GlcN) metabolites as a substrate or allosteric effectors. The enzyme from the actinobacteria Kocuria rhizophila exhibited marked and distinctive sensitivity to allosteric activation by GlcN-6P when producing ADP-Glc from glucose-1-phosphate (Glc-1P) and ATP. This behavior is also seen in the enzyme from Rhodococcus spp., the only one known so far to portray this activation. GlcN-6P had a more modest effect on the enzyme from other Actinobacteria (Streptomyces coelicolor), Firmicutes (Ruminococcus albus), and Proteobacteria (Agrobacterium tumefaciens) groups. In addition, we studied the catalytic capacity of ADP-Glc PPases from the different sources using GlcN-1P as a substrate when assayed in the presence of their respective allosteric activators. In all cases, the catalytic efficiency of Glc-1P was 1-2 orders of magnitude higher than GlcN-1P, except for the unregulated heterotetrameric protein (GlgC/GgD) from Geobacillus stearothermophilus. The Glc-1P substrate preference is explained using a model of ADP-Glc PPase from A. tumefaciens based on the crystallographic structure of the enzyme from potato tuber. The substrate-binding domain localizes near the N-terminal of an α-helix, which has a partial positive charge, thus favoring the interaction with a hydroxyl rather than a charged primary amine group. Results support the scenario where the ability of ADP-Glc PPases to use GlcN-1P as an alternative occurred during evolution despite the enzyme being selected to use Glc-1P and ATP for α-glucans synthesis. As an associated consequence in such a process, certain bacteria could have improved their ability to metabolize GlcN. The work also provides insights in designing molecular tools for producing oligo and polysaccharides with amino moieties.

Structure of the anthrax protective antigen D425A dominant negative mutant reveals a stalled intermediate state of pore maturation

The tripartite protein complex produced by anthrax bacteria (Bacillus anthracis) is a member of the AB family of β-barrel pore-forming toxins. The protective antigen (PA) component forms an oligomeric prepore that assembles on the host cell surface and serves as a scaffold for binding of lethal and edema factors. Following endocytosis, the acidic environment of the late endosome triggers a pH-induced conformational rearrangement to promote maturation of the PA prepore to a functional, membrane spanning pore that facilitates delivery of lethal and edema factors to the cytosol of the infected host. Here, we show that the dominant-negative D425A mutant of PA stalls anthrax pore maturation in an intermediate state at acidic pH. Our 2.7 Å cryo-EM structure of the intermediate state reveals structural rearrangements that involve constriction of the oligomeric pore combined with an intramolecular dissociation of the pore-forming module. In addition to defining the early stages of anthrax pore maturation, the structure identifies asymmetric conformational changes in the oligomeric pore that are influenced by the precise configuration of adjacent protomers.

Structural and functional insights into ABHD5, a ligand-regulated lipase co-activator

Alpha/beta hydrolase domain-containing protein 5 (ABHD5) is a highly conserved protein that regulates various lipid metabolic pathways via interactions with members of the perilipin (PLIN) and Patatin-like phospholipase domain-containing protein (PNPLA) protein families. Loss of function mutations in ABHD5 result in Chanarin-Dorfman Syndrome (CDS), characterized by ectopic lipid accumulation in numerous cell types and severe ichthyosis. Recent data demonstrates that ABHD5 is the target of synthetic and endogenous ligands that might be therapeutic beneficial for treating metabolic diseases and cancers. However, the structural basis of ABHD5 functional activities, such as protein-protein interactions and ligand binding is presently unknown. To address this gap, we constructed theoretical structural models of ABHD5 by comparative modeling and topological shape analysis to assess the spatial patterns of ABHD5 conformations computed in protein dynamics. We identified functionally important residues on ABHD5 surface for lipolysis activation by PNPLA2, lipid droplet targeting and PLIN-binding. We validated the computational model by examining the effects of mutating key residues in ABHD5 on an array of functional assays. Our integrated computational and experimental findings provide new insights into the structural basis of the diverse functions of ABHD5 as well as pathological mutations that result in CDS.

Structural Basis of Antibody Conformation and Stability Modulation by Framework Somatic Hypermutation

Accumulation of somatic hypermutation (SHM) is the primary mechanism to enhance the binding affinity of antibodies to antigens in vivo. However, the structural basis of the effects of many SHMs remains elusive. Here, we integrated atomistic molecular dynamics (MD) simulation and data mining to build a high-throughput structural bioinformatics pipeline to study the effects of individual and combination SHMs on antibody conformation, flexibility, stability, and affinity. By applying this pipeline, we characterized a common mechanism of modulation of heavy-light pairing orientation by frequent SHMs at framework positions 39_H, 91_H, 38_L, and 87_L through disruption of a conserved hydrogen-bond network. Q39L_H alone and in combination with light chain framework 4 (FWR4_L) insertions further modulated the elbow angle between variable and constant domains of many antibodies, resulting in improved binding affinity for a subset of anti-HIV-1 antibodies. Q39L_H also alleviated aggregation induced by FWR4_L insertion, suggesting remote epistasis between these SHMs. Altogether, this study provides tools and insights for understanding antibody affinity maturation and for engineering functionally improved antibodies.

Computational modelling of potentially emerging SARS-CoV-2 spike protein RBDs mutations with higher binding affinity towards ACE2: A structural modelling study

The spike protein of SARS-CoV-2 and the host ACE2 receptor plays a vital role in the entry to the cell. Among which the hotspot residue 501 is continuously subjected to positive selection pressure and induces unusual virulence. Keeping in view the importance of the hot spot residue 501, we predicted the potentially emerging structural variants of 501 residue. We analyzed the binding pattern of wild type and mutants (Spike RBD) to the ACE2 receptor by deciphering variations in the amino acids' interaction networks by graph kernels along with evolutionary, network metrics, and energetic information. Our analysis revealed that N501I, N501T, and N501V increase the binding affinity and alter the intra and inter-residue bonding networks. The N501T has shown strong positive selection and fitness in other animals. Docking results and repeated simulations (three times) confirmed the structural stability and tighter binding of these three variants, correlated with the previous results following the global stability trend. Consequently, we reported three variants N501I, N501T, and N501V could worsen the situation further if they emerged. The relations between the viral fitness and binding affinity is a complicated game thus the emergence of high affinity mutations in the SARS-CoV-2 RBD brings up the question of whether or not positive selection favours these mutations or not?

Fundamental Tick Vaccinomic Approach to Evade Host Autoimmune Reaction

Ticks are obligate hematophagous ectoparasites that infect domestic animals, humans, and wildlife. Ticks can transmit a wide range of pathogens (viruses, rickettsia, bacteria, parasites, etc.), and some of those are of zoonotic importance. Tick-borne diseases have a negative economic impact in several tropical and subtropical countries. With climate change, tick distribution and tick-associated pathogens have increased. Currently, tick control procedures have more environmental drawbacks and there are pitfalls in vaccination process. Since vaccinations have helped to prevent several diseases and infections, several vaccination trials are ongoing to control ticks and tick-borne pathogens. However, autoimmune reactions to vaccinations are reported as an adverse reaction since vaccines were used to protect against disease in humans and animals. The antibodies against the vaccine antigen might harm similar antigen in the host. Therefore, in this chapter, we attempt to shed light on the importance of raising awareness of possible adverse events associated with vaccinations and the methods that should be used to address this problem. In silico and lab work should be performed ahead of the vaccination process to evaluate the vaccine candidates and avoid the vaccination opposing consequences.