Contribution of Cation−π Interactions to Protein Stability

Lipases for targeted industrial applications, focusing on the development of biotechnologically significant aspects: A comprehensive review of recent trends in protein engineering

Lipases are remarkable biocatalysts, adept at catalyzing the breakdown of diverse compounds into glycerol, fatty acids, and mono- and di-glycerides via hydrolysis. Beyond this, they facilitate esterification, transesterification, alcoholysis, acidolysis, and more, making them versatile in industrial applications. In industrial processes, lipases that exhibit high stability are favored as they can withstand harsh conditions. However, most native lipases are unable to endure adverse conditions, making them unsuitable for industrial use. Protein engineering proves to be a potent technology in the development of lipases that can function effectively under challenging conditions and fulfill criteria for various industrial processes. This review concentrated on new trends in protein engineering to enhance the diversity of lipase genes and employed in silico methods for predicting and comprehensively analyzing target mutations in lipases. Additionally, key molecular factors associated with industrial characteristics of lipases, including thermostability, solvent tolerance, catalytic activity, and substrate preference have been elucidated. The present review delved into how industrial traits can be enhanced through directed evolution (epPCR, gene shuffling), rational design (FRESCO, ASR), combined engineering strategies (i.e. CAST, ISM, and FRISM) as protein engineering methodologies in contexts of biodiesel production, food processing, and applications of detergent, pharmaceutics, and plastic degradation.

Supra-biological performance of immobilized enzymes enabled by chaperone-like specific non-covalent interactions

Designing complex synthetic materials for enzyme immobilization could unlock the utility of biocatalysis in extreme environments. Inspired by biology, we investigate the use of random copolymer brushes as dynamic immobilization supports that enable supra-biological catalytic performance of immobilized enzymes. This is demonstrated by immobilizing Bacillus subtilis Lipase A on brushes doped with aromatic moieties, which can interact with the lipase through multiple non-covalent interactions. Incorporation of aromatic groups leads to a 50 °C increase in the optimal temperature of lipase, as well as a 50-fold enhancement in enzyme activity. Single-molecule FRET studies reveal that these supports act as biomimetic chaperones by promoting enzyme refolding and stabilizing the enzyme's folded and catalytically active state. This effect is diminished when aromatic residues are mutated out, suggesting the importance of π-stacking and π-cation interactions for stabilization. Our results underscore how unexplored enzyme-support interactions may enable uncharted opportunities for using enzymes in industrial biotransformations.

Mediator structure and function in transcription initiation

Recent advances in cryo-electron microscopy have led to multiple structures of Mediator in complex with the RNA polymerase II (Pol II) transcription initiation machinery. As a result we now hold in hands near-complete structures of both yeast and human Mediator complexes and have a better understanding of their interactions with the Pol II pre-initiation complex (PIC). Herein, we provide a summary of recent achievements and discuss their implications for future studies of Mediator and its role in gene regulation.

Systematic elucidation of the second coordination sphere effect on the structure and properties of a blue copper protein, pseudoazurin

The rational structural and computational studies of a blue copper protein, pseudoazurin (PAz), and its Met16X (X = Phe, Leu, Val, Ile) variants gave clear functional meanings of the noncovalent interaction (NCI) through the second coordination sphere. The high-resolution X-ray crystal structures of Met16X PAz demonstrated that the active site geometry is significantly affected by the substitution of Met16, which is located within the NCI distance from the His81 imidazole ring at the copper active site. The computational chemistry calculations based on the crystal structure analyses confirmed that the NCI of S-π/CH-π (wild-type), π-π (Met16Phe), double CH-π (Met16Leu), and single CH-π (Met16Val and Met16Ile). The estimated interaction energies for the NCI demonstrated that the fine-tuning of the protein stability and Cu site properties form the second coordination sphere of PAz.

Controlling Doxorubicin Release from a Peptide Hydrogel through Fine-Tuning of Drug-Peptide Fiber Interactions

Hydrogels are versatile materials that have emerged in the last few decades as promising candidates for a range of applications in the biomedical field, from tissue engineering and regenerative medicine to controlled drug delivery. In the drug delivery field, in particular, they have been the subject of significant interest for the spatially and temporally controlled delivery of anticancer drugs and therapeutics. Self-assembling peptide-based hydrogels, in particular, have recently come to the fore as potential candidate vehicles for the delivery of a range of drugs. In order to explore how drug–peptide interactions influence doxorubicin (Dox) release, five β-sheet-forming self-assembling peptides with different physicochemical properties were used for the purpose of this study, namely: FEFKFEFK (F8), FKFEFKFK (FK), FEFEFKFE (FE), FEFKFEFKK (F8K), and KFEFKFEFKK (KF8K) (F: phenylalanine; E: glutamic acid; K: lysine). First, Dox-loaded hydrogels were characterized to ensure that the incorporation of the drug did not significantly affect the hydrogel properties. Subsequently, Dox diffusion out of the hydrogels was investigated using UV absorbance. The amount of drug retained in F8/FE composite hydrogels was found to be directly proportional to the amount of charge carried by the peptide fibers. When cation−π interactions were used, the position and number of end-lysine were found to play a key role in the retention of Dox. In this case, the amount of Dox retained in F8/KF8K composite hydrogels was linked to the amount of end-lysine introduced, and an end-lysine/Dox interaction stoichiometry of 3/1 was obtained. For pure FE and KF8K hydrogels, the maximum amount of Dox retained was also found to be related to the overall concentration of the hydrogels and, therefore, to the overall fiber surface area available for interaction with the drug. For 14 mM hydrogel, ∼170–200 μM Dox could be retained after 24 h. This set of peptides also showed a broad range of susceptibilities to enzymatic degradation opening the prospect of being able to control also the rate of degradation of these hydrogels. Finally, the Dox released from the hydrogel was shown to be active and affect 3T3 mouse fibroblasts viability in vitro. Our study clearly shows the potential of this peptide design as a platform for the formulation of injectable or sprayable hydrogels for controlled drug delivery.

3D Interaction Homology: Hydropathic Analyses of the "π-Cation" and "π-π" Interaction Motifs in Phenylalanine, Tyrosine, and Tryptophan Residues

Three-dimensional (3D) maps of the hydropathic environments of protein amino acid residues are information-rich descriptors of preferred conformations, interaction types and energetics, and solvent accessibility. The interactions made by each residue are the primary factor for rotamer selection and the secondary, tertiary, and even quaternary protein structure. Our evolving basis set of environmental data for each residue type can be used to understand the protein structure. This work focuses on the aromatic residues phenylalanine, tyrosine, and tryptophan and their structural roles. We calculated and analyzed side chain-to-environment 3D maps for over 70,000 residues of these three types that reveal, with respect to hydrophobic and polar interactions, the environment around each. After binning with backbone ϕ/ψ and side chain χ₁, we clustered each bin by 3D similarities between map-map pairs. For each of the three residue types, four bins were examined in detail: one in the β-pleat, two in the right-hand α-helix, and one in the left-hand α-helix regions of the Ramachandran plot. For high degrees of side chain burial, encapsulation of the side chain by hydrophobic interactions is ubiquitous. The more solvent-exposed side chains are more likely to be involved in polar interactions with their environments. Evidence for π-π interactions was observed in about half of the residues surveyed [phenylalanine (PHE): 53.3%, tyrosine (TYR): 34.1%, and tryptophan (TRP): 55.7%], but on an energy basis, this contributed to only ∼4% of the total. Evidence for π-cation interactions was observed in 14.1% of PHE, 8.3% of TYR, and 26.8% of TRP residues, but on an energy basis, this contributed to only ∼1%. This recognition of even these subtle interactions in the 3D hydropathic environment maps is key support for our interaction homology paradigm of protein structure elucidation and possibly prediction.

Phylogenetic, functional and structural characterization of a GH10 xylanase active at extreme conditions of temperature and alkalinity

Endoxylanases active under extreme conditions of temperature and alkalinity can replace the use of highly pollutant chemicals in the pulp and paper industry. Searching for enzymes with these properties, we carried out a comprehensive bioinformatics study of the GH10 family. The phylogenetic analysis allowed the construction of a radial cladogram in which protein sequences putatively ascribed as thermophilic and alkaliphilic appeared grouped in a well-defined region of the cladogram, designated TAK Cluster. One among five TAK sequences selected for experimental analysis (Xyn11) showed extraordinary xylanolytic activity under simultaneous conditions of high temperature (90 °C) and alkalinity (pH 10.5). Addition of a carbohydrate binding domain (CBM2) at the C-terminus of the protein sequence further improved the activity of the enzyme at high pH. Xyn11 structure, which has been solved at 1.8 Å resolution by X-ray crystallography, reveals an unusually high number of hydrophobic, ionic and hydrogen bond atomic interactions that could account for the enzyme’s extremophilic nature.

Engineering Hydrogel Adhesion for Biomedical Applications via Chemical Design of the Junction

Hydrogel adhesion inherently relies on engineering the contact surface at soft and hydrated interfaces. Upon contact, adhesion normally occurs through the formation of chemical or physical interactions between the disparate surfaces. The ability to form these adhesion junctions is challenging for hydrogels as the interfaces are wet and deformable and often contain low densities of functional groups. In this Review, we link the design of the binding chemistries or adhesion junctions, whether covalent, dynamic covalent, supramolecular, or physical, to the emergent adhesive properties of soft and hydrated interfaces. Wet adhesion is useful for bonding to or between tissues and implants for a range of biomedical applications. We highlight several recent and emerging adhesive hydrogels for use in biomedicine in the context of efficient junction design. The main focus is on engineering hydrogel adhesion through molecular design of the junctions to tailor the adhesion strength, reversibility, stability, and response to environmental stimuli.

Parallel molecular mechanisms for enzyme temperature adaptation

The mechanisms that underly the adaptation of enzyme activities and stabilities to temperature are fundamental to our understanding of molecular evolution and how enzymes work. Here, we investigate the molecular and evolutionary mechanisms of enzyme temperature adaption, combining deep mechanistic studies with comprehensive sequence analyses of thousands of enzymes. We show that temperature adaptation in ketosteroid isomerase (KSI) arises primarily from one residue change with limited, local epistasis, and we establish the underlying physical mechanisms. This residue change occurs in diverse KSI backgrounds, suggesting parallel adaptation to temperature. We identify residues associated with organismal growth temperature across 1005 diverse bacterial enzyme families, suggesting widespread parallel adaptation to temperature. We assess the residue properties, molecular interactions, and interaction networks that appear to underly temperature adaptation.

Discovery of Novel Cyclic Salt Bridge in Thermophilic Bacterial Protease and Study of its Sequence and Structure

The plausible explanation behind the stability of thermophilic protein is still yet to be defined more clearly. Here, an in silico study has been undertaken by investigating the sequence and structure of protease from thermophilic (tPro) bacteria and mesophilic (mPro) bacteria. Results showed that charged and uncharged polar residues have higher abundance in tPro. In extreme environment, the tPro is stabilized by high number of isolated and network salt bridges. A novel cyclic salt bridge is also found in a structure of tPro. High number of metal ion-binding site also helps in protein stabilization of thermophilic protease. Aromatic-aromatic interactions also play a crucial role in tPro stabilization. Formation of long network aromatic-aromatic interactions also first time reported here. Finally, the present study provides a major insight with a newly identified cyclic salt bridge in the stability of the enzyme, which may be helpful for protein engineering. It is also used in industrial applications for human welfare.

Accurate Description of Cation-π Interactions in Proteins with a Nonpolarizable Force Field at No Additional Cost

Cation-π interactions play a significant role in a host of processes eminently relevant to biology. However, polarization effects arising from the interaction of cations with aromatic moieties have long been recognized to be inadequately described by pairwise additive force fields. In the present work, we address this longstanding shortcoming through the nonbonded FIX (NBFIX) feature of the CHARMM36 force field, modifying pair-specific Lennard-Jones (LJ) parameters, while circumventing the limitations of the Lorentz-Berthelot combination rules. The potentials of mean force (PMFs) characterizing prototypical cation-π interactions in aqueous solutions are first determined using a hybrid quantum mechanical/molecular mechanics (QM/MM) strategy in conjunction with an importance-sampling algorithm. The LJ parameters describing the cation-π pairs are then optimized to match the QM/MM PMFs. The standard binding free energies of nine cation-π complexes, i.e., toluene, para-cresol, and 3-methyl-indole interacting with either ammonium, guanidinium, or tetramethylammonium, determined with this new set of parameters agree well with the experimental measurements. Additional simulations were carried out on three different classes of biological objects featuring cation-π interactions, including five individual proteins, three protein-ligand complexes, and two protein-protein complexes. Our results indicate that the description of cation-π interactions is overall improved using NBFIX corrections, compared with the standard pairwise additive force field. Moreover, an accurate binding free energy calculation for a protein-ligand complex containing cation-π interactions (2BOK) shows that using the new parameters, the experimental binding affinity can be reproduced quantitatively. Put together, the present work suggests that the NBFIX parameters optimized here can be broadly utilized in the simulation of proteins in an aqueous solution to enhance the representation of cation-π interactions, at no additional computational cost.

Comparative bioinformatic and structural analyses of pepsin and renin

Pepsin, the archetypal pepsin-like aspartic protease, is irreversibly denatured when exposed to neutral pH conditions whereas renin, a structural homologue of pepsin, is fully stable and optimally active in the same conditions despite sharing highly similar enzyme architecture. To gain insight into the structural determinants of differential aspartic protease pH stability, the present study used comparative bioinformatic and structural analyses. In pepsin, an abundance of polar and aspartic acid residues were identified, a common trait with other acid-stable enzymes. Conversely, renin was shown to have increased levels of basic amino acids. In both pepsin and renin, the solvent exposure of these charged groups was high. Having similar overall acidic residue content, the solvent-exposed basic residues may allow for extensive salt bridge formation in renin, whereas in pepsin, these residues are protonated and serve to form stabilizing hydrogen bonds at low pH. Relative differences in structure and sequence in the turn and joint regions of the β-barrel and ψ-loop in both the N- and C-terminal lobes were identified as regions of interest in defining divergent pH stability. Compared to the structural rigidity of renin, pepsin has more instability associated with the N-terminus, specifically the B/C connector. By contrast, renin exhibits greater C-terminal instability in turn and connector regions. Overall, flexibility differences in connector regions, and amino acid composition, particularly in turn and joint regions of the β-barrel and ψ-loops, likely play defining roles in determining pH stability for renin and pepsin.

A missense variant remote from the active site impairs stability of human phosphoglucomutase 1

Missense variants of human phosphoglucomutase 1 (PGM1) cause the inherited metabolic disease known as PGM1 deficiency. This condition is categorised as both a glycogen storage disease and a congenital disorder of glycosylation. Approximately 20 missense variants of PGM1 are linked to PGM1 deficiency, and biochemical studies have suggested that they fall into two general categories: those affecting the active site and catalytic efficiency, and those that appear to impair protein folding and/or stability. In this study, we characterise a novel variant of Arg422, a residue distal from the active site of PGM1 and the site of a previously identified disease-related variant (Arg422Trp). In prior studies, the R422W variant was found to produce insoluble protein in a recombinant expression system, precluding further in vitro characterisation. Here we investigate an alternative variant of this residue, Arg422Gln, which is amenable to experimental characterisation presumably due to its more conservative physicochemical substitution. Biochemical, crystallographic, and computational studies of R422Q establish that this variant causes only minor changes in catalytic efficiency and 3D structure, but is nonetheless dramatically reduced in stability. Unexpectedly, binding of a substrate analog is found to further destabilise the protein, in contrast to its stabilising effect on wild-type PGM1 and several other missense variants. This work establishes Arg422 as a lynchpin residue for the stability of PGM1 and supports the impairment of protein stability as a pathomechanism for variants that cause PGM1 deficiency. SYNOPSIS: Biochemical and structural studies of a missense variant far from the active site of human PGM1 identify a residue with a key role in enzyme stability.

Improved Modeling of Cation-π and Anion-Ring Interactions Using the Drude Polarizable Empirical Force Field for Proteins

Cation-π interactions are noncovalent interactions between a π-electron system and a positively charged ion that are regarded as a strong noncovalent interaction and are ubiquitous in biological systems. Similarly, though less studied, anion-ring interactions are present in proteins along with in-plane interactions of anions with aromatic rings. As these interactions are between a polarizing ion and a polarizable π system, the accuracy of the treatment of these interactions in molecular dynamics (MD) simulations using additive force fields (FFs) may be limited. In the present work, to allow for a better description of ion-π interactions in proteins in the Drude-2013 protein polarizable FF, we systematically optimized the parameters for these interactions targeting model compound quantum mechanical (QM) interaction energies with atom pair-specific Lennard-Jones parameters along with virtual particles as selected ring centroids introduced to target the QM interaction energies and geometries. Subsequently, MD simulations were performed on a series of protein structures where ion-π pairs occur to evaluate the optimized parameters in the context of the Drude-2013 FF. The resulting FF leads to a significant improvement in reproducing the ion-π pair distances observed in experimental protein structures, as well as a smaller root-mean-square differences and fluctuations of the overall protein structures from experimental structures. Accordingly, the optimized Drude-2013 protein polarizable FF is suggested for use in MD simulations of proteins where cation-π and anion-ring interactions are critical. © 2019 Wiley Periodicals, Inc.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

A review on peptide functionalized graphene derivatives as nanotools for biosensing

Peptides exhibit unique binding behavior with graphene and its derivatives by forming bonds on its edges and planes. This makes them useful for sensing and imaging applications. This review with (155 refs.) summarizes the advances made in the last decade in the field of peptide-GO bioconjugation, and the use of these conjugates in analytical sciences and imaging. The introduction emphasizes the need for understanding the biotic-abiotic interactions in order to construct controllable peptide-functionalized graphitic material-based nanotools. The next section covers covalent and non-covalent interactions between peptide and oxidized graphene derivatives along with a discussion of the adsorption events during interfacing. We then describe applications of peptide-graphene conjugates in bioassays, with subsections on (a) detection of cancer cells, (b) monitoring protease activity, (c) determination of environmental pollutants and (d) determination of pathogenic microorganisms. The concluding section describes the current status of peptide functionalized graphitic bioconjugates and addresses future perspectives. Graphical abstractSchematic representation depicting biosensing applications of peptide functionalized graphene oxide.

Insight into Improved Thermostability of Cold-Adapted Staphylococcal Lipase by Glycine to Cysteine Mutation

Thermostability remains one of the most desirable traits in many lipases. Numerous studies have revealed promising strategies to improve thermostability and random mutagenesis often leads to unexpected yet interesting findings in engineering stability. Previously, the thermostability of C-terminal truncated cold-adapted lipase from Staphylococcus epidermidis AT2 (rT-M386) was markedly enhanced by directed evolution. The newly evolved mutant, G210C, demonstrated an optimal temperature shift from 25 to 45 °C and stability up to 50 °C. Interestingly, a cysteine residue was randomly introduced on the loop connecting the two lids and accounted for the only cysteine found in the lipase. We further investigated the structural and mechanistic insights that could possibly cause the significant temperature shift. Both rT-M386 and G210C were modeled and simulated at 25 °C and 50 °C. The results clearly portrayed the effect of cysteine substitution primarily on the lid stability. Comparative molecular dynamics simulation analysis revealed that G210C exhibited greater stability than the wild-type at high temperature simulation. The compactness of the G210C lipase structure increased at 50 °C and resulted in enhanced rigidity hence stability. This observation is supported by the improved and stronger non-covalent interactions formed in the protein structure. Our findings suggest that the introduction of a single cysteine residue at the lid region of cold-adapted lipase may result in unexpected increased in thermostability, thus this approach could serve as one of the thermostabilization strategies in engineering lipase stability.

Tryptophan usage by Helicobacter pylori differs among strains

Because of its association with severe gastric pathologies, including gastric cancer, Helicobacter pylori has been subject of research for more than 30 years. Its capacity to adapt and survive in the human stomach can be attributed to its genetic flexibility. Its natural competence and its capacity to turn genes on and off allows H. pylori to adapt rapidly to the changing conditions of its host. Because of its genetic variability, it is difficult to establish the uniqueness of each strain obtained from a human host. The methods considered to-date to deliver the best result for differentiation of strains are Rapid Amplification of Polymorphic DNA (RAPD), Multilocus Sequence Typing (MLST) and Whole Genome Sequencing (WGS) analysis. While RAPD analysis is cost-effective, it requires a stable genome for its reliability. MLST and WGS are optimal for strain identification, however, they require analysis of data at the bioinformatics level. Using the StainFree method, which modifies tryptophan residues on proteins using 2, 2, 2, - trichloroethanol (TCE), we observed a strain specific pattern of tryptophan in 1D acrylamide gels. In order to establish the effectiveness of tryptophan fingerprinting for strain identification, we compared the graphic analysis of tryptophan-labelled bands in the gel images with MLST results. Based on this, we find that tryptophan banding patterns can be used as an alternative method for the differentiation of H. pylori strains. Furthermore, investigating the origin for these differences, we found that H. pylori strains alters the number and/or position of tryptophan present in several proteins at the genetic code level, with most exchanges taking place in membrane- and cation-binding proteins, which could be part of a novel response of H. pylori to host adaptation.

Folding of Gα Subunits: Implications for Disease States

G-proteins play a central role in signal transduction by fluctuating between "on" and "off" phases that are determined by a conformational change. cAMP is a secondary messenger whose formation is inhibited or stimulated by activated G_iα1 or G_sα subunit. We used tryptophan fluorescence, UV/vis spectrophotometry, and circular dichroism to probe distinct structural features within active and inactive conformations from wild-type and tryptophan mutants of G_iα1 and G_sα. For all proteins studied, we found that the active conformations were more stable than the inactive conformations, and upon refolding from higher temperatures, activated wild-type subunits recovered significantly more native structure. We also observed that the wild-type subunits partially regained the ability to bind nucleotide. The increased compactness observed upon activation was consistent with the calculated decrease in solvent accessible surface area for wild-type G_iα1. We found that as the temperature increased, G_α subunits, which are known to be rich in α-helices, converted to proteins with increased content of β-sheets and random coil. For active conformations from wild-type and tryptophan mutants of G_iα1, melting temperatures indicated that denaturation starts around hydrophobic tryptophan microenvironments and then radiates toward tyrosine residues at the surface, followed by alteration of the secondary structure. For G_sα, however, disruption of secondary structure preceded unfolding around tyrosine residues. In the active conformations, a π-cation interaction between essential arginine and tryptophan residues, which was characterized by a fluorescence-measured red shift and modeled by molecular dynamics, was also shown to be a contributor to the stability of G_α subunits. The folding properties of G_α subunits reported here are discussed in the context of diseases associated to G-proteins.

Influenza D virus M2 protein exhibits ion channel activity in Xenopus laevis oocytes

BACKGROUND

A new type of influenza virus, known as type D, has recently been identified in cattle and pigs. Influenza D virus infection in cattle is typically asymptomatic; however, its infection in swine can result in clinical disease. Swine can also be infected with all other types of influenza viruses, namely A, B, and C. Consequently, swine can serve as a "mixing vessel" for highly pathogenic influenza viruses, including those with zoonotic potential. Currently, the only antiviral drug available targets influenza M2 protein ion channel is not completely effective. Thus, it is necessary to develop an M2 ion channel blocker capable of suppressing the induction of resistance to the genetic shift. To provide a basis for developing novel ion channel-blocking compounds, we investigated the properties of influenza D virus M2 protein (DM2) as a drug target.

RESULTS

To test the ion channel activity of DM2, the DNA corresponding to DM2 with cMyc-tag conjugated to its carboxyl end was cloned into the shuttle vector pNCB1. The mRNA of the DM2-cMyc gene was synthesized and injected into Xenopus oocytes. The translation products of DM2-cMyc mRNA were confirmed by immunofluorescence and mass spectrometry analyses. The DM2-cMyc mRNA-injected oocytes were subjected to the two-electrode voltage-clamp (TEVC) method, and the induced inward current was observed. The midpoint (Vmid) values in Boltzmann modeling for oocytes injected with DM2-cMyc RNA or a buffer were -152 and -200 mV, respectively. Assuming the same expression level in the Xenopus oocytes, DM2 without tag and influenza C virus M2 protein (CM2) were subjected to the TEVC method. DM2 exhibited ion channel activity under the condition that CM2 ion channel activity was reproduced. The gating voltages represented by Vmid for CM2 and DM2 were -141 and -146 mV, respectively. The reversal potentials observed in ND96 for CM2 and DM2 were -21 and -22 mV, respectively. Compared with intact DM2, DM2 variants with mutation in the YxxxK motif, namely Y72A and K76A DM2, showed lower Vmid values while showing no change in reversal potential.

CONCLUSION

The M2 protein from newly isolated influenza D virus showed ion channel activity similar to that of CM2. The gating voltage was shown to be affected by the YxxxK motif and by the hydrophobicity and bulkiness of the carboxyl end of the molecule.

Cation-π Interactions between Quaternary Ammonium Ions and Amino Acid Aromatic Groups in Aqueous Solution

Cation-π interactions play important roles in the stabilization of protein structures and protein-ligand complexes. They contribute to the binding of quaternary ammonium ligands (mainly RNH₃⁺ and RN(CH₃)₃⁺) to various protein receptors and are likely involved in the blockage of potassium channels by tetramethylammonium (TMA⁺) and tetraethylammonium (TEA⁺). Polarizable molecular models are calibrated for NH₄⁺, TMA⁺, and TEA⁺ interacting with benzene, toluene, 4-methylphenol, and 3-methylindole (representing aromatic amino acid side chains) based on the ab initio MP2(full)/6-311++G(d,p) properties of the complexes. Whereas the gas-phase affinity of the ions with a given aromatic follows the trend NH₄⁺ > TMA⁺ > TEA⁺, molecular dynamics simulations using the polarizable models show a reverse trend in water, likely due to a contribution from the hydrophobic effect. This reversed trend follows the solubility of aromatic hydrocarbons in quaternary ammonium salt solutions, which suggests a role for cation-π interactions in the salting-in of aromatic compounds in solution. Simulations in water show that the complexes possess binding free energies ranging from -1.3 to -3.3 kcal/mol (compared to gas-phase binding energies between -8.5 and -25.0 kcal/mol). Interestingly, whereas the most stable complexes involve TEA⁺ (the largest ion), the most stable solvent-separated complexes involve TMA⁺ (the intermediate-size ion).

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.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.

High stability and biological activity of the copper(II) complexes of alloferon 1 analogues containing tryptophan

Copper(II) complex formation processes between the alloferon 1 (Allo1) (HGVSGHGQHGVHG) analogues where the tryptophan residue is introducing in the place His residue H1W, H6W, H9W and H12W have been studied by potentiometric, UV-visible, CD and EPR spectroscopic, and MS methods. For all analogues of alloferon 1 complex speciation have been obtained for a 1:1 metal-to-ligand molar ratio and 2:1 of H1W because of precipitation at higher (2:1, 3:1 and 4:1) ratios. At physiological pH7.4 and a 1:1 metal-to-ligand molar ratio the tryptophan analogues of alloferon 1 form the CuH_-1L and/or CuH_-2L complexes with the 4N binding mode. The introduction of tryptophan in place of histidine residues changes the distribution diagram of the complexes formed with the change of pH and their stability constants compared to the respective substituted alanine analogues of alloferon 1. The CuH_-1L, CuH_-2L and CuH_-3L complexes of the tryptophan analogues are more stable from 1 to 5 log units in comparison to those of the alanine analogues. This stabilization of the complexes may result from cation(Cu(II))-π and indole/imidazole ring interactions. The induction of apoptosis in vivo, in Tenebrio molitor cells by the ligands and their copper(II) complexes at pH7.4 was studied. The biological results show that copper(II) ions in vivo did not cause any apparent apoptotic features. The most active were the H12W peptide and Cu(II)-H12W complex formed at pH7.4.

Highly Dynamic Anion-Quadrupole Networks in Proteins

The dynamics of anion-quadrupole (or anion-π) interactions formed between negatively charged (Asp/Glu) and aromatic (Phe) side chains are for the first time computationally characterized in RmlC (Protein Data Bank entry 1EP0 ), a homodimeric epimerase. Empirical force field-based molecular dynamics simulations predict anion-quadrupole pairs and triplets (anion-anion-π and anion-π-π) are formed by the protein during the simulated trajectory, which suggests that the anion-quadrupole interactions may provide a significant contribution to the overall stability of the protein, with an average of -1.6 kcal/mol per pair. Some anion-π interactions are predicted to form during the trajectory, extending the number of anion-quadrupole interactions beyond those predicted from crystal structure analysis. At the same time, some anion-π pairs observed in the crystal structure exhibit marginal stability. Overall, most anion-π interactions alternate between an "on" state, with significantly stabilizing energies, and an "off" state, with marginal or null stabilizing energies. The way proteins possibly compensate for transient loss of anion-quadrupole interactions is characterized in the RmlC aspartate 84-phenylalanine 112 anion-quadrupole pair observed in the crystal structure. A double-mutant cycle analysis of the thermal stability suggests a possible loss of anion-π interactions compensated by variations of hydration of the residues and formation of compensating electrostatic interactions. These results suggest that near-planar anion-quadrupole pairs can exist, sometimes transiently, which may play a role in maintaining the structural stability and function of the protein, in an otherwise very dynamic interplay of a nonbonded interaction network as well as solvent effects.

Molecular Determinants of Mutant Phenotypes, Inferred from Saturation Mutagenesis Data

Understanding how mutations affect protein activity and organismal fitness is a major challenge. We used saturation mutagenesis combined with deep sequencing to determine mutational sensitivity scores for 1,664 single-site mutants of the 101 residue Escherichia coli cytotoxin, CcdB at seven different expression levels. Active-site residues could be distinguished from buried ones, based on their differential tolerance to aliphatic and charged amino acid substitutions. At nonactive-site positions, the average mutational tolerance correlated better with depth from the protein surface than with accessibility. Remarkably, similar results were observed for two other small proteins, PDZ domain (PSD95^pdz3) and IgG-binding domain of protein G (GB1). Mutational sensitivity data obtained with CcdB were used to derive a procedure for predicting functional effects of mutations. Results compared favorably with those of two widely used computational predictors. In vitro characterization of 80 single, nonactive-site mutants of CcdB showed that activity in vivo correlates moderately with thermal stability and solubility. The inability to refold reversibly, as well as a decreased folding rate in vitro, is associated with decreased activity in vivo. Upon probing the effect of modulating expression of various proteases and chaperones on mutant phenotypes, most deleterious mutants showed an increased in vivo activity and solubility only upon over-expression of either Trigger factor or SecB ATP-independent chaperones. Collectively, these data suggest that folding kinetics rather than protein stability is the primary determinant of activity in vivo. This study enhances our understanding of how mutations affect phenotype, as well as the ability to predict fitness effects of point mutations.

Redox Activity of Copper(II) Complexes with NSFRY Pentapeptide and Its Analogues

The influence of cation-π interactions on the electrochemical properties of copper(II) complexes with synthesized pentapeptide C-terminal fragment of Atrial Natriuretic Factor (ANF) hormone was studied in this work. Molecular modeling performed for Cu(II)-NSFRY-NH₂ complex indicated that the cation-π interactions between Tyr and Cu(II), and also between Phe-Arg led to specific conformation defined as peptide box, in which the metal cation is isolated from the solvent by peptide ligand. Voltammetry experiments enabled to compare the redox properties and stability of copper(II) complexes with NSFRY-NH₂ and its analogues (namely: NSFRA-NH₂, NSFRF-NH_2, NSAAY-NH₂, NSAAA-NH₂, AAAAA-NH₂) as well as to evaluate the contribution of individual amino acid residues to these properties. The obtained results led to the conclusion, that cation-π interactions play a crucial role in the effective stabilization of copper(II) complexes with the fragments of ANF peptide hormone and therefore could control the redox processes in other metalloproteins.

A Miniature Protein Stabilized by a Cation-π Interaction Network

The design of folded miniature proteins is predicated on establishing noncovalent interactions that direct the self-assembly of discrete thermostable tertiary structures. In this work, we describe how a network of cation-π interactions present in proteins containing "WSXWS motifs" can be emulated to stabilize the core of a miniature protein. This 19-residue protein sequence recapitulates a set of interdigitated arginine and tryptophan residues that stabilize a distinctive β-strand:loop:PPII-helix topology. Validation of the compact fold determined by NMR was carried out by mutagenesis of the cation-π network and by comparison to the corresponding disulfide-bridged structure. These results support the involvement of a coordinated set of cation-π interactions that stabilize the tertiary structure.

Impact of cation-π interactions on the cell voltage of carbon nanotube-based Li batteries

Carbon nanotube (CNT)-based Li batteries have attracted wide attention because of their high capacity, high cyclability and high energy density and are believed to be one of the most promising electrochemical energy storage systems. In CNT-based Li batteries, the main interaction between the Li(+) ions and the CNT is the cation-π interaction. However, up to now, it is still not clear how this interaction affects the storage characteristics of CNT-based Li batteries. Here, using density functional theory (DFT) calculations, we report a highly favorable impact of cation-π interactions on the cell voltage of CNT-based Li batteries. Considering both Li(+)-π interaction and Li-π interaction, we show that cell voltage enhances with the increase of the CNT diameter. In addition, when the Li(+) ion adsorbs on the external wall, the cell voltage is larger than that when it adsorbs on the internal wall. This suggests that CNTs with a large diameter and a low array density are more advantageous to enhance storage performance of CNT-based Li batteries. Compared with Li(+) ions on the (4,4) CNT internal wall, the cell voltage of Li(+) on the (10,10) CNT external wall is 0.55 V higher, which indicates an improvement of about 38%. These results will be helpful for the design of more efficient CNT-based Li batteries.

Behaviour of cation–pi interaction in presence of external electric field

External electric field effects cation–π interaction.

Induced variations in brassinosteroid genes define barley height and sturdiness, and expand the green revolution genetic toolkit

Reduced plant height and culm robustness are quantitative characteristics important for assuring cereal crop yield and quality under adverse weather conditions. A very limited number of short-culm mutant alleles were introduced into commercial crop cultivars during the Green Revolution. We identified phenotypic traits, including sturdy culm, specific for deficiencies in brassinosteroid biosynthesis and signaling in semidwarf mutants of barley (Hordeum vulgare). This set of characteristic traits was explored to perform a phenotypic screen of near-isogenic short-culm mutant lines from the brachytic, breviaristatum, dense spike, erectoides, semibrachytic, semidwarf, and slender dwarf mutant groups. In silico mapping of brassinosteroid-related genes in the barley genome in combination with sequencing of barley mutant lines assigned more than 20 historic mutants to three brassinosteroid-biosynthesis genes (BRASSINOSTEROID-6-OXIDASE, CONSTITUTIVE PHOTOMORPHOGENIC DWARF, and DIMINUTO) and one brassinosteroid-signaling gene (BRASSINOSTEROID-INSENSITIVE1 [HvBRI1]). Analyses of F2 and M2 populations, allelic crosses, and modeling of nonsynonymous amino acid exchanges in protein crystal structures gave a further understanding of the control of barley plant architecture and sturdiness by brassinosteroid-related genes. Alternatives to the widely used but highly temperature-sensitive uzu1.a allele of HvBRI1 represent potential genetic building blocks for breeding strategies with sturdy and climate-tolerant barley cultivars.

Structure-based inhibition of protein-protein interactions

Protein-protein interactions (PPIs) are emerging as attractive targets for drug design because of their central role in directing normal and aberrant cellular functions. These interactions were once considered "undruggable" because their large and dynamic interfaces make small molecule inhibitor design challenging. However, landmark advances in computational analysis, fragment screening and molecular design have enabled development of a host of promising strategies to address the fundamental molecular recognition challenge. An attractive approach for targeting PPIs involves mimicry of protein domains that are critical for complex formation. This approach recognizes that protein subdomains or protein secondary structures are often present at interfaces and serve as organized scaffolds for the presentation of side chain groups that engage the partner protein(s). Design of protein domain mimetics is in principle rather straightforward but is enabled by a host of computational strategies that provide predictions of important residues that should be mimicked. Herein we describe a workflow proceeding from interaction network analysis, to modeling a complex structure, to identifying a high-affinity sub-structure, to developing interaction inhibitors. We apply the design procedure to peptidomimetic inhibitors of Ras-mediated signaling.

Fosfomycin induced structural change in fosfomycin resistance kinases FomA: molecular dynamics and molecular docking studies

Fosfomycin resistance kinases FomA, one of the key enzymes responsible for bacterial resistances to fosfomycin, has gained much attention recently due to the raising public concern for multi-drug resistant bacteria. Using molecular docking followed by molecular dynamics simulations, our group illustrated the process of fosfomycin induced conformational change of FomA. The detailed roles of the catalytic residues (Lys18, His58 and Thr210) during the formation of the enzyme-substrate complex were shown in our research. The organization functions of Gly53, Gly54, Ile61 and Leu75 were also highlighted. Furthermore, the cation-π interaction between Arg62 and Trp207 was observed and speculated to play an auxiliary role in the conformation change process of the enzyme. This detailed molecular level illustration of the formation of FomA·ATP·Mg·Fosfomycin complex could provide insight for both anti-biotic discovery and improvement of fosfomycin in the future.

The N-terminal strand modulates immunoglobulin light chain fibrillogenesis

It has been suggested that the N-terminal strand of the light chain variable domain (V(L)) protects the molecule from aggregation by hindering spurious intermolecular contacts. We evaluated the impact of mutations in the N-terminal strand on the thermodynamic stability and kinetic of fibrillogenesis of the V(L) protein 6aJL2. Mutations in this strand destabilized the protein in a position-dependent manner, accelerating the fibrillogenesis by shortening the lag time; an effect that correlated with the extent of destabilization. In contrast, the effect on the kinetics of fibril elongation, as assessed in seeding experiments was of different nature, as it was not directly dependant on the degree of destabilization. This finding suggests different factors drive the nucleation-dependent and elongation phases of light chain fibrillogenesis. Finally, taking advantage of the dependence of the Trp fluorescence upon environment, four single Trp substitutions were made in the N-terminal strand, and changes in solvent exposure during aggregation were evaluated by acrylamide-quenching. The results suggest that the N-terminal strand is buried in the fibrillar state of 6aJL2 protein. This finding suggest a possible explanation for the modulating effect exerted by the mutations in this strand on the aggregation behavior of 6aJL2 protein.

Emerging role of N- and C-terminal interactions in stabilizing (β/α)8 fold with special emphasis on Family 10 xylanases

Xylanases belong to an important class of industrial enzymes. Various xylanases have been purified and characterized from a plethora of organisms including bacteria, marine algae, plants, protozoans, insects, snails and crustaceans. Depending on the source, the enzymatic activity of xylanases varies considerably under various physico-chemical conditions such as temperature, pH, high salt and in the presence of proteases. Family 10 or glycosyl hydrolase 10 (GH10) xylanases are one of the well characterized and thoroughly studied classes of industrial enzymes. The TIM-barrel fold structure which is ubiquitous in nature is one of the characteristics of family 10 xylanases. Family 10 xylanases have been used as a “model system” due to their TIM-barrel fold to dissect and understand protein stability under various conditions. A better understanding of structure-stability-function relationships of family 10 xylanases allows one to apply these governing molecular rules to engineer other TIM-barrel fold proteins to improve their stability and retain function(s) under adverse conditions. In this review, we discuss the implications of N-and C-terminal interactions, observed in family 10 xylanases on protein stability under extreme conditions. The role of metal binding and aromatic clusters in protein stability is also discussed. Studying and understanding family 10 xylanase structure and function, can contribute to our protein engineering knowledge.