Strength and co-operativity of contributions of surface salt bridges to protein stability

Network of ion-pairs is responsible for the large differences in the thermal stability of two structurally similar aminopeptidases

Aminopeptidases are an important class of enzymes for protein metabolism. Leucyl aminopeptidase (PepL) preferably removes leucine from the N-terminus of small peptides. PepL of Lacticaseibacillus casei was observed to be thermally unstable, while a structurally similar aminopeptidases T (AmpT) of Thermus thermophilus is highly stable. To understand the molecular interaction responsible for large difference in their stability, molecular dynamics simulations were carried out to study thermal stability of PepL and AmpT in 300 K to 450 K temperature range over 100 ns. PepL sampled a larger conformational space with a rugged free-energy landscape, while AmpT navigated a smoother energy landscape to reach global minimum. The RMSD, RMSF, radius of gyration and principal component analysis suggested large movements in PepL than in AmpT with increase in temperature. Analysis of residue-interaction network revealed AmpT possessing a greater number of low, medium and high energy contacts in comparison to PepL. AmpT showed a higher abundance of ion-pair clusters and ionic residues per cluster compared to PepL, moreover retained a greater number of high energy contacts at elevated temperatures. These findings showed that the inherently lower stability of PepL originates from a comparatively smaller number of contacts and can be pivotal in engineering PepL for higher stability.

The S2 Pocket Governs the Genus-Specific Substrate Selectivity of Coronavirus 3C-Like Protease

Coronavirus 3C-like protease (CoV 3CLpro) is essential for viral replication, providing an attractive target for monitoring the evolution of CoV and developing anti-CoV drugs. Here, the substrate-binding modes of 3CLpros from four CoV genera are analyzed and found that the S2 pocket in 3CLpro is highly conserved within each genus but differs between genera. Functionally, the S2 pocket, in conjunction with S4 and S1' pockets, governs the genus-specific substrate selectivity of 3CLpro. Resurrected ancestral 3CLpros from four CoV genera validate the genus-specific divergence of S2 pocket. Drawing upon the genus-specific S2 pocket as evolutionary marker, eight newly identified 3CLpros uncover the ancestral state of modern 3CLpro and elucidate the possible evolutionary process for CoV. It is also demonstrated that the S2 pocket is highly correlated with the genus-specific inhibitory potency of PF-07321332 (an FDA-approved drug against COVID-19) on different CoV 3CLpros. This study on 3CLpro provides novel insights to inform evolutionary mechanisms for CoV and develop genera-specific or broad-spectrum drugs against CoVs.

Methylation and phosphorylation of formin homology domain proteins (Fhod1 and Fhod3) by protein arginine methyltransferase 7 (PRMT7) and Rho Kinase (ROCK1)

Protein post translational modifications (PTMs) can regulate biological processes by altering an amino acid's bulkiness, charge, and hydrogen bonding interactions. Common modifications include phosphorylation, methylation, acetylation and ubiquitylation. Although a primary focus of studying PTMs is understanding the effects of a single amino acid modification, the possibility of additional modifications increases the complexity. For example, substrate recognition motifs for arginine methyltransferases and some serine/threonine kinases overlap, leading to potential enzymatic crosstalk. In this study we have shown that the human family of formin homology domain containing proteins (Fhods) contain a substrate recognition motif specific for human protein arginine methyltransferase 7 (PRMT7). In particular, PRMT7 methylates two arginine residues in the diaphanous autoinhibitory domain (DAD) of the family of Fhod proteins: R1588 and/or R1590 of Fhod3 isoform 4. Additionally, we confirmed that S1589 and S1595 in the DAD domain of Fhod3 can be phosphorylated by Rho/ROCK1 kinase. Significantly, we have determined that if S1589 is phosphorylated then PRMT7 cannot subsequently methylate R1588 or R1590. In contrast, if R1588 or R1590 of Fhod3 is methylated then ROCK1 phosphorylation activity is only slightly affected. Finally, we show that the interaction of the N-terminal DID domain can also inhibit the methylation of the DAD domain. Taken together these results suggest that the family of Fhod proteins, potential in vivo substrates for PRMT7, might be regulated by a combination of methylation and phosphorylation.

Structural and functional effects of phosphopriming and scaffolding in the kinase GSK-3β

Glycogen synthase kinase 3β (GSK-3β) targets specific signaling pathways in response to distinct upstream signals. We used structural and functional studies to dissect how an upstream phosphorylation step primes the Wnt signaling component β-catenin for phosphorylation by GSK-3β and how scaffolding interactions contribute to this reaction. Our crystal structure of GSK-3β bound to a phosphoprimed β-catenin peptide confirmed the expected binding mode of the phosphoprimed residue adjacent to the catalytic site. An aspartate phosphomimic in the priming site of β-catenin adopted an indistinguishable structure but reacted approximately 1000-fold slower than the native phosphoprimed substrate. This result suggests that substrate positioning alone is not sufficient for catalysis and that native phosphopriming interactions are necessary. We also obtained a structure of GSK-3β with an extended peptide from the scaffold protein Axin that bound with greater affinity than that of previously crystallized Axin fragments. This structure neither revealed additional contacts that produce the higher affinity nor explained how substrate interactions in the GSK-3β active site are modulated by remote Axin binding. Together, our findings suggest that phosphopriming and scaffolding produce small conformational changes or allosteric effects, not captured in the crystal structures, that activate GSK-3β and facilitate β-catenin phosphorylation. These results highlight limitations in our ability to predict catalytic activity from structure and have potential implications for the role of natural phosphomimic mutations in kinase regulation and phosphosite evolution.

Contribution of a surface salt bridge to the protein stability of deep-sea Shewanella benthica cytochrome c'

Two homologous cytochromes c', SBCP and SVCP, from deep-sea Shewanella benthica and Shewanella violacea respectively exhibit only nine surface amino acid substitutions, along with one at the N-terminus. Despite the small sequence difference, SBCP is thermally more stable than SVCP. Here, we examined the thermal stability of SBCP variants, each containing one of the nine substituted residues in SVCP, and found that the SBCP K87V variant was the most destabilized. We then determined the X-ray crystal structure of the SBCP K87V variant at a resolution of 2.1 Å. The variant retains a four-helix bundle structure similar to the wild-type, but notable differences are observed in the hydration structure around the mutation site. Instead of forming of the intrahelical salt bridge between Lys-87 and Asp-91 in the wild-type, a clathrate-like hydration around Val-87 through a hydrogen bond network with the nearby amino acid residues is observed. This network potentially enhances the ordering of surrounding water molecules, leading to an entropic destabilization of the protein. These results suggest that the unfavorable hydrophobic hydration environment around Val-87 and the inability to form the Asp-91-mediated salt bridge contribute to the observed difference in stability between SBCP and SVCP. These findings will be useful in future protein engineering for controlling protein stability through the manipulation of surface intrahelical salt bridges.

Stabilization of Synthetic Collagen Triple Helices: Charge Pairs and Covalent Capture

Collagen mimetic peptides are composed of triple helices. Triple helical formation frequently utilizes charge pair interactions to direct protein assembly. The design of synthetic triple helices is challenging due to the large number of competing species and the overall fragile nature of collagen mimetics. A successfully designed triple helix incorporates both positive and negative criteria to achieve maximum specificity of the supramolecular assembly. Intrahelical charge pair interactions, particularly those involved in lysine-aspartate and lysine-glutamate pairs, have been especially successful both in driving helix specificity and for subsequent stabilization by covalent capture. Despite this progress, the important sequential and geometric relationships of charged residues in a triple helical context have not been fully explored for either supramolecular assembly or covalent capture stabilization. In this study, we compare the eight canonical axial and lateral charge pairs of lysine and arginine with glutamate and aspartate to their noncanonical, reversed charge pairs. These findings are put into the context of collagen triple helical design and synthesis.

Understanding the Role of Intramolecular Ion-Pair Interactions in Conformational Stability Using an Ab Initio Thermodynamic Cycle

Intramolecular ion-pair interactions yield shape and functionality to many molecules. With proper orientation, these interactions overcome steric factors and are responsible for the compact structures of several peptides. In this study, we present a thermodynamic cycle based on isoelectronic and alchemical mutation to estimate the intramolecular ion-pair interaction energy. We determine these energies for 26 benchmark molecules with common ion-pair combinations and compare them with results obtained using intramolecular symmetry-adapted perturbation theory. For systems with long linkers, the ion-pair energies evaluated using both approaches deviate by less than 2.5% in the vacuum phase. The thermodynamic cycle based on density functional theory facilitates calculations of salt-bridge interactions in model tripeptides with continuum/microsolvation modeling and four large peptides: 1EJG (crambin), 1BDK (bradykinin), 1L2Y (a mini-protein with a tryptophan cage), and 1SCO (a toxin from the scorpion venom).

A widely distributed family of eukaryotic and bacterial deubiquitinases related to herpesviral large tegument proteins

Distinct families of eukaryotic deubiquitinases (DUBs) are regulators of ubiquitin signaling. Here, we report on the presence of an additional DUB class broadly distributed in eukaryotes and several bacteria. The only described members of this family are the large tegument proteins of herpesviruses, which are attached to the outside of the viral capsid. By using a bioinformatics screen, we have identified distant homologs of this VTD (Viral tegument-like DUB) family in vertebrate transposons, fungi, insects, nematodes, cnidaria, protists and bacteria. While some VTD activities resemble viral tegument DUBs in that they favor K48-linked ubiquitin chains, other members are highly specific for K6- or K63-linked ubiquitin chains. The crystal structures of K48- and K6-specific members reveal considerable differences in ubiquitin recognition. The VTD family likely evolved from non-DUB proteases and spread through transposons, many of which became 'domesticated', giving rise to the Drosophila male sterile (3)76Ca gene and several nematode genes with male-specific expression.

The Hidden Intricacies of Aquaporins: Remarkable Details in a Common Structural Scaffold

Evolution turned aquaporins (AQPs) into the most efficient facilitators of passive water flow through cell membranes at no expense of solute discrimination. In spite of a plethora of solved AQP structures, many structural details remain hidden. Here, by combining extensive sequence- and structural-based analysis of a unique set of 20 non-redundant high-resolution structures and molecular dynamics simulations of four representatives, key aspects of AQP stability, gating, selectivity, pore geometry, and oligomerization, with a potential impact on channel functionality, are identified. The general view of AQPs possessing a continuous open water pore is challenged and it is depicted that AQPs' selectivity is not exclusively shaped by pore-lining residues but also by the relative arrangement of transmembrane helices. Moreover, this analysis reveals that hydrophobic interactions constitute the main determinant of protein thermal stability. Finally, a numbering scheme of the conserved AQP scaffold is established, facilitating direct comparison of, for example, disease-causing mutations and prediction of potential structural consequences. Additionally, the results pave the way for the design of optimized AQP water channels to be utilized in biotechnological applications.

Native Salt Bridges Are a Key Regulator of Ubiquitin's Mechanical Stability

Although it is known that various intramolecular interactions determine protein mechanical stability, a detailed molecular-level understanding of the key regulators of protein mechanical stability is still lacking. Here, we present evidence for salt bridges in ubiquitin as important intramolecular interactions that can affect protein mechanical stability. Ubiquitin has two salt bridges: one relatively surface-exposed (SB1:K11-E34) and the other relatively buried (SB2:K27-D52). Ubiquitin is a reversible post-translational modifier and is stable mechanically (F_avg^u = 185 pN). On breaking SB1, the mechanical stability of ubiquitin is slightly enhanced (F_avg^u = 193 pN). In contrast, the mechanical stability significantly decreased upon breaking SB2 (F_avg^u = 158 pN). These results suggest that SB1 are SB2 are regulators of the mechanical stability of ubiquitin. Interestingly, the mechanical stability decreased further (F_avg^u = 145 pN) for the double salt bridge (DB) null variant. Monte Carlo simulations elucidate that the main regulating factor is the spontaneous unfolding rate constant (k_u⁰), being the highest for the DB null variant followed by the SB2 null variant, and it remains unaltered for the SB1 null variant, while the native-to-transition-state distance (x_u) remains unchanged. Our study provides mechanistic understanding on how two native salt bridges can independently regulate the mechanical stability in a protein, which has implications in designing protein-based robust biomaterials in the future.

A newly introduced salt bridge cluster improves structural and biophysical properties of de novo TIM barrels

Protein stability can be fine‐tuned by modifying different structural features such as hydrogen‐bond networks, salt bridges, hydrophobic cores, or disulfide bridges. Among these, stabilization by salt bridges is a major challenge in protein design and engineering since their stabilizing effects show a high dependence on the structural environment in the protein, and therefore are difficult to predict and model. In this work, we explore the effects on structure and stability of an introduced salt bridge cluster in the context of three different de novo TIM barrels. The salt bridge variants exhibit similar thermostability in comparison with their parental designs but important differences in the conformational stability at 25°C can be observed such as a highly stabilizing effect for two of the proteins but a destabilizing effect to the third. Analysis of the formed geometries of the salt bridge cluster in the crystal structures show either highly ordered salt bridge clusters or only single salt bridges. Rosetta modeling of the salt bridge clusters results in a good prediction of the tendency on stability changes but not the geometries observed in the three‐dimensional structures. The results show that despite the similarities in protein fold, the salt bridge clusters differently influence the structural and stability properties of the de novo TIM barrel variants depending on the structural background where they are introduced.

PDB Code(s): 7OSU, 7OT7, 7OSV, 7OT8 and 7P12;

Elucidation of Charge Contribution in Iridium-Chelated Hydrogen-Bonding Systems

We present two iridium complexes 1H⁺ and 2H⁺ that contain cationic ligands to extend the knowledge of charge-assisted hydrogen bonding (CAHB), which counts among the strongest non-covalent bonding interactions. Upon protonation, both complexes were converted into new hydrogen-bonding arrays with various selectivity for respective H-bonding partners. This study compares the association strengths of four hydrogen-bonding co-systems, emphasizing the roles of CAHB in supramolecular systems. We determined that the cationic charge in these systems contributed up to 2.7 kJ mol⁻¹ in the H-bonding complexation processes.

Double Mutant Cycles as a Tool to Address Folding, Binding, and Allostery

Quantitative measurement of intramolecular and intermolecular interactions in protein structure is an elusive task, not easy to address experimentally. The phenomenon denoted 'energetic coupling' describes short- and long-range interactions between two residues in a protein system. A powerful method to identify and quantitatively characterize long-range interactions and allosteric networks in proteins or protein-ligand complexes is called double-mutant cycles analysis. In this review we describe the thermodynamic principles and basic equations that underlie the double mutant cycle methodology, its fields of application and latest employments, and caveats and pitfalls that the experimentalists must consider. In particular, we show how double mutant cycles can be a powerful tool to investigate allosteric mechanisms in protein binding reactions as well as elusive states in protein folding pathways.

In Silico Antibody Mutagenesis for Optimizing Its Binding to Spike Protein of Severe Acute Respiratory Syndrome Coronavirus 2

Coronavirus disease 2019 (COVID-19) is an ongoing global pandemic, and there are currently no FDA-approved medicines for treatment or prevention. Inspired by promising outcomes for convalescent plasma treatment, the development of antibody drugs (biologics) to block SARS-CoV-2 infection has been the focus of drug discovery, along with tremendous efforts in repurposing small-molecule drugs. In the past several months, experimentally, many human neutralizing monoclonal antibodies (mAbs) were successfully extracted from plasma of recovered COVID-19 patients. Currently, several mAbs targeting the SARS-CoV-2's spike glycoprotein (S-protein) are in clinical trials. With known atomic structures of the mAb and S-protein complex, it becomes possible to investigate in silico the molecular mechanism of mAb's binding with S-protein and to design more potent mAbs through protein mutagenesis studies, complementary to existing experimental efforts. Leveraging today's superb computing power, we propose a fully automated in silico protocol for quickly identifying possible mutations in a mAb (e.g., CB6) to enhance its binding affinity for S-protein for the design of more efficacious therapeutic mAbs.

Activation Loop Dynamics Are Coupled to Core Motions in Extracellular Signal-Regulated Kinase-2

The activation loop segment in protein kinases is a common site for regulatory phosphorylation. In extracellular signal-regulated kinase 2 (ERK2), dual phosphorylation and conformational rearrangement of the activation loop accompany enzyme activation. X-ray structures show the active conformation to be stabilized by multiple ion pair interactions between phosphorylated threonine and tyrosine residues in the loop and six arginine residues in the kinase core. Despite the extensive salt bridge network, nuclear magnetic resonance Carr-Purcell-Meiboom-Gill relaxation dispersion experiments show that the phosphorylated activation loop is conformationally mobile on a microsecond to millisecond time scale. The dynamics of the loop match those of previously reported global exchange within the kinase core region and surrounding the catalytic site that have been found to facilitate productive nucleotide binding. Mutations in the core region that alter these global motions also alter the dynamics of the activation loop. Conversely, mutations in the activation loop perturb the global exchange within the kinase core. Together, these findings provide evidence for coupling between motions in the activation loop and those surrounding the catalytic site in the active state of the kinase. Thus, the activation loop segment in dual-phosphorylated ERK2 is not held statically in the active X-ray conformation but instead undergoes exchange between conformers separated by a small energetic barrier, serving as part of a dynamic allosteric network controlling nucleotide binding and catalytic function.

Understanding Thermostability Factors of Barley Limit Dextrinase by Molecular Dynamics Simulations

Limit dextrinase (LD) is the only endogenous starch-debranching enzyme in barley (Hordeum vulgare, Hv), which is the key factor affecting the production of a high degree of fermentation. Free LD will lose its activity in the mashing process at high temperature in beer production. However, there remains a lack of understanding on the factor affecting the themostability of HvLD at the atomic level. In this work, the molecular dynamics simulations were carried out for HvLD to explore the key factors affecting the thermal stability of LD. The higher value of root mean square deviation (RMSD), radius of gyration (R_g), and surface accessibility (SASA) suggests the instability of HvLD at high temperatures. Intra-protein hydrogen bonds and hydrogen bonds between protein and water decrease at high temperature. Long-lived hydrogen bonds, salt bridges, and hydrophobic contacts are lost at high temperature. The salt bridge interaction analysis suggests that these salt bridges are important for the thermostability of HvLD, including E568–R875, D317–R378, D803–R884, D457–R214, D468–R395, D456–R452, D399–R471, and D541–R542. Root mean square fluctuation (RMSF) analysis identified the thermal-sensitive regions of HvLD, which will facilitate enzyme engineering of HvLD for enhanced themostability.

The unfolding transition state of ubiquitin with charged residues has higher energy than that with hydrophobic residues

The native-state structure and folding pathways of a protein are encoded in its amino acid sequence.

Structural insight into conformational change in prion protein by breakage of electrostatic network around H187 due to its protonation

A conformational change from normal prion protein(PrP^C) to abnormal prion protein(PrP^SC) induces fatal neurodegenerative diseases. Acidic pH is well-known factors involved in the conformational change. Because the protonation of H187 is strongly linked to the change in PrP stability, we examined the charged residues R156, E196, and D202 around H187. Interestingly, there have been reports on pathological mutants, such as H187R, E196A, and D202N. In this study, we focused on how an acidic pH and pathological mutants disrupt this electrostatic network and how this broken network destabilizes PrP structure. To do so, we performed a temperature-based replica-exchange molecular dynamics (T-REMD) simulation using a cumulative 252 μs simulation time. We measured the distance between amino acids comprising four salt bridges (R156–E196/D202 and H187–E196/D202). Our results showed that the spatial configuration of the electrostatic network was significantly altered by an acidic pH and mutations. The structural alteration in the electrostatic network increased the RMSF value around the first helix (H1). Thus, the structural stability of H1, which is anchored to the H2–H3 bundle, was decreased. It induces separation of R156 from the electrostatic network. Analysis of the anchoring energy also shows that two salt-bridges (R156-E196/D202) are critical for PrP stability.

Electrostatically Driven Guanidinium Interaction Domains that Control Hydrogel-Mediated Protein Delivery In Vivo

Protein biologics are an important class of drugs, but the necessity for frequent parenteral administration is a major limitation. Drug-delivery materials offer a potential solution, but protein-material adsorption can cause denaturation, which reduces their effectiveness. Here, we describe a new protein delivery platform that limits direct contact between globular protein domains and material matrix, yet from a single subcutaneous administration can be tuned for long-term drug release. The strategy utilizes complementary electrostatic interactions made between a suite of designed interaction domains (IDs), installed onto the terminus of a protein of interest, and a negatively charged self-assembled fibrillar hydrogel. These intermolecular interactions can be easily modulated by choice of ID to control material interaction and desorption energies, which allows regulation of protein release kinetics to fit desired release profiles. Molecular dynamics studies provided a molecular-level understanding of the mechanisms that govern release and identified optimal binding zones on the gel fibrils that facilitate strong ID-material interactions, which are crucial for sustained release of protein. This delivery platform can be easily loaded with cargo, is shear-thin syringe implantable, provides improved protein stability, is capable of a diverse range of in vitro release rates, and most importantly, can accomplish long-term control over in vivo protein delivery.

Design of a Helical-Stabilized, Cyclic, and Nontoxic Analogue of the Peptide Cm-p5 with Improved Antifungal Activity

Following the information obtained by a rational design study, a cyclic and helical-stabilized analogue of the peptide Cm-p5 was synthetized. The cyclic monomer showed an increased activity in vitro against Candida albicans and Candida parapsilosis, compared to Cm-p5. Initially, 14 mutants of Cm-p5 were synthesized following a rational design to improve the antifungal activity and pharmacological properties. Antimicrobial testing showed that the activity was lost in each of these 14 analogues, suggesting, as a main conclusion, that a Glu-His salt bridge could stabilize Cm-p5 helical conformation during the interaction with the plasma membrane. A derivative, obtained by substitution of Glu and His for Cys, was synthesized and oxidized with the generation of a cyclic monomer with improved antifungal activity. In addition, two dimers were generated during the oxidation procedure, a parallel and antiparallel one. The dimers showed a helical secondary structure in water, whereas the cyclic monomer only showed this conformation in SDS. Molecular dynamic simulations confirmed the helical stabilizations for all of them, therefore indicating the possible essential role of the Glu-His salt bridge. In addition, the antiparallel dimer showed a moderate activity against Pseudomonas aeruginosa and a significant activity against Listeria monocytogenes. Neither the cyclic monomer nor the dimers were toxic against macrophages or THP-1 human cells. Due to its increased capacity for fungal control compared to fluconazole, its low cytotoxicity, together with a stabilized α-helix and disulfide bridges, that may advance its metabolic stability, and in vivo activity, the new cyclic Cm-p5 monomer represents a potential systemic antifungal therapeutic candidate.

Bridging solvent molecules mediate RNase A - Ligand binding

Due to its high catalytic activity and readily available supply, ribonuclease A (RNase A) has become a pivotal enzyme in the history of protein science. Moreover, this great interest has carried over to computational chemistry and molecular dynamics, where RNase A has become a model system for various types of studies, all the while being an important drug design target in its own right. Here, we present a detailed molecular dynamics study of RNase–ligand binding involving 22 compounds, spanning nearly five orders of magnitude in affinity, and totaling 8.8 μs of sampling with the standard Amber parameters and an additional 8.8 μs of sampling with a modified potential. We show that short-lived, solvent-mediated bridging interactions are crucial to RNase–ligand binding. We characterize the behavior of bridging solvent molecules, uncovering a power-law dependence between the lifetime of a solvent bridge and the probability of its occurrence. We also demonstrate that from an energetic perspective, bridging solvent in RNase A–ligand binding behaves like part of the enzyme, rather than the ligands. Moreover, we describe an automated pipeline for the detection and processing of bridging interactions, and offer an independent assessment of the performance of the state-of-the-art fixed-charge force fields. Thus, our work has broad implications for drug design and computational chemistry in general.

The electrostatic core of the outer membrane protein X from E. coli

Electrostatic side chain contacts can contribute substantial interaction energy terms to the stability of proteins. The impact of electrostatic interactions on the structure and architecture of outer membrane proteins is however not well studied compared to soluble proteins. Here, we report the results of a systematic study of all charged side chains of the E. coli outer membrane protein X (OmpX). The data identify three distinct salt-bridge clusters in the core of OmpX that contribute significantly to protein stability in dodecylphosphocholine detergent micelles. The three clusters form an "electrostatic core" of the membrane protein OmpX, corresponding in its architectural role to the hydrophobic core of soluble proteins. This article is part of a Special Issue entitled: Molecular biophysics of membranes and membrane proteins.

Preservation of Protein Zwitterionic States in the Transition from Solution to Gas Phase Revealed by Sodium Adduction Mass Spectrometry

The structural characterization of proteins and their interaction network mapping in the gas phase highlights the need to preserve their most nativelike conformers in the transition from the solution to gas phase. Zwitterionic interactions in a protein are weak bonds between oppositely charged residues, which make an important contribution to protein stability. However, it is still not clear whether the native zwitterionic states of proteins can be retained or not when it is transferred from the solution to gas phase. Using the nonspecific Na⁺ adduction as a novel signature, here we show that the zwitterionic states of proteins can be preserved when a moderated droplet desolvation condition (temperature <30 °C) is used in native electrospray ionization mass spectrometry. The very low-level nonspecific metal adduction to proteins under such conditions also enables rapid and direct determination of the binding states of metal-binding proteins and sensitive detection of proteins from solutions containing highly concentrated involatile salts (e.g., 50 mM NaCl). We believe that our findings can be instructive for performing mass spectrometric analysis of proteins and useful for protein ions desalting which simply involves altering the temperature and flow rate of drying gas in the desolvation region.

Side chain electrostatic interactions and pH-dependent expansion of the intrinsically disordered, highly acidic carboxyl-terminus of γ-tubulin

Intramolecular electrostatic attraction and repulsion strongly influence the conformational sampling of intrinsically disordered proteins and domains (IDPs). In order to better understand this complex relationship, we have used nuclear magnetic resonance to measure side chain pK_a values and pH-dependent translational diffusion coefficients for the unstructured and highly acidic carboxyl-terminus of γ-tubulin (γ-CT), providing insight into how the net charge of an IDP relates to overall expansion or collapse of the conformational ensemble. Many of the pK_a values in the γ-CT are shifted upward by 0.3-0.4 units and exhibit negatively cooperative ionization pH profiles, likely due to the large net negative charge that accumulates on the molecule as the pH is raised. pK_a shifts of this magnitude correspond to electrostatic interaction energies between the affected residues and the rest of the charged molecule that are each on the order of 1 kcal mol^-1 . Diffusion of the γ-CT slowed with increasing net charge, indicative of an expanding hydrodynamic radius (r_H ). The degree of expansion agreed quantitatively with what has been seen from comparisons of IDPs with different charge content, yielding the general trend that every 0.1 increase in relative charge (|Q|/res) produces a roughly 5% increase in r_H . While γ-CT pH titration data followed this trend nearly perfectly, there were substantially larger deviations for the database of different IDP sequences. This suggests that other aspects of an IDP's primary amino acid sequence beyond net charge influence the sensitivity of r_H to electrostatic interactions.

The pre-mRNA splicing and transcription factor Tat-SF1 is a functional partner of the spliceosome SF3b1 subunit via a U2AF homology motif interface

The transcription elongation and pre-mRNA splicing factor Tat-SF1 associates with the U2 small nuclear ribonucleoprotein (snRNP) of the spliceosome. However, the direct binding partner and underlying interactions mediating the Tat-SF1-U2 snRNP association remain unknown. Here, we identified SF3b1 as a Tat-SF1-interacting subunit of the U2 snRNP. Our 1.1 Å resolution crystal structure revealed that Tat-SF1 contains a U2AF homology motif (UHM) protein-protein interaction module. We demonstrated that Tat-SF1 preferentially and directly binds the SF3b1 subunit compared with other U2AF ligand motif (ULM)-containing splicing factors, and further established that SF3b1 association depends on the integrity of the Tat-SF1 UHM. We next compared the Tat-SF1-binding affinities for each of the five known SF3b1 ULMs and then determined the structures of representative high- and low-affinity SF3b1 ULM complexes with the Tat-SF1 UHM at 1.9 Å and 2.1 Å resolutions, respectively. These structures revealed a canonical UHM-ULM interface, comprising a Tat-SF1 binding pocket for a ULM tryptophan (SF3b1 Trp³³⁸) and electrostatic interactions with a basic ULM tail. Importantly, we found that SF3b1 regulates Tat-SF1 levels and that these two factors influence expression of overlapping representative transcripts, consistent with a functional partnership of Tat-SF1 and SF3b1. Altogether, these results define a new molecular interface of the Tat-SF1-U2 snRNP complex for gene regulation.

Salt Bridge in Ligand-Protein Complexes-Systematic Theoretical and Statistical Investigations

Although the salt bridge is the strongest among all known noncovalent molecular interactions, no comprehensive studies have been conducted to date to examine its role and significance in drug design. Thus, a systematic study of the salt bridge in biological systems is reported herein, with a broad analysis of publicly available data from Protein Data Bank, DrugBank, ChEMBL, and GPCRdb. The results revealed the distance and angular preferences as well as privileged molecular motifs of salt bridges in ligand-receptor complexes, which could be used to design the strongest interactions. Moreover, using quantum chemical calculations at the MP2 level, the energetic, directionality, and spatial variabilities of salt bridges were investigated using simple model systems mimicking salt bridges in a biological environment. Additionally, natural orbitals for chemical valence (NOCV) combined with the extended-transition-state (ETS) bond-energy decomposition method (ETS-NOCV) were analyzed and indicated a strong covalent contribution to the salt bridge interaction. The present results could be useful for implementation in rational drug design protocols.

RankProt: A multi criteria-ranking platform to attain protein thermostabilizing mutations and its in vitro applications - Attribute based prediction method on the principles of Analytical Hierarchical Process

Attaining recombinant thermostable proteins is still a challenge for protein engineering. The complexity is the length of time and enormous efforts required to achieve the desired results. Present work proposes a novel and economic strategy of attaining protein thermostability by predicting site-specific mutations at the shortest possible time. The success of the approach can be attributed to Analytical Hierarchical Process and the outcome was a rationalized thermostable mutation(s) prediction tool- RankProt. Briefly the method involved ranking of 17 biophysical protein features as class predictors, derived from 127 pairs of thermostable and mesostable proteins. Among the 17 predictors, ionic interactions and main-chain to main-chain hydrogen bonds were the highest ranked features with eigen value of 0.091. The success of the tool was judged by multi-fold in silico validation tests and it achieved the prediction accuracy of 91% with AUC 0.927. Further, in vitro validation was carried out by predicting thermostabilizing mutations for mesostable Bacillus subtilis lipase and performing the predicted mutations by multi-site directed mutagenesis. The rationalized method was successful to render the lipase thermostable with optimum temperature stability and T_m increase by 20°C and 7°C respectively. Conclusively it can be said that it was the minimum number of mutations in comparison to the number of mutations incorporated to render Bacillus subtilis lipase thermostable, by directed evolution techniques. The present work shows that protein stabilizing mutations can be rationally designed by balancing the biophysical pleiotropy of proteins, in accordance to the selection pressure.

pH dependent membrane binding of the Solanum tuberosum plant specific insert: An in silico study

The Solanum tuberosum plant-specific insert (StPSI) has been shown to possess potent antimicrobial activity against both human and plant pathogens. Furthermore, in vitro, the StPSI is capable of fusing phospholipid vesicles, provided the conditions of net anionic vesicle charge and acidic pH are met. Constant pH replica-exchange simulations indicate several acidic residues on the dimer have highly perturbed pK_as (<3.0; E15, D28, E85 & E100) due to involvement in salt bridges. After setting the pH of the system to either 3.0 or 7.4, all-atom simulations provided details of the effect of pH on secondary structural elements, particularly in the previously unresolved crystallographic structure of the loop section. Coarse-grained dimer-bilayer simulations demonstrated that at pH 7.4, the dimer had no affinity for neutral or anionic membranes over the course of 1 μs simulations. Conversely, at pH 3.0 two binding modes were observed. Mode 1 is mediated primarily via strong N-terminal interactions on one monomer only, whereas in mode 2, N- and C-terminal residues of one monomer and numerous polar and basic residues on the second monomer, particularly in the third helix, participate in membrane interactions. Mode 2 was accompanied by re-orientation of the dimer to a more vertical position with respect to helices 1 and 4, positioning the dimer for membrane interactions. These results offer the first examination at near-atomic resolution of residues mediating the StPSI-membrane interactions, and allow for the postulation of a possible fusion mechanism.

Principles of Protein Stability and Their Application in Computational Design

Proteins are increasingly used in basic and applied biomedical research. Many proteins, however, are only marginally stable and can be expressed in limited amounts, thus hampering research and applications. Research has revealed the thermodynamic, cellular, and evolutionary principles and mechanisms that underlie marginal stability. With this growing understanding, computational stability design methods have advanced over the past two decades starting from methods that selectively addressed only some aspects of marginal stability. Current methods are more general and, by combining phylogenetic analysis with atomistic design, have shown drastic improvements in solubility, thermal stability, and aggregation resistance while maintaining the protein's primary molecular activity. Stability design is opening the way to rational engineering of improved enzymes, therapeutics, and vaccines and to the application of protein design methodology to large proteins and molecular activities that have proven challenging in the past.

Measuring inter-protein pairwise interaction energies from a single native mass spectrum by double-mutant cycle analysis

The strength and specificity of protein complex formation is crucial for most life processes and is determined by interactions between residues in the binding partners. Double-mutant cycle analysis provides a strategy for studying the energetic coupling between amino acids at the interfaces of such complexes. Here we show that these pairwise interaction energies can be determined from a single high-resolution native mass spectrum by measuring the intensities of the complexes formed by the two wild-type proteins, the complex of each wild-type protein with a mutant protein, and the complex of the two mutant proteins. This native mass spectrometry approach, which obviates the need for error-prone measurements of binding constants, can provide information regarding multiple interactions in a single spectrum much like nuclear Overhauser effects (NOEs) in nuclear magnetic resonance. Importantly, our results show that specific inter-protein contacts in solution are maintained in the gas phase.

Double mutant cycle (DMC) analyses can provide the interaction energies between amino acids at the interface of protein complexes. Here, the authors determine pairwise interaction energies using high-resolution native mass spectroscopy, offering a straightforward route for the DMC methodology.

Influence of Glu/Arg, Asp/Arg, and Glu/Lys Salt Bridges on α-Helical Stability and Folding Kinetics

Using a combination of ultraviolet circular dichroism, temperature-jump transient-infrared spectroscopy, and molecular dynamics simulations, we investigate the effect of salt bridges between different types of charged amino-acid residue pairs on α-helix folding. We determine the stability and the folding and unfolding rates of 12 alanine-based α-helical peptides, each of which has a nearly identical composition containing three pairs of positively and negatively charged residues (either Glu(-)/Arg(+), Asp(-)/Arg(+), or Glu(-)/Lys(+)). Within each set of peptides, the distance and order of the oppositely charged residues in the peptide sequence differ, such that they have different capabilities of forming salt bridges. Our results indicate that stabilizing salt bridges (in which the interacting residues are spaced and ordered such that they favor helix formation) speed up α-helix formation by up to 50% and slow down the unfolding of the α-helix, whereas salt bridges with an unfavorable geometry have the opposite effect. Comparing the peptides with different types of charge pairs, we observe that salt bridges between side chains of Glu(-) and Arg(+) are most favorable for the speed of folding, probably because of the larger conformational space of the salt-bridging Glu(-)/Arg(+) rotamer pairs compared to Asp(-)/Arg(+) and Glu(-)/Lys(+). We speculate that the observed impact of salt bridges on the folding kinetics might explain why some proteins contain salt bridges that do not stabilize the final, folded conformation.

A novel self-activation mechanism of Candida antarctica lipase B

The lipase CALB might have two competing lid-opening mechanisms: self-activation and surface-activation.

Smoothing a rugged protein folding landscape by sequence-based redesign

The rugged folding landscapes of functional proteins puts them at risk of misfolding and aggregation. Serine protease inhibitors, or serpins, are paradigms for this delicate balance between function and misfolding. Serpins exist in a metastable state that undergoes a major conformational change in order to inhibit proteases. However, conformational labiality of the native serpin fold renders them susceptible to misfolding, which underlies misfolding diseases such as α₁-antitrypsin deficiency. To investigate how serpins balance function and folding, we used consensus design to create conserpin, a synthetic serpin that folds reversibly, is functional, thermostable, and polymerization resistant. Characterization of its structure, folding and dynamics suggest that consensus design has remodeled the folding landscape to reconcile competing requirements for stability and function. This approach may offer general benefits for engineering functional proteins that have risky folding landscapes, including the removal of aggregation-prone intermediates, and modifying scaffolds for use as protein therapeutics.

Invertase Stabilization by Chemical Modification of Sugar Chains with Carboxymethylcellulose

Invertase from Saccharomyces cerevisiae was activated by periodate treatment and further reacted with ethylenediamine/sodium borohydride. Carboxymethylcellulose was then attached to ethylenediamine-modified invertase using 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide as coupling agent. The modified enzyme contained 3.5 mol of polysaccharide per mol of holoenzyme, and retained about 56% of the initial invertase activity. The thermostability of invertase increased from 64 to 70° C, the thermal inactivation at different temperatures ranging from 60 to 70° C was markedly reduced for polymer-modified enzyme. An increase of 9.1 kJ mol−1 in activation free energy of inactivation was determined for invertase after modification. Functional stability was increased for carboxymethylcellulose-invertase complex in the range of pH between 2.0 and 12.0. The conjugate was also more resistant to denaturation by 6 M urea solution.

Oxidation increases the strength of the methionine-aromatic interaction

Oxidation of methionine disrupts the structure and function of a range of proteins, but little is understood about the chemistry that underlies these perturbations. Using quantum mechanical calculations, we found that oxidation increased the strength of the methionine-aromatic interaction motif, a driving force for protein folding and protein-protein interaction, by 0.5-1.4 kcal/mol. We found that non-hydrogen-bonded interactions between dimethyl sulfoxide (a methionine analog) and aromatic groups were enriched in both the Protein Data Bank and Cambridge Structural Database. Thermal denaturation and NMR spectroscopy experiments on model peptides demonstrated that oxidation of methionine stabilized the interaction by 0.5-0.6 kcal/mol. We confirmed the biological relevance of these findings through a combination of cell biology, electron paramagnetic resonance spectroscopy and molecular dynamics simulations on (i) calmodulin structure and dynamics, and (ii) lymphotoxin-α binding toTNFR1. Thus, the methionine-aromatic motif was a determinant of protein structural and functional sensitivity to oxidative stress.

Tyrosine Residues from the S4-S5 Linker of Kv11.1 Channels Are Critical for Slow Deactivation

Slow deactivation of Kv11.1 channels is critical for its function in the heart. The S4-S5 linker, which joins the voltage sensor and pore domains, plays a critical role in this slow deactivation gating. Here, we use NMR spectroscopy to identify the membrane-bound surface of the S4S5 linker, and we show that two highly conserved tyrosine residues within the KCNH subfamily of channels are membrane-associated. Site-directed mutagenesis and electrophysiological analysis indicates that Tyr-542 interacts with both the pore domain and voltage sensor residues to stabilize activated conformations of the channel, whereas Tyr-545 contributes to the slow kinetics of deactivation by primarily stabilizing the transition state between the activated and closed states. Thus, the two tyrosine residues in the Kv11.1 S4S5 linker play critical but distinct roles in the slow deactivation phenotype, which is a hallmark of Kv11.1 channels.

A Soluble, Folded Protein without Charged Amino Acid Residues

Charges are considered an integral part of protein structure and function, enhancing solubility and providing specificity in molecular interactions. We wished to investigate whether charged amino acids are indeed required for protein biogenesis and whether a protein completely free of titratable side chains can maintain solubility, stability, and function. As a model, we used a cellulose-binding domain from Cellulomonas fimi, which, among proteins of more than 100 amino acids, presently is the least charged in the Protein Data Bank, with a total of only four titratable residues. We find that the protein shows a surprising resilience toward extremes of pH, demonstrating stability and function (cellulose binding) in the pH range from 2 to 11. To ask whether the four charged residues present were required for these properties of this protein, we altered them to nontitratable ones. Remarkably, this chargeless protein is produced reasonably well in Escherichia coli, retains its stable three-dimensional structure, and is still capable of strong cellulose binding. To further deprive this protein of charges, we removed the N-terminal charge by acetylation and studied the protein at pH 2, where the C-terminus is effectively protonated. Under these conditions, the protein retains its function and proved to be both soluble and have a reversible folding-unfolding transition. To the best of our knowledge, this is the first time a soluble, functional protein with no titratable side chains has been produced.

A comparison of X-ray and calculated structures of the enzyme MTH1

Modern computational chemistry methods provide a powerful tool for use in refining the geometry of proteins determined by X-ray crystallography. Specifically, computational methods can be used to correctly place hydrogen atoms unresolved by this experimental method and improve bond geometry accuracy. Using the semiempirical method PM7, the structure of the nucleotide-sanitizing enzyme MTH1, complete with hydrolyzed substrate 8-oxo-dGMP, was optimized and the resulting geometry compared with the original X-ray structure of MTH1. After determining hydrogen atom placement and the identification of ionized sites, the charge distribution in the binding site was explored. Where comparison was possible, all the theoretical predictions were in good agreement with experimental observations. However, when these were combined with additional predictions for which experimental observations were not available, the result was a new and alternative description of the substrate-binding site interaction. An estimate was made of the strengths and weaknesses of the PM7 method for modeling proteins on varying scales, ranging from overall structure to individual interatomic distances. An attempt to correct a known fault in PM7, the under-estimation of steric repulsion, is also described. This work sheds light on the specificity of the enzyme MTH1 toward the substrate 8-oxo-dGTP; information that would facilitate drug development involving MTH1. Graphical Abstract Overlay of the backbone traces of the two MTH1 protein chains (green and orange respectively) in PDB 3ZR0 and the equivalent PM7 structures (magenta and cyan respectively) each optimized separately.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

Using Kinetic Network Models To Probe Non-Native Salt-Bridge Effects on α-Helix Folding

Salt-bridge interactions play an important role in stabilizing many protein structures, and have been shown to be designable features for protein design. In this work, we study the effects of non-native salt bridges on the folding of a soluble alanine-based peptide (Fs peptide) using extensive all-atom molecular dynamics simulations performed on the Folding@home distributed computing platform. Using Markov State Models, we show how non-native salt-bridges affect the folding kinetics of Fs peptide by perturbing specific conformational states. Furthermore, we present methods for the automatic detection and analysis of such states. These results provide insight into helix folding mechanisms and useful information to guide simulation-based computational protein design.

Structural and dynamic properties that govern the stability of an engineered fibronectin type III domain

Consensus protein design is a rapid and reliable technique for the improvement of protein stability, which relies on the use of homologous protein sequences. To enhance the stability of a fibronectin type III (FN3) domain, consensus design was employed using an alignment of 2123 sequences. The resulting FN3 domain, FN3con, has unprecedented stability, with a melting temperature >100°C, a ΔG_D−N of 15.5 kcal mol⁻¹ and a greatly reduced unfolding rate compared with wild-type. To determine the underlying molecular basis for stability, an X-ray crystal structure of FN3con was determined to 2.0 Å and compared with other FN3 domains of varying stabilities. The structure of FN3con reveals significantly increased salt bridge interactions that are cooperatively networked, and a highly optimized hydrophobic core. Molecular dynamics simulations of FN3con and comparison structures show the cooperative power of electrostatic and hydrophobic networks in improving FN3con stability. Taken together, our data reveal that FN3con stability does not result from a single mechanism, but rather the combination of several features and the removal of non-conserved, unfavorable interactions. The large number of sequences employed in this study has most likely enhanced the robustness of the consensus design, which is now possible due to the increased sequence availability in the post-genomic era. These studies increase our knowledge of the molecular mechanisms that govern stability and demonstrate the rising potential for enhancing stability via the consensus method.