Salt Bridge Research Articles

The unique salt bridge network in GlacPETase: a key to its stability.

The extensive accumulation of polyethylene terephthalate (PET) has become a critical environmental issue. PET hydrolases can break down PET into its building blocks. Recently, we identified a glacial PET hydrolase GlacPETase sharing less than 31% amino acid identity with any known PET hydrolases. In this study, the crystal structure of GlacPETase was determined at 1.8 Å resolution, revealing unique structural features including a distinctive N-terminal disulfide bond and a specific salt bridge network. Site-directed mutagenesis demonstrated that the disruption of the N-terminal disulfide bond did not reduce GlacPETase's thermostability or its catalytic activity on PET. However, mutations in the salt bridges resulted in changes in melting temperature ranging from -8°C to +2°C and the activity on PET ranging from 17.5% to 145.5% compared to the wild type. Molecular dynamics simulations revealed that these salt bridges stabilized the GlacPETase's structure by maintaining their surrounding structure. Phylogenetic analysis indicated that GlacPETase represented a distinct branch within PET hydrolases-like proteins, with the salt bridges and disulfide bonds in this branch being relatively conserved. This research contributed to the improvement of our comprehension of the structural mechanisms that dictate the thermostability of PET hydrolases, highlighting the diverse characteristics and adaptability observed within PET hydrolases.IMPORTANCEThe pervasive problem of polyethylene terephthalate (PET) pollution in various terrestrial and marine environments is widely acknowledged and continues to escalate. PET hydrolases, such as GlacPETase in this study, offered a solution for breaking down PET. Its unique origin and less than 31% identity with any known PET hydrolases have driven us to resolve its structure. Here, we report the correlation between its unique structure and biochemical properties, focusing on an N-terminal disulfide bond and specific salt bridges. Through site-directed mutagenesis experiments and molecular dynamics simulations, the roles of the N-terminal disulfide bond and salt bridges were elucidated in GlacPETase. This research enhanced our understanding of the role of salt bridges in the thermostability of PET hydrolases, providing a valuable reference for the future engineering of PET hydrolases.

Effects of anthocyanidins on the conformational transition of Aβ(1-42) peptide: Insights from molecular docking and molecular dynamics simulations

Recent evidence from in vitro and in vivo studies has shown that anthocyanins and anthocyanidins can reduce and inhibit the amyloid beta (Aβ) species, one of the hallmarks of Alzheimer's disease (AD). However, their inhibition mechanisms on Aβ species at molecular details remain elusive. Therefore, in the present study, molecular modelling methods were employed to investigate their inhibitory mechanisms on Aβ(1-42) peptide. The results highlighted that anthocyanidins effectively inhibited the conformational transitions of helices into beta-sheet (β-sheet) conformation within Aβ(1-42) peptide by two different mechanisms: 1) the obstruction of two terminals from coming into contact due to the binding of anthocyanidins with residues of N- and second hydrophobic core (SHC)-C-terminals, and 2) the prevention of the folding process due to the binding of anthocyanidin with the central polar (Asp23 and Lys28) and native helix (Asp23, Lys28, and Leu34) residues. These new findings on the inhibition of β-sheet formation by targeting both N- and SHC-C-terminals, and the long-established target, D23-K28 salt bridge residues, not with the conventional central hydrophobic core (CHC) as reported in the literature, might aid in designing more potent inhibitors for AD treatment.

Darcy-Forchheimer flow of Ag–ZnO–CoFe2O4/H2O Casson ternary hybrid nanofluid induced by a rotatory disk with EMHD

This study comprehensively inspects the steady, three-dimensional electromagnetohydrodynamics Darcy-Forchheimer slip flow upon a Casson ternary hybrid nanofluid through a porous rotatory disk. The study plunges into the repercussions of key factors namely thermophoresis, chemical reaction, viscous dissipation, Brownian movement, thermal radiation, and activation energy attributes of the ternary hybrid nanofluid. The current endeavour is inspired by the flourishing requirement in several energy-related and industrial applications. A mathematical configuration is formed to express the fluid flow system and heat and nanoparticle mass transpiration processes. Then partial differential equations administrating our mathematical model, are reformed into ordinary differential equations with the assistance of specific similarity conversions and eventually interpreted numerically using the noteworthy bvp4c scheme. The consequences of discrete flow factors upon temperature, concentration, and velocity profiles and Nusselt number, skin friction’s coefficient, and Sherwood number are investigated. The achieved outcomes disclose the electric parameter aids in upgrading the tangential velocity and temperature profiles. The outcomes reveal that for Casson ternary hybrid nanofluid, the heat and mass transpiration rates elevate by 19.23% and 51.28%, respectively, compared to Casson nanofluid whereas skin friction coefficient is retrogressed by approximately 8.96% for Casson ternary hybrid nanofluid compared to ternary hybrid nanofluid.

Atomic Force Microscopy and Molecular Dynamic Simulation of Adsorption of Polyacrylamide with Different Chemistries onto Calcium Carbonate.

This study investigates the interaction of polyacrylamide (PAM) of different functional groups (sulfonate vs. carboxylate) and charge density (30% hydrolysed vs. 10% hydrolysed) with calcium carbonate (CaCO3) via atomic force microscopy (AFM) and partly via molecular dynamic (MD) simulations. The PAM used were F3330 (30% hydrolysed), AN125 (25% sulfonated), and AN910 (% hydrolysed). A total of 100 ppm of PAMs was prepared in 0.1% NaCl, 3% NaCl, and 4.36% NaNO3 to be employed in AFM experiments, while oligomeric models (30 repeating units) of hydrolysed polyacrylamide (HPAM), sulfonated polyacrylamide (SPAM), and neutral PAM (NPAM) were studied on a model calcite surface on MD simulations. AFM analysis indicated that F3330 has a higher average adhesion and interaction energy with CaCO3 than AN125 due to the bulky sulfonate side group of AN125 interfering with SPAM adsorption. Steric repulsion of both PAMs was similar due to their comparable molecular weights and densities of the charged group. In contrast, AN910 showed lower average adhesion and interaction energy, along with slightly longer steric repulsion with calcite than F3330, suggesting AN910 adopts more loops and tails than the slightly flatter F3330 configuration. An increase in salt concentration from 0.1% to 3% NaCl saw a reduction in adhesion and interaction energy for F3330 and AN125 due to charge screening, while AN910 saw an increase, and these values increased further at 4.36% NaNO3. MD simulations revealed that the salt ions in the system formed salt bridges between PAM and calcite, indicating that the adhesion and interaction energy observed from AFM are likely to be the net balance between PAM charged group screening and salt bridging by the salt ions present. Salt ions with larger bare radii and smaller hydrated radii were shown to form stronger salt bridges.

An evolutionarily conserved phosphoserine-arginine salt bridge in the interface between ribosomal proteins uS4 and uS5 regulates translational accuracy in Saccharomyces cerevisiae.

Protein-protein and protein-rRNA interactions at the interface between ribosomal proteins uS4 and uS5 are thought to maintain the accuracy of protein synthesis by increasing selection of cognate aminoacyl-tRNAs. Selection involves a major conformational change-domain closure-that stabilizes aminoacyl-tRNA in the ribosomal acceptor (A) site. This has been thought a constitutive function of the ribosome ensuring consistent accuracy. Recently, the Saccharomyces cerevisiae Ctk1 cyclin-dependent kinase was demonstrated to ensure translational accuracy and Ser238 of uS5 proposed as its target. Surprisingly, Ser238 is outside the uS4-uS5 interface and no obvious mechanism has been proposed to explain its role. We show that the true target of Ctk1 regulation is another uS5 residue, Ser176, which lies in the interface opposite to Arg57 of uS4. Based on site specific mutagenesis, we propose that phospho-Ser176 forms a salt bridge with Arg57, which should increase selectivity by strengthening the interface. Genetic data show that Ctk1 regulates accuracy indirectly;the data suggest that the kinase Ypk2 directly phosphorylates Ser176. A second kinase pathway involving TORC1 and Pkc1 can inhibit this effect. The level of accuracy appears to depend on competitive action of these two pathways to regulate the level of Ser176 phosphorylation.

Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey.

The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop's flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI's binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.

Tumor-acquired somatic mutation affects conformation to abolish ABCG2-mediated drug resistance

ABCG2 is an important ATP-binding cassette transporter impacting the absorption and distribution of over 200 chemical toxins and drugs. ABCG2 also reduces the cellular accumulation of diverse chemotherapeutic agents. Acquired somatic mutations in the phylogenetically conserved amino acids of ABCG2 might provide unique insights into its molecular mechanisms of transport. Here, we identify a tumor-derived somatic mutation (Q393K) that occurs in a highly conserved amino acid across mammalian species. This ABCG2 mutant seems incapable of providing ABCG2-mediated drug resistance. This was perplexing because it is localized properly and retained interaction with substrates and nucleotides. Using a conformationally sensitive antibody, we show that this mutant appears “locked” in a non-functional conformation. Structural modeling and molecular dynamics simulations based on ABCG2 cryo-EM structures suggested that the Q393K interacts with the E446 to create a strong salt bridge. The salt bridge is proposed to stabilize the inward-facing conformation, resulting in an impaired transporter that lacks the flexibility to readily change conformation, thereby disrupting the necessary communication between substrate binding and transport.

Computational study of the piperidine and FtsZ interaction in Salmonella Typhi: implications for disrupting cell division machinery

FtsZ, a bacterial cell division protein, is essential for assembling the contractile Z-ring crucial in bacterial cytokinesis. Consequently, inhibiting FtsZ could impede proto-filaments, disrupting FtsZ and other associated proteins vital for cell division machinery. Conduct an in-silico drug interaction study to identify novel drug candidates that inhibit the FtsZ protein, aiming to prevent Multi-Drug Resistant (MDR) Salmonella Typhi. Data mining was performed based on piperidine compounds, which were subsequently screened for safe pharmacokinetic profiles. Compounds that met favorable drug-likeness criteria underwent virtual screening against the FtsZ drug target. Two compounds were chosen for molecular docking and molecular dynamic simulation to verify the binding affinity and stability between the target protein and the potential compounds. The 400 isoforms of piperidine analogues were curated, among them potent compound ZINC000000005416 found to possess high binding affinity (-8.49 kcal/mol) and low dissociation constant (0.597 µM). The highest binding affinity shown by ZINC000000005416 was validated by hydrogen bonds, hydrophobic interaction, and salt bridges with the functional domain of the cell division regulatory protein. Docking profiles, when correlated with molecular dynamic simulation (MDS) depicted stable trajectories and compatible conformational changes in the FtsZ-ZINC000000005416 complex. The stable simulated trajectories were validated through free-energy calculations using the Molecular Mechanics-Poisson Boltzmann Surface Area (MM/PBSA) module. Low energy conformations, although the simulation trajectory confirmed the stable ZINC000000005416-FtsZ interaction, which encouraged experimental validations. This study encourages further exploration of the compound ZINC000000005416 as a drug candidate inhibiting FtsZ protein against MDR Salmonella Typhi. Communicated by Ramaswamy H. Sarma

Natural Polymorphisms D60E and I62V Stabilize a Closed Conformation in HIV-1 Protease in the Absence of an Inhibitor or Substrate.

HIV infection remains a global health issue plagued by drug resistance and virological failure. Natural polymorphisms (NPs) contained within several African and Brazilian protease (PR) variants have been shown to induce a conformational landscape of more closed conformations compared to the sequence of subtype B prevalent in North America and Western Europe. Here we demonstrate through experimental pulsed EPR distance measurements and molecular dynamic (MD) simulations that the two common NPs D60E and I62V found within subtypes F and H can induce a closed conformation when introduced into HIV-1PR subtype B. Specifically, D60E alters the conformation in subtype B through the formation of a salt bridge with residue K43 contained within the nexus between the flap and hinge region of the HIV-1 PR fold. On the other hand, I62V modulates the packing of the hydrophobic cluster of the cantilever and fulcrum, also resulting in a more closed conformation.

The NMR structure of the Ea22 lysogenic developmental protein from lambda bacteriophage

The ea22 gene resides in a relatively uncharacterized region of the lambda bacteriophage genome between the exo and xis genes and is among the earliest genes transcribed upon infection. In lambda and Shiga toxin-producing phages found in enterohemorrhagic E. coli (EHEC) associated with food poisoning, Ea22 favors a lysogenic over lytic developmental state. The Ea22 protein may be considered in terms of three domains: a short amino-terminal domain, a coiled-coiled domain, and a carboxy-terminal domain (CTD). While the full-length protein is tetrameric, the CTD is dimeric when expressed individually. Here, we report the NMR solution structure of the Ea22 CTD that is described by a mixed alpha–beta fold with a dimer interface reinforced by salt bridges. A conserved mobile loop may serve as a ligand for an unknown host protein that works with Ea22 to promote bacterial survival and the formation of new lysogens. From sequence and structural comparisons, the CTD distinguishes lambda Ea22 from homologs encoded by Shiga toxin-producing bacteriophages.

How to evict HP1 from H3: Use your biophysical-force toolbox and try a complex salt bridge

How to evict HP1 from H3: Use your biophysical-force toolbox and try a complex salt bridge

Pharmacogenomics of poor drug metabolism in greyhounds: Canine P450 oxidoreductase genetic variation, breed heterogeneity, and functional characterization.

Greyhounds metabolize cytochrome P450 (CYP) 2B11 substrates more slowly than other dog breeds. However, CYP2B11 gene variants associated with decreased CYP2B11 expression do not fully explain reduced CYP2B11 activity in this breed. P450 oxidoreductase (POR) is an essential redox partner for all CYPs. POR protein variants can enhance or repress CYP enzyme function in a CYP isoform and substrate dependent manner. The study objectives were to identify POR protein variants in greyhounds and determine their effect on coexpressed CYP2B11 and CYP2D15 enzyme function. Gene sequencing identified two missense variants (Glu315Gln and Asp570Glu) forming four alleles, POR-H1 (reference), POR-H2 (570Glu), POR-H3 (315Gln, 570Glu) and POR-H4 (315Gln). Out of 68 dog breeds surveyed, POR-H2 was widely distributed across multiple breeds, while POR-H3 was largely restricted to greyhounds and Scottish deerhounds (35% allele frequencies), and POR-H4 was rare. Three-dimensional protein structure modelling indicated significant effects of Glu315Gln (but not Asp570Glu) on protein flexibility through loss of a salt bridge between Glu315 and Arg519. Recombinant POR-H1 (reference) and each POR variant (H2-H4) were expressed alone or with CYP2B11 or CYP2D15 in insect cells. No substantial effects on POR protein expression or enzyme activity (cytochrome c reduction) were observed for any POR variant (versus POR-H1) when expressed alone or with CYP2B11 or CYP2D15. Furthermore, there were no effects on CYP2B11 or CYP2D15 protein expression, or on CYP2D15 enzyme kinetics by any POR variant (versus POR-H1). However, Vmax values for 7-benzyloxyresorufin, propofol and bupropion oxidation by CYP2B11 were significantly reduced by coexpression with POR-H3 (by 34-37%) and POR-H4 (by 65-72%) compared with POR-H1. Km values were unaffected. Our results indicate that the Glu315Gln mutation (common to POR-H3 and POR-H4) reduces CYP2B11 enzyme function without affecting at least one other major canine hepatic P450 (CYP2D15). Additional in vivo studies are warranted to confirm these findings.

Force-induced salt bridge formation enhances HP1 homodimer interaction, implying augmented chromatin crosslinking and phase separation

Force-induced salt bridge formation enhances HP1 homodimer interaction, implying augmented chromatin crosslinking and phase separation

Biochemical and in silico structural characterization of a cold-active arginase from the psychrophilic yeast, Glaciozyma antarctica PI12.

Glaciozyma antarctica PI12 is a psychrophilic yeast isolated from Antarctica. In this work, we describe the heterologous production, biochemical properties and in silico structure analysis of an arginase from this yeast (GaArg). GaArg is a metalloenzyme that catalyses thehydrolysis ofL-arginine to L-ornithine and urea. The cDNA of GaArg was reversed transcribed, cloned, expressed and purified as a recombinant protein in Escherichia coli. The purified protein was active against L-arginine as its substrate in a reaction at 20°C, pH 9. At 10-35°C and pH 7-9, the catalytic activity of the protein was still present around 50%. Mn2+, Ni2+, Co2+ and K+ were able to enhance the enzyme activity more than two-fold, while GaArg is most sensitive to SDS, EDTA and DTT. The predicted structure model of GaArg showed a very similar overall fold with other known arginases. GaArg possesses predominantly smaller and uncharged amino acids, fewer salt bridges, hydrogen bonds and hydrophobic interactions compared to the other counterparts. GaArg is the first reported arginase that is cold-active, facilitated by unique structural characteristics for its adaptation of catalytic functions at low-temperature environments. The structure and function of cold-active GaArg provide insights into the potentiality of new applications in various biotechnology and pharmaceutical industries.

Testing charge-scaling protocols for more realistic description of salt bridges in molecular dynamics simulations of insulin analog

Testing charge-scaling protocols for more realistic description of salt bridges in molecular dynamics simulations of insulin analog

Zn2+ Binding Increases Parallel Structure in the Aβ(16-22) Oligomer by Disrupting Salt Bridge in Antiparallel Structure.

The aggregation of monomeric amyloid β protein (Aβ) into oligomers and amyloid plaque in the brain is associated with Alzheimer's disease. The hydrophobic central core Aβ16-22 has been widely studied due to its essential role in the fibrillization of full-length Aβ peptides. Compared to the homogeneous antiparallel structure of Aβ16-22 at the late stage, the early-stage prefibrillar aggregates contain varying proportions of different β structures. In this work, we studied the appearance probabilities of various self-assembly structures of Aβ16-22 and the effects of Zn2+ on these probabilities by replica exchange molecular dynamics simulations. It was found that at room temperature, Aβ16-22 can readily form assembled β-sheet structures in pure water, where a typical antiparallel arrangement dominates (24.8% of all sampled trimer structures). The addition of Zn2+ to the Aβ16-22 solution will dramatically decrease the appearance probability of antiparallel trimer structures to 12.5% by disrupting the formation of the Lys16-Glu22 salt bridge. Meanwhile, the probabilities of hybrid antiparallel/parallel structures increase. Our simulation results not only reveal the competition between antiparallel and parallel structures in the Aβ16-22 oligomers but also show that Zn2+ can affect the oligomer structures. The results also provide insights into the role of metal ions in the self-assembly of short peptides.

Molecular Characterization of Alpha-Amylase of Aspergillus flavus and Aspergillus oryzae: An In Silico Study

Alpha amylases (EC 3.2.1.1.) hydrolyze polysaccharides like starch and glycogen to create oligosaccharides of different sizes by random cleavage of internal 1, 4-glucosidic linkages. These enzymes are produced by many different bacteria and fungi. In this investigation, α-amylase from Aspergillus flavus and Aspergillus oryzae was taken into consideration. Amino acid sequences and protein structure were retrieved from databases. Sequences of Aspergillus oryzae possess a high number of charged residues whereas Aspergillus flavus have a higher abundance of uncharged polar and hydrophobic amino acid residues. 3 common superfamilies were observed between both species. The alignment of sequences showed their similar conservative nature. Non-covalent intra-protein interactions like the salt bridge, aromatic-aromatic interactions, aromatic-sulfur interactions and cation-pi interactions helped to increase the stability of α-amylase of Aspergillus oryzae. The tunnels and cavities also help to increase the functionality and catalytic activity of α-amylase. This study will also be helpful for protein engineering.

THE EFFECT OF THE D185 MUTATION ON THE STABILITY AND FUNCTIONALITY OF RUBISCO LIKE PROTEIN (RLP) FROM CHROMOHALOBACTER SALEXIGEN BKL 5

THE EFFECT OF THE D185 MUTATION ON THE STABILITY AND FUNCTIONALITY OF RUBISCO LIKE PROTEIN (RLP) FROM CHROMOHALOBACTER SALEXIGEN BKL 5

How Chromophore Labels Shape the Structure and Dynamics of a Peptide Hydrogel.

Biocompatible and functionalizable hydrogels have a wide range of (potential) medicinal applications. The hydrogelation process, particularly for systems with very low polymer weight percentages (<1 wt %), remains poorly understood, making it challenging to predict the self-assembly of a given molecular building block into a hydrogel. This severely hinders the rational design of self-assembled hydrogels. In this study, we demonstrate the impact of an N-terminal group on the self-assembly and rheology of the peptide hydrogel hFF03 (hydrogelating, fibril forming peptide 03) using molecular dynamics simulations, oscillatory shear rheology, and circular dichroism spectroscopy. We find that the chromophore and even its specific regioisomers have a significant influence on the microscopic structure and dynamics of the self-assembled fibril, and on the macroscopic mechanical properties. This is because the chromophore influences the possible salt bridges, which form and stabilize the fibril formation. Furthermore, we find that the solvation shell fibrils by itself cannot explain the viscoelasticity of hFF03 hydrogels. Our atomistic model of the hFF03 fibril formation enables a more rational design of these hydrogels. In particular, altering the N-terminal chromophore emerges as a design strategy to tune the mechanic properties of these self-assembled peptide hydrogels.

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.