Search

Lactic acid bacterial exopolysaccharides (EPS) are used in the food industry to improve the stability and rheological properties of fermented dairy products. β-Lactoglobulin (BLG), the dominant whey protein in bovine milk, is well known to bind small molecules such as fatty acids, vitamins, and flavors, and to interact with neutral and anionic polysaccharides used in food and pharmaceuticals. While sparse data are available on the affinity of EPS-milk protein interactions, structural information on BLG-EPS complexes, including the EPS binding sites, is completely lacking. Here, binding sites on BLG variant A (BLGA), for oligosaccharides prepared by mild acid hydrolysis of two EPS produced by Streptococcus thermophilus LY03 and Lactobacillus delbrueckii ssp. bulgaricus CNRZ 1187, respectively, are identified by NMR spectroscopy and supplemented by isothermal titration calorimetry (ITC) and molecular docking of complexes. Evidence of two binding sites (site 1 and site 2) on the surface of BLGA is achieved for both oligosaccharides (LY03-OS and 1187-OS) through NMR chemical shift perturbations, revealing multivalency of BLGA for EPS. The affinities of LY03-OS and 1187-OS for BLGA gave K D values in the mM range obtained by both NMR (pH 2.65) and ITC (pH 4.0). Molecular docking suggested that the BLGA and EPS complexes depend on hydrogen bonds and hydrophobic interactions. The findings provide insights into how BLGA engages structurally different EPS-derived oligosaccharides, which may facilitate the design of BLG-EPS complexation, of relevance for formulation of dairy products and improve understanding of BLGA coacervation.

Conformational stability of human serum transferrin (Tf) at varying pH, salt, and excipient concentrations were investigated using molecular dynamics (MD) simulations and the results are compared with previously published small angle X-ray scattering (SAXS) experiments. SAXS study showed that at pH 5, Tf is predominantly present in partially open (PO) form, and the factions of PO differ based on the physicochemical condition and drifts towards closed form (HO) as the pH increases. Tf is a bilobal glycoprotein that is composed of homologous halves termed as N-lobe and C-lobe. The current study shows that the protonation of Y188 and K206 at pH 5 is the primary conformational drive into PO, which shifts towards the closed (HO) conformer as the pH increases. Furthermore, at pH 6.5, PO is unfavorable due to negative charge-charge repulsion at the N/C-lobe interface linker region causing increased hinge distance when compared to HO, which has favorable attractive electrostatics. Subsequently, the effect of salt concentration at 70 and 140 mM NaCl was studied. At 70 mM NaCl and pH 5, chloride ions bind strongly in the N-lobe iron-binding site, whereas these interactions are weak at pH 6.5. With increasing salt concentration at pH 5, regions surrounding the N-lobe iron-binding site are saturated and as a consequence sodium and chloride ions accumulate into the bulk. Additionally, protein-excipient interactions were investigated. At pH 5, excipients interact in similar loop regions, E89-T93, D416-D420, located in the C-lobe and N-lobe of the HO conformer, respectively. It is anticipated that interactions of additives in these two loop regions cause conformational changes that lead to iron coordinating residues in the N-lobe to drift away from iron and thus drive HO to PO conversion. Furthermore, at pH 6.5 and 140 mM histidine or phosphate, these interactions are negligible leading to the stabilization of HO.

Water play an important role in many structural and physicochemical properties of lyophilized proteins. Molecular dynamics simulations were employed to study the explicit water distributions on four structurally diversed proteins: insulin-like growth factor 1 (IGF1), immunoglobin G1 (IgG1), human serum albumin (HSA), and collagen. The MD simulations were combined with the literature data on water vapor sorption isotherms. To account for the heterogeneity of protein surface, the water molecules were classified into different groups according to the binding strengths. A mechanistic mathematical model was built to describe the type-II vapor sorption isotherms and successfully applied to all four model protein systems. Although commonly used Brunauer-Emmett-Teller (BET) theory has a good fitting to the experimental vapor sorption isotherms, the basic "monolayer" concept is not consistent with reality - covering too limited protein surface. Experimentally, several physicochemical properties did show a break point near the BET "monolayer" level. This study demonstrates that the water content threshold or BET "monolayer" is consistent with the onset of water cluster (n≥3) formation. Based on water distributions at different amino acid sidechains as well as the backbones, a simple formula was derived based on primary sequence and fractions of ordered secondary structures (i.e. alpha helix and beta sheet) to predict the BET "monolayer". We find that proteins with helical structural elements are more stable upon changes in water content compared to other protein architectures.

European journal of pharmaceutics and biopharmaceutics 171, 1 - 10 (2022). doi:10.1016/j.ejpb.2021.08.018, High throughput screening for measuring the stability of industrially relevant proteins and their variants is necessary for quality assessment in the development process. Advances in automation, measurement time and sample consumption for many techniques allow rapid measurements with minimal amount of protein. However, many methods include automated data analysis, potentially neglecting important aspects of the protein's behavior in certain conditions. In this study we implement small angle X-ray scattering (SAXS), typically not used to assess protein behavior in industrial screening, in a high throughput screening workflow to address problems of contradicting results and reproducibility among different high throughput methods. As a case study we use the lipases of Thermomyces lanuginosus and Rhizomucor miehei, widely used industrial biocatalysts. We show that even the initial analysis of the SAXS data without performing any time-consuming modelling provide valuable information on interparticle interactions. We conclude that recent advances in automation and data processing, have enabled SAXS to be used more widely as a tool to gain in-depth knowledge highly useful for protein formulation development. This is especially relevant in light of increasing accessibility to SAXS due to the commercial availability of benchtop instruments., Published by Elsevier, New York, NY [u.a.]

Self-assembly and fibril formation play important roles in protein behavior. Amyloid fibrils formation is well-studied due to its role in neurodegenerative diseases and characterized by refolding of the protein into predominant β-sheet form. However, much less is known about the assembly of proteins into other types of supramolecular structures. Using cryo-electron microscopy at a resolution of 1.97 Å, we show that a triple-mutant of the anti-microbial peptide plectasin assembles reversibly into helical non-amyloid fibrils. Plectasin contains a cysteine-stabilized α-helix-β-sheets structure, which remains intact upon fibril formation. Two fibrils form a right-handed superstructure with each fibril consisting of double helical, left-handed structures. The fibril formation is reversible and follows sigmoidal kinetics with a pH-dependent equilibrium between soluble monomer and protein fibril. The anti-microbial activity does not appear compromised by fibril formation. This is the first high-resolution structure of this type of α/β protein fibrils.

The pterin-dependent nonheme iron enzymes hydroxylate aromatic amino acids to perform the biosynthesis of neurotransmitters to maintain proper brain function. These enzymes activate oxygen using a pterin cofactor and an aromatic amino acid substrate bound to the Fe(II) active site to form a highly reactive Fe(IV) = O species that initiates substrate oxidation. In this study, using tryptophan hydroxylase, we have kinetically generated a pre-Fe(IV) = O intermediate and characterized its structure as a Fe(II)-peroxy-pterin species using absorption, Mössbauer, resonance Raman, and nuclear resonance vibrational spectroscopies. From parallel characterization of the pterin cofactor and tryptophan substrate–bound ternary Fe(II) active site before the O(2) reaction (including magnetic circular dichroism spectroscopy), these studies both experimentally define the mechanism of Fe(IV) = O formation and demonstrate that the carbonyl functional group on the pterin is directly coordinated to the Fe(II) site in both the ternary complex and the peroxo intermediate. Reaction coordinate calculations predict a 14 kcal/mol reduction in the oxygen activation barrier due to the direct binding of the pterin carbonyl to the Fe(II) site, as this interaction provides an orbital pathway for efficient electron transfer from the pterin cofactor to the iron center. This direct coordination of the pterin cofactor enables the biological function of the pterin-dependent hydroxylases and demonstrates a unified mechanism for oxygen activation by the cofactor-dependent nonheme iron enzymes.

While heterogeneous enzyme reactions play an essential role in both nature and green industries, computational predictions of their catalytic properties remain scarce. Recent experimental work demonstrated the applicability of the Sabatier principle for heterogeneous biocatalysis. This provides a simple relationship between binding strength and the catalytic rate and potentially opens a new way for inexpensive computational determination of kinetic parameters. However, broader implementation of this approach will require fast and reliable prediction of binding free energies of complex two-phase systems, and computational procedures for this are still elusive. Here, we propose a new framework for the assessment of the binding strengths of multidomain proteins, in general, and interfacial enzymes, in particular, based on an extended linear interaction energy (LIE) method. This two-domain LIE (2D-LIE) approach was successfully applied to predict binding and activation free energies of a diverse set of cellulases and resulted in robust models with high accuracy. Overall, our method provides a fast computational screening tool for cellulases that have not been experimentally characterized, and we posit that it may also be applicable to other heterogeneously acting biocatalysts.

Development of peptide therapeutics generally involves screening of excipients that inhibit peptide-peptide interactions, hence aggregation, and improve peptide stability. We used the therapeutic peptide plectasin to develop a fast screening method that combines microscale thermophoresis titration assays and molecular dynamics simulations to relatively rank the excipients with respect to binding affinity and to study key peptide-excipient interaction hotspots on a molecular level, respectively. Additionally, 1H-13C-HSQC NMR titration experiments were performed to validate the fast screening approach. The NMR results are in qualitative agreement with results from the fast screening method demonstrating that this approach can be reliably applied to other peptides and proteins as a fast screening method to relatively rank excipients and predict possible excipient binding sites.

Background Cellobiohydrolase from glycoside hydrolase family 7 is a major component of commercial enzymatic mixtures for lignocellulosic biomass degradation. For many years, Trichoderma reesei Cel7A (TrCel7A) has served as a model to understand structure–function relationships of processive cellobiohydrolases. The architecture of TrCel7A includes an N-glycosylated catalytic domain, which is connected to a carbohydrate-binding module through a flexible, O-glycosylated linker. Depending on the fungal expression host, glycosylation can vary not only in glycoforms, but also in site occupancy, leading to a complex pattern of glycans, which can affect the enzyme’s stability and kinetics. Results Two expression hosts, Aspergillus oryzae and Trichoderma reesei, were utilized to successfully express wild-types TrCel7A (WTAo and WTTr) and the triple N-glycosylation site deficient mutants TrCel7A N45Q, N270Q, N384Q (ΔN-glycAo and ΔN-glycTr). Also, we expressed single N-glycosylation site deficient mutants TrCel7A (N45QAo, N270QAo, N384QAo). The TrCel7A enzymes were studied by steady-state kinetics under both substrate- and enzyme-saturating conditions using different cellulosic substrates. The Michaelis constant (KM) was consistently found to be lowered for the variants with reduced N-glycosylation content, and for the triple deficient mutants, it was less than half of the WTs’ value on some substrates. The ability of the enzyme to combine productively with sites on the cellulose surface followed a similar pattern on all tested substrates. Thus, site density (number of sites per gram cellulose) was 30–60% higher for the single deficient variants compared to the WT, and about twofold larger for the triple deficient enzyme. Molecular dynamic simulation of the N-glycan mutants TrCel7A revealed higher number of contacts between CD and cellulose crystal upon removal of glycans at position N45 and N384. Conclusions The kinetic changes of TrCel7A imposed by removal of N-linked glycans reflected modifications of substrate accessibility. The presence of N-glycans with extended structures increased KM and decreased attack site density of TrCel7A likely due to steric hindrance effect and distance between the enzyme and the cellulose surface, preventing the enzyme from achieving optimal conformation. This knowledge could be applied to modify enzyme glycosylation to engineer enzyme with higher activity on the insoluble substrates.

Scientific reports 10(1), 10089 (2020). doi:10.1038/s41598-020-67002-9, Fusion technology is widely used in protein-drug development to increase activity, stability, and bioavailability of protein therapeutics. Fusion proteins, like any other type of biopharmaceuticals, need to remain stable during production and storage. Due to the high complexity and additional intramolecular interactions, it is not possible to predict the behavior of fusion proteins based on the behavior the individual proteins. Therefore, understanding the stability of fusion proteins on the molecular level is crucial for the development of biopharmaceuticals. The current study on the albumin-neprilysin (HSA-NEP) fusion protein uses a combination of thermal and chemical unfolding with small angle X-ray scattering and molecular dynamics simulations to show a correlation between decreasing stability and increasing repulsive interactions, which is unusual for most biopharmaceuticals. It is also seen that HSA-NEP is not fully flexible: it is present in both compact and extended conformations. Additionally, the volume fraction of each conformation changes with pH. Finally, the presence of NaCl and arginine increases stability at pH 6.5, but decreases stability at pH 5.0., Published by Macmillan Publishers Limited, part of Springer Nature, [London]

Secretory phospholipases A2 (sPLA2s) are a subclass of enzymes that catalyze the hydrolysis at the sn-2 position of glycerophospholipids, producing free fatty acids and lysophospholipids. In this study, different phospholipids with structural modifications close to the scissile sn-2 ester bond were studied to determine the effect of the structural changes on the formation of the Michaelis-Menten complex and the water entry/exit pathways using molecular dynamics simulations and the computational tracking tool AQUA-DUCT. Structural modifications include methylation, dehydrogenation, and polarization close to the sn-2 scissile bond. We found that all water molecules reaching the active site of sPLA2-IIA pass by the aromatic residues Phe5 and Tyr51 and enter the active site through an active-site cleft. The relative amount of water available for the enzymatic reaction of the different phospholipid-sPLA2 complexes was determined together with the distance between key atoms in the catalytic machinery. The results showed that (Z)-unsaturated phospholipid is a good substrate for sPLA2-IIA. The computational results are in good agreement with previously reported experimental data on the ability of sPLA2-IIA to hydrolyze liposomes made from the different phospholipids, and the results provide new insights into the necessary active-site solvation of the Michaelis-Menten complex and can pave the road for rational design in engineering applications.

Enzymes active on glycosidic bonds are defined according to the stereochemistry of both substrates and products of the reactions they catalyze. The CAZy classification further assigns these enzymes into sequence-based families sharing a common stereochemistry for substrates (either α- or β-) and products (i.e., inverting or retaining mechanism). Here we describe the N-acetylgalactosaminidases AmGH109A and AmGH109B (i.e., GH109: glycoside hydrolase family 109) from the human gut symbiont Akkermansia muciniphila. Notably, AmGH109A displays α-retaining and β-inverting N-acetylgalactosaminidase activities with comparable efficiencies on natural disaccharides. This dual specificity could provide an advantage in targeting a broader range of host-derived glycans. We rationalize this discovery through bioinformatics, structural, mutational, and computational studies, unveiling a histidine residing in a conserved GGHGG motif as the elusive catalytic acid-base of the GH109 family. (Less)

Therapeutic protein candidates should exhibit favorable properties that render them suitable to become drugs. Nevertheless, there are no well-established guidelines for the efficient selection of proteinaceous molecules with desired features during early stage development. Such guidelines can emerge only from a large body of published research that employs orthogonal techniques to characterize therapeutic proteins in different formulations. In this work, we share a study on a diverse group of proteins, including their primary sequences, purity data, and computational and biophysical characterization at different pH and ionic strength. We report weak linear correlations between many of the biophysical parameters. We suggest that a stability comparison of diverse therapeutic protein candidates should be based on a computational and biophysical characterization in multiple formulation conditions, as the latter can largely determine whether a protein is above or below a certain stability threshold. We use the presented data set to calculate several stability risk scores obtained with an increasing level of analytical effort and show how they correlate with protein aggregation during storage. Our work highlights the importance of developing combined risk scores that can be used for early stage developability assessment. We suggest that such scores can have high prediction accuracy only when they are based on protein stability characterization in different solution conditions.

We recently reported the discovery of a novel protein stabilizing dipeptide, glycyl-D-asparagine, through a structure-based approach. As the starting hypothesis leading to the discovery, we postulated a stabilizing effect achieved by binding of the dipeptide to an aggregation prone region on the protein's surface. Here we present a detailed study of the interaction mechanism between the dipeptide and Interferon-alpha-2A (IFN) through the construction of a Markov state model from molecular dynamics trajectories. We identify multiple binding sites and compare these to aggregation prone regions. Additionally, we calculate the lifetime of the protein-excipient complex. If the excipient remained bound to IFN after administration, it could alter the protein's therapeutic efficacy. We establish that the lifetime of the complex between IFN and glycyl-D-asparagine is extremely short. Under these circumstances, stabilization by stoichiometric binding is consequently no impediment for a safe use of an excipient.

Water has a critical role in the stability of the higher-order structure of proteins. In addition, it is considered to be a major destabilization factor for the physical and chemical stability of freeze-dried proteins and peptides. Physical and chemical aspects of protein/water relationships are commonly studied with the use of water vapor sorption isotherms for amorphous lyophilized proteins, which, in turn, are commonly analyzed using the Brunauer-Emmett-Teller (BET) equation to obtain the parameters, Wm and CB. The parameter Wm is generally referred to as the "monolayer limit of adsorption" and has a narrow range of 6-8% for most proteins. In this study, the water distribution on an IgG1 surface is investigated by molecular dynamics (MD) simulations at different water contents. The monolayer of water molecules was found to have limited coverage of the protein surface, and the true monolayer coverage of the protein globule actually occurs at a hydration level above 30%. The distribution of water molecules on the IgG1 surface is also highly heterogeneous, and the heterogeneity is not considered in the BET theory. In this study, a mechanistic model has been developed to describe the water vapor sorption isotherm. This model is based on the analysis of the hydrogen bonding network extracted from the MD simulations. The model is consistent with the experimental Type-II isotherm, which is usually observed for proteins. The physical meaning of the BET monolayer was redefined as the onset of water cluster formation. A simple model to calculate the onset water level, Wm, is proposed based on the hydration of different amino acids, as determined from the MD simulations.

We describe the synthesis of tetrahydroisoquinolines and tetrahydroisoquinolinium salts together with their pharmacological properties at various nicotinic acetylcholine receptors. In general, the compounds were α4β2 nAChR antagonists, with the tetrahydroisoquinolinium salts being more potent than the parent tetrahydroisoquinoline derivatives. The most potent α4β2 antagonist, 6c, exhibited submicromolar binding Ki and functional IC50 values and high selectivity for this receptor over the α4β4 and α3β4 nAChRs. Whereas the (S)-6c enantiomer was essentially inactive at α4β2, (R)-6c was a slightly more potent α4β2 antagonist than the reference β2-nAChR antagonist DHβΕ. The observation that the α4β2 activity resided exclusively in the (R)-enantiomer was in full agreement with docking studies. Several of tetrahydroisoquinolinium salts also displayed agonist activity at the α7 nAChR. Preliminary in vivo evaluation revealed antidepressant-like effects of both (R)-5c and (R)-6c in the mouse forced swim test, supporting the therapeutic potential of α4β2 nAChR antagonists for this indication.

The molecular mechanisms underlying the drug-solubilizing properties of γ-cyclodextrin were explored using hydrocortisone as a model drug. The BS-type phase-solubility diagram of hydrocortisone with Υ-cyclodextrin was thoroughly characterized by measuring the concentrations of hydrocortisone and Υ-cyclodextrin in solution and the solid phase. The drug-solubilizer interaction was also studied by isothermal titration calorimetry from which a precise value of the 1:1 binding constant (K 11 =4.01mM-1 at 20°C) was obtained. The formation of water-soluble 1:1 complexes is responsible for the initial increase in hydrocortisone solubility while the precipitation of entities with a 3:2 ratio of Υ-cyclodextrin:hydrocortisone is responsible for the plateau and the ensuing strong decrease in solubility once all solid hydrocortisone is used up. The complete phase-solubility diagram is well accounted for by a model employing the 1:1 binding constant and the solubility product of the precipitating 3:2 entity (K32S=5.51 mM5). For such systems, a small surplus of Υ-cyclodextrin above the optimum concentration may result in a significant decrease in drug solubility, and the implications for drug formulations are briefly discussed.

Phospholipases are ubiquitous in nature and the target of significant research aiming at both their physiological roles and technical applications in e.g. the food industry. In the search for sensitive and selective phospholipase assays, we have focused on synthetic FRET (Forster resonance energy transfer) substrates. This has led to the development of a facile, easily scalable and low cost synthesis of fluorogenic phospholipids featuring the dansyl/dabcyl fluorophore/quencher-pair on the fatty acid ω-position and on the phosphatidylethanolamine head group, respectively. Hence, the two substrates lyso-(dansyl-FA)-GPE-dabcyl (6) and (dansyl-FA)2-GPE-dabcyl (7) were synthesized by a chemoenzymatic strategy, in which preparation of (6) further included a novel selective enzymatic esterification step. As proof of concept, activity of a handful of phospholipases, one from each of the PLA1, PLA2, PLC and PLD classes, were assayed using substrates (6) and (7), and the kinetic parameter kcat/KM was determined. The PLA1 (Lecitase Ultra™) was found to be highly active on both substrates, whereas the PLD (from white cabbage) had no activity, presumably due to steric effects associated with the dabcyl-functionalization of the head group. It was further substantiated that the substrates are specific towards phospholipase activity as the tested lipase (Lipolase™) showed close to zero activity.

Graphical abstract, Highlights • Combination of SAXS and MD simulations can explain the molecular background for the stability studies. • Destabilizing excipients in general induce opening of the N-lobe by binding to specific residues in the protein. • The fully open conformation of transferrin is not seen without aggregation., Transferrin is an attractive candidate for drug delivery due to its ability to cross the blood brain barrier. However, in order to be able to use it for therapeutic purposes, it is important to investigate how its stability depends on different formulation conditions. Combining high-throughput thermal and chemical denaturation studies with small angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations, it was possible to connect the stability of transferrin with its conformational changes. Lowering pH induces opening of the transferrin N-lobe, which results in a negative effect on the stability. Presence of NaCl or arginine at low pH enhances the opening and has a negative impact on the overall protein stability. Statement of Significance Protein-based therapeutics have become an essential part of medical treatment. They are highly specific, have high affinity and fewer off-target effects. However, stabilization of proteins is critical, time-consuming, and expensive, and it is not yet possible to predict the behavior of proteins under different conditions. The current work is focused on a molecular understanding of the stability of human serum transferrin; a protein which is abundant in blood serum, may pass the blood brain barrier and therefore with high potential in drug delivery. Combination of high throughput unfolding techniques and structural studies, using small angle X-ray scattering and molecular dynamic simulations, allows us to understand the behavior of transferrin on a molecular level.

Plectasin is a small, cysteine-rich peptide antibiotic which belongs to the class of antimicrobial peptides and has potential antibacterial activity against various Gram-positive bacteria. In the current study, the effect of pH and ionic strength (NaCl) on the conformational stability of plectasin variants has been investigated. At all physiochemical conditions, peptide secondary structures are intact throughout simulations. However, flexibility increases with pH because of the change in electrostatics around the distinct anionic tetrapeptide (9-12) stretch. Furthermore, plectasin interactions with NaCl were measured by determining the preferential interaction coefficients, Γ23. Generally, wild-type plectasin has higher preference for sodium ions as 9ASP is mutated in other variants. Overall, the Γ23 trend with pH for the two salt conditions remain the same for all variants predominately having accumulation of sodium ions around 10GLU and 12ASP. Insignificant changes in the overall peptide conformational stability are in agreement with the fact that plectasin has three cystines. Thermodynamic integration molecular dynamics simulations supplemented with nuclear magnetic resonance were employed to determine the degree of involvement of three different cystines to the overall structural integrity of the peptide. Both methods show the same order of cystine reduction and complete unfolding is observed only upon reduction of all cystines.

An acoustically levitated droplet has been used to collect synchrotron SAXS data on human serum albumin protein solutions up to a protein concentration of 400 mg ml−1. A careful selection of experiments allows for fast data collection of a large amount of data, spanning a protein concentration/solvent concentration space with limited sample consumption (down to 3 µL per experiment) and few measurements. The data analysis shows data of high quality that are reproducible and comparable with data from standard flow-through capillary-based experiments. Furthermore, using this methodology, it is possible to achieve concentrations that would not be accessible by conventional cells. The protein concentration and ionic strength parameter space diagram may be covered easily and the amount of protein sample is significantly reduced (by a factor of 100 in this work). Used in routine measurements, the benefits in terms of protein cost and time spent are very significant.

Alginate is a linear polysaccharide from brown algae consisting of 1,4-linked β-D-mannuronic acid (M) and α-L-guluronic acid (G) arranged in M, G, and mixed MG blocks. Alginate was assumed to be indigestible in humans, but bacteria isolated from fecal samples can utilize alginate. Moreover, genomes of some human gut microbiome–associated bacteria encode putative alginate-degrading enzymes. Here, we genome-mined a polysaccharide lyase family 6 alginate lyase from the gut bacterium Bacteroides cellulosilyticus (BcelPL6). The structure of recombinant BcelPL6 was solved by X-ray crystallography to 1.3 Å resolution, revealing a single-domain, monomeric parallel β-helix containing a 10-step asparagine ladder characteristic of alginate-converting parallel β-helix enzymes. Substitutions of the conserved catalytic site residues Lys-249, Arg-270, and His-271 resulted in activity loss. However, imidazole restored the activity of BcelPL6-H271N to 2.5% that of the native enzyme. Molecular docking oriented tetra-mannuronic acid for syn attack correlated with M specificity. Using biochemical analyses, we found that BcelPL6 initially releases unsaturated oligosaccharides of a degree of polymerization of 2–7 from alginate and polyM, which were further degraded to di- and trisaccharides. Unlike other PL6 members, BcelPL6 had low activity on polyMG and none on polyG. Surprisingly, polyG increased BcelPL6 activity on alginate 7-fold. LC–electrospray ionization–MS quantification of products and lack of activity on NaBH4-reduced octa-mannuronic acid indicated that BcelPL6 is an endolyase that further degrades the oligosaccharide products with an intact reducing end. We anticipate that our results advance predictions of the specificity and mode of action of PL6 enzymes. Author’s Choice—Final version open access under the terms of the Creative Commons CC-BY license

The development of a new protein drug typically starts with the design, expression and biophysical characterization of many different protein constructs. The initially high number of constructs is radically reduced to a few candidates that exhibit the desired biological and physicochemical properties. This process of protein expression and characterization to find the most promising molecules is both expensive and time-consuming. Consequently, many companies adopt and implement philosophies, e.g. platforms for protein expression and formulation, computational approaches, machine learning, to save resources and facilitate protein drug development. Inspired by this, we propose the use of interpretable artificial neuronal networks (ANNs) to predict biophysical properties of therapeutic monoclonal antibodies i.e. melting temperature Tm, aggregation onset temperature Tagg, interaction parameter kD as a function of pH and salt concentration from the amino acid composition. Our ANNs were trained with typical early-stage screening datasets achieving high prediction accuracy. By only using the amino acid composition, we could keep the ANNs simple which allows for high general applicability, robustness and interpretability. Finally, we propose a novel "knowledge transfer" approach, which can be readily applied due to the simple algorithm design, to understand how our ANNs come to their conclusions.

Reducing the aggregation of proteins is of utmost interest to the pharmaceutical industry. Aggregated proteins are often less active and can cause severe immune reactions in the patient upon administration. At the same time the biopharmaceutical market is pushing for high concentration formulations and products that do not require refrigerated storage conditions. For a given protein, the only solution pH, ionic strength and concentration of a very limited number of excipients are the only parameters that can be varied to obtain a stable formulation. In this work, we present a structure-based approach to discover new molecules that successfully reduce the aggregation of proteins and apply the approach to the model protein Interferon-alpha-2a.

The ability of enzymes to operate in organic solvent is now widely accepted and is the basis for extensive research in enzymology. The challenge is to select the solvent media that allows the modulation of enzyme activity. For a rational selection of a solvent, it is necessary to understand the effect of organic solvent molecules on enzyme structure and the enzymatic reaction on a molecular level. To gain such insight, we combined experimental kinetics studies with full atomic molecular dynamics simulations and found a correlation between the activity of Candida antarctica lipase B (CALB) [for the esterification reaction between butyric acid and ethanol at a fixed water activity] and the binding of the solvent/substrate molecules in the active site region of CALB. We have investigated the influence of four organic solvents—hexane (HEX), methyl tertiary butyl ether (MTBE), acetonitrile (ACN), and tertiary butanol (TBU)—on the catalytic activity of CALB for the esterification reaction. The solvents have bee...

The lipidation of peptide drugs is one strategy to obtain extended half-lives, enabling once-daily or even less frequent injections for patients. The half-life extension results from a combination of self-association and association with human serum albumin (albumin). The self-association and association with albumin of two insulin analogues, insulin detemir and insulin degludec, were investigated by small-angle X-ray scattering (SAXS) and dynamic light scattering (DLS) in phenolic buffers. Detemir shows concentration-dependent self-association, with an equilibrium between hexamer, dihexamer, trihexamer and larger species, while degludec appears as a dihexamer independent of concentration. The solution structure of the detemir trihexamer has a bent shape. The stoichiometry of the association with albumin was studied using DLS. For albumin–detemir the molar stoichiometry was determined to be 1:6 (albumin:detemir ratio) and for albumin–degludec it was between 1:6 and 1:12 (albumin:degludec ratio). Batch SAXS measurements of a 1:6 albumin:detemir concentration series revealed a concentration dependence of complex formation. The data allowed the modelling of a complex between albumin and a detemir hexamer and a complex consisting of two albumins binding to opposite ends of a detemir dihexamer. Measurements of size-exclusion chromatography coupled to SAXS revealed a complex between a degludec dihexamer and albumin. Based on the results, equilibria for the albumin–detemir and albumin–degludec mixtures are proposed.

Despite a very large number of bacterial exopolysaccharides have been reported, detailed knowledge on their molecular structures and associative interactions with proteins is lacking. Small-angle X-ray scattering, dynamic light scattering and analytical ultracentrifugation (AUC) were used to characterize the interactions of six lactic acid bacterial heteroexopolysaccharides (HePS-1–HePS-6) with β-lactoglobulin (BLG). Compared to free HePSs, a large increase in the X-ray radius of gyration RG, maximum length L and hydrodynamic diameter dH of HePS-1–HePS-4 mixed with BLG revealed strong aggregation, the extent of which depended on the compact conformation and degree of branching of these HePSs. No significant effects were observed with HePS-5 and HePS-6. Turbidity and AUC analyses showed that both soluble and insoluble BLG–HePS complexes were formed. The findings provide new insights into the role of molecular structures in associative interactions between HePSs and BLG which has relevance for various industrial applications.

In the present study, small-angle X-ray scattering (SAXS) and static light scattering (SLS) have been used to study the solution properties and self-interaction of recombinant human serum albumin (rHSA) molecules in three pharmaceutically relevant buffer systems. Measurements are carried out up to high protein concentrations and as a function of ionic strength by adding sodium chloride to probe the role of electrostatic interactions. The effective structure factors (Seff) as a function of the scattering vector magnitude q have been extracted from the scattering profiles and fit to the solution of the Ornstein-Zernike equation using a screened Yukawa potential to describe the double-layer force. Although only a limited q range is used, accurate fits required including an electrostatic repulsion element in the model at low ionic strength, while only a hard sphere model with a tunable diameter is necessary for fitting to high-ionic-strength data. The fit values of net charge agree with available data from potentiometric titrations. Osmotic compressibility data obtained by extrapolating the SAXS profiles or directly from SLS measurements has been fit to a 10-term virial expansion for hard spheres and an equation of state for hard biaxial ellipsoids. We show that modeling rHSA as an ellipsoid, rather than a sphere, provides a much more accurate fit for the thermodynamic data over the entire concentration range. Osmotic virial coefficient data, derived at low protein concentration, can be used to parameterize the model for predicting the behavior up to concentrations as high as 450 g/L. The findings are especially important for the biopharmaceutical sector, which require approaches for predicting concentrated protein solution behavior using minimal sample consumption.

A recently developed membrane-mimetic model was applied to study membrane interaction and binding of the two anchoring C2-like discoidin domains of human coagulation factor VIIIa (FVIIIa), the C1 and C2 domains. Both individual domains, FVIII C1 and FVIII C2, were observed to bind the phospholipid membrane by partial or full insertion of their extruding loops (the spikes). However, the two domains adopted different molecular orientations in their membrane-bound states; FVIII C2 roughly was positioned normal to the membrane plane, while FVIII C1 displayed a multitude of tilted orientations. The results indicate that FVIII C1 may be important in modulating the orientation of the FVIIIa molecule to optimize the interaction with FIXa, which is anchored to the membrane via its γ-carboxyglutamic acid-rich (Gla) domain. Additionally, a structural change was observed in FVIII C1 in the coiled main chain leading the first spike. A tight interaction with one lipid per domain, similar to what has been suggested for the homologous FVa C2, is characterized. Finally, we rationalize known FVIII antibody epitopes and the scarcity of documented hemophilic missense mutations related to improper membrane binding of FVIIIa, based on the prevalent nonspecificity of ionic interactions in the simulated membrane-bound states of FVIII C1 and FVIII C2.

The glucagon-like peptide 1 (GLP-1) analog, liraglutide, is a GLP-1 agonist and is used in the treatment of type-2 diabetes mellitus and obesity. From a pharmaceutical perspective, it is important to know the oligomerization state of liraglutide with respect to stability. Compared to GLP-1, liraglutide has an added fatty acid (FA) moiety that causes oligomerization of liraglutide as suggested by small-angle x-ray scattering (SAXS) and multiangle static light scattering (MALS) results. SAXS data suggested a global shape of a hollow elliptical cylinder of size hexa-, hepta-, or octamer, whereas MALS data indicate a hexamer. To elaborate further on the stability of these oligomers and the role of the FA chains, a series of molecular-dynamics simulations were carried out on 11 different hexa-, hepta-, and octameric systems. Our results indicate that interactions of the fatty acid chains contribute noticeably to the stabilization. The simulation results indicate that the heptamer with paired FA chains is the most stable oligomer when compared to the 10 other investigated structures. Theoretical SAXS curves extracted from the simulations qualitatively agree with the experimentally determined SAXS curves supporting the view that liraglutide forms heptamers in solution. In agreement with the SAXS data, the heptamer forms a water-filled oligomer of elliptical cylindrical shape.

We have performed a series of molecular dynamics simulations to study the interactions between the neurotransmitters (NTs) γ-aminobutyrate (GABA), glycine (GLY), acetylcholine (ACH) and glutamate (GLU) as well as the amidated/acetylated γ-aminobutyrate (GABA(neu)) and the osmolyte molecule glycerol (GOL) with a dipalmitoylphosphatidylcholine (DPPC) bilayer. In agreement with previously published experimental data, we found the lowest membrane affinity for the charged molecules and a moderate affinity for zwitterionic and polar molecules. The affinity can be ranked as follows: ACH-GLU<

Various cellulases consist of a catalytic domain connected to a carbohydrate-binding module (CBM) by a flexible linker peptide. The linker if often strongly O-glycosylated and typically has a length of 20-50 amino acid residues. Functional roles, other than connecting the two folded domains, of the linker and its glycans, have been widely discussed, but experimental evidence remains sparse. One of the most studied cellulose degrading enzymes is the multi-domain cellobiohydrolase Cel7A from Hypocrea jecorina. Here, we designed variants of Cel7A with mutations in the linker region to elucidate the role of the linker. We found that moderate modification of the linker could result in significant changes in substrate affinity and catalytic efficacy. These changes were quite different for different linker variants. Thus, deletion of six residues near the catalytic domain had essentially no effects on enzyme function. Conversely, a substitution of four glycosylation sites near the middle of the linker reduced substrate affinity and increased maximal turnover. The observation of weaker binding provides some support of recent suggestions that linker glycans may be directly involved in substrate interactions. However, a variant with several inserted glycosylation sites near the CBM also showed lower affinity for the substrate compared to the wild-type, and we suggest that substrate interactions of the glycans depend on their exact location as well as other factors such as changes in structure and dynamics of the linker peptide.

Molecular structures of exopolysaccharides are required to understand their functions and the relationships between the structure and physical and rheological properties. Small-angle X-ray scattering and dynamic light scattering were used in conjunction with molecular modeling to characterize solution structures of three lactic acid bacterial heteroexopolysaccharides (HePS-1, HePS-2, and HePS-3). Values of radius of gyration RG, cross-sectional radius of gyration Rxs, approximate length L, and hydrodynamic diameter were not directly proportional to the molar mass and indicated the HePSs adopted a compact coil-like rather than an extended conformation. Constrained molecular modeling of 15000 randomized HePS-1 conformers resulted in five best-fit structures with R factor of 3.9-4.6% revealing random coil-like structure. φ and ψ angle analysis of glycosidic linkages in HePS-1 structures suggests Galf residues significantly influence the conformation. Ab initio scattering modeling of HePS-2 and HePS-3 gave excellent curve fittings with χ2 of 0.43 and 0.34 for best-fit models, respectively, compatible with coil-like conformation. The findings disclose solution behavior of HePS relevant for their interactions with biomacromolecules, for example, milk proteins.

Tryptophan hydroxylase (TPH) catalyzes the initial and rate-limiting step in the biosynthesis of serotonin, which is associated with a variety of disorders such as depression and irritable bowel syndrome. TPH exists in two isoforms: TPH1 and TPH2. TPH1 catalyzes the initial step in the synthesis of serotonin in the peripheral tissues, while TPH2 catalyzes this step in the brain. In this study, the steady-state kinetic mechanism for the catalytic domain of human TPH1 has been determined. Varying substrate tryptophan (Trp) and tetrahydrobiopterin (BH4) results in a hybrid Ping Pong-ordered mechanism in which the reaction can either occur through a Ping Pong or a sequential mechanism depending on the concentration of tryptophan. The catalytic domain of TPH1 shares a sequence identity of 81% with TPH2. Despite the high sequence identity, differences in the kinetic parameters of the isoforms have been identified; i.e., only TPH1 displays substrate tryptophan inhibition. This study demonstrates that the difference can be traced to an active site loop which displays different properties in the TPH isoforms. Steady-state kinetic results of the isoforms, and variants with point mutations in a loop lining the active site, show that the kinetic parameters of only TPH1 are significantly changed upon mutations. Mutations in the active site loop of TPH1 result in an increase in the substrate inhibition constant, Ki, and therefore turnover rate. Molecular dynamics simulations reveal that this substrate inhibition mechanism occurs through a closure of the cosubstrate, BH4, binding pocket, which is induced by Trp binding.

The glucagon-like peptide-1 receptor (GLP-1R) belongs to family B of the G-protein coupled receptors (GPCRs), and has become a promising target for the treatment of type 2 diabetes. Here we describe the development and characterization of a fully functional cysteine-deprived and C-terminally truncated GLP-1R. Single cysteines were initially substituted with alanine, and functionally redundant cysteines were subsequently changed simultaneously. Our results indicate that Cys(174), Cys(226), Cys(296) and Cys(403) are important for the GLP-1-mediated response, whereas Cys(236), Cys(329), Cys(341), Cys(347), Cys(438), Cys(458) and Cys(462) are not. Extensive deletions were made in the C-terminal tail of GLP-1R in order to determine the limit for truncation. As for other family B GPCRs, we observed a direct correlation between the length of the C-terminal tail and specific binding of (125)I-GLP-1, indicating that the membrane proximal part of the C-terminal is involved in receptor expression at the cell surface. The results show that seven cysteines and more than half of the C-terminal tail can be removed from GLP-1R without compromising GLP-1 binding or function.

The formation of inclusion complexes of hydroxypropylated β-cyclodextrins (CDs) with three bile salts are investigated to shed light on the role played by the hydroxypropyl (HP) substituents. The HP-chains are situated at the rim of the CD and may thus extend the hydrophobic cavity of the CD. Calorimetric titrations in a broad temperature range and molecular dynamics simulations confirm previous speculations that the HP-chains cause an increase in dehydrated nonpolar surface area upon formation of the complexes. This additional burial of nonpolar surface area, 12–16 A2 per HP-chain according to the MD simulations, results in more negative values of ΔC p °, which are in quantitative agreement with what is expected for hydrophobic dehydration. Although these observations support the picture of an extended hydrophobic cavity, HPβCD complexes were less stable than their unsubstituted counterparts. This indicates that increased hydrophobic contacts are not always accompanied by increased binding strength. The linear dependence of ΔC p °, ΔH° and ΔS° on the number of HP-chains give rise to isoentropic and isoenthalpic temperatures at which ΔH° and ΔS° are independent of the number of HP-chains on the host CD (but depend on the type of bile salt). Interestingly, these convergence temperatures are close to what is observed for unfolding of proteins and may be a common feature of hydrophobic dehydration.

The aim of the present work was to investigate the interaction between bile salts present in the intestine of man, dog and rat with the negatively charged cyclodextrin (CD), sulfobutylether-β-cyclodextrin (SBEβCD). The interactions between bile salts and CDs are of importance for the release of CD-complexed drugs upon oral administration. This makes a good understanding of this particular interaction important for rational drug formulation. SBEβCD is a modified CD, which has attracted particular interest in formulation science. It is unique in the sense that it carries approximately seven negatively charged side chains, which can potentially interact electrostatically with the guest molecule. Bile salts are negatively charged at physiological pH, and the concomitant repulsion from SBEβCD could potentially reduce their affinity for this CD and hence their ability to expel drugs delivered as SBEβCD complexes. However, this study has demonstrated that the interaction, between a bile salts and SBEβCD is only slightly weaker than the corresponding interactions with natural βCD. Significant differences between the thermodynamics of bile salt complexes with respectively HPβCD and SBEβCD were found, when comparing the same degree of substitution. This underscores the importance of the substituents on the interactions of modified CDs with bile salts.

The bio-protective properties of monosaccharaides, namely mannose, fructose and fucose, on the stability and dynamical properties of the NMR determined hen egg-white lysozyme structure have been investigated by means of molecular dynamics simulations at room temperature in aqueous solution and in 7 and 13 wt % concentrations of the three sugars. Results are discussed in the framework of the bio-protective phenomena. The three sugars show similar bio-protective behaviours at room temperature (300 K) in the concentration range studied as shown by the small RMSDs of the resulting MD structures from that of starting NMR structure. The effects of sugars on protein conformation are found to be relatively strong in that the conformation of lysozyme is stable after an initial 9 ns equilibration for fucose and mannose and 12 ns equilibration for fructose, respectively, at high concentrations. For mannose the final RMSD is significantly smaller than that of fucose and fructose at the higher concentration, while at the lower concentration the RMSD are essentially the same. The radial distribution function of the water and sugars around lysozyme was used to monitor the preferential hydration. Analysis of the solvent and sugar distributions around lysozyme was used to investigate the interfacial solvent and sugar structure near the protein surface.

Fluorinated phospholipid analogues are investigated as potential substrates for phospholipase A2 (PLA2) using classical molecular dynamics simulations and quantum mechanics/density functional theory calculations. The fluorinated phospholipid analogues are α-fluoro (HF-ProAEL) and α,α-difluoro (F2-ProAEL) conjugates of (R)-1-O-hexadecyl-2-palmitoyl-sn-glycero-3-phoshocholineglycerol (ProAEL). Our results provide a theoretical assessment of the potential usefulness of these fluorinated lipids in the rational design of liposomal drug-delivery systems. The α-fluorine-substituted phospholipid analogues are found to be substrates for secretory PLA2, with sufficient accessibility of water to the active site to allow for enzymatic hydrolysis. Because of the inherently less stable nature of HF-ProAEL and F2-ProAEL when compared to that of ProAEL, the hydrolytic reaction is predicted to occur at a progressively faster rate; the more electronegative substituent at the α-position effectively lowers the energy barrier for hydrolysis. We conclude that the partially fluorinated phospholipid analogues facilitate rational design of liposomal vesicles of phospholipid mixtures with desirable physicochemical properties and that are still subjects for important and pharmaceutically proven drug-delivery mechanisms.

Tryptophan hydroxylase 2 (TPH2) catalyses the initial and rate-limiting step in the biosynthesis of serotonin, which is associated with a variety of disorders such as depression, obsessive compulsive disorder, and schizophrenia. Full-length TPH2 is poorly characterized due to low purification quantities caused by its inherent instability. Three truncated variants of human TPH2 (rch TPH2; regulatory and catalytic domain, NΔ47-rch TPH2; truncation of 47 residues in the N terminus of rch TPH2, and ch TPH2; catalytic domain) were expressed, purified, and examined for changes in transition temperature, inactivation rate, and oligomeric state. ch TPH2 displayed 14- and 11-fold higher half-lives compared to rch TPH2 and NΔ47-rch TPH2, respectively. Differential scanning calorimetry experiments demonstrated that this is caused by premature unfolding of the less stable regulatory domain. By differential scanning fluorimetry, the unfolding transitions of rch TPH2 and NΔ47-rch TPH2 are found to shift from polyphasic to apparent two-state by the addition of l-Trp or l-Phe. Analytical gel filtration revealed that rch TPH2 and NΔ47-rch TPH2 reside in a monomer-dimer equilibrium which is significantly shifted toward dimer in the presence of l-Phe. The dimerizing effect induced by l-Phe is accompanied by a stabilizing effect, which resulted in a threefold increase in half-lives of rch TPH2 and NΔ47-rch TPH2. Addition of l-Phe to the purification buffer significantly increases the purification yields, which will facilitate characterization of hTPH2.

Magnesium, the second and fourth most abundant cation in the intracellular compartment and whole body, respectively, is of great physiologic importance. Magnesium exists as bound and free ionized forms depending on temperature, pH, ionic strength, and competing ions. Free magnesium participates in many biochemical processes and is most commonly measured by ion-selective electrode. This analytical approach is problematic because complete selectivity is not possible due to competition with other ions, i.e., calcium, and pH interference. Unfortunately, many studies have focused on measurement of total magnesium rather than its free bioactive form making it difficult to correlate to disease states. This systematic literature review presents current analytical challenges in obtaining accurate and reproducible test results for magnesium.

The stability of cyclic peptide assemblies (CPs) forming a macromolecular nanotube structure was investigated in solvents of different polarity using computational methods. The stability and structure of the complexes were studied using traditional molecular dynamics (MD). Energy of dissociation was estimated from steered MD in combination with umbrella sampling simulations. A cyclic peptide nanotube (CPNT) was constructed by stacking of eight cyclo[(d-Trp-l-Gln-d-Trp-l-Glu)2], and hereafter is referred to as (WQWE)8. Its dissociation was studied by pulling 1, 2, or 3 subunits from the nanotube. The crucial point in the dissociation event of the CP subunit(s) is the breaking of backbone–backbone hydrogen bonds and consecutive annihilation of side chain interactions. Gibbs free energy calculations to estimate the binding affinity of CP subunit(s) reveal that the (WQWE)8 nanotube is significantly more stable in non-polar environments than in polar environments. The presently investigated nanotube, (WQWE)8, displays a higher stability in polar solvent than the previously studied nanotube, (QAEA)8. It appears that tryptophan contributes favorable to the improved stability by forming side chain–side chain hydrogen bonds.

The possibility for obtaining derivative properties for mixtures from integrals of spatial total and direct correlation functions obtained from molecular dynamics simulations is explored. Theoretically well-supported methods are examined to extend simulation radial distribution functions to long range so the integrals can converge. A previously published method developed for pure atomic fluids is here extended to handle simulations of molecular mixtures using all-atom force fields. We first test the method on simulations of Lennard-Jones/Stockmayer mixtures and show that that the results are consistent with an excess Helmholtz energy model fitted to available simulations. In addition, simulations of water/methanol and water/t-butanol mixtures have been carried out. The method yields results for partial molar volumes, activity coefficient derivatives, and individual correlation function integrals in reasonable agreement with smoothed experimental data. The proposed method for obtaining correlation function integrals is shown to perform at least as well as or better than two previously published approaches.

Structural derivatives of 4-MTA, an illegal amphetamine analogue have been previously shown to have anticancer effects in vitro. In this study we report the synthesis of a series of novel 1,3-bis(aryl)-2-nitro-1-propene derivatives related in structure to 4-MTA. A number of these compounds containing a classic nitrostyrene structure are shown to have antiproliferative activities in vitro in a range of malignant cell lines, particularly against Burkitt’s lymphoma derived cell lines, whilst having no effect on ‘normal’ peripheral blood mononuclear cells. Such effects appear to be independent of the serotonin transporter, a high affinity target for amphetamines and independent of protein tyrosine phosphatases and tubulin dynamics both of which have been previously associated with nitrostyrene-induced cell death. We demonstrate that a number of these compounds induce caspase activation, PARP cleavage, chromatin condensation and membrane blebbing in a Burkitt’s lymphoma derived cell line, consistent with these compounds inducing apoptosis in vitro. Although no specific target has yet been identified for the action of these compounds, the cell death elicited is potent, selective and worthy of further investigation.