A simple clustering algorithm can be accurate enough for use in calculations of pKs in macromolecules

Computer-aided engineering of stabilized fibroblast growth factor 21

FGF21 is an endocrine signaling protein belonging to the family of fibroblast growth factors (FGFs). It has emerged as a molecule of interest for treating various metabolic diseases due to its role in regulating glucogenesis and ketogenesis in the liver. However, FGF21 is prone to heat, proteolytic, and acid-mediated degradation, and its low molecular weight makes it susceptible to kidney clearance, significantly reducing its therapeutic potential. Protein engineering studies addressing these challenges have generally shown that increasing the thermostability of FGF21 led to improved pharmacokinetics. Here, we describe the computer-aided design and experimental characterization of FGF21 variants with enhanced melting temperature up to 15 °C, uncompromised efficacy at activation of MAPK/ERK signaling in Hep G2 cell culture, and ability to stimulate proliferation of Hep G2 and NIH 3T3 fibroblasts cells comparable with FGF21-WT. We propose that stabilizing the FGF21 molecule by rational design should be combined with other reported stabilization strategies to maximize the pharmaceutical potential of FGF21.

Exploring the selectivity of cytochrome P450 for enhanced novel anticancer agent synthesis

Cytochrome P450 monooxygenases are an extensive and unique class of enzymes, which can regio- and stereo-selectively functionalise hydrocarbons by way of oxidation reactions. These enzymes are naturally occurring but have also been extensively applied in a synthesis context, where they are used as efficient biocatalysts. Recently, a biosynthetic pathway where a cytochrome P450 monooxygenase catalyses a critical step of the pathway was uncovered, leading to the production of a number of products that display high antitumour potency. In this work, we use computational techniques to gain insight into the factors that determine the relative yields of the different products. We use conformational search algorithms to understand the substrate stereochemistry. On a machine-learned 3D protein structure, we use molecular docking to obtain a library of favourable poses for substrate-protein interaction. With molecular dynamics, we investigate the most favourable poses for reactivity on a molecular level, allowing us to investigate which protein-substrate interactions favour a given product and thus gain insight into the product selectivity.

Mechanism of 7H-Dibenzo[c,g]carbazole metabolism in cytochrome P450 1A1: Insights from computational studies

7H-Dibenzo[c,g]carbazole (DBC) is a prevalent environmental contaminant that induces tumorigenesis in several experimental animals. Recently, it has been utilized to develop high-performance solar cells and organic phosphorescent materials. It is imperative to strengthen investigations of DBC metabolism to understand its potential risks to human health. In this study, human CYP1A1 was employed as the metabolic enzyme to investigate the metabolic mechanism of DBC by molecular docking, molecular dynamics (MD) simulation, and quantum mechanical (QM) calculation. The results indicate that DBC binds to CYP1A1 in two modes (mode 1 and mode 2) mainly through nonpolar solvation energies (ΔGnonpolar). The formation of the two binding modes is attributed to the anchoring effect of the hydrogen bond formed by DBC with Asp320 (mode 1) or Ser116 (mode 2). Mode 1 is a "reactive" conformation, while mode 2 is not considered a "reactive" conformation. C5 is identified as the dominant site, and the pyrrole nitrogen cannot participate in the metabolism. DBC is metabolized mainly by a distinct electrophilic addition-rearrangement mechanism, with an energy barrier of 21.74 kcal/mol. The results provide meaningful insights into the biometabolic process of DBC and contribute to understanding its environmental effects and health risks.

Proton-Coupled Electron Transfer upon Oxidation of Tyrosine in a De Novo Protein: Analysis of Proton Acceptor Candidates

Redox-active residues, such as tyrosine and tryptophan, play important roles in a wide range of biological processes. The α3Y de novo protein, which is composed of three α helices and a tyrosine residue Y32, provides a platform for investigating the redox properties of tyrosine in a well-defined protein environment. Herein, the proton-coupled electron transfer (PCET) reaction that occurs upon oxidation of tyrosine in this model protein by a ruthenium photosensitizer is studied by using a vibronically nonadiabatic PCET theory that includes hydrogen tunneling and excited vibronic states. The input quantities to the analytical nonadiabatic rate constant expression, such as the diabatic proton potential energy curves and associated proton vibrational wave functions, reorganization energy, and proton donor-acceptor distribution functions, are obtained from density functional theory calculations on model systems and molecular dynamics simulations of the solvated α3Y protein. Two possible proton acceptors, namely, water or a glutamate residue in the protein scaffold, are explored. The PCET rate constant is greater when glutamate is the proton acceptor, mainly due to the more favorable driving force and shorter equilibrium proton donor-acceptor distance, although contributions from excited vibronic states mitigate these effects. Nevertheless, water could be the dominant proton acceptor if its equilibrium constant associated with hydrogen bond formation is significantly greater than that for glutamate. Although these calculations do not definitively identify the proton acceptor for this PCET reaction, they elucidate the conditions under which each proton acceptor can be favored. These insights have implications for tyrosine-based PCET in a wide variety of biochemical processes.

Structures of BlEst2 from Bacillus licheniformis in its propeptide and mature forms reveal autoinhibitory effects of the C-terminal domain

Carboxylesterases comprise a major class of α/β-fold hydrolases responsible for the cleavage and formation of ester bonds. Found ubiquitously in nature, these enzymes are crucial for the metabolism of both endogenous and exogenous carboxyl esters in animals, plants and microorganisms. Beyond their essential physiological roles, carboxylesterases stand out as one of the important classes of biocatalysts for biotechnology. BlEst2, an enzyme previously classified as Bacillus licheniformis esterase, remains largely uncharacterized. In the present study, we elucidate the structural biology, molecular dynamics and biochemical features of BlEst2. Our findings reveal a canonical α/β-hydrolase fold similar to the ESTHER block L of lipases, further augmented by two additional accessory C-terminal domains. Notably, the catalytic domain demonstrates two insertions, which occupy conserved locations in α/β-hydrolase proteins and commonly form the lid domain in lipase structures. Intriguingly, our in vitro cleavage of C-terminal domains revealed the structure of the active form of BlEst2. Upon activation, BlEst2 showed a markedly elevated hydrolytic activity. This observation implies that the intramolecular C-terminal domain serves as a regulatory intramolecular inhibitor. Interestingly, despite exhibiting esterase-like activity, BlEst2 structural characteristics align more closely with lipases. This suggests that BlEst2 could potentially represent a previously unrecognized subgroup within the realm of carboxyl ester hydrolases.

Lytic transglycosylase Slt of Pseudomonas aeruginosa as a periplasmic hub protein

Peptidoglycan is a major constituent of the bacterial cell wall. Its integrity as a polymeric edifice is critical for bacterial survival and, as such, it is a preeminent target for antibiotics. The peptidoglycan is a dynamic crosslinked polymer that undergoes constant biosynthesis and turnover. The soluble lytic transglycosylase (Slt) of Pseudomonas aeruginosa is a periplasmic enzyme involved in this dynamic turnover. Using amber-codon-suppression methodology in live bacteria, we incorporated a fluorescent chromophore into the structure of Slt. Fluorescent microscopy shows that Slt populates the length of the periplasmic space and concentrates at the sites of septation in daughter cells. This concentration persists after separation of the cells. Amber-codon-suppression methodology was also used to incorporate a photoaffinity amino acid for the capture of partner proteins. Mass-spectrometry-based proteomics identified 12 partners for Slt in vivo. These proteomics experiments were complemented with in vitro pulldown analyses. Twenty additional partners were identified. We cloned the genes and purified to homogeneity 22 identified partners. Biophysical characterization confirmed all as bona fide Slt binders. The identities of the protein partners of Slt span disparate periplasmic protein families, inclusive of several proteins known to be present in the divisome. Notable periplasmic partners (KD < 0.5 μM) include PBPs (PBP1a, KD = 0.07 μM; PBP5 = 0.4 μM); other lytic transglycosylases (SltB2, KD = 0.09 μM; RlpA, KD = 0.4 μM); a type VI secretion system effector (Tse5, KD = 0.3 μM); and a regulatory protease for alginate biosynthesis (AlgO, KD < 0.4 μM). In light of the functional breadth of its interactome, Slt is conceptualized as a hub protein within the periplasm.

De novo design of proteins housing excitonically coupled chlorophyll special pairs

Natural photosystems couple light harvesting to charge separation using a 'special pair' of chlorophyll molecules that accepts excitation energy from the antenna and initiates an electron-transfer cascade. To investigate the photophysics of special pairs independently of the complexities of native photosynthetic proteins, and as a first step toward creating synthetic photosystems for new energy conversion technologies, we designed C2-symmetric proteins that hold two chlorophyll molecules in closely juxtaposed arrangements. X-ray crystallography confirmed that one designed protein binds two chlorophylls in the same orientation as native special pairs, whereas a second designed protein positions them in a previously unseen geometry. Spectroscopy revealed that the chlorophylls are excitonically coupled, and fluorescence lifetime imaging demonstrated energy transfer. The cryo-electron microscopy structure of a designed 24-chlorophyll octahedral nanocage with a special pair on each edge closely matched the design model. The results suggest that the de novo design of artificial photosynthetic systems is within reach of current computational methods.

Nanocellulose-Bovine Serum Albumin Interactions in an Aqueous Medium: Investigations Using In Situ Nanocolloidal Probe Microscopy and Reactive Molecular Dynamics Simulations

As a versatile nanomaterial derived from renewable sources, nanocellulose has attracted considerable attention for its potential applications in various sectors, especially those focused on water treatment and remediation. Here, we have combined atomic force microscopy (AFM) and reactive molecular dynamics (RMD) simulations to characterize the interactions between cellulose nanofibers modified with carboxylate or phosphate groups and the protein foulant model bovine serum albumin (BSA) at pH 3.92, which is close to the isoelectric point of BSA. Colloidal probes were prepared by modification of the AFM probes with the nanofibers, and the nanofiber coating on the AFM tip was for the first time confirmed through fluorescence labeling and confocal optical sectioning. We have found that the wet-state normalized adhesion force is approximately 17.87 ± 8.58 pN/nm for the carboxylated cellulose nanofibers (TOCNF) and about 11.70 ± 2.97 pN/nm for the phosphorylated ones (PCNF) at the studied pH. Moreover, the adsorbed protein partially unfolded at the cellulose interface due to the secondary structure's loss of intramolecular hydrogen bonds. We demonstrate that nanocellulose colloidal probes can be used as a sensitive tool to reveal interactions with BSA at nano and molecular scales and under in situ conditions. RMD simulations helped to gain a molecular- and atomistic-level understanding of the differences between these findings. In the case of PCNF, partially solvated metal ions, preferentially bound to the phosphates, reduced the direct protein-cellulose connections. This understanding can lead to significant advancements in the development of cellulose-based antifouling surfaces and provide crucial insights for expanding the pH range of use and suggesting appropriate recalibrations.

Atomic Insight into the Enzymatic Selectivity of Acetyltransferase for Endogenous Polyamines

The N1-Spermidine/spermine acetyltransferase (SSAT) serves as the rate-limiting enzyme in the polyamine metabolism pathway, specifically catalyzing the acetylation of spermidine, spermine, and other specific polyamines. The source of its enzymatic selectivity remains elusive. Here, we used quantum mechanics and molecular mechanics simulations combined with various technologies to explore the enzymatic mechanism of SSAT for endogenous polyamines from an atomic perspective. The static binding and chemical transformation were considered. The binding affinity was identified to be dependent on protonated state of polyamine. The order of the binding affinity for Spm, Spd, and Put is consistent with the experimental results, which is also verified by the dynamic separation of polyamine and SSAT. Hydrogen bond interactions and salt bridges contribute most, and the common hot residues were identified. In addition, the transfer of acetyl and proton between polyamine and AcCoA was discovered to follow a concert mechanism, and thermodynamic properties are responsible for the catalytic efficiency of SSAT. This work may be helpful for development of polyamine derivatives based on catalysis to regulate polyamine metabolism.

Allosteric Modulation of the YAP/TAZ-TEAD Interaction by Palmitoylation and Small-Molecule Inhibitors

The Hippo signaling pathway is a highly conserved signaling network that plays a central role in regulating cellular growth, proliferation, and organ size. This pathway consists of a kinase cascade that integrates various upstream signals to control the activation or inactivation of YAP/TAZ proteins. Phosphorylated YAP/TAZ is sequestered in the cytoplasm; however, when the Hippo pathway is deactivated, it translocates into the nucleus, where it associates with TEAD transcription factors. This partnership is instrumental in regulating the transcription of progrowth and antiapoptotic genes. Thus, in many cancers, aberrantly hyperactivated YAP/TAZ promotes oncogenesis by contributing to cancer cell proliferation, metastasis, and therapy resistance. Because YAP and TAZ exert their oncogenic effects by binding with TEAD, it is critical to understand this key interaction to develop cancer therapeutics. Previous research has indicated that TEAD undergoes autopalmitoylation at a conserved cysteine, and small molecules that inhibit TEAD palmitoylation disrupt effective YAP/TAZ binding. However, how exactly palmitoylation contributes to YAP/TAZ-TEAD interactions and how the TEAD palmitoylation inhibitors disrupt this interaction remains unknown. Utilizing molecular dynamics simulations, our investigation not only provides detailed atomistic insight into the YAP/TAZ-TEAD dynamics but also unveils that the inhibitor studied influences the binding of YAP and TAZ to TEAD in distinct manners. This discovery has significant implications for the design and deployment of future molecular interventions targeting this interaction.

Role of UDP-N-acetylmuramic acid in the regulation of MurA activity revealed by molecular dynamics simulations

The peptidoglycan biosynthesis pathway plays a vital role in bacterial cells, and facilitates peptidoglycan layer formation, a fundamental structural component of the bacterial cell wall. The enzymes in this pathway are candidates for antibiotic development, as most do not have mammalian homologues. The UDP-N-acetylglucosamine (UNAG) enolpyruvyl transferase enzyme (MurA) in the peptidoglycan pathway cytoplasmic step is responsible for the phosphoenolpyruvate (PEP)-UNAG catalytic reaction, forming UNAG enolpyruvate and inorganic phosphate. Reportedly, UDP-N-acetylmuramic acid (UNAM) binds tightly to MurA forming a dormant UNAM-PEP-MurA complex and acting as a MurA feedback inhibitor. MurA inhibitors are complex, owing to competitive binding interactions with PEP, UNAM, and UNAG at the MurA active site. We used computational methods to explore UNAM and UNAG binding. UNAM showed stronger hydrogen-bond interactions with the Arg120 and Arg91 residues, which help to stabilize the closed conformation of MurA, than UNAG. Binding free energy calculations using end-point computational methods showed that UNAM has a higher binding affinity than UNAG, when PEP is attached to Cys115. The unbinding process, simulated using τ-random acceleration molecular dynamics, showed that UNAM has a longer relative residence time than UNAG, which is related to several complex dissociation pathways, each with multiple intermediate metastable states. This prevents the loop from opening and exposing the Arg120 residue to accommodate UNAG and potential new ligands. Moreover, we demonstrate the importance of Cys115-linked PEP in closed-state loop stabilization. We provide a basis for evaluating novel UNAM analogues as potential MurA inhibitors. PUBLIC SIGNIFICANCE: MurA is a critical enzyme involved in bacterial cell wall biosynthesis and is involved in antibiotic resistance development. UNAM can remain in the target protein's active site for an extended time compared to its natural substrate, UNAG. The prolonged interaction of this highly stable complex known as the 'dormant complex' comprises UNAM-PEP-MurA and offers insights into antibiotic development, providing potential options against drug-resistant bacteria and advancing our understanding of microbial biology.

Myrtucommulones and Related Acylphloroglucinols from Myrtaceae as a Promising Source of Multitarget SARS-CoV-2 Cycle Inhibitors

The LABEXTRACT plant extract bank, featuring diverse members of the Myrtaceae family from Brazilian hot spot regions, provides a promising avenue for bioprospection. Given the pivotal roles of the Spike protein and 3CLpro and PLpro proteases in SARS-CoV-2 infection, this study delves into the correlations between the Myrtaceae species from the Atlantic Forest and these targets, as well as an antiviral activity through both in vitro and in silico analyses. The results uncovered notable inhibitory effects, with Eugenia prasina and E. mosenii standing out, while E. mosenii proved to be multitarget, presenting inhibition values above 72% in the three targets analyzed. All extracts inhibited viral replication in Calu-3 cells (EC50 was lower than 8.3 µg·mL-1). Chemometric analyses, through LC-MS/MS, encompassing prediction models and molecular networking, identified potential active compounds, such as myrtucommulones, described in the literature for their antiviral activity. Docking analyses showed that one undescribed myrtucommulone (m/z 841 [M - H]-) had a higher fitness score when interacting with the targets of this study, including ACE2, Spike, PLpro and 3CLpro of SARS-CoV-2. Also, the study concludes that Myrtaceae extracts, particularly from E. mosenii and E. prasina, exhibit promising inhibitory effects against crucial stages in SARS-CoV-2 infection. Compounds like myrtucommulones emerge as potential anti-SARS-CoV-2 agents, warranting further exploration.

MD investigation on the differences in the dynamic interactions between the specific ligand azamulin and two CYP3A isoforms, 3A4 and 3A5

The unmarked potential drug molecule azamulin has been found to be a specific inhibitor of CYP3A4 and CYP3A5 in recent years, but this molecule also shows different binding ability and affinity to the two CYP3A isoforms. In order to explore the microscopic mechanism, conventional molecular dynamics (MD) simulation methods were performed to study the dynamic interactions between two isoforms and azamulin. The simulation results show that the binding of the ligand leads to different structural properties of two CYP3A proteins. First of all, compared with apo-CYP3A4, the binding of the ligand azamulin can lead to changes in the structural flexibility of CYP3A4, i.e., holo-CYP3A4 is more flexible than apo-CYP3A4. The structural changes of CYP3A5 are just the opposite. The ligand binding increases the rigidity of CYP3A5. Furthermore, the representative structures of the production phase in the MD simulation were in details analyzed to obtain the microscopic interactions between the ligand azamulin and two CYP3A isoforms at the atomic level. It is speculated that the difference of composition and interaction of the active sites is the fundamental cause of the change of structural properties of the two proteins.Communicated by Ramaswamy H. Sarma.

Learnt representations of proteins can be used for accurate prediction of small molecule binding sites on experimentally determined and predicted protein structures

Protein-ligand binding site prediction is a useful tool for understanding the functional behaviour and potential drug-target interactions of a novel protein of interest. However, most binding site prediction methods are tested by providing crystallised ligand-bound (holo) structures as input. This testing regime is insufficient to understand the performance on novel protein targets where experimental structures are not available. An alternative option is to provide computationally predicted protein structures, but this is not commonly tested. However, due to the training data used, computationally-predicted protein structures tend to be extremely accurate, and are often biased toward a holo conformation. In this study we describe and benchmark IF-SitePred, a protein-ligand binding site prediction method which is based on the labelling of ESM-IF1 protein language model embeddings combined with point cloud annotation and clustering. We show that not only is IF-SitePred competitive with state-of-the-art methods when predicting binding sites on experimental structures, but it performs better on proxies for novel proteins where low accuracy has been simulated by molecular dynamics. Finally, IF-SitePred outperforms other methods if ensembles of predicted protein structures are generated.

Revealing the binding mechanism of BACE1 inhibitors through molecular dynamics simulations

Alzheimer's disease is a debilitating neurodegenerative disorder, and the Beta-Site Amyloid Precursor Protein Cleaving Enzyme 1 (BACE1) is a key therapeutic target in its treatment. This study employs molecular dynamics simulations and binding energy analysis to investigate the binding interactions between BACE1 and four selected small molecules: CNP520, D9W, NB641, and NB360. The binding model analysis indicates that the binding of BACE1 with four molecules are stable, except the loop regions show significant fluctuation. The binding free energy analyses reveal that NB360 exhibits the highest binding affinity with BACE1, surpassing other molecules (CNP520, D9W, and NB641). Detailed energy component assessments highlight the critical roles of electrostatic interactions and van der Waals forces in the binding process. Furthermore, residue contribution analysis identifies key amino acids influencing the binding process across all systems. Hydrogen bond analysis reveals a limited number of bonds between BACE1 and each small molecule, highlighting the importance of structural modifications to enable more stable hydrogen bonds. This research provides valuable insights into the molecular mechanisms of potential Alzheimer's disease therapeutics, guiding the way for improved drug design and the development of effective treatments targeting BACE1.Communicated by Ramaswamy H. Sarma.

Computational approaches to investigate fluoride binding, selectivity and transport across the membrane

The use of molecular dynamics (MD) simulations to study biomolecular systems has proven reliable in elucidating atomic-level details of structure and function. In this chapter, MD simulations were used to uncover new insights into two phylogenetically unrelated bacterial fluoride (F-) exporters: the CLCF F-/H+ antiporter and the Fluc F- channel. The CLCF antiporter, a member of the broader CLC family, has previously revealed unique stoichiometry, anion-coordinating residues, and the absence of an internal glutamate crucial for proton import in the CLCs. Through MD simulations enhanced with umbrella sampling, we provide insights into the energetics and mechanism of the CLCF transport process, including its selectivity for F- over HF. In contrast, the Fluc F- channel presents a novel architecture as a dual topology dimer, featuring two pores for F- export and a central non-transported sodium ion. Using computational electrophysiology, we simulate the electrochemical gradient necessary for F- export in Fluc and reveal details about the coordination and hydration of both F- and the central sodium ion. The procedures described here delineate the specifics of these advanced techniques and can also be adapted to investigate other membrane protein systems.

Suitability of Adenosine Derivatives in Improving the Activity and Stability of Cytochrome c under Stress: Insights into the Effect of Phosphate Groups

It is well known that adenosine and its phosphate derivatives play a crucial role in biological phenomena such as apoptosis and cell signaling and act as the energy currency of the cell. Although their interactions with various proteins and enzymes have been described, the focus of this work is to demonstrate the effect of the phosphate group on the activity and stability of the native heme metalloprotein cytochrome c (Cyt c), which is important from both biological and industrial aspects. In situ and in silico characterizations are used to correlate the relationship between the binding affinity of adenosine and its phosphate groups with unfolding behavior, corresponding peroxidase activities, and stability factors. Interaction of adenosine (ADN), adenosine monophosphate (AMP), adenosine 5'-diphosphate (ADP), and adenosine 5'-triphosphate (ATP) with Cyt c increases peroxidase-like activity by up to 1.8-6.5-fold compared to native Cyt c. This activity is significantly maintained even after multiple stress conditions such as oxidative stress and the presence of a chaotropic agent such as guanidine hydrochloride (GuHCl). With binding affinities on the order of ADN < AMP < ADP < ATP, adenosine derivatives were found to stabilize Cyt c by varying the secondary structural features of the protein. Thus, in addition to being a fundamental study, the current work also proposes a way of stabilizing protein systems to be used for real-time biocatalytic applications.

Ligand 1 H NMR Chemical Shifts as Accurate Reporters for Protein-Ligand Binding Interfaces in Solution

The availability of high-resolution 3D structural information is crucial for investigating guest-host systems across a wide range of fields. In the context of drug discovery, the information is routinely used to establish and validate structure-activity relationships, grow initial hits from screening campaigns, and to guide molecular docking. For the generation of protein-ligand complex structural information, X-ray crystallography is the experimental method of choice, however, with limited information on protein flexibility. An experimentally verified structural model of the binding interface in the native solution-state would support medicinal chemists in their molecular design decisions. Here we demonstrate that protein-bound ligand 1 H NMR chemical shifts are highly sensitive and accurate probes for the immediate chemical environment of protein-ligand interfaces. By comparing the experimental ligand 1 H chemical shift values with those computed from the X-ray structure using quantum mechanics methodology, we identify significant disagreements for parts of the ligand between the two experimental techniques. We show that quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) ensembles can be used to refine initial X-ray co-crystal structures resulting in a better agreement with experimental 1 H ligand chemical shift values. Overall, our findings highlight the usefulness of ligand 1 H NMR chemical shift information in combination with a QM/MM MD workflow for generating protein-ligand ensembles that accurately reproduce solution structural data.

Structural insights in the permeation mechanism of an activated GIRK2 channel

G protein-gated inwardly rectifying potassium (GIRK) channels play a significant role in physiopathology by the regulation of cell excitability. This regulation depends on the K+ ion conduction induced by structural constrictions: the selectivity filters (SFs), helix bundle crossings (HBCs), and G-loop gates. To explore why no permeation occurred when the constrictions were kept in the open state, a 4-K+-related occupancy mechanism was proposed. Unfortunately, this hypothesis was neither assessed, nor was the energetic characteristics presented. To identify the permeation mechanism on an atomic level, all-atom molecular dynamic (MD) simulations and a coupled quantum mechanics and molecular mechanics (QM/MM) method were used for the GIRK2 mutant R201A. It was found that the R201A had a moderate conductive capability in the presence of PIP2. Furthermore, the 4-K+ group of ions was found to dominate the conduction through the activated HBC gate. This shielding-like mechanism was assessed by the potential energy barrier along the conduction pathway. Mutation studies did further support the assumption that E152 was responsible for the mechanism. Moreover, E152 was most probably facilitating the inflow of ions from the SF to the cavity. On the contrary, N184 had no remarkable effect on this mechanism, except for the conduction efficiency. These findings highlighted the necessity of a multi-ion distribution for the conduction to take place, and indicated that the K+ migration was not only determined by the channel conductive state in the GIRK channel. The here presented multi-ion permeation mechanism may help to provide an effective way to regulate the channelopathies.

Biophysical characterization of the Plasmodium falciparum circumsporozoite protein's N-terminal domain

The circumsporozoite protein (CSP) is the main surface antigen of the Plasmodium sporozoite (SPZ) and forms the basis of the currently only licensed anti-malarial vaccine (RTS,S/AS01). CSP uniformly coats the SPZ and plays a pivotal role in its immunobiology, in both the insect and the vertebrate hosts. Although CSP's N-terminal domain (CSPN ) has been reported to play an important role in multiple CSP functions, a thorough biophysical and structural characterization of CSPN is currently lacking. Here, we present an alternative method for the recombinant production and purification of CSPN from Plasmodium falciparum (PfCSPN ), which provides pure, high-quality protein preparations with high yields. Through an interdisciplinary approach combining in-solution experimental methods and in silico analyses, we provide strong evidence that PfCSPN is an intrinsically disordered region displaying some degree of compaction.

Candida antarctica lipase B performance in organic solvent at varying water activities studied by molecular dynamics simulations

Applications of lipases in low-water environments are found across a broad range of industries, including the pharmaceutical and oleochemical sectors. This includes condensation reactions in organic solvents where the enzyme activity has been found to depend strongly on both the solvent and the water activity (aw). Despite several experimental and computational studies, knowledge is largely empirical, and a general predictive approach is much needed. To close this gap, we chose native Candida antarctica lipase B (CALB) and two mutants thereof and used molecular dynamics (MD) simulations to gain a molecular understanding of the effect of aw on the specific activity of CALB in hexane. Based on the simulations, we propose four criteria to understand the performance of CALB in organic media, which is supported by enzyme kinetics experiments. First, the lipase must be stable in the organic solvent, which was the case for native CALB and the two mutants studied here. Secondly, water clusters that form and grow close to the active site must not block the path of substrate molecules into the active site. Thirdly, the lipase's lid must not cover the active site. Finally, mutations and changes in aw must not disrupt the geometry of the active site. We show that mutating specific residues close to the active site can hinder water cluster formation and growth, making the lipase resistant to changes in aw. Our computational screening criteria could potentially be used to screen in-silico designed variants, so only promising candidates could be pushed forward to characterisation.

Revealing Drivers for Carboxy-S-adenosyl-l-methionine Use by Neomorphic Variants of a DNA Methyltransferase

Methylation of DNA plays a key role in diverse biological processes spanning from bacteria to mammals. DNA methyltransferases (MTases) typically employ S-adenosyl-l-methionine (SAM) as a critical cosubstrate and the relevant methyl donor for modification of the C5 position of cytosine. Recently, work on the CpG-specific bacterial MTase, M.MpeI, has shown that a single N374K point mutation can confer the enzyme with the neomorphic ability to use the sparse, naturally occurring metabolite carboxy-S-adenosyl-l-methionine (CxSAM) in order to generate the unnatural DNA modification, 5-carboxymethylcytosine (5cxmC). Here, we aimed to investigate the mechanistic basis for this DNA carboxymethyltransferase (CxMTase) activity by employing a combination of computational modeling and in vitro characterization. Modeling of substrate interactions with the enzyme variant allowed us to identify a favorable salt bridge between CxSAM and N374K that helps to rationalize selectivity of the CxMTase. Unexpectedly, we also discovered a potential role for a key active site E45 residue that makes a bidentate interaction with the ribosyl sugar of CxSAM, located on the opposite face of the CxMTase active site. Prompted by these modeling results, we further explored the space-opening E45D mutation and found that the E45D/N374K double mutant in fact inverts selectivity, preferring CxSAM over SAM in biochemical assays. These findings provide new insight into CxMTase active site architecture and may offer broader utility given the numerous opportunities offered by using SAM analogs for selective molecular labeling in concert with nucleic acid or even protein-modifying MTases.

In silico analysis of SOD1 aggregation inhibition modes of tertiary amine pyrazolone and pyrano coumarin ferulate as ALS drug candidates

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease, the familial form (fALS) of which is often cognate to mutations in the antioxidant enzyme Cu/Zn superoxide dismutase 1 (SOD1) leading to misfolding and aggregation. Two small molecules, a tertiary amine pyrazolone (TAP) and a pyrano coumarin ferulate (PCF) were suggested to be ALS drug candidates following experimental observation of their ability to inhibit SOD1 protein misfolding and aggregation. The present work aims at computational investigation of these experimentally proposed drug candidates to gain insight into their mechanism of SOD1 misfolding and aggregation inhibition. On the basis of molecular docking, molecular dynamics simulation, MM-PBSA and per-residue energy decomposition analysis, we examined the specific interactions of TAP and PCF with three probable binding sites of SOD1, namely, dimeric interface cavity, W32 and, UMP binding sites. Results suggest that the binding of TAP at W32 and at UMP sites are least probable due to absence of any favorable interaction. The binding of TAP to dimeric cavity is also unstable due to strong unfavorable interactions. In case of PCF, binding at the UMP site is least probable while binding at dimeric cavity is accompanied by unfavorable interactions. PCF, however, exhibits stable binding with the W32 binding site of SOD1 by stabilizing the solvent accessible hydrophobic residues, which otherwise would have acted as contact points for aggregation. Thus the results imply that compound PCF functions as an inhibitior of SOD1 misfolding/aggregation through direct interaction with the protein SOD1 at the W32 binding site. However, TAP is likely to act as an inhibitor through a different mechanism rather than direct interaction with the protein SOD1. These results apart from reinforcing previous experimental findings, shed light on the probable mechanism of action of the proposed drug candidates.

Vibrio cholerae's ToxRS bile sensing system

The seventh pandemic of the diarrheal cholera disease, which began in 1960, is caused by the Gram-negative bacterium Vibrio cholerae. Its environmental persistence provoking recurring sudden outbreaks is enabled by V. cholerae's rapid adaption to changing environments involving sensory proteins like ToxR and ToxS. Located at the inner membrane, ToxR and ToxS react to environmental stimuli like bile acid, thereby inducing survival strategies for example bile resistance and virulence regulation. The presented crystal structure of the sensory domains of ToxR and ToxS in combination with multiple bile acid interaction studies, reveals that a bile binding pocket of ToxS is only properly folded upon binding to ToxR. Our data proposes an interdependent functionality between ToxR transcriptional activity and ToxS sensory function. These findings support the previously suggested link between ToxRS and VtrAC-like co-component systems. Besides VtrAC, ToxRS is now the only experimentally determined structure within this recently defined superfamily, further emphasizing its significance. In-depth analysis of the ToxRS complex reveals its remarkable conservation across various Vibrio species, underlining the significance of conserved residues in the ToxS barrel and the more diverse ToxR sensory domain. Unravelling the intricate mechanisms governing ToxRS's environmental sensing capabilities, provides a promising tool for disruption of this vital interaction, ultimately inhibiting Vibrio's survival and virulence. Our findings hold far-reaching implications for all Vibrio strains that rely on the ToxRS system as a shared sensory cornerstone for adapting to their surroundings.

Reparameterization of Non-Bonded Parameters for Copper Ions in Plastocyanin: An Adaptive Force Matching Study

Molecular mechanics rely on existing experimental and theoretical inputs to confidently calculate the trajectories of molecular systems. These calculations, however, are often hindered by missing force field parameters. A notable subject of this problem is metal centers of proteins. This study parameterized, through an adaptive force matching (AFM) workflow, the copper cofactor of plastocyanin in its two oxidation states. New 12-6 Lennard-Jones (LJ) parameters and atomic partial charges were generated to complete the non-bonded description of the copper site. Our models show uniform distorted tetrahedral structures for reduced plastocyanin, Cu(I), and oxidized plastocyanin, Cu(II). These structures align with the QM/MM MD results and existing crystallography studies. TD-DFT calculations, meanwhile, showed that conformations with elongated axial Cu-SMet and shortened equatorial Cu-SCys bonds retain the experimental UV-Vis profile of blue copper (BC) proteins, thus signifying the importance of Cu-S interactions on BC proteins' unique spectroscopic properties.

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

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

Effect of temperature and nucleotide on the binding of BiP chaperone to a protein substrate

BiP (immunoglobulin heavy-chain binding protein) is a Hsp70 monomeric ATPase motor that plays broad and crucial roles in maintaining proteostasis inside the cell. Structurally, BiP is formed by two domains, a nucleotide-binding domain (NBD) with ATPase activity connected by a flexible hydrophobic linker to the substrate-binding domain. While the ATPase and substrate binding activities of BiP are allosterically coupled, the latter is also dependent on nucleotide binding. Recent structural studies have provided new insights into BiP's allostery; however, the influence of temperature on the coupling between substrate and nucleotide binding to BiP remains unexplored. Here, we study BiP's binding to its substrate at the single molecule level using thermo-regulated optical tweezers which allows us to mechanically unfold the client protein and explore the effect of temperature and different nucleotides on BiP binding. Our results confirm that the affinity of BiP for its protein substrate relies on nucleotide binding, by mainly regulating the binding kinetics between BiP and its substrate. Interestingly, our findings also showed that the apparent affinity of BiP for its protein substrate in the presence of nucleotides remains invariable over a wide range of temperatures, suggesting that BiP may interact with its client proteins with similar affinities even when the temperature is not optimal. Thus, BiP could play a role as a "thermal buffer" in proteostasis.

A Theoretical Study of Organotin Binding in Aromatase

The widely used organotin compounds are notorious for their acute toxicity. Experiments revealed that organotin might cause reproductive toxicity by reversibly inhibiting animal aromatase functioning. However, the inhibition mechanism is obscure, especially at the molecular level. Compared to experimental methods, theoretical approaches via computational simulations can help to gain a microscopic view of the mechanism. Here, in an initial attempt to uncover the mechanism, we combined molecular docking and classical molecular dynamics to investigate the binding between organotins and aromatase. The energetics analysis indicated that the van der Waals interaction is the primary driving force of binding the organic tail of organotin and the aromatase center. The hydrogen bond linkage trajectory analysis revealed that water plays a significant role in linking the ligand-water-protein triangle network. As an initial step in studying the mechanism of organotin inhibiting aromatase, this work provides an in-depth understanding of the binding mechanism of organotin. Further, our study will help to develop effective and environmentally friendly methods to treat animals that have already been contaminated by organotin, as well as sustainable solutions for organotin degradation.

AI-Driven Discovery of SARS-CoV-2 Main Protease Fragment-like Inhibitors with Antiviral Activity In Vitro

SARS-CoV-2 is the causative agent of COVID-19 and is responsible for the current global pandemic. The viral genome contains 5 major open reading frames of which the largest ORF1ab codes for two polyproteins, pp1ab and pp1a, which are subsequently cleaved into 16 nonstructural proteins (nsp) by two viral cysteine proteases encoded within the polyproteins. The main protease (Mpro, nsp5) cleaves the majority of the nsp's, making it essential for viral replication and has been successfully targeted for the development of antivirals. The first oral Mpro inhibitor, nirmatrelvir, was approved for treatment of COVID-19 in late December 2021 in combination with ritonavir as Paxlovid. Increasing the arsenal of antivirals and development of protease inhibitors and other antivirals with a varied mode of action remains a priority to reduce the likelihood for resistance emerging. Here, we report results from an artificial intelligence-driven approach followed by in vitro validation, allowing the identification of five fragment-like Mpro inhibitors with IC50 values ranging from 1.5 to 241 μM. The three most potent molecules (compounds 818, 737, and 183) were tested against SARS-CoV-2 by in vitro replication in Vero E6 and Calu-3 cells. Compound 818 was active in both cell models with an EC50 value comparable to its measured IC50 value. On the other hand, compounds 737 and 183 were only active in Calu-3, a preclinical model of respiratory cells, showing selective indexes twice as high as those for compound 818. We also show that our in silico methodology was successful in identifying both reversible and covalent inhibitors. For instance, compound 818 is a reversible chloromethylamide analogue of 8-methyl-γ-carboline, while compound 737 is an N-pyridyl-isatin that covalently inhibits Mpro. Given the small molecular weights of these fragments, their high binding efficiency in vitro and efficacy in blocking viral replication, these compounds represent good starting points for the development of potent lead molecules targeting the Mpro of SARS-CoV-2.

Uncovering the Mechanism of the Proton-Coupled Fluoride Transport in the CLCF Antiporter

Fluoride is a natural antibiotic abundantly present in the environment and, in micromolar concentrations, is able to inhibit enzymes necessary for bacteria to survive. However, as is the case with many antibiotics, bacteria have evolved resistance methods, including through the use of recently discovered membrane proteins. One such protein is the CLC^F F^-/H⁺ antiporter protein, a member of the CLC superfamily of anion-transport proteins. Though previous studies have examined this F^- transporter, many questions are still left unanswered. To reveal details of the transport mechanism used by CLC^F, we have employed molecular dynamics simulations and umbrella sampling calculations. Our results have led to several discoveries, including the mechanism of proton import and how it is able to aid in the fluoride export. Additionally, we have determined the role of the previously identified residues Glu118, Glu318, Met79, and Tyr396. This work is among the first studies of the CLC^F F^-/H⁺ antiporter and is the first computational investigation to model the full transport process, proposing a mechanism which couples the F^- export with the H⁺ import.

A Computational Biology Study on the Structure and Dynamics Determinants of Thermal Stability of the Chitosanase from Aspergillus fumigatus

Environmentally friendly and efficient biodegradation with chitosanase for degrading chitosan to oligosaccharide has been gaining more importance. Here, we studied a chitosanase from Aspergillus fumigatus with potential for production, but does not have the ideal thermal stability. The structure predicted by the Alphafold2 model, especially the binding site and two catalytic residues, has been found to have a high similarity with the experimental structure of the chitosanase V-CSN from the same family. The effects of temperature on structure and function were studied by dynamic simulation and the results showed that the binding site had high flexibility. After heating up from 300 K to 350 K, the RMSD and RMSF of the binding site increased significantly, in particular, the downward shift of loop6 closed the binding site, resulting in the spatial hindrance of binding. The time proportions of important hydrogen bonds at the binding site decreased sharply, indicating that serious disruption of hydrogen bonds should be the main interaction factor for conformational changes. The residues contributing energetically to binding were also revealed to be in the highly flexible region, which inevitably leads to the decrease in the activity stability at high temperature. These findings provide directions for the modification of thermal stability and perspectives on the research of proteins without experimental structures.

Simplified Derivatives of Tetrandrine as Potent and Specific P-gp Inhibitors to Reverse Multidrug Resistance in Cancer Chemotherapy

Targeted inhibition of a drug efflux transporter P-glycoprotein (P-gp) is an important strategy to reverse multidrug resistance in cancer chemotherapy. In this study, a rationally structural simplification to natural tetrandrine was performed based on molecular dynamics simulation and fragment growth, leading to an easily prepared, novel, and simplified compound OY-101 with high reversal activity and low cytotoxicity. Its excellent synergistic anti-cancer effect with vincristine (VCR) against drug-resistant cells Eca109/VCR was confirmed by reversal activity assay, flow cytometry, plate clone formation assay, and drug synergism analysis (IC50 = 9.9 nM, RF = 690). Further mechanism study confirmed that the OY-101 was a specific and efficient P-gp inhibitor. Importantly, OY-101 increased VCR sensitization in vivo without obvious toxicity. Overall, our findings may provide an alternative strategy for the design of novel specific P-gp inhibitor as an anti-tumor chemotherapy sensitizer.

Using High-Throughput Molecular Dynamics Simulation to Enhance the Computational Design of Kemp Elimination Enzymes

Computational enzyme design has been successfully applied to identify new alternatives to natural enzymes for the biosynthesis of important compounds. However, the moderate catalytic activities of de novo designed enzymes indicate that the modeling accuracy of current computational enzyme design methods should be improved. Here, high-throughput molecular dynamics simulations were used to enhance computational enzyme design, thus allowing the identification of variants with higher activities in silico. Different time schemes of high-throughput molecular dynamics simulations were tested to identify the catalytic features of evolved Kemp eliminases. The 20 × 1 ns molecular dynamics simulation scheme was sufficiently accurate and computationally viable to screen the computationally designed massive variants of Kemp elimination enzymes. The developed hybrid computational strategy was used to redesign the most active Kemp eliminase, HG3.17, and five variants were generated and experimentally confirmed to afford higher catalytic efficiencies than that of HG3.17, with one double variant (D52Q/A53S) exhibiting a 55% increase. The hybrid computational enzyme design strategy is general and computationally economical, with which we anticipate the efficient creation of practical enzymes for industrial biocatalysis.

Enhanced Grand Canonical Sampling of Occluded Water Sites Using Nonequilibrium Candidate Monte Carlo

Water molecules play a key role in many biomolecular systems, particularly when bound at protein-ligand interfaces. However, molecular simulation studies on such systems are hampered by the relatively long time scales over which water exchange between a protein and solvent takes place. Grand canonical Monte Carlo (GCMC) is a simulation technique that avoids this issue by attempting the insertion and deletion of water molecules within a given structure. The approach is constrained by low acceptance probabilities for insertions in congested systems, however. To address this issue, here, we combine GCMC with nonequilibium candidate Monte Carlo (NCMC) to yield a method that we refer to as grand canonical nonequilibrium candidate Monte Carlo (GCNCMC), in which the water insertions and deletions are carried out in a gradual, nonequilibrium fashion. We validate this new approach by comparing GCNCMC and GCMC simulations of bulk water and three protein binding sites. We find that not only is the efficiency of the water sampling improved by GCNCMC but that it also results in increased sampling of ligand conformations in a protein binding site, revealing new water-mediated ligand-binding geometries that are not observed using alternative enhanced sampling techniques.

Emergence of a proton exchange-based isomerization and lactonization mechanism in the plant coumarin synthase COSY

Plants contain rapidly evolving specialized enzymes that support the biosynthesis of functionally diverse natural products. In coumarin biosynthesis, a BAHD acyltransferase-family enzyme COSY was recently discovered to accelerate coumarin formation as the only known BAHD enzyme to catalyze an intramolecular acyl transfer reaction. Here we investigate the structural and mechanistic basis for COSY's coumarin synthase activity. Our structural analyses reveal an unconventional active-site configuration adapted to COSY's specialized activity. Through mutagenesis studies and deuterium exchange experiments, we identify a unique proton exchange mechanism at the α-carbon of the o-hydroxylated trans-hydroxycinnamoyl-CoA substrates during the catalytic cycle of COSY. Quantum mechanical cluster modeling and molecular dynamics further support this key mechanism for lowering the activation energy of the rate-limiting trans-to-cis isomerization step in coumarin production. This study unveils an unconventional catalytic mechanism mediated by a BAHD-family enzyme, and sheds light on COSY's evolutionary origin and its recruitment to coumarin biosynthesis in eudicots.

Thermostabilizing mechanisms of canonical single amino acid substitutions at a GH1 β-glucosidase probed by multiple MD and computational approaches

β-glucosidases play a pivotal role in second-generation biofuel (2G-biofuel) production. For this application, thermostable enzymes are essential due to the denaturing conditions on the bioreactors. Random amino acid substitutions have originated new thermostable β-glucosidases, but without a clear understanding of their molecular mechanisms. Here, we probe by different molecular dynamics simulation approaches with distinct force fields and submitting the results to various computational analyses, the molecular bases of the thermostabilization of the Paenibacillus polymyxa GH1 β-glucosidase by two-point mutations E96K (TR1) and M416I (TR2). Equilibrium molecular dynamic simulations (eMD) at different temperatures, principal component analysis (PCA), virtual docking, metadynamics (MetaDy), accelerated molecular dynamics (aMD), Poisson-Boltzmann surface analysis, grid inhomogeneous solvation theory and colony method estimation of conformational entropy allow to converge to the idea that the stabilization carried by both substitutions depend on different contributions of three classic mechanisms: (i) electrostatic surface stabilization; (ii) efficient isolation of the hydrophobic core from the solvent, with energetic advantages at the solvation cap; (iii) higher distribution of the protein dynamics at the mobile active site loops than at the protein core, with functional and entropic advantages. Mechanisms i and ii predominate for TR1, while in TR2, mechanism iii is dominant. Loop A integrity and loops A, C, D, and E dynamics play critical roles in such mechanisms. Comparison of the dynamic and topological changes observed between the thermostable mutants and the wildtype protein with amino acid co-evolutive networks and thermostabilizing hotspots from the literature allow inferring that the mechanisms here recovered can be related to the thermostability obtained by different substitutions along the whole family GH1. We hope the results and insights discussed here can be helpful for future rational approaches to the engineering of optimized β-glucosidases for 2G-biofuel production for industry, biotechnology, and science.

The pursuit of new alternative ways to eradicate Helicobacter pylori continues: Detailed characterization of interactions in the adenylosuccinate synthetase active site

Purine nucleotide synthesis is realised only through the salvage pathway in pathogenic bacterium Helicobacter pylori. Therefore, the enzymes of this pathway, among them also the adenylosuccinate synthetase (AdSS), present potential new drug targets. This paper describes characterization of His₆-tagged AdSS from H. pylori. Thorough analysis of 3D-structures of fully ligated AdSS (in a complex with guanosine diphosphate, 6-phosphoryl-inosine monophosphate, hadacidin and Mg²⁺) and AdSS in a complex with inosine monophosphate (IMP) only, enabled identification of active site interactions crucial for ligand binding and enzyme activity. Combination of experimental and molecular dynamics (MD) simulations data, particularly emphasized the importance of hydrogen bond Arg135-IMP for enzyme dimerization and active site formation. The synergistic effect of substrates (IMP and guanosine triphosphate) binding was suggested by MD simulations. Several flexible elements of the structure (loops) are stabilized by the presence of IMP alone, however loops comprising residues 287-293 and 40-44 occupy different positions in two solved H. pylori AdSS structures. MD simulations discovered the hydrogen bond network that stabilizes the closed conformation of the residues 40-50 loop, only in the presence of IMP. Presented findings provide a solid basis for the design of new AdSS inhibitors as potential drugs against H. pylori.

Cassava pullulanase and its synergistic debranching action with isoamylase 3 in starch catabolism

Pullulanase (EC 3.2.1.41, PUL), a debranching enzyme belonging to glycoside hydrolase family 13 subfamily 13, catalyses the cleavage of α-1,6 linkages of pullulan and β-limit dextrin. The present work studied PUL from cassava Manihot esculenta Crantz (MePUL) tubers, an important economic crop. The Mepul gene was successfully cloned and expressed in E. coli and rMePUL was biochemically characterised. MePUL was present as monomer and homodimer, as judged by apparent mass of ~ 84 - 197 kDa by gel permeation chromatography analysis. Optimal pH and temperature were at pH 6.0 and 50 °C, and enzyme activity was enhanced by the addition of Ca²⁺ ions. Pullulan is the most favourable substrate for rMePUL, followed by β-limit dextrin. Additionally, maltooligosaccharides were potential allosteric modulators of rMePUL. Interestingly, short-chain maltooligosaccharides (DP 2 - 4) were significantly revealed at a higher level when rMePUL was mixed with cassava isoamylase 3 (rMeISA3), compared to that of each single enzyme reaction. This suggests that MePUL and MeISA3 debranch β-limit dextrin in a synergistic manner, which represents a major starch catabolising process in dicots. Additionally, subcellular localisation suggested the involvement of MePUL in starch catabolism, which normally takes place in plastids.

Experimental-theoretical study of laccase as a detoxifier of aflatoxins

We investigate laccase-mediated detoxification of aflatoxins, fungal carcinogenic food contaminants. Our experimental comparison between two aflatoxins with similar structures (AFB₁ and AFG₂) shows significant differences in laccase-mediated detoxification. A multi-scale modeling approach (Docking, Molecular Dynamics, and Density Functional Theory) identifies the highly substrate-specific changes required to improve laccase detoxifying performance. We employ a large-scale density functional theory-based approach, involving more than 7000 atoms, to identify the amino acid residues that determine the affinity of laccase for aflatoxins. From this study we conclude: (1) AFB₁ is more challenging to degrade, to the point of complete degradation stalling; (2) AFG₂ is easier to degrade by laccase due to its lack of side products and favorable binding dynamics; and (3) ample opportunities to optimize laccase for aflatoxin degradation exist, especially via mutations leading to π-π stacking. This study identifies a way to optimize laccase for aflatoxin bioremediation and, more generally, contributes to the research efforts aimed at rational enzyme optimization.

Interactions between curcumin and human salt-induced kinase 3 elucidated from computational tools and experimental methods

Natural products are widely used for treating mitochondrial dysfunction-related diseases and cancers. Curcumin, a well-known natural product, can be potentially used to treat cancer. Human salt-induced kinase 3 (SIK3) is one of the target proteins for curcumin. However, the interactions between curcumin and human SIK3 have not yet been investigated in detail. In this study, we studied the binding models for the interactions between curcumin and human SIK3 using computational tools such as homology modeling, molecular docking, molecular dynamics simulations, and binding free energy calculations. The open activity loop conformation of SIK3 with the ketoenol form of curcumin was the optimal binding model. The I72, V80, A93, Y144, A145, and L195 residues played a key role for curcumin binding with human SIK3. The interactions between curcumin and human SIK3 were also investigated using the kinase assay. Moreover, curcumin exhibited an IC₅₀ (half-maximal inhibitory concentration) value of 131 nM, and it showed significant antiproliferative activities of 9.62 ± 0.33 µM and 72.37 ± 0.37 µM against the MCF-7 and MDA-MB-23 cell lines, respectively. This study provides detailed information on the binding of curcumin with human SIK3 and may facilitate the design of novel salt-inducible kinases inhibitors.

A theoretical study on binding and stabilization of galactose and novel galactose analogues to the human α-galactosidase A variant causing Fabry disease

α-galactosidase A (α-Gal A) catalyzes the hydrolysis of terminal α-galactosyl moieties from globotriaosylceramide, and mutations in this enzyme lead to the lipid metabolism disorder "Fabry disease". Mutation in α-Gal A possibly causes the protein misfolding, which reduces catalytic activity and stability of the enzyme. A recent study demonstrated that the binding of galactose on the α-Gal A catalytic site significantly increases its stability. Herein, the effect of mutation on secondary structure, structural energy, and galactose affinity of α-Gal A (wild type and A143T variant) was investigated using molecular dynamics simulations and free energy calculations based on MM/GBSA method. The results showed that A143T mutation caused the formation of unusual H-bonds that induced the change in secondary structure and binding affinities toward galactose. The amino acid residues involved in galactose binding were identified. The molecular binding mechanism obtained from this study could be helpful for optimizations and designs of new galactose analogs as pharmacological chaperones against Fabry disease.

Biologically relevant small variations of intra-cellular pH can have significant effect on stability of protein-DNA complexes, including the nucleosome

Stability of a protein-ligand complex may be sensitive to pH of its environment. Here we explore, computationally, stability of a set of protein-nucleic acid complexes using fundamental thermodynamic linkage relationship. The nucleosome, as well as an essentially random selection of 20 protein complexes with DNA or RNA, are included in the analysis. An increase in intra-cellular/intra-nuclear pH destabilizes most complexes, including the nucleosome. We propose to quantify the effect by ΔΔG^0.3-the change in the binding free energy due to pH increase of 0.3 units, corresponding to doubling of the H ⁺ activity; variations of pH of this amplitude can occur in living cells, including in the course of the cell cycle, and in cancer cells relative to normal ones. We suggest, based on relevant experimental findings, a threshold of biological significance of 1 2 k B T ( ∼ 0.3 k c a l / m o l ) for changes of stability of chromatin-related protein-DNA complexes: a change in the binding affinity above the threshold may have biological consequences. We find that for 70% of the examined complexes, Δ Δ G 0.3 > 1 2 k B T (for 10%, ΔΔG^0.3 is between 3 and 4 k _B T). Thus, small but relevant variations of intra-nuclear pH of 0.3 may have biological consequences for many protein-nucleic acid complexes. The binding affinity between the histone octamer and its DNA, which directly affects the DNA accessibility in the nucleosome, is predicted to be highly sensitive to intra-nuclear pH. A variation of 0.3 units results in ΔΔG^0.3 ∼ 10k _B T ( ∼ 6 k c a l / m o l ) ; for spontaneous unwrapping of 20 bp long entry/exit fragments of the nucleosomal DNA, ΔΔG^0.3 = 2.2k _B T; partial disassembly of the nucleosome into the tetrasome is characterized by ΔΔG^0.3 = 5.2k _B T. The predicted pH -induced modulations of the nucleosome stability are significant enough to suggest that they may have consequences relevant to the biological function of the nucleosome. Accessibility of the nucleosomal DNA is predicted to positively correlate with pH variations during the cell cycle; an increase in intra-cellular pH seen in cancer cells is predicted to lead to a more accessible nucleosomal DNA; a drop in pH associated with apoptosis is predicted to make nucleosomal DNA less accessible. We speculate that processes that depend on accessibility to the DNA in the nucleosomes, such as transcription or DNA replication, might become upregulated due to relatively small, but nevertheless realistic increases of intra-nuclear pH.

Molecular Dynamics Simulation of E412 Catalytic Residue Mutation of GOx-IPBCC

The enzyme glucose oxidase from Aspergillus niger has a homodimeric structure, consisting of two identical subunits with a molecular weight of 150,000 Daltons. In this study, we used the structure of the enzyme glucose oxidase from Aspergillus niger IPBCC.08.610 (GOx-IPBCC), this enzyme had a total activity of 92.87 U (μmol/min) and a Michaelis-Menten constant (Km) of 2.9 mM (millimolar). This study was conducted to predict the molecular dynamics of E412 (Glu412) residue catalytic mutation belonging to the GOx-IPBCC enzyme was determine the effect of changes in the catalytic residue on substrate binding (β-D-glucose). The results of molecular docking of 19 mutant structures, six E412 mutant homologous structures were selected (E412C, E412K, E412Q, E412T, E412, E412V, and E412W), which were evaluated using molecular dynamics simulation for 50 ns. The results showed a decrease in ∆G values in two mutant structures is E412C and E412T, and there is one mutant structure that increased ∆G values, namely E412W, these three mutant structures showed the best stability, bond interaction, and salt bridge profile according to molecular dynamics simulation.

RBD spatial orientation of the spike protein and its binding to ACE2: insight into the high infectivity of the SARS-CoV-2 Delta variant from MD simulations

The spike glycoprotein on the surface of the SARS-CoV-2 envelope plays an important role in its invasion into host cells. The binding of the spike glycoprotein RBD to the angiotensin-converting enzyme 2 (ACE2) receptor as a critical step in the spread of the virus has been explored intensively since the outbreak of COVID-19, but the high transmissibility of the virus such as the Delta variant is still not fully understood. Here, molecular simulations on the binding interactions of the wild-type spike protein and its four variants (Beta, Kappa, Delta, and Mu) with ACE2 and the antibody were performed, and the present results reveal that the residue mutations will not strengthen the binding affinity of the variant for ACE2, but remarkably influences the spatial orientation of the spike protein. Only the up-right conformational receptor binding domain (RBD) can bind ACE2, which is stabilized by the nearby RBDs in the down state, revealing that the RBD bears dual functional characteristics. The present results provide new insights into plausible mechanisms for high infectivity of the virus variants and their immune escape.

Role of the Triplet State and Protein Dynamics in the Formation and Stability of the Tryptophan Radical in an Apoazurin Mutant

The protein, azurin, has enabled the study of the tryptophan radical. Upon UV excitation of tyrosine-deficient apoazurin and in the presence of a Co(III) electron acceptor, the neutral radical (W48•) is formed. The lifetime of W48• in apoazurin is 41 s, which is shorter than the lifetime of several hours in Zn-substituted azurin. Molecular dynamics simulations revealed enhanced fluctuations of apoazurin which likely destabilize W48•. The photophysics of W48 was investigated to probe the precursor state for ET. The phosphorescence intensity was eliminated in the presence of an electron acceptor while the fluorescence was unchanged; this quenching of the phosphorescence is attributed to ET. The kinetics associated with W48• were examined with a model that incorporates intersystem crossing, ET, deprotonation, and decay of the cation radical. The estimated rate constants for ET (6 × 10⁶ s^–1) and deprotonation (3 × 10⁵ s^–1) are in agreement with a photoinduced mechanism where W48• is derived from the triplet state. The triplet as the precursor state for ET was supported by photolysis of apoazurin with 280 nm in the absence and presence of triplet-absorbing 405 nm light. Absorption bands from the neutral radical were observed only in the presence of blue light.

Structural features of Dnase1L3 responsible for serum antigen clearance

Autoimmunity develops when extracellular DNA released from dying cells is not cleared from serum. While serum DNA is primarily digested by Dnase1 and Dnase1L3, Dnase1 cannot rescue autoimmunity arising from Dnase1L3 deficiencies. Dnase1L3 uniquely degrades antigenic forms of cell-free DNA, including DNA complexed with lipids and proteins. The distinct activity of Dnase1L3 relies on its unique C-terminal Domain (CTD), but the mechanism is unknown. We used multiple biophysical techniques and functional assays to study the interplay between the core catalytic domain and the CTD. While the core domain resembles Dnase1, there are key structural differences between the two enzymes. First, Dnase1L3 is not inhibited by actin due to multiple differences in the actin recognition site. Second, the CTD augments the ability of the core to bind DNA, thereby facilitating the degradation of complexed DNA. Together, these structural insights will inform the development of Dnase1L3-based therapies for autoimmunity.

Streptomyces coelicolor macrodomain hydrolase SCO6735 cleaves thymidine-linked ADP-ribosylation of DNA

ADP-ribosylation is an ancient, highly conserved, and reversible covalent modification critical for a variety of endogenous processes in both prokaryotes and eukaryotes. ADP-ribosylation targets proteins, nucleic acids, and small molecules (including antibiotics). ADP-ribosylation signalling involves enzymes that add ADP-ribose to the target molecule, the (ADP-ribosyl)transferases; and those that remove it, the (ADP-ribosyl)hydrolases. Recently, the toxin/antitoxin pair DarT/DarG composed of a DNA ADP-ribosylating toxin, DarT, and (ADP-ribosyl)hydrolase antitoxin, DarG, was described. DarT modifies thymidine in single-stranded DNA in a sequence-specific manner while DarG reverses this modification, thereby rescuing cells from DarT toxicity. We studied the DarG homologue SCO6735 which is highly conserved in all Streptomyces species and known to be associated with antibiotic production in the bacterium S. coelicolor. SCO6735 shares a high structural similarity with the bacterial DarG and human TARG1. Like DarG and TARG1, SCO6735 can also readily reverse thymidine-linked ADP-ribosylation catalysed by DarT in vitro and in cells. SCO6735 active site analysis including molecular dynamic simulations of its complex with ADP-ribosylated thymidine suggests a novel catalytic mechanism of DNA-(ADP-ribose) hydrolysis. Moreover, a comparison of SCO6735 structure with ALC1-like homologues revealed an evolutionarily conserved feature characteristic for this subclass of macrodomain hydrolases.

Exploration of Diazaspiro Cores as Piperazine Bioisosteres in the Development of σ2 Receptor Ligands

A series of σ2R compounds containing benzimidazolone and diazacycloalkane cores was synthesized and evaluated in radioligand binding assays. Replacing the piperazine moiety in a lead compound with diazaspiroalkanes and the fused octahydropyrrolo[3,4-b] pyrrole ring system resulted in a loss in affinity for the σ2R. On the other hand, the bridged 2,5-diazabicyclo[2.2.1]heptane, 1,4-diazepine, and a 3-aminoazetidine analog possessed nanomolar affinities for the σ2R. Computational chemistry studies were also conducted with the recently published crystal structure of the σ2R/TMEM97 and revealed that hydrogen bond interactions with ASP29 and π-stacking interactions with TYR150 were largely responsible for the high binding affinity of small molecules to this protein.