Bioprocess in-line monitoring and control using Raman spectroscopy and Indirect Hard Modeling (IHM).

Process in-line monitoring and control are crucial to optimize the productivity of bioprocesses. A frequently applied Process Analytical Technology (PAT) tool for bioprocess in-line monitoring is Raman spectroscopy. However, evaluating bioprocess Raman spectra is complex and calibrating state-of-the-art statistical evaluation models is effortful. To overcome this challenge, we developed an Indirect Hard Modeling (IHM) prediction model in a previous study. The combination of Raman spectroscopy and the IHM prediction model enables non-invasive in-line monitoring of glucose and ethanol mass fractions during yeast fermentations with significantly less calibration effort than comparable approaches based on statistical models. In this study, we advance this IHM-based approach and successfully demonstrate that the combination of Raman spectroscopy and IHM is capable of not only bioprocess monitoring but also bioprocess control. For this purpose, we used this combination's in-line information as input of a simple on-off glucose controller to control the glucose mass fraction in Saccharomyces cerevisiae fermentations. When we performed two of these fermentations with different predefined glucose set points, we achieved similar process control quality as approaches using statistical models, despite considerably smaller calibration effort. Therefore, this study reaffirms that the combination of Raman spectroscopy and IHM is a powerful PAT tool for bioprocesses.

Interdisciplinary development of an overall process concept from glucose to 4,5-dimethyl-1,3-dioxolane via 2,3-butanediol.

To reduce carbon dioxide emissions, carbon-neutral fuels have recently gained renewed attention. Here we show the development and evaluation of process routes for the production of such a fuel, the cyclic acetal 4,5-dimethyl-1,3-dioxolane, from glucose via 2,3-butanediol. The selected process routes are based on the sequential use of microbes, enzymes and chemo-catalysts in order to exploit the full potential of the different catalyst systems through a tailor-made combination. The catalysts (microbes, enzymes, chemo-catalysts) and the reaction medium selected for each conversion step are key factors in the development of the respective production methods. The production of the intermediate 2,3-butanediol by combined microbial and enzyme catalysis is compared to the conventional microbial route from glucose in terms of specific energy demand and overall yield, with the conventional route remaining more efficient. In order to be competitive with current 2,3-butanediol production, the key performance indicator, enzyme stability to high aldehyde concentrations, needs to be increased. The target value for the enzyme stability is an acetaldehyde concentration of 600 mM, which is higher than the current maximum concentration (200 mM) by a factor of three.

The Good pH probe: non-invasive pH in-line monitoring using Good buffers and Raman spectroscopy.

In bioprocesses, the pH value is a critical process parameter that requires monitoring and control. For pH monitoring, potentiometric methods such as pH electrodes are state of the art. However, they are invasive and show measurement value drift. Spectroscopic pH monitoring is a non-invasive alternative to potentiometric methods avoiding this measurement value drift. In this study, we developed the Good pH probe, which is an approach for spectroscopic pH monitoring in bioprocesses with an effective working range between pH 6 and pH 8 that does not require the estimation of activity coefficients. The Good pH probe combines for the first time the Good buffer 3-(N-morpholino)propanesulfonic acid (MOPS) as pH indicator with Raman spectroscopy as spectroscopic technique, and Indirect Hard Modeling (IHM) for the spectral evaluation. During a detailed characterization, we proved that the Good pH probe is reversible, exhibits no temperature dependence between 15 and 40 °C, has low sensitivity to the ionic strength up to 1100 mM, and is applicable in more complex systems, in which other components significantly superimpose the spectral features of MOPS. Finally, the Good pH probe was successfully used for non-invasive pH in-line monitoring during an industrially relevant enzyme-catalyzed reaction with a root mean square error of prediction (RMSEP) of 0.04 pH levels. Thus, the Good pH probe extends the list of critical process parameters monitorable using Raman spectroscopy and IHM by the pH value.

Bioprocess in-line monitoring using Raman spectroscopy and Indirect Hard Modeling (IHM): A simple calibration yields a robust model.

To increase the process productivity and product quality of bioprocesses, the in-line monitoring of critical process parameters is highly important. For monitoring substrate, metabolite, and product concentrations, Raman spectroscopy is a commonly used Process Analytical Technology (PAT) tool that can be applied in-situ and non-invasively. However, evaluating bioprocess Raman spectra with a robust state-of-the-art statistical model requires effortful model calibration. In the present study, we in-line monitored a glucose to ethanol fermentation by Saccharomyces cerevisiae (S. cerevisiae) using Raman spectroscopy in combination with the physics-based Indirect Hard Modeling (IHM) and showed successfully that IHM is an alternative to statistical models with significantly lower calibration effort. The IHM prediction model was developed and calibrated with only 16 Raman spectra in total, which did not include any process spectra. Nevertheless, IHM's root mean square errors of prediction (RMSEPs) for glucose (3.68 g/L) and ethanol (1.69 g/L) were comparable to the prediction quality of similar studies that used statistical models calibrated with several calibration batches. Despite our simple calibration, we succeeded in developing a robust model for evaluating bioprocess Raman spectra.

Site-based description of R1ρ relaxation in local reference frames.

Nuclear magnetic spin relaxation in the presence of an applied radiofrequency field depends critically on chemical exchange processes that transfer nuclear spins between chemical or conformational environments with distinct resonance frequencies. Characterization of chemical exchange processes in R1ρ relaxation dispersion, CEST, and DEST experiments provides powerful insights into chemical and conformational kinetics of biological macromolecules. The present work reformulates expressions for magnetization evolution and the R1ρ relaxation rate constant by focussing on the orientations of the tilted rotating frames of reference for magnetization components in different sites, by treating the spin-locking field strength as a perturbation to free-precession evolution, and by applying the Homotopy Analysis and Laplace transform methods to approximate solutions to the Bloch-McConnell equations. The results provide an expression for R1ρ that is invariant to the topology of the kinetic scheme, an approximate equation for evolution of spin-locked z-magnetization, and an approach for effective simplification of chemical exchange topologies for 2- and N-site chemical exchange processes. The theoretical approach also provides an accurate approximation for relaxation during a constant-amplitude radiofrequency field in the absence of exchange.

Completely Computational Model Setup for Spectroscopic Techniques: The Ab Initio Molecular Dynamics Indirect Hard Modeling Approach.

The spectroscopic quantification of mixture compositions usually requires pure compounds and mixtures of known compositions for calibration. Since they are not always available, methods to fill such gaps have evolved, which are, however, not generally applicable. Therefore, calibration can be extremely challenging, especially when multiple unstable species, for example, intermediates, exist in a system. This study presents a new calibration approach that uses ab initio molecular dynamics (AIMD)-simulated spectra to set up and calibrate models for the physics-based spectral analysis method indirect hard modeling (IHM). To demonstrate our approach called AIMD-IHM, we analyze Raman spectra of ternary hydrogen-bonding mixtures of acetone, methanol, and ethanol. The derived AIMD-IHM pure-component models and calibration coefficients are in good agreement with conventionally generated experimental results. The method yields compositions with prediction errors of less than 5% without any experimental calibration input. Our approach can be extended, in principle, to infrared and NMR spectroscopy and allows for the analysis of systems that were hitherto inaccessible to quantitative spectroscopic analysis.

Receptor tyrosine kinases regulate signal transduction through a liquid-liquid phase separated state.

The recruitment of signaling proteins into activated receptor tyrosine kinases (RTKs) to produce rapid, high-fidelity downstream response is exposed to the ambiguity of random diffusion to the target site. Liquid-liquid phase separation (LLPS) overcomes this by providing elevated, localized concentrations of the required proteins while impeding competitor ligands. Here, we show a subset of phosphorylation-dependent RTK-mediated LLPS states. We then investigate the formation of phase-separated droplets comprising a ternary complex including the RTK, (FGFR2); the phosphatase, SHP2; and the phospholipase, PLCγ1, which assembles in response to receptor phosphorylation. SHP2 and activated PLCγ1 interact through their tandem SH2 domains via a previously undescribed interface. The complex of FGFR2 and SHP2 combines kinase and phosphatase activities to control the phosphorylation state of the assembly while providing a scaffold for active PLCγ1 to facilitate access to its plasma membrane substrate. Thus, LLPS modulates RTK signaling, with potential consequences for therapeutic intervention.

Dimerization of Cadherin-11 involves multi-site coupled unfolding and strand swapping.

Cadherin extracellular domain 1 (EC1) mediates homophilic dimerization in adherens junctions. Conserved Trp2 and Trp4 residues in type II cadherins anchor the EC1 A strand intermolecularly in strand-swapped dimers. Herein, NMR spectroscopy is used to elucidate the roles of Trp2 and Trp4 in Cadherin-11 dimerization. The monomeric state, with the A strand and Trp side chains packed intramolecularly, is in equilibrium with sparsely populated partially and fully A-strand-exposed states, in which Trp2 (and Trp4, respectively) side-chain packing is disrupted. Exchange kinetics between the major state and the partially (fully) A-strand-exposed state is slow-intermediate (intermediate-fast). A separate very fast process exchanges ordered and random-coil BC-loop conformations with populations dependent on A-strand exposure and dimerization status. In addition, very slow processes connect the folded A-strand-exposed conformation to partially unfolded states, which may represent additional domain-swapping intermediates. The dimerization mechanism of type II cadherins is revealed as coupled folding and strand swapping.

Design, synthesis, and evaluation of peptide-imidazo[1,2-a]pyrazine bioconjugates as potential bivalent inhibitors of the VirB11 ATPase HP0525.

Helicobacter pylori (H. pylori) infections have been implicated in the development of gastric ulcers and various cancers: however, the success of current therapies is compromised by rising antibiotic resistance. The virulence and pathogenicity of H. pylori is mediated by the type IV secretion system (T4SS), a multiprotein macromolecular nanomachine that transfers toxic bacterial factors and plasmid DNA between bacterial cells, thus contributing to the spread of antibiotic resistance. A key component of the T4SS is the VirB11 ATPase HP0525, which is a hexameric protein assembly. We have previously reported the design and synthesis of a series of novel 8-amino imidazo[1,2-a]pyrazine derivatives as inhibitors of HP0525. In order to improve their selectivity, and potentially develop these compounds as tools for probing the assembly of the HP0525 hexamer, we have explored the design and synthesis of potential bivalent inhibitors. We used the structural details of the subunit-subunit interactions within the HP0525 hexamer to design peptide recognition moieties of the subunit interface. Different methods (cross metathesis, click chemistry, and cysteine-malemide) for bioconjugation to selected 8-amino imidazo[1,2-a]pyrazines were explored, as well as peptides spanning larger or smaller regions of the interface. The IC50 values of the resulting linker-8-amino imidazo[1,2-a]pyrazine derivatives, and the bivalent inhibitors, were related to docking studies with the HP0525 crystal structure and to molecular dynamics simulations of the peptide recognition moieties.

Algebraic expressions for Carr-Purcell-Meiboom-Gill relaxation dispersion for N-site chemical exchange.

The Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiment measures the effective relaxation rate constant during a train of spin-echo pulse sequence elements as a function of the echo time. The CPMG experiment is a powerful method for characterizing chemical and conformational dynamic processes, termed chemical and conformational exchange, on µs-ms time scales, comparable to the experimentally accessible echo times. Approximate theoretical expressions for the effective relaxation rate constant for N-site chemical exchange have been reported (H. Koss, M. Rance, and A. G. Palmer, Biochemistry 57, 4753-4763 (2018)). Expressions for the effective relaxation rate constant have been improved by using the Cayley-Hamilton theorem to obtain simple and accurate approximations of the average Liouvillian for the CPMG experiment. The improved accuracy of the results allows efficient analyses of experimental data. In addition, the relationship is clarified between the approach of Koss and coworkers and that of Jen (1978).

Chemical Exchange.

The phenomenon of chemical or conformational exchange in NMR spectroscopy has enabled detailed characterization of time-dependent aspects of biomolecular function, including folding, molecular recognition, allostery, and catalysis, on timescales from microsecond to second. Importantly, NMR methods based on a variety of spin relaxation parameters have been developed that provide quantitative information on interconversion kinetics, thermodynamic properties, and structural features of molecular states populated to a fraction of a percent at equilibrium and otherwise unobservable by other NMR approaches. The ongoing development of more sophisticated experimental techniques and the necessity to apply these methods to larger and more complex molecular systems engenders a corresponding need for theoretical advances describing such techniques and facilitating data analysis in applications. This review surveys current aspects of the theory of chemical exchange, as utilized in ZZ-exchange; Hahn and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo; and R1ρ, chemical exchange saturation transfer (CEST), and dark state saturation transfer (DEST) spin-locking experiments. The review emphasizes theoretical results for kinetic topologies with more than two interconverting states, both to obtain compact analytical forms suitable for data analysis and to establish conditions for distinguishability between alternative kinetic schemes.

General Expressions for Carr-Purcell-Meiboom-Gill Relaxation Dispersion for N-Site Chemical Exchange.

The Carr-Purcell-Meiboom-Gill (CPMG) nuclear magnetic resonance experiment is widely used to characterize chemical exchange phenomena in biological macromolecules. Theoretical expressions for the nuclear spin relaxation rate constant for two-site chemical exchange during CPMG pulse trains valid for all time scales are well-known as are descriptions of N-site exchange in the fast limit. We have obtained theoretical expressions for N-site exchange outside of the fast limit by using approximations to an average Liouvillian describing the decay of magnetization during a CPMG pulse train. We obtain general expressions for CPMG experiments for any N-site scheme and all experimentally accessible time scales. For sufficiently slow chemical exchange, we obtain closed-form expressions for the relaxation rate constant and a general characteristic polynomial for arbitrary kinetic schemes. Furthermore, we highlight features that qualitatively characterize CPMG curves obtained for various N-site kinetic topologies, quantitatively characterize CPMG curves obtained from systems in various N-site exchange situations, and test distinguishability of kinetic models.

Dynamic Allostery in PLCγ1 and Its Modulation by a Cancer Mutation Revealed by MD Simulation and NMR.

Phosphatidylinositol phospholipase Cγ (PLCγ) is an intracellular membrane-associated second-messenger signaling protein activated by tyrosine kinases such as fibroblast growth factor receptor 1. PLCγ contains the regulatory γ-specific array (γSA) comprising a tandem Src homology 2 (SH2) pair, an SH3 domain, and a split pleckstrin homology domain. Binding of an activated growth factor receptor to γSA leads to Tyr783 phosphorylation and consequent PLCγ activation. Several disease-relevant mutations in γSA have been identified; all lead to elevated phospholipase activity. In this work, we describe an allosteric mechanism that connects the Tyr783 phosphorylation site to the nSH2-cSH2 junction and involves dynamic interactions between the cSH2-SH3 linker and cSH2. Molecular dynamics simulations of the tandem SH2 protein suggest that Tyr783 phosphorylation is communicated to the nSH2-cSH2 junction by modulating cSH2 binding to sections of the cSH2-SH3 linker. NMR chemical shift perturbation analyses for designed tandem SH2 constructs reveal combined fast and slow dynamic processes that can be attributed to allosteric communication involving these regions of the protein, establishing an example in which complex N-site exchange can be directly inferred from 1H,15N-HSQC spectra. Furthermore, in tandem SH2 and γSA constructs, molecular dynamics and NMR results show that the Arg687Trp mutant in PLCγ1 (equivalent to the cancer mutation Arg665Trp in PLCγ2) perturbs the dynamic allosteric pathway. This combined experimental and computational study reveals a rare example of multistate kinetics involved in a dynamic allosteric process that is modulated in the context of a disease-relevant mutation. The allosteric influences and the weakened binding of the cSH2-SH3 linker to cSH2 should be taken into account in any more holistic investigation of PLCγ regulation.

NMR backbone assignments of the tyrosine kinase domain of human fibroblast growth factor receptor 3 in apo state and in complex with inhibitor PD173074.

Fibroblast growth factors receptors (FGFR) are transmembrane protein tyrosine kinases involved in many cellular process, including growth, differentiation and angiogenesis. Dysregulation of FGFR enzymatic activity is associated with developmental disorders and cancers; therefore FGFRs have become attractive targets for drug discovery, with a number of agents in late-stage clinical trials. Here, we present the backbone resonance assignments of FGFR3 tyrosine kinase domain in the ligand-free form and in complex with the canonical FGFR kinase inhibitor PD173074. Analysis of chemical shift changes upon inhibitor binding highlights a characteristic pattern of allosteric network perturbations that is of relevance for future drug discovery activities aimed at development of conformationally-selective FGFR inhibitors.

Multicomponent diffusion coefficients from microfluidics using Raman microspectroscopy.

Diffusion is slow. Thus, diffusion experiments are intrinsically time-consuming and laborious. Additionally, the experimental effort is multiplied for multicomponent systems as the determination of multicomponent diffusion coefficients typically requires several experiments. To reduce the experimental effort, we present the first microfluidic diffusion measurement method for multicomponent liquid systems. The measurement setup combines a microfluidic chip with Raman microspectroscopy. Excellent agreement between experimental results and literature data is achieved for the binary system cyclohexane + toluene and the ternary system 1-propanol + 1-chlorobutane + heptane. The Fick diffusion coefficients are obtained from fitting a multicomponent convection-diffusion model to the mole fractions measured in experiments. Ternary diffusion coefficients can be obtained from a single experiment; high accuracy is already obtained from two experiments. Advantages of the presented measurement method are thus short measurement times, reduced sample consumption, and less experiments for the determination of a multicomponent diffusion coefficient.

General expressions for R1ρ relaxation for N-site chemical exchange and the special case of linear chains.

Exploration of dynamic processes in proteins and nucleic acids by spin-locking NMR experiments has been facilitated by the development of theoretical expressions for the R1ρ relaxation rate constant covering a variety of kinetic situations. Herein, we present a generalized approximation to the chemical exchange, Rex, component of R1ρ for arbitrary kinetic schemes, assuming the presence of a dominant major site population, derived from the negative reciprocal trace of the inverse Bloch-McConnell evolution matrix. This approximation is equivalent to first-order truncation of the characteristic polynomial derived from the Bloch-McConnell evolution matrix. For three- and four-site chemical exchange, the first-order approximations are sufficient to distinguish different kinetic schemes. We also introduce an approach to calculate R1ρ for linear N-site schemes, using the matrix determinant lemma to reduce the corresponding 3N×3N Bloch-McConnell evolution matrix to a 3×3 matrix. The first- and second order-expansions of the determinant of this 3×3 matrix are closely related to previously derived equations for two-site exchange. The second-order approximations for linear N-site schemes can be used to obtain more accurate approximations for non-linear N-site schemes, such as triangular three-site or star four-site topologies. The expressions presented herein provide powerful means for the estimation of Rex contributions for both low (CEST-limit) and high (R1ρ-limit) radiofrequency field strengths, provided that the population of one state is dominant. The general nature of the new expressions allows for consideration of complex kinetic situations in the analysis of NMR spin relaxation data.

The Effect of Mutations on Drug Sensitivity and Kinase Activity of Fibroblast Growth Factor Receptors: A Combined Experimental and Theoretical Study.

Fibroblast growth factor receptors (FGFRs) are recognized therapeutic targets in cancer. We here describe insights underpinning the impact of mutations on FGFR1 and FGFR3 kinase activity and drug efficacy, using a combination of computational calculations and experimental approaches including cellular studies, X-ray crystallography and biophysical and biochemical measurements. Our findings reveal that some of the tested compounds, in particular TKI258, could provide therapeutic opportunity not only for patients with primary alterations in FGFR but also for acquired resistance due to the gatekeeper mutation. The accuracy of the computational methodologies applied here shows a potential for their wider application in studies of drug binding and in assessments of functional and mechanistic impacts of mutations, thus assisting efforts in precision medicine.

2- and 3-substituted imidazo[1,2-a]pyrazines as inhibitors of bacterial type IV secretion.

A novel series of 8-amino imidazo[1,2-a]pyrazine derivatives has been developed as inhibitors of the VirB11 ATPase HP0525, a key component of the bacterial type IV secretion system. A flexible synthetic route to both 2- and 3-aryl substituted regioisomers has been developed. The resulting series of imidazo[1,2-a]pyrazines has been used to probe the structure-activity relationships of these inhibitors, which show potential as antibacterial agents.

Dysfunction of phospholipase Cγ in immune disorders and cancer.

The surge in genetic and genomic investigations over the past 5 years has resulted in many discoveries of causative variants relevant to disease pathophysiology. Although phospholipase C (PLC) enzymes have long been recognized as important components in intracellular signal transmission, it is only recently that this approach highlighted their role in disease development through gain-of-function mutations. In this review we describe the new findings that link the PLCγ family to immune disorders and cancer, and illustrate further efforts to elucidate the molecular mechanisms that underpin their dysfunction.