Forces stabilizing proteins

Relationships between conformational and vibrational features of tryptophan characteristic Raman markers

Tryptophan (Trp) residue provides characteristic vibrational markers to the middle wavenumber spectral region of the Raman spectra recorded from peptides and proteins. In this report, we were particularly interested in eight Trp Raman markers, referred to as Wi (i = 1,…,8). All responsible for pronounced Raman lines, these markers originate from indole moiety, a bicyclic conjugated segment involved in the Trp structure. Numerous investigations have previously attempted to relate the variations observed in the spectral features of these markers to the environmental changes of Trp residues. To emphasize the most important points we can mention (i) the variations in the Raman profile of W4 (∼1360 cm-1) and W5 (∼1340 cm-1), frequently observed as a doublet with variable intensity ratio. These two markers were thought to result from a Fermi-resonance effect between certain planar and nonplanar modes; (ii) the changes observed in the wavenumbers and relative intensities of W4, W7 (∼880 cm-1) and W8 (∼760 cm-1) were supposed to be related to the accessibility of Trp to surrounding water molecules; and (iii) the wavenumber fluctuations of W3 (∼1550 cm-1), taken as a Trp side chain orientational marker. However, some ambiguities still exist regarding the interpretation of these markers, needing further clarification. Herein, upon a joint experimental and theoretical analysis based on a multiconformational approach, attention was paid to the relationships between structural and vibrational features of three indole-containing compounds with increasing structural complexity, i.e., skatole (3-methylindole), tryptophan, and tripeptide Gly-Trp-Gly. This study clearly shows that the existing assignments given to certain Trp Raman markers should be reconsidered, especially those based on the Fermi-resonance origin of W4-W5 (∼1360-1340 cm-1) doublet, as well as the purely environmental dependence of W7 and W8 markers.

Influence of different pH milieu on the structure and function of human Aurora kinase B protein (AURK-B): Amalgamation of both spectroscopic and computational approach

Aurora kinase B (AURK-B) is a serine/threonine kinase protein that plays an essential role in chromosomal separation during the cell cycle event. AURK-B is highly expressed in various types of cancer such as human seminoma, thyroid carcinoma, non-small cell lung carcinoma (NSCLC), oral carcinoma, and gastric cancer. Hence, it is a potential therapeutic target in the treatment of various cancers. The structure of AURK-B in complex with one of its substrate inner centromeric protein (INCENP) is present, but the structural and functional characterization of native AURK-B at different pH environment is still unexplored.This study determines the effect of different pH milieu on the structure and function of AURK-B protein wherein the influence of pH on the protein conformation was probed using Circular dichroism (CD) and fluorescence spectroscopy. The structural studies were further combined with functional activity assay to observe the change in kinase activity at various pH milieu (2.0-11.0). CD and fluorescence spectroscopy experiments dictate that at high acidic conditions (pH 2.0 - 5.0), the secondary and tertiary structures of AURK-B become distorted, leading to diminished activity. The protein, however, was observed to stabilize towards pH 7.0 - 8.0 with minimal structure alteration over the basic pH range (pH 9.0 -11.0). The measured spectroscopic structural features were found to be in-line with obtained experimental kinase activity assays. Further, in-vitro experiments indicate that the enzyme is maximally active at pH 8.0. More ordered conformation and compact structure was observed at this pH (pH 8.0) as compared to other pH values through molecular dynamics simulation studies (MDS). As AURK-B localizes itself in the intracellular compartment, this study may provide a clue about the role of different pH environments in enhancing cancer growth, proliferation, and invasion.

Protein Thermodynamic Properties, Crystallisation, and the Hofmeister Series

The Hofmeister series is a series of ions ordered according to their ability to precipitate proteins. It has also been linked to a host of (bio)chemical phenomena. Several attempts over the years to correlate the series to the varying success of different salts in crystallising proteins have been largely inconclusive. A correlation, based on published data and crystallisation conditions for several proteins, is proposed here between some thermodynamic properties of proteins and the position in the Hofmeister series of the salts from which they preferentially crystallise. Namely, a high ratio between the entropic or enthalpic protein-solvent interactions contribution to thermodynamic stability and the total thermodynamic stability of a given protein, indicate the protein's high propensity to crystallise in solutions of highly kosmotropic salts. Low such ratios on the other hand, indicate that chaotropic salts can be equally successful, i. e. that the protein in question is rather indifferent to the Hofmeister character of the salt. Testing various model proteins for crystallisation against screens containing salts found at different points on the Hofmeister series, as well as further bibliographic analysis, have yielded results that appear to largely corroborate this hypothesis. These conclusions may conceivably be used as a crystallisation predictive tool.

Interpretation of Hydrogen/Deuterium Exchange Mass Spectrometry

This paper sheds light on the meaning of hydrogen/deuterium exchange-mass spectrometry (HDX-MS) data. HDX-MS data provide not structural information but dynamic information on an analyte protein. First, the reaction mechanism of backbone amide HDX reaction is considered and the correlation between the parameters from an X-ray crystal structure and the protection factors of HDX reactions of cytochrome c is evaluated. The presence of H-bonds in a protein structure has a strong influence on HDX rates which represent protein dynamics, while the solvent accessibility only weakly affects the HDX rates. Second, the energy diagrams of the HDX reaction at each residue in the presence and absence of perturbation are described. Whereas the free energy change upon mutation can be directly measured by the HDX rates, the free energy change upon ligand binding may be complicated due to the presence of unbound analyte protein in the protein-ligand mixture. Third, the meanings of HDX and other biophysical techniques are explained using a hypothetical protein folding well. The shape of the protein folding well describes the protein dynamics and provides Boltzmann distribution of open and closed states which yield HDX protection factors, while a protein's crystal structure represents a snapshot near the bottom of the well. All biophysical data should be consistent yet provide different information because they monitor different parts of the same protein folding well.

Inquiry lipaseoring the mechanism of pancreatic lipase inhibition by isovitexin based on multispectral method and enzyme inhibition assay

Isovitexin is a main natural flavonoid component in various plants. Currently, the inhibitory effect of isovitexin on pancreatic lipase (PL) and its mechanism have not been elucidated yet. In the present study, we investigated the inhibitory effect of isovitexin on PL, as well as its interaction mechanism, using enzyme inhibition methods, spectroscopic analysis, and molecular simulations. Results showed that isovitexin possessed significant PL inhibitory activity, with IC50 values of 0.26 ± 0.02 mM. The interaction between isovitexin and PL was dominated by static quenching, and mainly through hydrogen bonding and hydrophobic interaction forces. Analysis of fluorescence spectroscopy confirmed that isovitexin binding altered the conformation of the PL. Circular dichroism (CD) spectrum indicated that isovitexin altered the secondary structure of PL by decreasing the α-helix content and increasing the β-fold content. Molecular simulations further characterize the conformational changes produced by the interaction between isovitexin with PL. The performed study may provide a new insight into the inhibitory mechanism of isovitexin as a novel PL inhibitor.

Identification of PS1/gamma-secretase and glutamate transporter GLT-1 interaction sites

The recently discovered interaction between Presenilin 1 (PS1), a catalytic subunit of γ-secretase responsible for generating amyloid-β peptides, and GLT-1, a major glutamate transporter in the brain (EAAT2), provides a mechanistic link between these two key factors involved in Alzheimer's disease (AD) pathology. Modulating this interaction can be crucial to understand the consequence of such crosstalk in AD context and beyond. However, the interaction sites between these two proteins are unknown. Herein, we utilized an alanine scanning approach coupled with FRET-based fluorescence lifetime imaging microscopy to identify the interaction sites between PS1 and GLT-1 in their native environment within intact cells. We found that GLT-1 residues at position 276 to 279 (TM5) and PS1 residues at position 249 to 252 (TM6) are crucial for GLT-1-PS1 interaction. These results have been cross validated using AlphaFold Multimer prediction. To further investigate whether this interaction of endogenously expressed GLT-1 and PS1 can be prevented in primary neurons, we designed PS1/GLT-1 cell-permeable peptides (CPPs) targeting the PS1 or GLT-1 binding site. We used HIV TAT domain to allow for cell penetration which was assayed in neurons. First, we assessed the toxicity and penetration of CPPs by confocal microscopy. Next, to ensure the efficiency of CPPs, we monitored the modulation of GLT-1-PS1 interaction in intact neurons by fluorescence lifetime imaging microscopy. We saw significantly less interaction between PS1 and GLT-1 with both CPPs. Our study establishes a new tool to study the functional aspect of GLT-1-PS1 interaction and its relevance in normal physiology and AD models.

Polysaccharides screening for pulmonary mucus penetration by molecular dynamics simulation and in vitro verification

Mucus penetration is one of the physiologic barriers of inhalation and nanocarriers can effectively facilitate the permeation of drugs. The interactions between the nanocarriers and mucin are crucial for penetration across the mucus layer on the respiratory tract. In this study, we proposed a molecular dynamics (MD) simulation method for the screening of polysaccharides that acted as the surface modification materials for inhalable nano-preparations to facilitate mucus penetration. MD revealed all-atom interactions between the monomers of polysaccharides, including dextran (DEX)/hyaluronic acid (HA)/carboxymethyl chitosan (CMCS) and the human mucin protein MUC5AC (hMUC5AC). The obtained data showed that DEX formed stronger non-covalent bonds with hMUC5AC compared to HA and CMCS, which suggested that HA and CMCS had better mucus permeability than DEX. For the in vitro verification, HA/CMCS-coated liposomes and DEX/PEG-inserted liposomes were prepared. The results of mucin interactions and mucus penetration studies confirmed that HA and CMCS possessed the weakest interactions with mucin and facilitated the mucus penetration, which was in consistent with the data from MD simulation. This work may shed light on the MD simulation-based screening of surface modification materials for inhalable nano-preparations to facilitate mucus penetration.

Multidomain chimeric enzymes as a promising alternative for biocatalysts improvement: a minireview

Searching for new and better biocatalysts is an area of study in constant development. In nature, mechanisms generally occurring in evolution, such as genetic duplication, recombination, and natural selection processes, produce various enzymes with different architectures and properties. The recombination of genes that code proteins produces multidomain chimeric enzymes that contain two or more domains that sometimes enhance their catalytic properties. Protein engineering has mimicked this process to enhance catalytic activity and the global stability of enzymes, searching for new and better biocatalysts. Here, we present and discuss examples from both natural and synthetic multidomain chimeric enzymes and how additional domains heighten their stability and catalytic activity. Moreover, we also describe progress in developing new biocatalysts using synthetic fusion enzymes and revise some methodological strategies to improve their biological fitness.

Stepwise introduction of stabilizing mutations reveals nonlinear additive effects in de novo TIM barrels

Over the past decades, the TIM-barrel fold has served as a model system for the exploration of how changes in protein sequences affect their structural, stability, and functional characteristics, and moreover, how this information can be leveraged to design proteins from the ground up. After numerous attempts to design de novo proteins with this specific fold, sTIM11 was the first validated de novo design of an idealized four-fold symmetric TIM barrel. Subsequent efforts to enhance the stability of this initial design resulted in the development of DeNovoTIMs, a family of de novo TIM barrels with various stabilizing mutations. In this study, we present an investigation into the biophysical and thermodynamic effects upon introducing a varying number of stabilizing mutations per quarter along the sequence of a four-fold symmetric TIM barrel. We compared the base design DeNovoTIM0 without any stabilizing mutations with variants containing mutations in one, two, three, and all four quarters-designated TIM1q, TIM2q, TIM3q, and DeNovoTIM6, respectively. This analysis revealed a stepwise and nonlinear change in the thermodynamic properties that correlated with the number of mutated quarters, suggesting positive nonadditive effects. To shed light on the significance of the location of stabilized quarters, we engineered two variants of TIM2q which contain the same number of mutations but positioned in different quarter locations. Characterization of these TIM2q variants revealed that the mutations exhibit varying effects on the overall protein stability, contingent upon the specific region in which they are introduced. These findings emphasize that the amount and location of stabilized interfaces among the four quarters play a crucial role in shaping the conformational stability of these four-fold symmetric TIM barrels. Analysis of de novo proteins, as described in this study, enhances our understanding of how sequence variations can finely modulate stability in both naturally occurring and computationally designed proteins.

Living in trinity of extremes: Genomic and proteomic signatures of halophilic, thermophilic, and pH adaptation

Since nucleic acids and proteins of unicellular prokaryotes are directly exposed to extreme environmental conditions, it is possible to explore the genomic-proteomic compositional determinants of molecular mechanisms of adaptation developed by them in response to harsh environmental conditions. Using a wealth of currently available complete genomes/proteomes we were able to explore signatures of adaptation to three environmental factors, pH, salinity, and temperature, observing major trends in compositions of their nucleic acids and proteins. We derived predictors of thermostability, halophilic, and pH adaptations and complemented them by the principal components analysis. We observed a clear difference between thermophilic and salinity/pH adaptations, whereas latter invoke seemingly overlapping mechanisms. The genome-proteome compositional trade-off reveals an intricate balance between the work of base paring and base stacking in stabilization of coding DNA and r/tRNAs, and, at the same time, universal requirements for the stability and foldability of proteins regardless of the nucleotide biases. Nevertheless, we still found hidden fingerprints of ancient evolutionary connections between the nucleotide and amino acid compositions indicating their emergence, mutual evolution, and adjustment. The evolutionary perspective on the adaptation mechanisms is further studied here by means of the comparative analysis of genomic/proteomic traits of archaeal and bacterial species. The overall picture of genomic/proteomic signals of adaptation obtained here provides a foundation for future engineering and design of functional biomolecules resistant to harsh environments.

Protein Misfolding in Pregnancy: Current Insights, Potential Mechanisms, and Implications for the Pathogenesis of Preeclampsia

Protein misfolding disorders are a group of diseases characterized by supra-physiologic accumulation and aggregation of pathogenic proteoforms resulting from improper protein folding and/or insufficiency in clearance mechanisms. Although these processes have been historically linked to neurodegenerative disorders, such as Alzheimer's disease, evidence linking protein misfolding to other pathologies continues to emerge. Indeed, the deposition of toxic protein aggregates in the form of oligomers or large amyloid fibrils has been linked to type 2 diabetes, various types of cancer, and, in more recent years, to preeclampsia, a life-threatening pregnancy-specific disorder. While extensive physiological mechanisms are in place to maintain proteostasis, processes, such as aging, genetic factors, or environmental stress in the form of hypoxia, nutrient deprivation or xenobiotic exposures can induce failure in these systems. As such, pregnancy, a natural physical state that already places the maternal body under significant physiological stress, creates an environment with a lower threshold for aberrant aggregation. In this review, we set out to discuss current evidence of protein misfolding in pregnancy and potential mechanisms supporting a key role for this process in preeclampsia pathogenesis. Improving our understanding of this emerging pathophysiological process in preeclampsia can lead to vital discoveries that can be harnessed to create better diagnoses and treatment modalities for the disorder.

Structural characterization of solvent-based food preparation of jellyfish

Jellyfish as a potential sustainable food material has recently gained increasing interest. However, with their soft gel-like texture and easy spoilage, it remains challenging to achieve desirable edible structures from jellyfish. The culinary preparation of jellyfish is a complex process and extends beyond conventional cooking methods. In this study, we investigate the transformation of jellyfish into crispy-like structures by manipulating their microstructural and mechanical properties through a solvent-based preparation. The study focuses on the use of "poor solvents", namely ethanol and acetone, and employs rheology measurements and quantitative microscopy techniques to analyze the effects of these solvents on the mechanical properties and microstructure of jellyfish. Our findings reveal that both ethanol and acetone lead to a significant increase in jellyfish hardness and deswelling. Notably, a micro-scale network is formed within the jellyfish matrix, and this network is then mechanically reinforced before a crispy-like texture can be obtained. Our study points to solvent polarity as also being a crucial factor for creating these effects and determines an upper polarity limit in the range of 12.2-12.9 MPa1/2 for added solvents, corresponding to approximately 60% of added ethanol or 70% of added acetone. Our study highlights that solvent-based preparation serves as a "reverse cooking" technique, where mechanical modification rather than traditional softening mechanisms are employed to stabilize and strengthen the microstructures and fibers of jellyfish. By elucidating the underlying mechanisms of solvent-induced stabilization, our findings may facilitate the development of innovative and sustainable culinary practices, paving the way for broader applications of jellyfish and other soft edible materials in the gastronomic landscape.

Computational and Biophysical Approaches to Identify Cell Wall-Associated Modulators in Salmonella enterica serovar Typhi

An increase in the number of antibiotic-resistant bacterial pathogens, in recent times, has posed a great challenge for treating the affected patients. This has paved the way for the development and design of antibiotics against the previously less explored newer targets. Among these, peptidoglycan (PG) biosynthesis serves as a promising target for the design and development of novel drugs. The peptidoglycan cell wall synthesis in bacteria is essential for its viability. The enzyme class, Mur ligases, plays a key role in PG biosynthesis. Therefore, compounds with the ability to inhibit these enzymes (Mur ligase) can serve as potential candidates for developing small modulators. The enzyme, UDP-N-acetyl pyruvyl-glucosamine reductase (MurB), is essential for PG biosynthesis, a crucial part of the bacterial cell wall. The development of novel drugs to treat infections may thus focus on inhibiting MurB function. Understanding the mechanism of action of Mur B is central to developing efficient inhibitors. For the treatment of S. typhi infections, it is also critical to find therapeutic drugs that specifically target MurB. The enzyme Mur B from Salmonella enterica serovar Typhi (stMurB) was expressed and purified for biophysical characterization to gauge the molecular interactions and estimate thermodynamic stability, for determining attributes for possible therapeutic intervention. The thermal melting profile of MurB was monitored by circular dichroism (CD) and validated by performing differential scanning calorimetry (DSC). An in silico virtual screening of various natural inhibitors was conducted with modelled stMurB structure. The three top hits (quercetin, berberine, and scopoletin) obtained from in silico screening were validated for complex stability through molecular dynamics (MD) simulation. Further, fluorescence binding studies were undertaken for the selected natural inhibitors with stMurB alone and with its NADPH-bound form. The natural inhibitors, scopoletin and berberine, displayed lesser binding to stMurB compared to quercetin. Also, a stronger binding affinity was exhibited between quercetin and stMurB compared to NADPH and stMurB. Based on the above two findings, quercetin can be developed as an inhibitor of stMurB enzyme.

Stabilization challenges and aggregation in protein-based therapeutics in the pharmaceutical industry

Protein-based therapeutics have revolutionized the pharmaceutical industry and become vital components in the development of future therapeutics. They offer several advantages over traditional small molecule drugs, including high affinity, potency and specificity, while demonstrating low toxicity and minimal adverse effects. However, the development and manufacturing processes of protein-based therapeutics presents challenges related to protein folding, purification, stability and immunogenicity that should be addressed. These proteins, like other biological molecules, are prone to chemical and physical instabilities. The stability of protein-based drugs throughout the entire manufacturing, storage and delivery process is essential. The occurrence of structural instability resulting from misfolding, unfolding, and modifications, as well as aggregation, poses a significant risk to the efficacy of these drugs, overshadowing their promising attributes. Gaining insight into structural alterations caused by aggregation and their impact on immunogenicity is vital for the advancement and refinement of protein therapeutics. Hence, in this review, we have discussed some features of protein aggregation during production, formulation and storage as well as stabilization strategies in protein engineering and computational methods to prevent aggregation.

The Monomer-Dimer Equilibrium of Triscatechol Titanium(IV)-Based Hierarchical Helicates as a Tool for the Development of Molecular Balances and Molecular Switches

ConspectusHierarchical helicates are formed by noncovalent connection of two or more monomeric metal complex units, e.g., by bridging metal cations. A unique kind of hierarchical helicate is obtained from 3-carbonyl substituted catechol ligands with titanium(IV) ions in the presence of lithium cations. This kind of supramolecular complex shows in solution a "monomer-dimer" equilibrium. There are different possibilities (solvent, countercation, substituents at carbonyl unit, etc.) to shift this equilibrium to either the monomer or the dimer side. Thus, the lithium-bridged catecholate-based hierarchical helicates resemble a molecular switch. In this Account, different aspects are discussed of how this unique behavior of the dimeric titanium catecholates can be used for application.Thorough investigation of the energetics of the monomer-dimer equilibrium leads to a deeper understanding of the thermodynamic and kinetic effects of the dimerization (or dissociation) process. In this context, even weak interaction of substituents in the periphery of the complexes can be observed. Hereby on the one hand, solvent effects have an important influence and can be easily evaluated. The thorough understanding of the behavior of the monomer-dimer equilibrium allows one to develop some novel applications. In this respect, the use of the hierarchical helicate-based switch as a platform for reaction control and catalysis is described. Decent enantioselectivities up to ee = 58% can be found in Diels-Alder reactions in the periphery of the dimers, while switching to the monomer as a reaction platform still allows the cycloaddition reaction but turns the selectivity off. Additionally, it is described that catalytically important units can be introduced and hydrogenation reactions as well as Michael-type reactions are catalyzed at the helicates.Covalent connection of two catechol ester units leads to classical helicates. Depending on the alkaline metal cation, those can be switched from a compressed to an expanded form or vice versa. Hereby the monomer-dimer equilibrium is transformed into a structural switch. The switching process can be initiated by removal or addition of lithium cations (e.g., by addition of [2.1.1]cryptand). Alternative switching possibilities are based in the case of glycol bridged helicates on cation translocation isomerism and with thioester derivatives it occurs spontaneously in DMSO. Introduction of chiral tethers results in a three state switch allowing expansion/compression as well as switching of the helicity.

Quantifying Amide-Aromatic Interactions at Molecular and Atomic Levels: Experimentally Determined Enthalpic and Entropic Contributions to Interactions of Amide sp2O, N, C and sp3C Unified Atoms with Naphthalene sp2C Atoms in Water

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp2 atoms of amides, aromatics, and other groups occur in protein self-assembly processes including folding, oligomerization, and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp2O with amide sp2N unified atoms (presumably C═O···H-N hydrogen bonds) and amide/aromatic sp2C (lone pair π, n-π*) are particularly favorable. Sp3C-sp3C (hydrophobic), sp3C-sp2C (hydrophobic, CH-π), sp2C-sp2C (hydrophobic, π-π), and sp3C-sp2N interactions are favorable, sp2C-sp2N interactions are neutral, while sp2O-sp2O and sp2N-sp2N self-interactions and sp2O-sp3C interactions are unfavorable. Here, from determinations of favorable effects of 14 amides on naphthalene solubility at 10, 25, and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp2O, sp2N, sp2C, and sp3C unified atoms with aromatic sp2C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp2O-aromatic sp2C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g., lone pair-π), while amide sp3C- and sp2C-aromatic sp2C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp2 atoms in protein processes.

Identification of hits as anti-obesity agents against human pancreatic lipase via docking, drug-likeness, in-silico ADME(T), pharmacophore, DFT, molecular dynamics, and MM/PB(GB)SA analysis

Obesity, characterized by excessive fat accumulation, is a major health concern. Inhibition of human pancreatic lipase, an enzyme involved in fat digestion, offers a potential strategy for weight loss and obesity treatment. This study aimed to identify polyphenols capable of forming stable complexes with human pancreatic lipase to block its activity. Molecular docking, density functional theory (DFT), molecular dynamics (MD) simulations, and MMPBGBSA calculations were employed to evaluate ligand binding, stability, and energy profiles. Pharmacophore modeling was also performed to identify key structural features for effective inhibition. Virtual screening identified ZINC000015120539, ZINC000000899200, ZINC000001531702, and ZINC000013340267 as potential candidates, exhibiting favorable binding and stable interactions over 100 ns MD simulations. These findings provide insights into the inhibitory potential of selected polyphenols on human pancreatic lipase and support further experimental investigations for obesity treatment.Communicated by Ramaswamy H. Sarma.

Investigation of thermal stability characteristic in family A DNA polymerase - A theoretical study

DNA polymerases create complementary DNA strands in living cells and are crucial to genome transmission and maintenance. These enzymes possess similar human right-handed folds which contain thumb, fingers, and palm subdomains and contribute to polymerization activities. These enzymes are classified into seven evolutionary families, A, B, C, D, X, Y, and RT, based on amino acid sequence analysis and biochemical characteristics. Family A DNA polymerases exist in an extended range of organisms including mesophilic, thermophilic, and hyper-thermophilic bacteria, participate in DNA replication and repair, and have a broad application in molecular biology and biotechnology. In this study, we attempted to detect factors that play a role in the thermostability properties of this family member despite their remarkable similarities in structure and function. For this purpose, similarities and differences in amino acid sequences, structure, and dynamics of these enzymes have been inspected. Our results demonstrated that thermophilic and hyper-thermophilic enzymes have more charged, aromatic, and polar residues than mesophilic ones and consequently show further electrostatic and cation-pi interactions. In addition, in thermophilic enzymes, aliphatic residues tend to position in buried states more than mesophilic enzymes. These residues within their aliphatic parts increase hydrophobic core packing and therefore enhance the thermostability of these enzymes. Furthermore, a decrease in thermophilic cavities volumes assists in the protein compactness enhancement. Moreover, molecular dynamic simulation results revealed that increasing temperature impacts mesophilic enzymes further than thermophilic ones that reflect on polar and aliphatic residues surface area and hydrogen bonds changes.

Structural Insights into the Stability and Recognition Mechanism of the Antiquinalphos Nanobody for the Detection of Quinalphos in Foods

Nanobodies (Nbs) have great potential in immunoassays due to their exceptional physicochemical properties. With the immortal nature of Nbs and the ability to manipulate their structures using protein engineering, it will become increasingly valuable to understand what structural features of Nbs drive high stability, affinity, and selectivity. Here, we employed an anti-quinalphos Nb as a model to illustrate the structural basis of Nbs' distinctive physicochemical properties and the recognition mechanism. The results indicated that the Nb-11A-ligand complexes exhibit a "tunnel" binding mode formed by CDR1, CDR2, and FR3. The orientation and hydrophobicity of small ligands are the primary determinants of their diverse affinities to Nb-11A. In addition, the primary factors contributing to Nb-11A's limited stability at high temperatures and in organic solvents are the rearrangement of the hydrogen bonding network and the enlargement of the binding cavity. Importantly, Ala 97 and Ala 34 at the active cavity's bottom and Arg 29 and Leu 73 at its entrance play vital roles in hapten recognition, which were further confirmed by mutant Nb-F3. Thus, our findings contribute to a deeper understanding of the recognition and stability mechanisms of anti-hapten Nbs and shed new light on the rational design of novel haptens and directed evolution to produce high-performance antibodies.

Quantifying Amide-Aromatic Interactions at Molecular and Atomic Levels: Experimentally-determined Enthalpic and Entropic Contributions to Interactions of Amide sp 2 O, N, C and sp 3 C Unified Atoms with Naphthalene sp 2 C Atoms in Water

In addition to amide hydrogen bonds and the hydrophobic effect, interactions involving π-bonded sp 2 atoms of amides, aromatics and other groups occur in protein self-assembly processes including folding, oligomerization and condensate formation. These interactions also occur in aqueous solutions of amide and aromatic compounds, where they can be quantified. Previous analysis of thermodynamic coefficients quantifying net-favorable interactions of amide compounds with other amides and aromatics revealed that interactions of amide sp 2 O with amide sp 2 N unified atoms (presumably C=O···H-N hydrogen bonds) and amide/aromatic sp 2 C (lone pair-π, n-π * ) are particularly favorable. Sp 3 C-sp 3 C (hydrophobic), sp 3 C-sp 2 C (hydrophobic, CH-π), sp 2 C-sp 2 C (hydrophobic, π-π) and sp 3 C-sp 2 N interactions are favorable, sp 2 C-sp 2 N interactions are neutral, while sp 2 O-sp 2 O and sp 2 N-sp 2 N self-interactions and sp 2 O-sp 3 C interactions are unfavorable. Here, from determinations of favorable effects of fourteen amides on naphthalene solubility at 10, 25 and 45 °C, we dissect amide-aromatic interaction free energies into enthalpic and entropic contributions and find these vary systematically with amide composition. Analysis of these results yields enthalpic and entropic contributions to intrinsic strengths of interactions of amide sp 2 O, sp 2 N, sp 2 C and sp 3 C unified atoms with aromatic sp 2 C atoms. For each interaction, enthalpic and entropic contributions have the same sign and are much larger in magnitude than the interaction free energy itself. The amide sp 2 O-aromatic sp 2 C interaction is enthalpy-driven and entropically unfavorable, consistent with direct chemical interaction (e.g. lone pair-π) while amide sp 3 C- and sp 2 C-aromatic sp 2 C interactions are entropy-driven and enthalpically unfavorable, consistent with hydrophobic effects. These findings are relevant for interactions involving π-bonded sp 2 atoms in protein processes.

Table of Contents Graphic

The "Beacon" Structural Model of Protein Folding: Application for Trp-Cage in Water

Protein folding is a process in which a polypeptide must undergo folding process to obtain its three-dimensional structure. Thermodynamically, it is a process of enthalpy to overcome the loss of conformational entropy in folding. Folding is primarily related to hydrophobic interactions and intramolecular hydrogen bondings. During folding, hydrophobic interactions are regarded to be the driving forces, especially in the initial structural collapse of a protein. Additionally, folding is guided by the strong interactions within proteins, such as intramolecular hydrogen bondings related to the α-helices and β-sheets of proteins. Therefore, a protein is divided into the folding key (FK) regions related to intramolecular hydrogen bondings and the non-folding key (non-FK) regions. Various conformations are expected for FK and non-FK regions. Different from non-FK regions, it is necessary for FK regions to form the specific conformations in folding, which are regarded as the necessary folding pathways (or "beacons"). Additionally, sequential folding is expected for the FK regions, and the intermediate state is found during folding. They are reflected on the local basins in the free energy landscape (FEL) of folding. To demonstrate the structural model, molecular dynamics (MD) simulations are conducted on the folding pathway of the TRP-cage in water.

Ordered-domain unfolding of thermophilic isolated β subunit ATP synthase

The flexibility of the ATP synthase's β subunit promotes its role in the ATP synthase rotational mechanism, but its domains stability remains unknown. A reversible thermal unfolding of the isolated β subunit (Tβ) of the ATP synthase from Bacillus thermophilus PS3, tracked through circular dichroism and molecular dynamics, indicated that Tβ shape transits from an ellipsoid to a molten globule through an ordered unfolding of its domains, preserving the β-sheet residual structure at high temperature. We determined that part of the stability origin of Tβ is due to a transversal hydrophobic array that crosses the β-barrel formed at the N-terminal domain and the Rossman fold of the nucleotide-binding domain (NBD), while the helix bundle of the C-terminal domain is the less stable due to the lack of hydrophobic residues, and thus the more flexible to trigger the rotational mechanism of the ATP synthase.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Identification of PS1/gamma-secretase and glutamate transporter GLT-1 interaction sites

The recently discovered interaction between Presenilin 1 (PS1), a catalytic subunit of γ-secretase responsible for generating amyloid-β (Aβ) peptides, and GLT-1, a major glutamate transporter in the brain (EAAT2) provides a mechanistic link between these two key factors involved in Alzheimer's disease (AD) pathology. Modulating this interaction can be crucial to understand the consequence of such crosstalk in AD context and beyond. However, the interaction sites between these two proteins are unknown. Herein, we utilized an alanine scanning approach coupled with FRET-based fluorescence lifetime imaging microscopy (FLIM) to identify the interaction sites between PS1 and GLT-1 in their native environment within intact cells. We found that GLT-1 residues at position 276 to 279 (TM5) and PS1 residues at position 249 to 252 (TM6) are crucial for GLT-1/PS1 interaction. These results have been cross validated using AlphaFold Multimer prediction. To further investigate whether this interaction of endogenously expressed GLT-1 and PS1 can be prevented in primary neurons, we designed PS1/GLT-1 cell-permeable peptides (CPPs) targeting the PS1 or GLT-1 binding site. We used HIV TAT domain to allow for cell penetration which was assayed in neurons. First, we assessed the toxicity and penetration of CPPs by confocal microscopy. Next, to ensure the efficiency of CPPs, we monitored the modulation of GLT-1/PS1 interaction in intact neurons by FLIM. We saw significantly less interaction between PS1 and GLT-1 with both CPPs. Our study establishes a new tool to study the functional aspect of GLT-1/PS1 interaction and its relevance in normal physiology and AD models.

Human Hemoglobin and Antipsychotics Clozapine, Ziprasidone and Sertindole: Friends or Foes?

Packed with hemoglobin, an essential protein for oxygen transport, human erythrocytes are a suitable model system for testing the pleiotropic effects of lipophilic drugs. Our study investigated the interaction between antipsychotic drugs clozapine, ziprasidone, sertindole, and human hemoglobin under simulated physiological conditions. Analysis of protein fluorescence quenching at different temperatures and data obtained from the van't Hoff diagram and molecular docking indicate that the interactions are static and that the tetrameric human hemoglobin has one binding site for all drugs in the central cavity near αβ interfaces and is dominantly mediated through hydrophobic forces. The association constants were lower-moderate strength (~10⁴ M^-1), the highest observed for clozapine (2.2 × 10⁴ M^-1 at 25 °C). The clozapine binding showed "friendly" effects: increased α-helical content, a higher melting point, and protein protection from free radical-mediated oxidation. On the other hand, bound ziprasidone and sertindole had a slightly pro-oxidative effect, increasing ferrihemoglobin content, a possible "foe". Since the interaction of proteins with drugs plays a vital role in their pharmacokinetic and pharmacodynamic properties, the physiological significance of the obtained findings is briefly discussed.

Lipid Peroxidation Products Influence Calpain-1 Functionality In Vitro by Covalent Binding

The objective of the current study was to evaluate the effects of lipid peroxidation products, malondialdehyde (MDA), hexenal, and 4-hydroxynonenal (HNE), on calpain-1 function, and liquid chromatography and tandem mass spectrometry (LC-MS/MS) identification of adducts on calpain-1. Calpain-1 activity slightly increased after incubation with 100 μM MDA but not with 500 and 1000 μM MDA. However, calpain-1 activity was lowered by hexenal and HNE at 100, 500, and 1000 μM. No difference in calpain-1 autolysis was observed between the control and 1000 μM MDA. However, 1000 μM hexenal and HNE treatments slowed the calpain-1 autolysis. Adducts of MDA were detected on glutamine, arginine, lysine, histidine, and asparagine residues via Schiff base formation, while HNE adducts were detected on histidine, lysine, glutamine, and asparagine residues via Michael addition. These results are the first to demonstrate that lipid peroxidation products can impact calpain-1 activity in a concentration-dependent manner and may impact the development of meat tenderness postmortem.

In silico study of potential SARS-CoV-2 antagonist from Clitoria ternatea

Objectives

In this study, we implemented a structure-based virtual screening protocol in search of natural bioactive compounds in Clitoria ternatea that could inhibit the viral Mpro.

Methods

A library of twelve main bioactive compounds in C. ternatea was created from PubChem database by minimizing ligand structure in PyRx software to increase the ligand flexibility. Molecular docking studies were performed by targeting Mpro (PDB ID: 6lu7) via Discovery Studio Visualiser and PyRx platforms. Top hits compounds were then selected to study their Adsorption, distribution, metabolism, excretion, and toxicity (ADMET) and drug likeness properties through pkCSM pharmacokinetics tool to understand the stability, interaction, conformational changes, and pharmaceutical relevant parameters.

Results

This investigation found that, in the molecular docking simulation, four bioactive compounds (procyanidin A2 [-9.3 kcal/mol], quercetin-3-rutinoside [-8.9 kcal/mol], delphinidin-3-O-glucoside [-8.3 kcal/mol], and ellagic acid [-7.4 kcal/mol]) showed producing the strongest binding affinity to the Mpro of severe acute respiratory syndrome coronavirus 2, as compared to positive control (N3 inhibitor) (-7.5 kcal/mol). These binding energies were found to be favorable for an efficient docking and resultant. In addition, the stability of quercetin-3-rutinoside and ellagic acid is higher without any unfavorable bond. The ADMET and drug likeness of these two compounds were found that they are considered an effective and safe coronavirus disease 2019 (COVID-19) inhibitors through Lipinski's Rule, absorption, distribution, metabolism, and toxicity properties.

Conclusion

From these results, it was concluded that C. ternatea possess potential therapeutic properties against COVID-19.

Sulfur-mediated chalcogen versus hydrogen bonds in proteins: a see-saw effect in the conformational space

Divalent sulfur (S) forms a chalcogen bond (Ch-bond) via its σ-holes and a hydrogen bond (H-bond) via its lone pairs. The relevance of these interactions and their interplay for protein structure and function is unclear. Based on the analyses of the crystal structures of small organic/organometallic molecules and proteins and their molecular electrostatic surface potential, we show that the reciprocity of the substituent-dependent strength of the σ-holes and lone pairs correlates with the formation of either Ch-bond or H-bond. In proteins, cystines preferentially form Ch-bonds, metal-chelated cysteines form H-bonds, while methionines form either of them with comparable frequencies. This has implications for the positioning of these residues and their role in protein structure and function. Computational analyses reveal that the S-mediated interactions stabilise protein secondary structures by mechanisms such as helix capping and protecting free β-sheet edges by negative design. The study highlights the importance of S-mediated Ch-bond and H-bond for understanding protein folding and function, the development of improved strategies for protein/peptide structure prediction and design and structure-based drug discovery.

In silico repurposing of CNS drugs for multiple sclerosis

Multiple sclerosis (MS) is an autoimmune neurodegenerative disease affecting numerous people worldwide. While the relapsing subtypes of MS are to some extent treatable, the disease remains incurable leading to progressive disability. Limited efficacy of current small molecule drugs necessitates development of efficient and safe MS medications. Accordingly, drug repurposing is an invaluable strategy that recognizes new targets for known drugs especially in the field of poorly addressed therapeutic areas. Drug discovery largely depends on the identification of potential binding molecules to the intended biomolecular target(s). In this regard, current study was devoted to in silico repurposing of 263 small molecule CNS drugs to achieve superior binders to some MS-related targets. On the basis of molecular docking scores, thioxanthene and benzisothiazole-based antipsychotics could be identified as potential binders to sphingosine-1-phosphate lyase (S1PL) and cyclophilin D (CypD). Tightest interaction modes were observed for zuclopenthixol-S1PL (ΔG_b -7.96 kcal/mol) and lurasidone-CypD (ΔG_b -8.84 kcal/mol) complexes. Molecular dynamics (MD) simulations proved the appropriate and stable accommodation of top-ranked drugs inside enzyme binding sites during 100 ns. Hydroxyethyl piperazine of zuclopenthixol and benzisothiazole of lurasidone flipped inside the binding pocket to interact with adjacent polar and apolar residues. Solvent accessible surface area (SASA) fluctuations confirmed the results of binding trajectory analysis and showed that non-polar hydrophobic interactions played significant roles in acquired stabilities. Our results on lurasidone binding pattern were interestingly in accordance with previous reports on X-ray structures of other norbornane maleimide derivatives as CypD inhibitors. According to this, Asn144, Phe102 and Phe155 served as important residues in providing stable binding pose of lurasidone through both exo and endo conformations. Although experimental results are necessary to be achieved, the outcomes of this study proposed the potentiality of some thioxanthene and benzisothiazole-based antipsychotics for binding to S1PL and CypD, respectively, as MS-related targets.

Identification of natural small molecule modulators of MurB from Salmonella enterica serovar Typhi Ty2 strain using computational and biophysical approaches

The increase of antibiotic-resistant bacterial pathogens has created challenges in treatment and warranted the design of antibiotics against comparatively less exploited targets. The peptidoglycan (PG) biosynthesis delineates unique pathways for the design and development of a novel class of drugs. Mur ligases are an essential component of bacterial cell wall synthesis that play a pivotal role in PG biosynthesis to maintain internal osmotic pressure and cell shape. Inhibition of these enzymes can interrupt bacterial replication and hence, form attractive targets for drug discovery. In the present work, we focused on the PG biosynthesis pathway enzyme, UDP-N-acetylpyruvylglucosamine reductase, from Salmonella enterica serovar Typhi (stMurB). Biophysical characterization of purified StMurB was performed to gauge the molecular interactions and estimate thermodynamic stability for determination of attributes for possible therapeutic intervention. The thermal melting profile of MurB was monitored by circular dichroism and validated through differential scanning calorimetry experiment. Frequently used chemical denaturants, GdmCl and urea, were employed to study the chemical-induced denaturation of stMurB. In the search for natural compound-based inhibitors, against this important drug target, an in silico virtual screening based investigation was conducted with modeled stMurB structure. The three top hits (quercetin, berberine, and scopoletin) returned were validated for complex stability through molecular dynamics simulation. Further, fluorescence binding studies were undertaken for the selected natural compounds with stMurB alone and with NADPH bound form. The compounds scopoletin and berberine, displayed lesser binding to stMurB whereas quercetin exhibited stronger binding affinity than NADPH. This study suggests that quercetin can be evolved as an inhibitor of stMurB enzyme.

Urea induced unfolding of rai seed cystatin: Influence of glycerol as a chemical chaperone

Cystatin superfamily members, by virtue of their thiol protease regulatory properties, show involvement in myriad physiological processes important for survival and well-being. The current study involves urea-induced denaturation of a novel variant of the cystatin superfamily, rai seed cystatin (RSC), employing a variety of biophysical assays in order to characterize different folding intermediates generated on unfolding. Urea as a denaturant presented the passage of RSC through a series of events resulting in the loss of RSC functional capability, accompanied by changes in the archetype at secondary and tertiary structural levels, as evident from protease inhibitory, UV absorption, and intrinsic fluorescence assays, respectively. ANS fluorescence also revealed routing of RSC through discrete multiple sub-states thus presenting the generation of intermediate states somewhat close to the pre-molten globule and/or molten globule forms of RSC. Furthermore, far-UV circular dichroism analysis revealed a concentration-dependent gradual loss in typical -helical RSC peaks, indicating a nearly 50 % loss in secondary structural elements around 5 M urea treatment. The study also reports the possible role of glycerol in the refolding and/or reactivation of the urea unfolded RSC form. Glycerol presented itself as a potent structural stabilizer as it assisted in the refolding and reactivation of the unfolded RSC in a dosage-dependent manner, concomitantly paving the way for unravelling the mechanistic approach involved in the phenomenon, which can facilitate future studies.

Caseins: Versatility of Their Micellar Organization in Relation to the Functional and Nutritional Properties of Milk

The milk of mammals is a complex fluid mixture of various proteins, minerals, lipids, and other micronutrients that play a critical role in providing nutrition and immunity to newborns. Casein proteins together with calcium phosphate form large colloidal particles, called casein micelles. Caseins and their micelles have received great scientific interest, but their versatility and role in the functional and nutritional properties of milk from different animal species are not fully understood. Caseins belong to a class of proteins that exhibit open and flexible conformations. Here, we discuss the key features that maintain the structures of the protein sequences in four selected animal species: cow, camel, human, and African elephant. The primary sequences of these proteins and their posttranslational modifications (phosphorylation and glycosylation) that determine their secondary structures have distinctively evolved in these different animal species, leading to differences in their structural, functional, and nutritional properties. The variability in the structures of milk caseins influence the properties of their dairy products, such as cheese and yogurt, as well as their digestibility and allergic properties. Such differences are beneficial to the development of different functionally improved casein molecules with variable biological and industrial utilities.

Mutation of the polyadenylation complex subunit CstF77 reveals that mRNA 3' end formation and HSP101 levels are critical for a robust heat stress response

Heat shock protein 101 (HSP101) in plants, and bacterial and yeast orthologs, is essential for thermotolerance. To investigate thermotolerance mechanisms involving HSP101, we performed a suppressor screen in Arabidopsis thaliana of a missense HSP101 allele (hot1-4). hot1-4 plants are sensitive to acclimation heat treatments that are otherwise permissive for HSP101 null mutants, indicating that the hot1-4 protein is toxic. We report one suppressor (shot2, suppressor of hot1-4 2) has a missense mutation of a conserved residue in CLEAVAGE STIMULATION FACTOR77 (CstF77), a subunit of the polyadenylation complex critical for mRNA 3' end maturation. We performed ribosomal RNA depletion RNA-Seq and captured transcriptional readthrough with a custom bioinformatics pipeline. Acclimation heat treatment caused transcriptional readthrough in hot1-4 shot2, with more readthrough in heat-induced genes, reducing the levels of toxic hot1-4 protein and suppressing hot1-4 heat sensitivity. Although shot2 mutants develop like the wild type in the absence of stress and survive mild heat stress, reduction of heat-induced genes and decreased HSP accumulation makes shot2 in HSP101 null and wild-type backgrounds sensitive to severe heat stress. Our study reveals the critical function of CstF77 for 3' end formation of mRNA and the dominant role of HSP101 in dictating the outcome of severe heat stress.

Immunoinformatics Study: Multi-Epitope Based Vaccine Design from SARS-CoV-2 Spike Glycoprotein

The coronavirus disease 2019 outbreak has become a huge challenge in the human sector for the past two years. The coronavirus is capable of mutating at a higher rate than other viruses. Thus, an approach for creating an effective vaccine is still needed to induce antibodies against multiple variants with lower side effects. Currently, there is a lack of research on designing a multiepitope of the COVID-19 spike protein for the Indonesian population with comprehensive immunoinformatic analysis. Therefore, this study aimed to design a multiepitope-based vaccine for the Indonesian population using an immunoinformatic approach. This study was conducted using the SARS-CoV-2 spike glycoprotein sequences from Indonesia that were retrieved from the GISAID database. Three SARS-CoV-2 sequences, with IDs of EIJK-61453, UGM0002, and B.1.1.7 were selected. The CD8+ cytotoxic T-cell lymphocyte (CTL) epitope, CD4+ helper T lymphocyte (HTL) epitope, B-cell epitope, and IFN-γ production were predicted. After modeling the vaccines, molecular docking, molecular dynamics, in silico immune simulations, and plasmid vector design were performed. The designed vaccine is antigenic, non-allergenic, non-toxic, capable of inducing IFN-γ with a population reach of 86.29% in Indonesia, and has good stability during molecular dynamics and immune simulation. Hence, this vaccine model is recommended to be investigated for further study.

Engineering of a T7 Bacteriophage Endolysin Variant with Enhanced Amidase Activity

The therapeutic use of bacteriophage-encoded endolysins as enzybiotics has increased significantly in recent years due to the emergence of antibiotic resistant bacteria. Phage endolysins lyse the bacteria by targeting their cell wall. Various engineering strategies are commonly used to modulate or enhance the utility of therapeutic enzymes. This study employed a structure-guided mutagenesis approach to engineer a T7 bacteriophage endolysin (T7L) with enhanced amidase activity and lysis potency via replacement of a noncatalytic gating residue (His 37). Two H37 variants (H37A and H37K) were designed and characterized comprehensively using integrated biophysical and biochemical techniques to provide mechanistic insights into their structure-stability-dynamics-activity paradigms. Among the studied proteins, cell lysis data suggested that the obtained H37A variant exhibits amidase activity (∼35%) enhanced compared to that of wild-type T7 endolysin (T7L-WT). In contrast to this, the H37K variant is highly unstable, prone to aggregation, and less active. Comparison of the structure and dynamics of the H37A variant to those of T7L-WT evidenced that the alteration at the site of H37 resulted in long-range structural perturbations, attenuated the conformational heterogeneity, and quenched the microsecond to millisecond time scale motions. Stability analysis confirmed the altered stability of H37A compared to that of its WT counterpart. All of the obtained results established that the H37A variant enhances the lysis activity by regulating the stability-activity trade-off. This study provided deeper atomic level insights into the structure-function relationships of endolysin proteins, thus aiding researchers in the rational design of engineered endolysins with enhanced therapeutic properties.

Insights Into The Effects of Amino Acid Substitutions on The Stability of Reteplase Structure: A Molecular Dynamics Simulation Study

Background

Reteplase (recombinant plasminogen activator, r-PA) is a recombinant protein designed to imitate the endogenous tissue plasminogen activator and catalyze the plasmin production. It is known that the application of reteplase is limited by the complex production processes and protein's stability challenges. Computational redesign of proteins has gained momentum in recent years, particularly as a powerful tool for improving protein stability and consequently its production efficiency. Hence, in the current study, we implemented computational approaches to improve r-PA conformational stability, which fairly correlates with protein's resistance to proteolysis.

Objectives

The current study was developed in order to evaluate the effect of amino acid substitutions on the stability of reteplase structure using molecular dynamic simulations and computational predictions.

Materials and Methods

Several web servers designed for mutation analysis were utilized to select appropriate mutations. Additionally, the experimentally reported mutation, R103S, converting wild type r-PA into non-cleavable form, was also employed. Firstly, mutant collection, consisting of 15 structures, was constructed based on the combinations of four designated mutations. Then, 3D structures were generated using MODELLER. Finally, 17 independent 20-ns molecular dynamics (MD) simulations were conducted and different analysis were performed like root-mean-square deviation (RMSD), root-mean-square fluctuations (RMSF), secondary structure analysis, number of hydrogen bonds, principal components analysis (PCA), eigenvector projection, and density analysis.

Results

Predicted mutations successfully compensated the more flexible conformation caused by R103S substitution, so, improved conformational stability was analyzed from MD simulations. In particular, R103S/A286I/G322I indicated the best results and remarkably enhanced the protein stability.

Conclusion

The conformational stability conferred by these mutations will probably lead to more protection of r-PA in protease-rich environments in various recombinant systems and potentially enhance its production and expression level.

Impact of Cd(II) on the stability of human uracil DNA glycosylase enzyme; an implication of molecular dynamics trajectories on stability analysis

Uracil DNA glycosylase is a key enzyme that identifies and removes damaged bases from DNA in the base excision repair pathway. Experimentalists have identified the possibility of Cd(II) reducing the activity of human uracil DNA glycosylase (hUNG) by binding with the enzyme replacing the catalytic water molecule. The present study focus on the stability variation of the enzyme in the presence and absence of Cd(II) and confirms the reported results with the stability analysis done using molecular dynamic (MD) simulation trajectories. The CavityPlus web server identified seven cavities for the free enzyme as possible binding sites and a cavity containing the active site of the enzyme as the best binding cavity for a ligand. Based on the CavityPlus results and the previously reported work, a free hUNG system and two systems of the enzyme with Cd(II); one with Cd(II) replacing the catalytic water molecule in the active site of the enzyme and the other replacing a non-catalytic water molecule in the active site were generated for the simulation. The simulation trajectories were used for the structural stability analysis of the enzyme in all three systems. The binding free energy of the Cd(II) with the enzyme was calculated using molecular mechanics Poisson Boltzmann surface area method. The results showed that the enzyme achieves comparatively high stability with the removal of catalytic water of the enzyme by Cd(II). Therefore, this supports the previously reported idea that Cd(II) replaces catalytic water molecules and affects enzyme activity.

The use of pharmacological chaperones in rare diseases caused by reduced protein stability

Rare diseases are most often caused by inherited genetic disorders that, after translation, will result in a protein with altered function. Decreased protein stability is the most frequent mechanism associated with a congenital pathogenic missense mutation and it implies the destabilization of the folded conformation in favour of unfolded or misfolded states. In the cellular context and when experimental data is available, a mutant protein with altered thermodynamic stability often also results in impaired homeostasis, with the deleterious accumulation of protein aggregates, metabolites and/or metabolic by-products. In the last decades, a significant effort has enabled the characterization of rare diseases associated to protein stability defects and triggered the development of innovative therapeutic intervention lines, say, the use of pharmacological chaperones to correct the intracellular impaired homeostasis. Here, we review the current knowledge on rare diseases caused by reduced protein stability, paying special attention to the thermodynamic aspects of the protein destabilization, also focusing on some examples where pharmacological chaperones are being tested.

Dissecting the stability determinants of a challenging de novo protein fold using massively parallel design and experimentation

Designing entirely new protein structures remains challenging because we do not fully understand the biophysical determinants of folding stability. Yet, some protein folds are easier to design than others. Previous work identified the 43-residue ɑββɑ fold as especially challenging: The best designs had only a 2% success rate, compared to 39 to 87% success for other simple folds [G. J. Rocklin et al., Science 357, 168-175 (2017)]. This suggested the ɑββɑ fold would be a useful model system for gaining a deeper understanding of folding stability determinants and for testing new protein design methods. Here, we designed over 10,000 new ɑββɑ proteins and found over 3,000 of them to fold into stable structures using a high-throughput protease-based assay. NMR, hydrogen-deuterium exchange, circular dichroism, deep mutational scanning, and scrambled sequence control experiments indicated that our stable designs fold into their designed ɑββɑ structures with exceptional stability for their small size. Our large dataset enabled us to quantify the influence of universal stability determinants including nonpolar burial, helix capping, and buried unsatisfied polar atoms, as well as stability determinants unique to the ɑββɑ topology. Our work demonstrates how large-scale design and test cycles can solve challenging design problems while illuminating the biophysical determinants of folding.

Stability profile of the neuronal SNARE complex reflects its potency to drive fast membrane fusion

Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) form the SNARE complex to mediate most fusion events of the secretory pathway. The neuronal SNARE complex is featured by its high stability and half-zippered conformation required for driving robust and fast synaptic exocytosis. However, these two features seem to be thermodynamically mutually exclusive. In this study, we have employed temperature-dependent disassociation assays and single-molecule Förster resonance energy transfer (FRET) experiments to analyze the stability and conformation of the neuronal SNARE complex. We reclassified the amino acids of the SNARE motif into four sub-groups (core, core-side I and II, and non-contact). Our data showed that the core residues predominantly contribute to the complex stability to meet a basal requirement for SNARE-mediated membrane fusion, while the core-side residues exert an unbalanced effect on the N- and C-half bundle stability that determines the half-zippered conformation of the neuronal SNARE complex, which would accommodate essential regulations by complexins and synaptotagmins for fast Ca2+-triggered membrane fusion. Furthermore, our data confirmed a strong coupling of folding energy between the N- and C-half assembly of the neuronal SNARE complex, which rationalizes the strong potency of the half-zippered conformation to conduct robust and fast fusion. Overall, these results uncovered that the stability profile of the neuronal SNARE complex reflects its potency to drive fast and robust membrane fusion. Based on these results, we also developed a new parameter, the stability factor (Fs), to characterize the overall stability of the neuronal SNARE complex and resolved a linear correlation between the stability and inter-residue coulombic interactions of the neuronal SNARE complex, which would help rationally design artificial SNARE complexes and remold functional SNARE complexes with desirable stability.

Protein Albumin Manipulation and Electrical Quantification of Molecular Dielectrophoresis Responses for Biomedical Applications

Research relating to dielectrophoresis (DEP) has been progressing rapidly through time as it is a strong and controllable technique for manipulation, separation, preconcentration, and partitioning of protein. Extensive studies have been carried out on protein DEP, especially on Bovine Serum Albumin (BSA). However, these studies involve the usage of dye and fluorescent probes to observe DEP responses as the physical properties of protein albumin molecular structure are translucent. The use of dye and the fluorescent probe could later affect the protein's physiology. In this article, we review three methods of electrical quantification of DEP responses: electrochemical impedance spectroscopy (EIS), cyclic voltammetry (CV) and capacitance measurement for protein BSA DEP manipulation. The correlation of these methods with DEP responses is further discussed. Based on the observations on capacitance measurement, it can be deduced that the electrical quantifying method is reliable for identifying DEP responses. Further, the possibility of manipulating the protein and electrically quantifying DEP responses while retaining the original physiology of the protein and without the usage of dye or fluorescent probe is discussed.

Sustainable isomaltulose production in Corynebacterium glutamicum by engineering the thermostability of sucrose isomerase coupled with one-step simplified cell immobilization

Sucrose isomerase (SI), catalyzing sucrose to isomaltulose, has been widely used in isomaltulose production, but its poor thermostability is still resisted in sustainable batches production. Here, protein engineering and one-step immobilized cell strategy were simultaneously coupled to maintain steady state for long-term operational stabilities. First, rational design of Pantoea dispersa SI (PdSI) for improving its thermostability by predicting and substituting the unstable amino acid residues was investigated using computational analysis. After screening mutagenesis library, two single mutants (PdSIV280L and PdSIS499F) displayed favorable characteristics on thermostability, and further study found that the double mutant PdSIV280L/S499F could stabilize PdSIWT better. Compared with PdSIWT, PdSIV280L/S499F displayed a 3.2°C-higher T m , and showed a ninefold prolonged half-life at 45°C. Subsequently, a one-step simplified immobilization method was developed for encapsulation of PdSIV280L/S499F in food-grade Corynebacterium glutamicum cells to further enhance the recyclability of isomaltulose production. Recombinant cells expressing combinatorial mutant (RCSI2) were successfully immobilized in 2.5% sodium alginate without prior permeabilization. The immobilized RCSI2 showed that the maximum yield of isomaltulose by batch conversion reached to 453.0 g/L isomaltulose with a productivity of 41.2 g/l/h from 500.0 g/L sucrose solution, and the conversion rate remained 83.2% after 26 repeated batches.

Protein stability [determination] problems

Human health depends on the correct folding of proteins, for misfolding and aggregation lead to diseases. An unfolded (denatured) protein can refold to its original folded state. How does this occur is known as the protein folding problem. One of several related questions to this problem is that how much more stable is the folded state than the unfolded state. There are several measures of protein stability. In this article, protein stability is given a thermodynamic definition and is measured by Gibbs free energy change (ΔGD0) associated with the equilibrium, native (N) conformation ↔ denatured (D) conformation under the physiological condition usually taken as dilute buffer (or water) at 25 °C. We show that this thermodynamic quantity (ΔGD0), where subscript D represents transition between N and D states, and superscript 0 (zero) represents the fact that the transition occurs in the absence of denaturant, can be neither measured nor predicted under physiological conditions. However, ΔGD can be measured in the presence of strong chemical denaturants such as guanidinium chloride and urea which are shown to destroy all noncovalent interactions responsible for maintaining the folded structure. A problem with this measurement is that the estimate of ΔGD0 comes from the analysis of the plot of ΔGDversus denaturant concentration, which requires a long extrapolation of values of ΔGD, and all the three methods of extrapolation give three different values of ΔGD0 for a protein. Thus, our confidence in the authentic value of ΔGD0 is eroded. Another problem with this in vitro measurement of ΔGD0 is that it is done on the pure protein sample in dilute buffer which is a very large extrapolation of the in vivo conditions, for the crowding effect on protein stability is ignored.

The Flexible, Extended Coil of the PDZ-Binding Motif of the Three Deadly Human Coronavirus E Proteins Plays a Role in Pathogenicity

The less virulent human (h) coronaviruses (CoVs) 229E, NL63, OC43, and HKU1 cause mild, self-limiting respiratory tract infections, while the more virulent SARS-CoV-1, MERS-CoV, and SARS-CoV-2 have caused severe outbreaks. The CoV envelope (E) protein, an important contributor to the pathogenesis of severe hCoV infections, may provide insight into this disparate severity of the disease. We, therefore, generated full-length E protein models for SARS-CoV-1 and -2, MERS-CoV, HCoV-229E, and HCoV-NL63 and docked C-terminal peptides of each model to the PDZ domain of the human PALS1 protein. The PDZ-binding motif (PBM) of the SARS-CoV-1 and -2 and MERS-CoV models adopted a more flexible, extended coil, while the HCoV-229E and HCoV-NL63 models adopted a less flexible alpha helix. All the E peptides docked to PALS1 occupied the same binding site and the more virulent hCoV E peptides generally interacted more stably with PALS1 than the less virulent ones. We hypothesize that the increased flexibility of the PBM in the more virulent hCoVs facilitates more stable binding to various host proteins, thereby contributing to more severe disease. This is the first paper to model full-length 3D structures for both the more virulent and less virulent hCoV E proteins, providing novel insights for possible drug and/or vaccine development.

Investigation of the causal etiology in a patient with T-B+NK+ immunodeficiency

Newborn screening for severe combined immunodeficiency (SCID) has not only accelerated diagnosis and improved treatment for affected infants, but also led to identification of novel genes required for human T cell development. A male proband had SCID newborn screening showing very low T cell receptor excision circles (TRECs), a biomarker for thymic output of nascent T cells. He had persistent profound T lymphopenia, but normal numbers of B and natural killer (NK) cells. Despite an allogeneic hematopoietic stem cell transplant from his brother, he failed to develop normal T cells. Targeted resequencing excluded known SCID genes; however, whole exome sequencing (WES) of the proband and parents revealed a maternally inherited X-linked missense mutation in MED14 (MED14V763A), a component of the mediator complex. Morpholino (MO)-mediated loss of MED14 function attenuated T cell development in zebrafish. Moreover, this arrest was rescued by ectopic expression of cDNA encoding the wild type human MED14 ortholog, but not by MED14V763A , suggesting that the variant impaired MED14 function. Modeling of the equivalent mutation in mouse (Med14V769A) did not disrupt T cell development at baseline. However, repopulation of peripheral T cells upon competitive bone marrow transplantation was compromised, consistent with the incomplete T cell reconstitution experienced by the proband upon transplantation with bone marrow from his healthy male sibling, who was found to have the same MED14V763A variant. Suspecting that the variable phenotypic expression between the siblings was influenced by further mutation(s), we sought to identify genetic variants present only in the affected proband. Indeed, WES revealed a mutation in the L1 cell adhesion molecule (L1CAMQ498H); however, introducing that mutation in vivo in mice did not disrupt T cell development. Consequently, immunodeficiency in the proband may depend upon additional, unidentified gene variants.

A ricin-based peptide BRIP from Hordeum vulgare inhibits Mpro of SARS-CoV-2

COVID-19 pandemic caused by SARS-CoV-2 led to the research aiming to find the inhibitors of this virus. Towards this world problem, an attempt was made to identify SARS-CoV-2 main protease (M^pro) inhibitory peptides from ricin domains. The ricin-based peptide from barley (BRIP) was able to inhibit M^pro in vitro with an IC₅₀ of 0.52 nM. Its low and no cytotoxicity upto 50 µM suggested its therapeutic potential against SARS-CoV-2. The most favorable binding site on M^pro was identified by molecular docking and steered molecular dynamics (MD) simulations. The M^pro-BRIP interactions were further investigated by evaluating the trajectories for microsecond timescale MD simulations. The structural parameters of M^pro-BRIP complex were stable, and the presence of oppositely charged surfaces on the binding interface of BRIP and M^pro complex further contributed to the overall stability of the protein-peptide complex. Among the components of thermodynamic binding free energy, Van der Waals and electrostatic contributions were most favorable for complex formation. Our findings provide novel insight into the area of inhibitor development against COVID-19.

PSMC1 variant causes a novel neurological syndrome

Proteasome 26S, the eukaryotic proteasome, serves as the machinery for cellular protein degradation. It is composed of the 20S core particle and one or two 19S regulatory particles, composed of a base and a lid. To date, several human diseases have been associated with mutations within the 26S proteasome subunits; only one of them affects a base subunit. We now delineate an autosomal recessive syndrome of failure to thrive, severe developmental delay and intellectual disability, spastic tetraplegia with central hypotonia, chorea, hearing loss, micropenis and undescended testes, as well as mild elevation of liver enzymes. None of the affected individuals achieved verbal communication or ambulation. Ventriculomegaly was evident on MRI. Homozygosity mapping combined with exome sequencing revealed a disease‐associated p.I328T PSMC1 variant. Protein modeling demonstrated that the PSMC1 variant is located at the highly conserved putative ATP binding and hydrolysis domain, and is suggested to interrupt a hydrophobic core within the protein. Fruit flies in which we silenced the Drosophila ortholog Rpt2 specifically in the eye exhibited an apparent phenotype that was highly rescued by the human wild‐type PSMC1, yet only partly by the mutant PSMC1, proving the functional effect of the p.I328T disease‐causing variant.

Endophytes in Lignin Valorization: A Novel Approach

Lignin, one of the essential components of lignocellulosic biomass, comprises an abundant renewable aromatic resource on the planet earth. Although 15%––40% of lignocellulose pertains to lignin, its annual valorization rate is less than 2% which raises the concern to harness and/or develop effective technologies for its valorization. The basic hindrance lies in the structural heterogeneity, complexity, and stability of lignin that collectively makes it difficult to depolymerize and yield common products. Recently, microbial delignification, an eco-friendly and cheaper technique, has attracted the attention due to the diverse metabolisms of microbes that can channelize multiple lignin-based products into specific target compounds. Also, endophytes, a fascinating group of microbes residing asymptomatically within the plant tissues, exhibit marvellous lignin deconstruction potential. Apart from novel sources for potent and stable ligninases, endophytes share immense ability of depolymerizing lignin into desired valuable products. Despite their efficacy, ligninolytic studies on endophytes are meagre with incomplete understanding of the pathways involved at the molecular level. In the recent years, improvement of thermochemical methods has received much attention, however, we lagged in exploring the novel microbial groups for their delignification efficiency and optimization of this ability. This review summarizes the currently available knowledge about endophytic delignification potential with special emphasis on underlying mechanism of biological funnelling for the production of valuable products. It also highlights the recent advancements in developing the most intriguing methods to depolymerize lignin. Comparative account of thermochemical and biological techniques is accentuated with special emphasis on biological/microbial degradation. Exploring potent biological agents for delignification and focussing on the basic challenges in enhancing lignin valorization and overcoming them could make this renewable resource a promising tool to accomplish Sustainable Development Goals (SDG’s) which are supposed to be achieved by 2030.

Engineered Design of the E-Helix Structure on Ferritin Nanoparticles

Insertion of an immunogenic epitope at the C-terminus of ferritin has shown the potential to produce a stable and efficacious vaccine. There is however limited understanding of how C-terminus insertion affects ferritin protein stability. The E-helix at the C-terminus has attracted interest because there are contradictory reports as to whether it has a role in protein stabilization. Here, we report, for the first time, combining molecular dynamics simulation (MDS) with experiment to engineer the design of the E-helix at the C-terminus of engineered human ferritin heavy chain (F₁) inserted with Epstein-Barr nuclear antigen 1 (EBNA1, E₁) and flexible linker (L₃) residues (to afford F₁L₃E₁). Hot spots on the E-helix of the C-terminus were predicted by MDS at aa 167 (Glu) and aa 171 (Asp). Five (5) variants of F₁L₃E₁ were constructed by considering hot spots and alteration of electrostatic or hydrophobic interfaces, namely, (1) C1, hot spots substituted with noncharged residue Gln; (2) C2, hot spots substituted with positively charged residue Arg; (3) C3, hydrophobic residues substituted with the most hydrophobic residues Val and Ile; (4) C4, hydrophobic residues substituted with the most hydrophilic residues Gln and Asn; and (5) C5, a heptad repeat structure in the E-helix disrupted by substituting "a" and "d" heptad residues with noncharged polar residue Gln. It was found that the E-helix is essential to maintain integrated protein stability and that changing the hydrophobic interface (C3 and C4) had more significant effects on protein folding and stability than changing the electrostatic interface (C1 and C2). It was confirmed by both MDS and experiment that variants C1, C2, and C5 were able to fold to form stable conformational structures with protein surface hydrophobicity similar to that of F₁L₃E₁. However, they are less thermally stable than F₁L₃E₁. Significant changes in hydrophobicity drove significant protein aggregation for variants C3 and C4. It is concluded that the molecular design of the C-terminus in engineered ferritin, especially the E-helix, is important to ensure the epitope-based chimeric vaccine is safe (aggregate free) and efficacious.