Very fast empirical prediction and rationalization of protein pKa values

Comparison of two peroxidases with high potential for biotechnology applications - HRP vs. APEX2

Peroxidases are essential elements in many biotechnological applications. An especially interesting concept involves split enzymes, where the enzyme is separated into two smaller and inactive proteins that can dimerize into a fully active enzyme. Such split forms were developed for the horseradish peroxidase (HRP) and ascorbate peroxidase (APX) already. Both peroxidases have a high potential for biotechnology applications. In the present study, we performed biophysical comparisons of these two peroxidases and their split analogues. The active site availability is similar for all four structures. The split enzymes are comparable in stability with their native analogues, meaning that they can be used for further biotechnology applications. Also, the tertiary structures of the two peroxidases are similar. However, differences that might help in choosing one system over another for biotechnology applications were noticed. The main difference between the two systems is glycosylation which is not present in the case of APX/sAPEX2, while it has a high impact on the HRP/sHRP stability. Further differences are calcium ions and cysteine bridges that are present only in the case of HRP/sHRP. Finally, computational results identified sAPEX2 as the systems with the smallest structural variations during molecular dynamics simulations showing its dominant stability comparing to other simulated proteins. Taken all together, the sAPEX2 system has a high potential for biotechnological applications due to the lack of glycans and cysteines, as well as due to high stability.

Differential modulation of mutant CALR and JAK2 V617F-driven oncogenesis by HLA genotype in myeloproliferative neoplasms

It has been demonstrated previously that human leukocyte antigen class I (HLA-I) and class II (HLA-II) alleles may modulate JAK2 V617F and CALR mutation (CALRmut)-associated oncogenesis in myeloproliferative neoplasms (MPNs). However, the role of immunogenetic factors in MPNs remains underexplored. We aimed to investigate the potential involvement of HLA genes in CALRmut+ MPNs. High-resolution genotyping of HLA-I and -II loci was conducted in 42 CALRmut+ and 158 JAK2 V617F+ MPN patients and 1,083 healthy controls. A global analysis of the diversity of HLA-I genotypes revealed no significant differences between CALRmut+ patients and controls. However, one HLA-I allele (C*06:02) showed an inverse correlation with presence of CALR mutation. A meta-analysis across independent cohorts and healthy individuals from the 1000 Genomes Project confirmed an inverse correlation between the presentation capabilities of the HLA-I loci for JAK2 V617F and CALRmut-derived peptides in both patients and healthy individuals. scRNA-Seq analysis revealed low expression of TAP1 and CIITA genes in CALRmut+ hematopoietic stem and progenitor cells. In conclusion, the HLA-I genotype differentially restricts JAK2 V617F and CALRmut-driven oncogenesis potentially explaining the mutual exclusivity of the two mutations and differences in their presentation latency. These findings have practical implications for the development of neoantigen-based vaccines in MPNs.

Assessment of the in vitro anti-diabetic activity with molecular dynamic simulations of limonoids isolated from Adalia lemon peels

Limonoids are important constituents of citrus that have a significant impact on promoting human health. Therefore, the primary focus of this research was to assess the overall limonoid content and isolate limonoids from Adalia lemon (Citrus limon L.) peels for their potential use as antioxidants and anti-diabetic agents. The levels of limonoid aglycones in the C. limon peel extract were quantified through a colorimetric assay, revealing a concentration of 16.53 ± 0.93 mg/L limonin equivalent. Furthermore, the total concentration of limonoid glucosides was determined to be 54.38 ± 1.02 mg/L. The study successfully identified five isolated limonoids, namely limonin, deacetylnomilin, nomilin, obacunone 17-O-β-D-glucopyranoside, and limonin 17-O-β-D-glucopyranoside, along with their respective yields. The efficacy of the limonoids-rich extract and the five isolated compounds was evaluated at three different concentrations (50, 100, and 200 µg/mL). It was found that both obacunone 17-O-β-D-glucopyranoside and limonin 17-O-β-D-glucopyranoside possessed the highest antioxidant, free radical scavenging, and anti-diabetic activities, followed by deacetylnomilin, and then the limonoids-rich extract. The molecular dynamic simulations were conducted to predict the behavior of the isolated compounds upon binding to the protein's active site, as well as their interaction and stability. The results revealed that limonin 17-O-β-D-glucopyranoside bound to the protein complex system exhibited a relatively more stable conformation than the Apo system. The analysis of Solvent Accessible Surface Area (SASA), in conjunction with the data obtained from Root-Mean-Square Deviation (RMSD), Root-Mean-Square Fluctuation (RMSF), and Radius of Gyration (ROG) computations, provided further evidence that the limonin 17-O-β-D-glucopyranoside complex system remained stable within the catalytic domain binding site of the human pancreatic alpha-amylase (HPA)-receptor. The research findings suggest that the limonoids found in Adalia lemon peels have the potential to be used as effective natural substances in creating innovative therapeutic treatments for conditions related to oxidative stress and disorders in carbohydrate metabolism.

Computational study of the mechanism of a polyurethane esterase A (PueA) from Pseudomonas chlororaphis

The effective management of plastic waste has become a global imperative, given our reliance on a linear model in which plastics are manufactured, used once, and then discarded. This has led to the pervasive accumulation of plastic debris in landfills and environmental contamination. Recognizing this issue, numerous initiatives are underway to address the environmental repercussions associated with plastic disposal. In this study, we investigate the possible molecular mechanism of polyurethane esterase A (PueA), which has been previously identified as responsible for the degradation of a polyester polyurethane (PU) sample in Pseudomonas chlororaphis, as an effort to develop enzymatic biodegradation solutions. After generating the unsolved 3D structure of the protein by AlphaFold2 from its known genome, the enzymatic hydrolysis of the same model PU compound previously used in experiments has been explored employing QM/MM molecular dynamics simulations. This required a preliminary analysis of the 3D structure of the apo-enzyme, identifying the putative active site, and the search for the optimal protein-substrate binding site. Finally, the resulting free energy landscape indicates that wild-type PueA can degrade PU chains, although with low-level activity. The reaction takes place by a characteristic four-step path of the serine hydrolases, involving an acylation followed by a diacylation step. Energetics and structural analysis of the evolution of the active site along the reaction suggests that PueA can be considered a promising protein scaffold for further development to achieve efficient biodegradation of PU.

High terephthalic acid purity: Effective polyethylene terephthalate degradation process based on pH regulation with dual-function hydrolase

Eco-friendly enzymatic-recycling has been widely utilized in tackling plastic pollution. However, the limited activity on the polyethylene terephthalate (PET) degradation product mono-hydroxyethyl terephthalate (MHET) leads to the formation of heterogeneous hydrolysis products, resulting in PET downcycling. Herein, by applying a dual-function PET hydrolase IsPETasePA with balanced PET and MHET degradation efficiency, an effective PET hydrolysis process was developed to enhance the terephthalic acid (TPA) product purity. Firstly, the impact of pH on the catalytic activity of IsPETasePA revealed that the pH reduction caused by TPA generation hindered the complete conversion of MHET to TPA. Further investigation of the catalytic mechanism showed that the pH-induced protonation of His208 in the catalytic triad destabilized the interaction between IsPETasePA and MHET. Thus, by introducing pH regulation strategy on the bifunctional IsPETasePA, the single-enzyme process could achieve high-purity TPA recovery (>99 %). Overall, this work ensured the high-quality PET enzymatic-recycling for effectively addressing plastic pollution.

Blocking the PD-1 signal transduction by occupying the phosphorylated ITSM recognition site of SHP-2

Targeting the PD-1/PD-L1 axis with small-molecular inhibitors is a promising approach for immunotherapy. Here, we identify a natural pentacyclic triterpenoid, Pygenic Acid A (PA), as a PD-1 signaling inhibitor. PA exerts anti-tumor activity in hPD-1 knock-in C57BL/6 mice and enhances effector functions of T cells to promote immune responses by disrupting the PD-1 signaling transduction. Furthermore, we identify SHP-2 as the direct molecular target of PA for inhibiting the PD-1 signaling transduction. Subsequently, mechanistic studies suggest that PA binds to a new druggable site in the phosphorylated PD-1 ITSM recognition site of SHP-2, inhibiting the recruitment of SHP-2 by PD-1. Taken together, our findings demonstrate that PA has a potential application in cancer immunotherapy and occupying the phosphorylated ITSM recognition site of SHP-2 may serve as an alternative strategy to develop PD-1 signaling inhibitors. In addition, our success in target recognition provides a paradigm of target identification and confirmation for natural products.

Magnesium induced structural reorganization in the active site of adenylate kinase

Phosphoryl transfer is a fundamental reaction in cellular signaling and metabolism that requires Mg2+ as an essential cofactor. While the primary function of Mg2+ is electrostatic activation of substrates, such as ATP, the full spectrum of catalytic mechanisms exerted by Mg2+ is not known. In this study, we integrate structural biology methods, molecular dynamic (MD) simulations, phylogeny, and enzymology assays to provide molecular insights into Mg2+-dependent structural reorganization in the active site of the metabolic enzyme adenylate kinase. Our results demonstrate that Mg2+ induces a conformational rearrangement of the substrates (ATP and ADP), resulting in a 30° adjustment of the angle essential for reversible phosphoryl transfer, thereby optimizing it for catalysis. MD simulations revealed transitions between conformational substates that link the fluctuation of the angle to large-scale enzyme dynamics. The findings contribute detailed insight into Mg2+ activation of enzymes and may be relevant for reversible and irreversible phosphoryl transfer reactions.

Deepening insights into cholinergic agents for intraocular pressure reduction: systems genetics, molecular modeling, and in vivo perspectives

Parasympathetic activation in the anterior eye segment regulates various physiological functions. This process, mediated by muscarinic acetylcholine receptors, also impacts intraocular pressure (IOP) through the trabecular meshwork. While FDA-approved M3 muscarinic receptor (M3R) agonists exist for IOP reduction, their systemic cholinergic adverse effects pose limitations in clinical use. Therefore, advancing our understanding of the cholinergic system in the anterior segment of the eye is crucial for developing additional IOP-reducing agents with improved safety profiles. Systems genetics analyses were utilized to explore correlations between IOP and the five major muscarinic receptor subtypes. Molecular docking and dynamics simulations were applied to human M3R homology model using a comprehensive set of human M3R ligands and 1,667 FDA-approved or investigational drugs. Lead compounds from the modeling studies were then tested for their IOP-lowering abilities in mice. Systems genetics analyses unveiled positive correlations in mRNA expressions among the five major muscarinic receptor subtypes, with a negative correlation observed only in M3R with IOP. Through modeling studies, rivastigmine and edrophonium emerged as the most optimally suited cholinergic drugs for reducing IOP via a potentially distinct mechanism from pilocarpine or physostigmine. Subsequent animal studies confirmed comparable IOP reductions among rivastigmine, edrophonium, and pilocarpine, with longer durations of action for rivastigmine and edrophonium. Mild cholinergic adverse effects were observed with pilocarpine and rivastigmine but absent with edrophonium. These findings advance ocular therapeutics, suggesting a more nuanced role of the parasympathetic system in the anterior eye segment for reducing IOP than previously thought.

Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations

Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-μs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.

Dual mechanisms contribute to enhanced voltage dependence of an electric fish potassium channel

The voltage dependence of different voltage-gated potassium channels, described by the voltage at which half of the channels are open (V1/2), varies over a range of 80 mV and is influenced by factors such as the number of positive gating charges and the identity of the hydrophobic amino acids in the channel's voltage sensor (S4). Here we explore by experimental manipulations and molecular dynamics simulation the contributions of two derived features of an electric fish potassium channel (Kv1.7a) that is among the most voltage-sensitive Shaker family potassium channels known. These are a patch of four contiguous negatively charged glutamates in the S3-S4 extracellular loop and a glutamate in the S3b helix. We find that these negative charges affect V1/2 by separate, complementary mechanisms. In the closed state, the S3-S4 linker negative patch reduces the membrane surface charge biasing the channel to enter the open state while, upon opening, the negative amino acid in the S3b helix faces the second (R2) gating charge of the voltage sensor electrostatically biasing the channel to remain in the open state. This work highlights two evolutionary novelties that illustrate the potential influence of negatively charged amino acids in extracellular loops and adjacent helices to voltage dependence.

Pepper catalase: a broad analysis of its modulation during fruit ripening and by nitric oxide

Catalase is a major antioxidant enzyme located in plant peroxisomes that catalyzes the decomposition of H2O2. Based on our previous transcriptomic (RNA-Seq) and proteomic (iTRAQ) data at different stages of pepper (Capsicum annuum L.) fruit ripening and after exposure to nitric oxide (NO) enriched atmosphere, a broad analysis has allowed us to characterize the functioning of this enzyme. Three genes were identified, and their expression was differentially modulated during ripening and by NO gas treatment. A dissimilar behavior was observed in the protein expression of the encoded protein catalases (CaCat1-CaCat3). Total catalase activity was down-regulated by 50% in ripe (red) fruits concerning immature green fruits. This was corroborated by non-denaturing polyacrylamide gel electrophoresis, where only a single catalase isozyme was identified. In vitro analyses of the recombinant CaCat3 protein exposed to peroxynitrite (ONOO-) confirmed, by immunoblot assay, that catalase underwent a nitration process. Mass spectrometric analysis identified that Tyr348 and Tyr360 were nitrated by ONOO-, occurring near the active center of catalase. The data indicate the complex regulation at gene and protein levels of catalase during the ripening of pepper fruits, with activity significantly down-regulated in ripe fruits. Nitration seems to play a key role in this down-regulation, favoring an increase in H2O2 content during ripening. This pattern can be reversed by the exogenous NO application. While plant catalases are generally reported to be tetrameric, the analysis of the protein structure supports that pepper catalase has a favored quaternary homodimer nature. Taken together, data show that pepper catalase is down-regulated during fruit ripening, becoming a target of tyrosine nitration, which provokes its inhibition.

Unrevealing the multitargeted potency of 3-1-BCMIYPPA against lung cancer structural maintenance and suppression proteins through pharmacokinetics, QM-DFT, and multiscale MD simulation studies

Lung cancer, a relentless and challenging disease, demands unwavering attention in drug design research. Single-target drugs have yielded limited success, unable to effectively address this malignancy's profound heterogeneity and often developed resistance. Consequently, the clarion call for lung cancer drug design echoes louder than ever, and multitargeted drug design emerges as an imperative approach in this landscape, which is done by concurrently targeting multiple proteins and pathways and offering a beacon of hope. This study is focused on the multitargeted drug designing approach by identifying drug candidates against human cyclin-dependent kinase-2, SRC-2 domains of C-ABL, epidermal growth factor and receptor extracellular domains, and insulin-like growth factor-1 receptor kinase. We performed the multitargeted molecular docking studies of Drug Bank compounds using HTVS, SP and XP algorithms and poses filter with MM\GBSA against all proteins and identified DB02504, namely [3-(1-Benzyl-3-Carbamoylmethyl-2-Methyl-1h-Indol-5-Yloxy)-Propyl-]-Phosphonic Acid (3-1-BCMIYPPA) as multitargeted lead with docking and MM\GBSA score range from -8.242 to -6.274 and -28.2 and -44.29 Kcal/mol, respectively. Further, the QikProp-based pharmacokinetic computations and QM-based DFT showed acceptance results against standard values, and interaction fingerprinting reveals that THR, MET, GLY, VAL, LEU, GLU and ASP were among the most interacting residues. The NPT ensemble-based 100ns MD simulation in a neutralised state with an SPC water model has also shown a stable performance and produced deviation and fluctuations <2Å with huge interactions, making it a promising multitargeted drug candidate-however, experimental studies are suggested.

Intermolecular Electrostatic Interactions in Cytochrome c Protein Monolayer on Montmorillonite Alumosilicate Surface: A Positive Cooperative Effect

Montmorillonite (MM) crystal nanoplates acquire anticancer properties when coated with the mitochondrial protein cytochrome c (cytC) due to the cancer cells' capability to phagocytize cytC-MM colloid particles. The introduced exogenous cytC initiates apoptosis: an irreversible cascade of biochemical reactions leading to cell death. In the present research, we investigate the organization of the cytC layer on the MM surface by employing physicochemical and computer methods-microelectrophoresis, static, and electric light scattering-to study cytC adsorption on the MM surface, and protein electrostatics and docking to calculate the local electric potential and Gibbs free energy of interacting protein globules. The found protein concentration dependence of the adsorbed cytC quantity is nonlinear, manifesting a positive cooperative effect that emerges when the adsorbed cytC globules occupy more than one-third of the MM surface. Computer analysis reveals that the cooperative effect is caused by the formation of protein associates in which the cytC globules are oriented with oppositely charged surfaces. The formation of dimers and trimers is accompanied by a strong reduction in the electrostatic component of the Gibbs free energy of protein association, while the van der Waals component plays a secondary role.

SilE-R and SilE-S-DABB Proteins Catalying Enantiospecific Hydrolysis of Organosilyl Ethers

Silyl ethers fulfil a fundamental role in synthetic organic chemistry as protecting groups and their selective cleavage is an important factor in their application. We present here for the first time two enzymes, SilE-R and SilE-S, which are able to hydrolyse silyl ethers. They belong to the stress-response dimeric A/B barrel domain (DABB) family and are able to cleave the Si-O bond with opposite enantiopreference. Silyl ethers containing aromatic, cyclic or aliphatic alcohols and, depending on the alcohol moiety, silyl functions as large as TBDMS are accepted. The X-ray crystal structure of SilE-R, determined to a resolution of 1.98 Å, in combination with mutational studies, revealed an active site featuring two histidine residues, H8 and H79, which likely act synergistically as nucleophile and Brønsted base in the hydrolytic mechanism, which has not previously been described for enzymes. Although the natural function of SilE-R and SilE-S is unknown, we propose that these 'silyl etherases' may have significant potential for synthetic applications.

Targeted removal of the FA2 site on human albumin prevents fatty acid-mediated inhibition of Zn2+ binding

Zinc is required for virtually all biological processes. In plasma, Zn2+ is predominantly transported by human serum albumin (HSA), which possesses two Zn2+-binding sites of differing affinities (sites A and B). Fatty acids (FAs) are also transported by HSA, with seven structurally characterized FA-binding sites (named FA1-FA7) known. FA binding inhibits Zn2+-HSA interactions, in a manner that can impact upon hemostasis and cellular zinc uptake, but the degree to which binding at specific FA sites contributes to this inhibition is unclear. Wild-type HSA and H9A, H67A, H247A, and Y150F/R257A/S287A (FA2-KO) mutant albumins were expressed in Pichia pastoris. Isothermal titration calorimetry studies revealed that the Zn2+-binding capacity at the high-affinity Zn2+ site (site A) was reduced in H67A and H247A mutants, with site B less affected. The H9A mutation decreased Zn2+ binding at the lower-affinity site, establishing His9 as a site B ligand. Zn2+ binding to HSA and H9A was compromised by palmitate, consistent with FA binding affecting site A. 13C-NMR experiments confirmed that the FA2-KO mutations prohibited FA binding at site FA2. Zn2+ binding to the FA2-KO mutant was unaffected by myristate, suggesting binding at FA2 is solely responsible for inhibition. Molecular dynamics studies identified the steric obstruction exerted by bound FA in site FA2, which impedes the conformational change from open (FA-loaded) to closed (FA-free) states, required for Zn2+ to bind at site A. The successful targeting of the FA2 site will aid functional studies exploring the interplay between circulating FA levels and plasma Zn2+ speciation in health and disease.

Human indoleamine-2,3-dioxygenase 2 cofactor lability and low substrate affinity explained by homology modeling, molecular dynamics and molecular docking

The human indoleamine-2,3-dioxygenase 2 (hIDO2) protein is growing of interest as it is increasingly implicated in multiple diseases (cancer, autoimmune diseases, COVID-19). However, it is only poorly reported in the literature. Its mode of action remains unknown because it does not seem to catalyze the reaction for which it is attributed: the degradation of the L-Tryptophan into N-formyl-kynurenine. This contrasts with its paralog, the human indoleamine-2,3-dioxygenase 1 (hIDO1), which has been extensively studied in the literature and for which several inhibitors are already in clinical trials. Yet, the recent failure of one of the most advanced hIDO1 inhibitors, the Epacadostat, could be caused by a still unknown interaction between hIDO1 and hIDO2. In order to better understand the mechanism of hIDO2, and in the absence of experimental structural data, a computational study mixing homology modeling, Molecular Dynamics, and molecular docking was conducted. The present article highlights an exacerbated lability of the cofactor as well as an inadequate positioning of the substrate in the active site of hIDO2, which might bring part of an answer to its lack of activity.Communicated by Ramaswamy H. Sarma.

Characterization of the Bottlenecks and Pathways for Inhibitor Dissociation from [NiFe] Hydrogenase

[NiFe] hydrogenases can act as efficient catalysts for hydrogen oxidation and biofuel production. However, some [NiFe] hydrogenases are inhibited by gas molecules present in the environment, such as O2 and CO. One strategy to engineer [NiFe] hydrogenases and achieve O2- and CO-tolerant enzymes is by introducing point mutations to block the access of inhibitors to the catalytic site. In this work, we characterized the unbinding pathways of CO in the complex with the wild-type and 10 different mutants of [NiFe] hydrogenase from Desulfovibrio fructosovorans using τ-random accelerated molecular dynamics (τRAMD) to enhance the sampling of unbinding events. The ranking provided by the relative residence times computed with τRAMD is in agreement with experiments. Extensive data analysis of the simulations revealed that from the two bottlenecks proposed in previous studies for the transit of gas molecules (residues 74 and 122 and residues 74 and 476), only one of them (residues 74 and 122) effectively modulates diffusion and residence times for CO. We also computed pathway probabilities for the unbinding of CO, O2, and H2 from the wild-type [NiFe] hydrogenase, and we observed that while the most probable pathways are the same, the secondary pathways are different. We propose that introducing mutations to block the most probable paths, in combination with mutations to open the main secondary path used by H2, can be a feasible strategy to achieve CO and O2 resistance in the [NiFe] hydrogenase from Desulfovibrio fructosovorans.

3-Chymotrypsin-like Protease (3CLpro) of SARS-CoV-2: Validation as a Molecular Target, Proposal of a Novel Catalytic Mechanism, and Inhibitors in Preclinical and Clinical Trials

Proteases represent common targets in combating infectious diseases, including COVID-19. The 3-chymotrypsin-like protease (3CLpro) is a validated molecular target for COVID-19, and it is key for developing potent and selective inhibitors for inhibiting viral replication of SARS-CoV-2. In this review, we discuss structural relationships and diverse subsites of 3CLpro, shedding light on the pivotal role of dimerization and active site architecture in substrate recognition and catalysis. Our analysis of bioinformatics and other published studies motivated us to investigate a novel catalytic mechanism for the SARS-CoV-2 polyprotein cleavage by 3CLpro, centering on the triad mechanism involving His41-Cys145-Asp187 and its indispensable role in viral replication. Our hypothesis is that Asp187 may participate in modulating the pKa of the His41, in which catalytic histidine may act as an acid and/or a base in the catalytic mechanism. Recognizing Asp187 as a crucial component in the catalytic process underscores its significance as a fundamental pharmacophoric element in drug design. Next, we provide an overview of both covalent and non-covalent inhibitors, elucidating advancements in drug development observed in preclinical and clinical trials. By highlighting various chemical classes and their pharmacokinetic profiles, our review aims to guide future research directions toward the development of highly selective inhibitors, underscore the significance of 3CLpro as a validated therapeutic target, and propel the progression of drug candidates through preclinical and clinical phases.

Structure and inhibition of the human lysosomal transporter Sialin

Sialin, a member of the solute carrier 17 (SLC17) transporter family, is unique in its ability to transport not only sialic acid using a pH-driven mechanism, but also transport mono and diacidic neurotransmitters, such as glutamate and N-acetylaspartylglutamate (NAAG), into synaptic vesicles via a membrane potential-driven mechanism. While most transporters utilize one of these mechanisms, the structural basis of how Sialin transports substrates using both remains unclear. Here, we present the cryogenic electron-microscopy structures of human Sialin: apo cytosol-open, apo lumen-open, NAAG-bound, and inhibitor-bound. Our structures show that a positively charged cytosol-open vestibule accommodates either NAAG or the Sialin inhibitor Fmoc-Leu-OH, while its luminal cavity potentially binds sialic acid. Moreover, functional analyses along with molecular dynamics simulations identify key residues in binding sialic acid and NAAG. Thus, our findings uncover the essential conformational states in NAAG and sialic acid transport, demonstrating a working model of SLC17 transporters.

Proton Coulomb Blockade Effect Involving Covalent Oxygen-Hydrogen Bond Switching

Instead of the canonical Grotthuss mechanism, we show that a knock-on proton transport process is preferred between organic functional groups (e.g., -COOH and -OH) and adjacent water molecules in biological proton channel and synthetic nanopores through comprehensive quantum and classical molecular dynamics simulations. The knock-on process is accomplished by the switching of covalent O─H bonds of the functional group under externally applied electric fields. The proton transport through the synthetic nanopore exhibits nonlinear current-voltage characteristics, suggesting an unprecedented proton Coulomb blockade effect. These findings not only enhance the understanding of proton transport in nanoconfined systems but also pave the way for the design of a variety of proton-based nanofluidic devices.

Facile and Novel Synthetic Approach, Molecular Docking, Molecular Dynamics, and Drug-Likeness Evaluation of 9-Substituted Acridine Derivatives as Dual Anticancer and Antimicrobial Agents

In the present study, numerous acridine derivatives A1-A20 were synthesized via aromatic nucleophilic substitution (SNAr) reaction of 9-chloroacridine with carbonyl hydrazides, amines, or phenolic derivatives depending upon facile, novel, and eco-friendly approaches (Microwave and ultrasonication assisted synthesis). The structures of the new compounds were elucidated using spectroscopic methods. The title products were assessed for their antimicrobial, antioxidant, and antiproliferative activities using numerous assays. Promisingly, the investigated compounds mainstream revealed promising antibacterial and anticancer activities. Thereafter, the investigated compounds' expected mode of action was debated by using an array of in silico studies. Compounds A2 and A3 were the most promising antimicrobial agents, while compounds A2, A5, and A7 revealed the most cytotoxic activities. Accordingly, RMSD, RMSF, Rg, and SASA analyses of compounds A2 and A3 were performed, and MMPBSA was calculated. Lastly, the ADMET (absorption, distribution, metabolism, excretion, and toxicity) analyses of the novel acridine derivatives were investigated. The tested compounds' existing screening results afford an inspiring basis leading to developing new compelling antimicrobial and anticancer agents based on the acridine scaffold.

Characterization of host receptor interaction with envelop protein of Kyasanur forest disease virus and predicting suitable epitopes for vaccine candidate

Kyasanur forest disease (KFD) is a zoonotic disease that is endemic to southern India and caused by KFD virus (KFDV) belonging to the family Flaviviridae. Humans are the dead-end host of the KFDV life cycle. The absence of effective treatment strategies against KFD can be attributed to a lack of studies on the mechanistic part of the spread of the disease. Hypothesizing molecular etiological similarity of KFDV to other well characterized flaviviruses, such as dengue virus (DENV), we focused on predicting the target receptor protein(s) in host and provided molecular basis of receptor-mediated recognition of the human host by KFDV envelop protein (EKFDV), drawing from the extant knowledge on the dengue counterpart, EDENV. Indeed, in silico approach helped to identify that the EKFDV structure closely resembles the EDENV structure and indicated DC-SIGN and/or Mannose receptors to be the plausible target host receptors. Immune-informatics approach aided in predicting 10 epitopes from E, NS1, NS2A, and NS2B proteins of the KFDV-P9605 genotype for vaccine design against KFDV. Further, molecular dynamics simulation (MDS) analyses of their complexes with human leukocyte antigens (HLAs) identified the epitopes DISLTCRVT and YAMEIRPVH as two high ranking candidates for vaccine design.Communicated by Ramaswamy H. Sarma.

Ali S, Efferth T, Hegazy MEF. Salvimulticanol from Salvia multicaulis suppresses LPS-induced inflammation in RAW264.7 macrophages: in vitro and in silico studies

Sustained inflammatory responses can badly affect several vital organs and lead to chronic inflammation-related disorders, such as atherosclerosis, pneumonia, rheumatoid arthritis, obesity, diabetes, Alzheimer's disease, and cancers. Salvia multicaulis is one of the widely distributed plants that contains several biologically active phytochemicals and diterpenoids with anti-inflammatory effects. Therefore, finding alternative and safer natural plant-extracted compounds with good curative anti-inflammatory efficiencies is an urgent need for the clinical treatment of inflammation-related diseases. In the current study, S. multicaulis Vahl was used to extract and isolate two compounds identified as salvimulticanol and candesalvone B methyl ester to examine their effects against inflammation in murine macrophage RAW264.7 cells that were induced by lipopolysaccharide (LPS). Accordingly, after culturing RAW264.7 cells and induction of inflammation by LPS (100 ng/ml), cells were exposed to different concentrations (9, 18, 37.5, 75, and 150 µM) of each compound. Then, Griess assay for detection of nitric oxide (NO) levels and western blotting for the determination of inducible nitric oxide synthase (iNOS) expression were performed. Molecular docking and molecular dynamics (MD) simulation studies were employed to investigate the anti-inflammatory mechanism. Our obtained results validated that the level of NO was significantly decreased in the macrophage cell suspensions as a response to salvimulticanol treatment in a dose-dependent manner (IC50: 25.1 ± 1.2 µM) as compared to the methyl ester of candesalvone B which exerted a weaker inhibition (IC50: 69.2 ± 3.0 µM). This decline in NO percentage was comparable with a down-regulation of iNOS expression by western blotting. Salvimulticanol strongly interacted with both the Toll-like receptor 4 (TLR4)/myeloid differentiation factor 2 (MD-2) complex and the inhibitor of nuclear factor kappa-B (NF-κB) kinase subunit beta (IKKβ) to disrupt their inflammatory activation due to the significant hydrogen bonds and effective interactions with amino acid residues present in the target proteins' active sites. S.multicaulis is a rich natural source of the aromatic abietane diterpenoid, salvimulticanol, which exerted a strong anti-inflammatory effect through targeting iNOS and diminishing NO production in LPS-induced RAW264.7 cells in a mechanism that is dependent on the inhibition of TLR4-MD-2 and IKKβ as activators of the classical NF-κB-mediated inflammatory pathway.

Automated preparation of nanoscopic structures: Graph-based sequence analysis, mismatch detection, and pH-consistent protonation with uncertainty estimates

Structure and function in nanoscale atomistic assemblies are tightly coupled, and every atom with its specific position and even every electron will have a decisive effect on the electronic structure, and hence, on the molecular properties. Molecular simulations of nanoscopic atomistic structures therefore require accurately resolved three-dimensional input structures. If extracted from experiment, these structures often suffer from severe uncertainties, of which the lack of information on hydrogen atoms is a prominent example. Hence, experimental structures require careful review and curation, which is a time-consuming and error-prone process. Here, we present a fast and robust protocol for the automated structure analysis and pH-consistent protonation, in short, ASAP. For biomolecules as a target, the ASAP protocol integrates sequence analysis and error assessment of a given input structure. ASAP allows for pK a  prediction from reference data through Gaussian process regression including uncertainty estimation and connects to system-focused atomistic modeling described in Brunken and Reiher (J. Chem. Theory Comput. 16, 2020, 1646). Although focused on biomolecules, ASAP can be extended to other nanoscopic objects, because most of its design elements rely on a general graph-based foundation guaranteeing transferability. The modular character of the underlying pipeline supports different degrees of automation, which allows for (i) efficient feedback loops for human-machine interaction with a low entrance barrier and for (ii) integration into autonomous procedures such as automated force field parametrizations. This facilitates fast switching of the pH-state through on-the-fly system-focused reparametrization during a molecular simulation at virtually no extra computational cost.

Origin of Enantioselectivity in Engineered Cytochrome c-Catalyzed Carbon-Radical FePP Hydrolysis Revealed Using QM/MM (ABEEM Polarizable Force Field) and MD Simulations

The origin of highly efficient asymmetric aminohydroxylation of styrene catalyzed by engineered cytochrome c is investigated by the developed Atom-Bond Electronegativity Equalization Method polarizable force field (ABEEM PFF), which is a combined outcome of electronic and steric effects. Model molecules were used to establish the charge parameters of the ABEEM PFF, for which the bond-stretching and angle-bending parameters were obtained by using a combination of modified Seminario and scan methods. The interactions between carbon-radical Fe-porphyrin (FePP) and waters are simulated by molecular dynamics, which shows a clear preference for the pre-R over the pre-S. This preference is attributed to the hydrogen-bond between the mutated 100S and 101P residues as well as van der Waals interactions, enforcing a specific conformation of the carbon-radical FePP complex within the binding pocket. Meanwhile, the hydrogen-bond between water and the nitrogen atom in the active intermediate dictates the stereochemical outcome. Quantum mechanics/molecular mechanics (QM/MM (ABEEM PFF)) and free-energy perturbation calculations elucidate that the 3RTS is characterized by sandwich-like structure among adjacent amino acid residues, which exhibits greater stability than crowed arrangement in 3STS and enables the R enantiomer to form more favorably. Thus, this study provides mechanistic insight into the catalytic reaction of hemoproteins.

Molecular Dynamics Insights into the Aggregation Behavior of N-Terminal β-Lactoglobulin Peptides

β-lactoglobulin (BLG) forms amyloid-like aggregates at high temperatures, low pH, and low ionic strengths. At a pH below 2, BLG undergoes hydrolysis into peptides, with N-terminal peptides 1-33 and 1-52 being prone to fibrillization, forming amyloid-like fibrils. Due to their good mechanical properties, BLG amyloids demonstrate great potential for diverse applications, including biosensors, nanocomposites, and catalysts. Consequently, further studies are essential to comprehensively understand the factors governing the formation of BLG amyloid-like morphologies. In this study, all-atom molecular dynamics simulations were employed to explore the aggregation of N-terminal 1-33 and 1-52 BLG peptides under conditions of pH 2 and at 10 mM NaCl concentration. The simulations revealed that the peptides spontaneously assembled into aggregates of varying sizes. The aggregation process was enabled by the low charge of peptides and the presence of hydrophobic residues within them. As the peptides associated into aggregates, there was a concurrent increase in β-sheet structures and the establishment of hydrogen bonds, enhancing the stability of the aggregates. Notably, on average, 1-33 peptides formed larger aggregates compared to their 1-52 counterparts, while the latter exhibited a slightly higher content of β-sheets and higher cluster orderliness. The applied approach facilitated insights into the early stages of amyloid-like aggregation and molecular-level insight into the formation of β-sheets, which serve as nucleation points for further fibril growth.

DeepKa Web Server: High-Throughput Protein pKa Prediction

DeepKa is a deep-learning-based protein pKa predictor proposed in our previous work. In this study, a web server was developed that enables online protein pKa prediction driven by DeepKa. The web server provides a user-friendly interface where a single step of entering a valid PDB code or uploading a PDB format file is required to submit a job. Two case studies have been attached in order to explain how pKa's calculated by the web server could be utilized by users. Finally, combining the web server with post processing as described in case studies, this work suggests a quick workflow of investigating the relationship between protein structure and function that are pH dependent. The web server of DeepKa is freely available at http://www.computbiophys.com/DeepKa/main.

Theoretical Analysis of Selectivity Differences in Ketoreductases toward Aldehyde and Ketone Carbonyl Groups

Lactobacillus kefir alcohol dehydrogenase (LkADH) and ketoreductase from Chryseobacterium sp. CA49 (ChKRED12) exhibit different chemoselectivity and stereoselectivity toward a substrate with both keto and aldehyde carbonyl groups. LkADH selectively reduces the keto carbonyl group while retaining the aldehyde carbonyl group, producing optically pure R-alcohols. In contrast, ChKRED12 selectively reduces the aldehyde group and exhibits low reactivity toward ketone carbonyls. This study investigated the structural basis for these differences and the role of specific residues in the active site. Molecular dynamics (MD) simulations and quantum chemical calculations were used to investigate the interactions between the substrate and the enzymes and the essential cause of this phenomenon. The present study has revealed that LkADH and ChKRED12 exhibit significant differences in the structure of their respective active pockets, which is a crucial determinant of their distinct chemoselectivity toward the same substrate. Moreover, residues N89, N113, and E144 within LkADH as well as Q151 and D190 within ChKRED12 have been identified as key contributors to substrate stabilization within the active pocket through electrostatic interactions and van der Waals forces, followed by hydride transfer utilizing the coenzyme NADPH. Furthermore, the enantioselectivity mechanism of LkADH has been elucidated using quantum chemical methods. Overall, these findings not only provide fundamental insights into the underlying reasons for the observed differences in selectivity but also offer a detailed mechanistic understanding of the catalytic reaction.

Chlordiazepoxide against signalling, receptor and regulatory proteins of breast cancer: a structure-based in-silico approach

Among the most prevalent forms of cancer are breast, lung, colon-rectum, and prostate cancers, and breast cancer is a major global health challenge, contributing to 2.26 million cases with approximately 685,000 deaths worldwide in 2020 alone, typically beginning in the milk ducts or lobules that produce and transport milk during lactation and it is becoming challenging to treat as the tissues are developing resistance, which makes urgent calls for new multitargeted drugs. The multitargeted drug design provides a better solution, simultaneously targeting multiple pathways, even when the drug resists one, it remains effective for others. In this study, we included four crucial proteins that perform signalling, receptor, and regulatory action, namely- NUDIX Hydrolases, Dihydrofolate Reductase, HER2/neu Kinase and EGFR and performed multitargeted molecular docking studies against human-approved drugs using HTVS, SP and extra precise algorithms and filtered the poses with MM\GBSA, suggested a benzodiazepine derivative chlordiazepoxide, used as an anxiolytic agent, can be a multitargeted inhibitor with docking and MM\GBSA score ranging from - 4.628 to - 7.877 and - 18.59 to - 135.86 kcal/mol, respectively, and the most interacted residues were 6ARG, 6GLU, 3TRP, and 3VAL. The QikProp-based ADMET and DFT computations showed the suitability and stability of the drug candidate followed by 100 ns MD simulation in water and MMGBSA on trajectories, resulting in stable performance and many intermolecular interactions to make the complexes stable, which favours that chlordiazepoxide can be a multitargeted breast cancer inhibitor. However, experimental validation is needed before its use.

Rational molecular design converting fascaplysin derivatives to potent broad-spectrum inhibitors against bacterial pathogens via targeting FtsZ

The filamentous temperature-sensitive mutant Z protein (FtsZ), a key player in bacterial cell division machinery, emerges as an attractive target to tackle the plight posed by the ever growing antibiotic resistance over the world. Therefore in this regard, agents with scaffold diversities and broad-spectrum antibacterial activity against Gram-positive and Gram-negative pathogens are highly needed. In this study, a new class of marine-derived fascaplysin derivatives has been designed and synthesized by Suzuki-Miyaura cross-coupling. Some compounds exhibited potent bactericidal activities against a panel of Gram-positive (MIC = 0.024-6.25 μg/mL) and Gram-negative (MIC = 1.56-12.5 μg/mL) bacteria including methicillin-resistant S. aureus (MRSA). They exerted their effects by dual action mechanism via disrupting the integrity of the bacterial cell membrane and targeting FtsZ protein. These compounds stimulated polymerization of FtsZ monomers and bundling of the polymers, and stabilized the resulting polymer network, thus leading to the dysfunction of FtsZ in cell division. In addition, these agents showed negligible hemolytic activity and low cytotoxicity to mammalian cells. The studies on docking and molecular dynamics simulations suggest that these inhibitors bind to the hydrophilic inter-domain cleft of FtsZ protein and the insights obtained in this study would facilitate the development of potential drugs with broad-spectrum bioactivities.

Functional convergence in gastric lysozymes of foregut-fermenting rodents, ruminants, and primates is not attributed to convergent molecular evolution

Convergent evolution is a widespread phenomenon. While there are many examples of convergent evolution at the phenotypic scale, convergence at the molecular level has been more difficult to identify. A classic example of convergent evolution across scales is that of the digestive lysozyme found in ruminants and Colobine monkeys. These herbivorous species rely on foregut fermentation, which has evolved to function more optimally under acidic conditions. Here, we explored if rodents with similar dietary strategies and digestive morphologies have convergently evolved a lysozyme with digestive functions. At the phenotypic level, we find that rodents with bilocular stomach morphologies exhibited a lysozyme that maintained higher relative activities at low pH values, similar to the lysozymes of ruminants and Colobine monkeys. Additionally, the lysozyme of Peromyscus leucopus shared a similar predicted protonation state as that observed in previously identified digestive lysozymes. However, we found limited evidence of positive selection acting on the lysozyme gene in foregut-fermenting species and did not identify patterns of convergent molecular evolution in this gene. This study emphasizes that phenotypic convergence need not be the result of convergent genetic modifications, and we encourage further exploration into the mechanisms regulating convergence across biological scales.

Recent Progress in the Endosomal Escape Mechanism and Chemical Structures of Polycations for Nucleic Acid Delivery

Nucleic acid-based therapies are seeing a spiralling surge. Stimuli-responsive polymers, especially pH-responsive ones, are gaining widespread attention because of their ability to efficiently deliver nucleic acids. These polymers can be synthesized and modified according to target requirements, such as delivery sites and the nature of nucleic acids. In this regard, the endosomal escape mechanism of polymer-nucleic acid complexes (polyplexes) remains a topic of considerable interest owing to various plausible escape mechanisms. This review describes current progress in the endosomal escape mechanism of polyplexes and state-of-the-art chemical designs for pH-responsive polymers. The importance is also discussed of the acid dissociation constant (i.e., pKa) in designing the new generation of pH-responsive polymers, along with assays to monitor and quantify the endosomal escape behavior. Further, the use of machine learning is addressed in pKa prediction and polymer design to find novel chemical structures for pH responsiveness. This review will facilitate the design of new pH-responsive polymers for advanced and efficient nucleic acid delivery.

Titratable residues that drive RND efflux: Insights from molecular simulations

The resistance-nodulation-division efflux machinery confers antimicrobial resistance to Gram-negative bacteria by actively pumping antibiotics out of the cell. The protein complex is powered by proton motive force; however, the proton transfer mechanism itself and indeed even its stoichiometry is still unclear. Here we review computational studies from the last decade that focus on elucidating the number of protons transferred per conformational cycle of the pump. Given the difficulties in studying proton movement using even state-of-the-art structural biology methods, the contributions from computational studies have been invaluable from a mechanistic perspective.

Disrupted protein interaction dynamics in a genetic neurodevelopmental disorder revealed by structural bioinformatics and genetic code expansion

Deciphering the structural effects of gene variants is essential for understanding the pathophysiological mechanisms of genetic diseases. Using a neurodevelopmental disorder called Bosch-Boonstra-Schaaf Optic Atrophy Syndrome (BBSOAS) as a genetic disease model, we applied structural bioinformatics and Genetic Code Expansion (GCE) strategies to assess the pathogenic impact of human NR2F1 variants and their binding with known and novel partners. While the computational analyses of the NR2F1 structure delineated the molecular basis of the impact of several variants on the isolated and complexed structures, the GCE enabled covalent and site-specific capture of transient supramolecular interactions in living cells. This revealed the variable quaternary conformations of NR2F1 variants and highlighted the disrupted interplay with dimeric partners and the newly identified co-factor, CRABP2. The disclosed consequence of the pathogenic mutations on the conformation, supramolecular interplay, and alterations in the cell cycle, viability, and sub-cellular localization of the different variants reflect the heterogeneous disease spectrum of BBSOAS and set up novel foundation for unveiling the complexity of neurodevelopmental diseases.

Identifying novel aryl hydrocarbon receptor (AhR) modulators from clinically approved drugs: In silico screening and In vitro validation

The aryl hydrocarbon receptor (AhR) functions as a vital ligand-activated transcription factor, governing both physiological and pathophysiological processes. Notably, it responds to xenobiotics, leading to a diverse array of outcomes. In the context of drug repurposing, we present here a combined approach of utilizing structure-based virtual screening and molecular dynamics simulations. This approach aims to identify potential AhR modulators from Drugbank repository of clinically approved drugs. By focusing on the AhR PAS-B binding pocket, our screening protocol included binding affinities calculations, complex stability, and interactions within the binding site as a filtering method. Comprehensive evaluations of all DrugBank small molecule database revealed ten promising hits. This included flibanserin, butoconazole, luliconazole, naftifine, triclabendazole, rosiglitazone, empagliflozin, benperidol, nebivolol, and zucapsaicin. Each exhibiting diverse binding behaviors and remarkably very low binding free energy. Experimental studies further illuminated their modulation of AhR signaling, and showing that they are consistently reducing AhR activity, except for luliconazole, which intriguingly enhances the AhR activity. This work demonstrates the possibility of using computational modelling as a quick screening tool to predict new AhR modulators from extensive drug libraries. Importantly, these findings hold immense therapeutic potential for addressing AhR-associated disorders. Consequently, it offers compelling prospects for innovative interventions through drug repurposing.

Structure-based multitargeted docking screening, pharmacokinetics, DFT, and dynamics simulation studies reveal mitoglitazone as a potent inhibitor of cellular survival and stress response proteins of lung cancer

Lung cancer is a disease in which lung cells grow abnormally and uncontrollably, and the cause of it is direct smoking, secondhand smoke, radon, asbestos, and certain chemicals. The worldwide leading cause of death is lung cancer, which is responsible for more than 1.8 million deaths yearly and is expected to rise to 2.2 million by 2030. The most common type of lung cancer is non-small cell lung cancer (NSCLC), which accounts for about 80% and small cell lung cancer (SCLC), which is more aggressive than NSCLC and is often diagnosed later and accounts for 20% of cases. The global concern for lung cancer demands efficient drugs with the slightest chance of developing resistance, and the idea of multitargeted drug designing came up with the solution. In this study, we have performed multitargeted molecular docking studies of Drug Bank compounds with HTVS, SP and XP algorithms followed by MM\GBSA against the four proteins of lung cancer cellular survival and stress responses, which revealed Mitoglitazone as a multitargeted inhibitor with a docking and MM\GBSA score ranging from - 5.784 to - 7.739 kcal/mol and - 25.81 to - 47.65kcal/mol, respectively. Moreover, we performed pharmacokinetics studies and QM-based DFT analysis, showing suitable candidate and interaction pattern analysis revealed the most count of interacting residues was 4GLY, 5PHE, 6ASP, 6GLU, 6LYS, and 6THR. Further, the results were validated with SPC water model-based MD simulation for 100ns in neutralised condition, showing the cumulative deviation and fluctuation < 2Å with many intermolecular interactions. The whole analysis has suggested that Mitoglitazone can be used as a multitargeted inhibitor against lung cancer-however, experimental studies are needed before human use.

Quantitative Assessment of Energetic Contributions of Residues in a SARS-CoV-2 Viral Enzyme/Nanobody Interface

The highly conserved protease enzyme from SARS-CoV-2 (MPro) is crucial for viral replication and is an attractive target for the design of novel inhibitory compounds. MPro is known to be conformationally flexible and has been stabilized in an extended conformation in a complex with a novel nanobody (NB2B4), which inhibits the dimerization of the enzyme via binding to an allosteric site. However, the energetic contributions of the nanobody residues stabilizing the MPro/nanobody interface remain unresolved. We probed these residues using all-atom MD simulations in combination with alchemical free energy calculations by studying the physical residue-residue interactions and discovered the role of hydrophobic and electrostatic interactions in stabilizing the complex. Specifically, we found via mutational analysis that three interfacial nanobody residues (Y59, R106, and L109) contributed significantly, two residues (L107 and P110) contributed moderately, and two residues (H112 and T113) contributed minimally to the overall binding affinity of the nanobody. We also discovered that the nanobody affinity could be enhanced via a charge-reversal mutation (D62R) that alters the local interfacial electrostatic environment of this residue in the complex. These findings are potentially useful in designing novel synthetic nanobodies as allosteric inhibitors of MPro.

Protein3D: Enabling analysis and extraction of metal-containing sites from the Protein Data Bank with molSimplify

Metalloenzymes catalyze a wide range of chemical transformations, with the active site residues playing a key role in modulating chemical reactivity and selectivity. Unlike smaller synthetic catalysts, a metalloenzyme active site is embedded in a larger protein, which makes interrogation of electronic properties and geometric features with quantum mechanical calculations challenging. Here we implement the ability to fetch crystallographic structures from the Protein Data Bank and analyze the metal binding sites in the program molSimplify. We show the usefulness of the newly created protein3D class to extract the local environment around non-heme iron enzymes containing a two histidine motif and prepare 372 structures for quantum mechanical calculations. Our implementation of protein3D serves to expand the range of systems molSimplify can be used to analyze and will enable high-throughput study of metal-containing active sites in proteins.

Cryo-EM of soft-landed β-galactosidase: Gas-phase and native structures are remarkably similar

Native mass spectrometry (MS) has become widely accepted in structural biology, providing information on stoichiometry, interactions, homogeneity, and shape of protein complexes. Yet, the fundamental assumption that proteins inside the mass spectrometer retain a structure faithful to native proteins in solution remains a matter of intense debate. Here, we reveal the gas-phase structure of β-galactosidase using single-particle cryo-electron microscopy (cryo-EM) down to 2.6-Å resolution, enabled by soft landing of mass-selected protein complexes onto cold transmission electron microscopy (TEM) grids followed by in situ ice coating. We find that large parts of the secondary and tertiary structure are retained from the solution. Dehydration-driven subunit reorientation leads to consistent compaction in the gas phase. By providing a direct link between high-resolution imaging and the capability to handle and select protein complexes that behave problematically in conventional sample preparation, the approach has the potential to expand the scope of both native mass spectrometry and cryo-EM.

Identification of Residues Potentially Involved in Optical Shifts in the Water-Soluble Chlorophyll a-Binding Protein through Molecular Dynamics Simulations

Reversible light and thermally induced spectral shifts are universally observed in a wide variety of pigment-protein complexes at temperatures ranging from cryogenic to ambient. In this paper, we employed large-scale molecular dynamics (MD) simulations of a prototypical pigment-protein complex to better understand these shifts at a molecular scale. Although multiple mechanisms have been proposed over the years, no verification of these proposals via MD simulations has thus far been performed; our work represents the first step in this direction. From simulations of the water-soluble chlorophyll-binding protein complex, we determined that rearrangements of long hydrogen bonds were unlikely to be the origin of the multiwell landscape features necessary to explain observed spectral shifts. We also assessed small motions of amino acid residues and identified side chain rotations of some of these residues as likely candidates for the origin of relevant multiwell landscape features. The protein free-energy landscapes associated with side chain rotations feature energy barriers of around 1100-1600 cm-1, in agreement with optical spectroscopy results, with the most promising residue type associated with experimental signatures being serine, which possesses a symmetric triple-well landscape and moment of inertia of a relevant magnitude.

Multiscale In Silico Study of the Mechanism of Activation of the RtcB Ligase by the PTP1B Phosphatase

Inositol-requiring enzyme 1 (IRE1) is a transmembrane sensor that is part of a trio of sensors responsible for controlling the unfolded protein response within the endoplasmic reticulum (ER). Upon the accumulation of unfolded or misfolded proteins in the ER, IRE1 becomes activated and initiates the cleavage of a 26-nucleotide intron from human X-box-containing protein 1 (XBP1). The cleavage is mediated by the RtcB ligase enzyme, which splices together two exons, resulting in the formation of the spliced isoform XBP1s. The XBP1s isoform activates the transcription of genes involved in ER-associated degradation to maintain cellular homeostasis. The catalytic activity of RtcB is controlled by the phosphorylation and dephosphorylation of three tyrosine residues (Y306, Y316, and Y475), which are regulated by the ABL1 tyrosine kinase and PTP1B phosphatase, respectively. This study focuses on investigating the mechanism by which the PTP1B phosphatase activates the RtcB ligase using a range of advanced in silico methods. Protein-protein docking identified key interacting residues between RtcB and PTP1B. Notably, the phosphorylated Tyr306 formed hydrogen bonds and salt bridge interactions with the "gatekeeper" residues Arg47 and Lys120 of the inactive PTP1B. Classical molecular dynamics simulation emphasized the crucial role of Asp181 in the activation of PTP1B, driving the conformational change from an open to a closed state of the WPD-loop. Furthermore, QM/MM-MD simulations provided insights into the free energy landscape of the dephosphorylation reaction mechanism of RtcB, which is mediated by the PTP1B phosphatase.

Three-Dimensional Structural Stability and Local Electrostatic Potential at Point Mutations in Spike Protein of SARS-CoV-2 Coronavirus

The contagiousness of SARS-CoV-2 β-coronavirus is determined by the virus-receptor electrostatic association of its positively charged spike (S) protein with the negatively charged angiotensin converting enzyme-2 (ACE2 receptor) of the epithelial cells. If some mutations occur, the electrostatic potential on the surface of the receptor-binding domain (RBD) could be altered, and the S-ACE2 association could become stronger or weaker. The aim of the current research is to investigate whether point mutations can noticeably alter the electrostatic potential on the RBD and the 3D stability of the S1-subunit of the S-protein. For this purpose, 15 mutants with different hydrophilicity and electric charge (positive, negative, or uncharged) of the substituted and substituting amino acid residues, located on the RBD at the S1-ACE2 interface, are selected, and the 3D structure of the S1-subunit is reconstructed on the base of the crystallographic structure of the S-protein of the wild-type strain and the amino acid sequence of the unfolded polypeptide chain of the mutants. Then, the Gibbs free energy of folding, isoelectric point, and pH-dependent surface electrostatic potential of the S1-subunit are computed using programs for protein electrostatics. The results show alterations in the local electrostatic potential in the vicinity of the mutant amino acid residue, which can influence the S-ACE2 association. This approach allows prediction of the relative infectivity, transmissibility, and contagiousness (at equal social immune status) of new SARS-CoV-2 mutants by reconstruction of the 3D structure of the S1-subunit and calculation of the surface electrostatic potential.

Structure-Based In Silico Approaches Reveal IRESSA as a Multitargeted Breast Cancer Regulatory, Signalling, and Receptor Protein Inhibitor

Breast cancer begins in the breast cells, mainly impacting women. It starts in the cells that line the milk ducts or lobules responsible for producing milk and can spread to nearby tissues and other body parts. In 2020, around 2.3 million women across the globe received a diagnosis, with an estimated 685,000 deaths. Additionally, 7.8 million women were living with breast cancer, making it the fifth leading cause of cancer-related deaths among women. The mutational changes, overexpression of drug efflux pumps, activation of alternative signalling pathways, tumour microenvironment, and cancer stem cells are causing higher levels of drug resistance, and one of the major solutions is to identify multitargeted drugs. In our research, we conducted a comprehensive screening using HTVS, SP, and XP, followed by an MM/GBSA computation of human-approved drugs targeting HER2/neu, BRCA1, PIK3CA, and ESR1. Our analysis pinpointed IRESSA (Gefitinib-DB00317) as a multitargeted inhibitor for these proteins, revealing docking scores ranging from -4.527 to -8.809 Kcal/mol and MM/GBSA scores between -49.09 and -61.74 Kcal/mol. We selected interacting residues as fingerprints, pinpointing 8LEU, 6VAL, 6LYS, 6ASN, 5ILE, and 5GLU as the most prevalent in interactions. Subsequently, we analysed the ADMET properties and compared them with the standard values of QikProp. We extended our study for DFT computations with Jaguar and plotted the electrostatic potential, HOMO and LUMO regions, and electron density, followed by a molecular dynamics simulation for 100 ns in water, showing an utterly stable performance, making it a suitable drug candidate. IRESSA is FDA-approved for lung cancer, which shares some pathways with breast cancers, clearing the hurdles of multitargeted drugs against breast and lung cancer. This has the potential to be groundbreaking; however, more studies are needed to concreate IRESSA's role.

Integration of persistent Laplacian and pre-trained transformer for protein solubility changes upon mutation

Protein mutations can significantly influence protein solubility, which results in altered protein functions and leads to various diseases. Despite tremendous effort, machine learning prediction of protein solubility changes upon mutation remains a challenging task as indicated by the poor scores of normalized Correct Prediction Ratio (CPR). Part of the challenge stems from the fact that there is no three-dimensional (3D) structures for the wild-type and mutant proteins. This work integrates persistent Laplacians and pre-trained Transformer for the task. The Transformer, pretrained with hundreds of millions of protein sequences, embeds wild-type and mutant sequences, while persistent Laplacians track the topological invariant change and homotopic shape evolution induced by mutations in 3D protein structures, which are rendered from AlphaFold2. The resulting machine learning model was trained on an extensive data set labeled with three solubility types. Our model outperforms all existing predictive methods and improves the state-of-the-art up to 15%.

In-silico guided design, screening, and molecular dynamic simulation studies for the identification of potential SARS-CoV-2 main protease inhibitors for the targeted treatment of COVID-19

COVID-19, the disease responsible for the recent pandemic, is caused by a novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The main protease (Mpro) of SARS-CoV-2 is an essential proteolytic enzyme that plays a number of important roles in the replication of the virus in human host cells. Blocking the function of SARS-CoV-2 Mpro offers a promising and targeted, therapeutic option for the treatment of the COVID-19 infection. Such an inhibitory strategy is currently successful in treating COVID-19 under FDA's emergency use authorization, although with limited benefit to the immunocompromised along with an unfortunate number of side effects and drug-drug interactions. Current COVID vaccines protect against severe disease and death but are mostly ineffective toward long COVID which has been seen in 5-36% of patients. SARS-CoV-2 is a rapidly mutating virus and is here to stay endemically. Hence, alternate therapeutics to treat SARS-CoV-2 infections are still needed. Moreover, because of the high degree of conservation of Mpro among different coronaviruses, any newly developed antiviral agents should better prepare us for potential future epidemics or pandemics. In this paper, we first describe the design and computational docking of a library of novel 188 first-generation peptidomimetic protease inhibitors using various electrophilic warheads with aza-peptide epoxides, α-ketoesters, and β-diketones identified as the most effective. Second-generation designs, 192 compounds in total, focused on aza-peptide epoxides with drug-like properties, incorporating dipeptidyl backbones and heterocyclic ring motifs such as proline, indole, and pyrrole groups, yielding 8 hit candidates. These novel and specific inhibitors for SARS-CoV-2 Mpro can ultimately serve as valuable alternate and broad-spectrum antivirals against COVID-19.Communicated by Ramaswamy H. Sarma.

Structural characterization and functional insights into the type II secretion system of the poly-extremophile Deinococcus radiodurans

The extremophile bacterium D. radiodurans boasts a distinctive cell envelope characterized by the regular arrangement of three protein complexes. Among these, the Type II Secretion System (T2SS) stands out as a pivotal structural component. We used cryo-electron microscopy to reveal unique features, such as an unconventional protein belt (DR_1364) around the main secretin (GspD), and a cap (DR_0940) found to be a separated subunit rather than integrated with GspD. Furthermore, a novel region at the N-terminus of the GspD constitutes an additional second gate, supplementing the one typically found in the outer membrane region. This T2SS was found to contribute to envelope integrity, while also playing a role in nucleic acid and nutrient trafficking. Studies on intact cell envelopes show a consistent T2SS structure repetition, highlighting its significance within the cellular framework.

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.

AMBERff at Scale: Multimillion-Atom Simulations with AMBER Force Fields in NAMD

All-atom molecular dynamics (MD) simulations are an essential structural biology technique with increasing application to multimillion-atom systems, including viruses and cellular machinery. Classical MD simulations rely on parameter sets, such as the AMBER family of force fields (AMBERff), to accurately describe molecular motion. Here, we present an implementation of AMBERff for use in NAMD that overcomes previous limitations to enable high-performance, massively parallel simulations encompassing up to two billion atoms. Single-point potential energy comparisons and case studies on model systems demonstrate that the implementation produces results that are as accurate as running AMBERff in its native engine.

Molecular basis of the reaction mechanism of the methyltransferase HENMT1

PIWI-interacting RNAs (piRNAs) are important for ensuring the integrity of the germline. 3'-terminal 2'-O-methylation is essential for piRNA maturation and to protect them from degradation. HENMT1 (HEN Methyltransferase 1) carries out the 2'-O-methylation, which is of key importance for piRNA stability and functionality. However, neither the structure nor the catalytic mechanism of mammalian HENMT1 have been studied. We have constructed a catalytic-competent HENMT1 complex using computational approaches, in which Mg2+ is primarily coordinated by four evolutionary conserved residues, and is further auxiliary coordinated by the 3'-O and 2'-O on the 3'-terminal nucleotide of the piRNA. Our study suggests that metal has limited effects on substrate and cofactor binding but is essential for catalysis. The reaction consists of deprotonation of the 2'-OH to 2'-O and a methyl transfer from SAM to the 2'-O. The methyl transfer is spontaneous and fast. Our in-depth analysis suggests that the 2'-OH may be deprotonated before entering the active site or it may be partially deprotonated at the active site by His800 and Asp859, which are in a special alignment that facilitates the proton transfer out of the active site. Furthermore, we have developed a detailed potential reaction scenario indicating that HENMT1 is Mg2+ utilizing but is not a Mg2+ dependent enzyme.

Effects of LDLR variants rs5928, rs750518671 and rs879254797 on protein structure and functional activity in HepG2 cells transfected with CRISPR/Cas9 constructs

Familial Hypercholesterolemia (FH) is a genetic disorder associated with premature atherosclerosis and increased risk of cardiovascular diseases. LDLR deleterious mutations are associated with FH, however the role of some missense variants in FH pathogenicity remains to be elucidated. This study explored the predictive impact of LDLR missense variants on protein structure and investigated their functional effects on LDLR expression in HepG2 cells transfected with CRISPR/Cas9 constructs. FH (n = 287) and non-FH patients (n = 45) were selected, and lipid profile was obtained from medical records. LDLR variants were identified using an exon-targeted gene sequencing strategy, considering its cost-effective to increase accuracy in the identification step of the most likely FH-related variants in a less laborious process. LDLR variants were selected based on conflicting pathogenicity results found in Clinvar, in silico prediction tools, affected LDLR domains, and less common variants considering minor allele frequency < 0.05. Molecular modeling studies were used to predict the effects of LDLR missense variants on protein structure. Recombinant LDLR variants were constructed using CRISPR/Cas9 system and were used to transfect HepG2 cells. Functional assays in transfected cells were performed to assess LDLR expression using flow cytometry and western blotting, and LDLR activity using flow cytometry and confocal microscopy. The variants rs121908039 (c.551G>A, p.C184Y), rs879254797 (c.1118G>A, p.G373D), rs28941776 (c.1646G>A, p.G549D), rs750518671 (c.2389G>C, p.V797L), rs5928 (c.2441G>A, p.R814Q) and rs137853964 (c.2479G>A, p.V827I) were selected for molecular docking analysis. The p.C184Y exhibited a favorable energy change for protein stability due to its interaction with EGF-A/EGF-B regions; p.G373D and p.G549D displayed intermediate energy changes; and p.R814Q and p.V827I showed smaller energy changes. The results of functional assays showed that p.G373D, p.V797L and p.R814Q reduced LDLR expression and activity (p < 0.05). Microscopic analysis of the p.V797L and p.G373D variants revealed altered lipid localization and accumulation in transfected HepG2 cells. Carriers of p.G549D, p.V797L and p.R814Q had higher LDL cholesterol levels than non-FH group, and (p < 0.05). p.G373D and p.G549D were associated with clinical manifestations of FH. In conclusion, the p.C184Y, p.G373D, p.G549D and p.R814Q variants alter protein stability and intramolecular interactions, while p.V797L has a minimal impact on protein stability, and p.V827I has no significant intramolecular interactions. p.G373D, p.V767L and p.R814Q are associated with impaired LDLR expression and activity.