Deconvolutions based on singular value decomposition and the pseudoinverse: a guide for beginners

Drugs Form Ternary Complexes with Human Liver Fatty Acid Binding Protein 1 (FABP1) and FABP1 Binding Alters Drug Metabolism

Liver fatty acid binding protein 1 (FABP1) binds diverse endogenous lipids and is highly expressed in the human liver. Binding to FABP1 alters the metabolism and homeostasis of endogenous lipids in the liver. Drugs have also been shown to bind to rat FABP1, but limited data are available for human FABP1 (hFABP1). FABP1 has a large binding pocket, and up to two fatty acids can bind to FABP1 simultaneously. We hypothesized that drug binding to hFABP1 results in formation of ternary complexes and that FABP1 binding alters drug metabolism. To test these hypotheses, native protein mass spectrometry (MS) and fluorescent 11-(dansylamino)undecanoic acid (DAUDA) displacement assays were used to characterize drug binding to hFABP1, and diclofenac oxidation by cytochrome P450 2C9 (CYP2C9) was studied in the presence and absence of hFABP1. DAUDA binding to hFABP1 involved high (Kd,1 = 0.2 μM) and low (Kd,2 > 10 μM) affinity binding sites. Nine drugs bound to hFABP1 with equilibrium dissociation constant (Kd) values ranging from 1 to 20 μM. None of the tested drugs completely displaced DAUDA from hFABP1, and fluorescence spectra showed evidence of ternary complex formation. Formation of DAUDA-hFABP1-diclofenac ternary complex was verified with native MS. Docking predicted diclofenac binding in the portal region of FABP1 with DAUDA in the binding cavity. The catalytic rate constant of diclofenac hydroxylation by CYP2C9 was decreased by ∼50% (P < 0.01) in the presence of FABP1. Together, these results suggest that drugs form ternary complexes with hFABP1 and that hFABP1 binding in the liver will alter drug metabolism and clearance. SIGNIFICANCE STATEMENT: Many commonly prescribed drugs bind fatty acid binding protein 1 (FABP1), forming ternary complexes with FABP1 and the fluorescent fatty acid 11-(dansylamino)undecanoic acid. These findings suggest that drugs will bind to apo-FABP1 and fatty acid-bound FABP1 in the human liver. The high expression of FABP1 in the liver, together with drug binding to FABP1, may alter drug disposition processes in vivo.

Drugs Form Ternary Complexes with Human Liver Fatty Acid Binding Protein (FABP1) and FABP1 Binding Alters Drug Metabolism

Liver fatty acid binding protein (FABP1) binds diverse endogenous lipids and is highly expressed in the human liver. Binding to FABP1 alters the metabolism and homeostasis of endogenous lipids in the liver. Drugs have also been shown to bind to rat FABP1, but limited data is available for human FABP1 (hFABP1). FABP1 has a large binding pocket and multiple fatty acids can bind to FABP1 simultaneously. We hypothesized that drug binding to hFABP1 results in formation of ternary complexes and that FABP1 binding alters drug metabolism. To test these hypotheses native protein mass spectrometry (MS) and fluorescent 11-(dansylamino)undecanoic acid (DAUDA) displacement assays were used to characterize drug binding to hFABP1 and diclofenac oxidation by cytochrome P450 2C9 (CYP2C9) was studied in the presence and absence of hFABP1. DAUDA binding to hFABP1 involved high (Kd,1=0.2 µM) and low affinity (Kd,2 >10 µM) binding sites. Nine drugs bound to hFABP1 with Kd values ranging from 1 to 20 µM. None of the tested drugs completely displaced DAUDA from hFABP1 and fluorescence spectra showed evidence of ternary complex formation. Formation of DAUDA-diclofenac-hFABP1 ternary complex was verified with native MS. Docking placed diclofenac in the portal region of FABP1 with DAUDA in the binding cavity. Presence of hFABP1 decreased the kcat and Km,u of diclofenac with CYP2C9 by ~50% suggesting that hFABP1 binding in the liver will alter drug metabolism and clearance. Together, these results suggest that drugs form ternary complexes with hFABP1 and that hFABP1 interacts with CYP2C9.

Nanoelectromechanical Infrared Spectroscopy with In Situ Separation by Thermal Desorption: NEMS-IR-TD

We present a novel method for the quantitative analysis of mixtures of semivolatile chemical compounds. For the first time, thermal desorption is integrated directly with nanoelectromechanical infrared spectroscopy (NEMS-IR-TD). In this new technique, an analyte mixture is deposited via nebulization on the surface of a NEMS sensor and subsequently desorbed using heating under vacuum. The desorption process is monitored in situ via infrared spectroscopy and thermogravimetric analysis. The resulting spectro-temporal maps allow for selective identification and analysis of the mixture. In addition, the corresponding thermogravimetric data allow for analysis of the desorption dynamics of the mixture components. As a demonstration, caffeine and theobromine were selectively identified and quantified from a mixture with a detection limit of less than 6 pg (about 30 fmol). With its exceptional sensitivity, NEMS-IR-TD allows for the analysis of low abundance and complex analytes with potential applications ranging from environmental sensing to life sciences.

MitoNEET's Reactivity of Lys55 toward Pyridoxal Phosphate Demonstrates its Activity as a Transaminase Enzyme

MitoNEET is a [2Fe-2S] redox active mitochondrial protein belonging to the CDGSH iron-sulfur domain (CISD) family of proteins. MitoNEET has been implicated as a potential target for drug development to treat various disorders, including type-2 diabetes, cancer, and Parkinson's disease. However, the specific cellular function(s) for mitoNEET still remains to be fully elucidated, and this presents a significant roadblock in rational drug development. Here, we show that mitoNEET binds the enzymatic cofactor pyridoxal phosphate (PLP) specifically at only one of its 11 lysine residues, Lys55. Lys55 is part of the soluble portion of the protein and is in a hydrogen-bonding network with the histidine residue that ligates the [2Fe-2S] cluster. In the presence of mitoNEET, PLP catalyzes the transamination reaction of the amino acid cysteine and the alpha-keto acid 2-oxoglutarate to form 3-mercaptopyruvate and glutamate. This work identifies, for the first time, mitoNEET as an enzyme with cysteine transaminase activity.

KiMoPack: A python Package for Kinetic Modeling of the Chemical Mechanism

Herein, we present KiMoPack, an analysis tool for the kinetic modeling of transient spectroscopic data. KiMoPack enables a state-of-the-art analysis routine including data preprocessing and standard fitting (global analysis), as well as fitting of complex (target) kinetic models, interactive viewing of (fit) results, and multiexperiment analysis via user accessible functions and a graphical user interface (GUI) enhanced interface. To facilitate its use, this paper guides the user through typical operations covering a wide range of analysis tasks, establishes a typical workflow and is bridging the gap between ease of use for less experienced users and introducing the advanced interfaces for experienced users. KiMoPack is open source and provides a comprehensive front-end for preprocessing, fitting and plotting of 2-dimensional data that simplifies the access to a powerful python-based data-processing system and forms the foundation for a well documented, reliable, and reproducible data analysis.

Reversible, β-sheet-dependent self-assembly of the phosphoprotein phosvitin is controlled by the concentration and valency of cations

The hyperphosphorylated protein phosvitin (PV) undergoes a pH-dependent transition between P_II^- and β-sheet secondary structures, a process deemed crucial for its role in the promotion of biogenic apatite formation. The transition occurs surprisingly slowly (minutes to hours). This is consistent with a slow aggregation process involving ionic interactions of charged groups on the protein surface. Herein, we determined the associated transition pK values and time constants through matrix least-squares (MLS) global fitting of a series of pH- and time-dependent circular dichroism (CD) spectra recorded in the presence of different mono-, bi- and trivalent cations. Supporting our results with dynamic light scattering data, we clearly identified a close correlation of β-sheet transition and the formation of small aggregates at low pH. This process is inhibited in the presence of all tested cations with the strongest effects for trivalent cations (Fe³⁺ and Al³⁺). In the presence of Ca²⁺ and Mg²⁺, larger higher-order particles are formed from PV in the β-sheet conformation, as identified from the interpretation of differential scattering observed in the CD spectra. Our observations are consistent with the existence of a multi-step equilibrium between aggregated and non-aggregated species of PV. The equilibrium is highly sensitive to the environment pH and salt concentration with exceptional behavior in the presence of divalent cations such as Ca²⁺ and Mg²⁺.

Controlling amphipathic peptide adsorption by smart switchable germanium interfaces

The in situ control of reversible protein adsorption to a surface is a critical step towards biofouling prevention and finds utilisation in bioanalytical applications. In this work, adsorption of peptides is controlled by employing the electrode potential induced, reversible change of germanium (100) surface termination between a hydrophobic, hydrogen terminated and a hydrophilic, hydroxyl terminated surface. This simple but effective 'smart' interface is used to direct adsorption of two peptides models, representing the naturally highly abundant structural motifs of amphipathic helices and coiled-coils. Their structural similarity coincides with their opposite overall charge and hence allows the examination of the influence of charge and hydrophobicity on adsorption. Polarized attenuated total reflection infrared (ATR-IR) spectroscopy at controlled electrode potential has been used to follow the adsorption process at physiological pH in deuterated buffer. The delicate balance of hydrophobic and electrostatic peptide/surface interactions leads to two different processes upon switching that are both observed in situ: reversible adsorption and reversible reorientation. Negatively charged peptide adsorption can be fully controlled by switching to the hydrophobic interface, while the same switch causes the positively charged, helical peptide to tilt down. This principle can be used for 'smart' adsorption control of a wider variety of proteins and peptides and hence find application, as e.g. a bioanalytical tool or functional biosensor.

Dynamic Structural Changes and Thermodynamics in Phase Separation Processes of an Intrinsically Disordered-Ordered Protein Model

Elastin-like proteins (ELPs) are biologically important proteins and models for intrinsically disordered proteins (IDPs) and dynamic structural transitions associated with coacervates and liquid-liquid phase transitions. However, the conformational status below and above coacervation temperature and its role in the phase separation process is still elusive. Employing matrix least-squares global Boltzmann fitting of the circular dichroism spectra of the ELPs (VPGVG)20 , (VPGVG)40 , and (VPGVG)60 , we found that coacervation occurs sharply when a certain number of repeat units has acquired β-turn conformation (in our sequence setting a threshold of approx. 20 repeat units). The character of the differential scattering of the coacervate suspensions indicated that this fraction of β-turn structure is still retained after polypeptide assembly. Such conformational thresholds may also have a role in other protein assembly processes with implications for the design of protein-based smart materials.

Control of distal lysine coordination in a monomeric hemoglobin: A role for heme peripheral interactions

THB1 is a monomeric truncated hemoglobin (TrHb) found in the cytoplasm of the green alga Chlamydomonas reinhardtii. The canonical heme coordination scheme in hemoglobins is a proximal histidine ligand and an open distal site. In THB1, the latter site is occupied by Lys53, which is likely to facilitate Fe(II)/Fe(III) redox cycling but hinders dioxygen binding, two features inherent to the NO dioxygenase activity of the protein. TrHb surveys show that a lysine at a position aligning with Lys53 is an insufficient determinant of coordination, and in this study, we sought to identify factors controlling lysine affinity for the heme iron. We solved the "Lys-off" X-ray structure of THB1, represented by the cyanide adduct of the Fe(III) protein, and hypothesized that interactions that differ between the known "Lys-on" structure and the Lys-off structure participate in the control of Lys53 affinity for the heme iron. We applied an experimental approach (site-directed mutagenesis, heme modification, pH titrations in the Fe(III) and Fe(II) states) and a computational approach (MD simulations in the Fe(II) state) to assess the role of heme propionate-protein interactions, distal helix capping, and the composition of the distal pocket. All THB1 modifications resulted in a weakening of lysine affinity and affected the coupling between Lys53 proton binding and heme redox potential. The results supported the importance of specific heme peripheral interactions for the pH stability of iron coordination and the ability of the protein to undergo redox reactions.

REGALS: a general method to deconvolve X-ray scattering data from evolving mixtures

Mixtures of biological macromolecules are inherently difficult to study using structural methods, as increasing complexity presents new challenges for data analysis. Recently, there has been growing interest in studying evolving mixtures using small-angle X-ray scattering (SAXS) in conjunction with time-resolved, high-throughput or chromatography-coupled setups. Deconvolution and interpretation of the resulting datasets, however, are nontrivial when neither the scattering components nor the way in which they evolve are known a priori. To address this issue, the REGALS method (regularized alternating least squares) is introduced, which incorporates simple expectations about the data as prior knowledge, and utilizes parameterization and regularization to provide robust deconvolution solutions. The restraints used by REGALS are general properties such as smoothness of profiles and maximum dimensions of species, making it well suited for exploring datasets with unknown species. Here, REGALS is applied to the analysis of experimental data from four types of SAXS experiment: anion-exchange (AEX) coupled SAXS, ligand titration, time-resolved mixing and time-resolved temperature jump. Based on its performance with these challenging datasets, it is anticipated that REGALS will be a valuable addition to the SAXS analysis toolkit and enable new experiments. The software is implemented in both MATLAB and Python and is available freely as an open-source software package.

Polymorphism of human telomeric quadruplexes with drugs: a multi-technique biophysical study

The human telomeric G-quadruplex structural motif of DNA has come to be known as a new and stimulating target for anticancer drug discovery. Small molecules that interact with G-quadruplex structures in a selective way have gained impressive interest in recent years as they may serve as potential therapeutic agents. Here, we show how circular dichroism, UV resonance Raman and small angle X-ray scattering spectroscopies can be effectively combined to provide insights into structural and molecular aspects of the interaction between human telomeric quadruplexes and ligands. This study focuses on the ability of berberine and palmatine to bind with human telomeric quadruplexes and provides analysis of the conformational landscape visited by the relevant complexes upon thermal unfolding. With increasing temperature, both free and bound G-quadruplexes undergo melting through a multi-state process, populating different intermediate states. Despite the structural similarity of the two ligands, valuable distinctive features characterising their interaction with the G-quadruplex emerged from our multi-technique approach.

Spatial confinement alters the ultrafast photoisomerization dynamics of azobenzenes

Ultrafast transient absorption spectroscopy reveals new excited-state dynamics following excitation of trans-azobenzene (t-Az) and several alkyl-substituted t-Az derivatives encapsulated in a water-soluble supramolecular host-guest complex. Encapsulation increases the excited-state lifetimes and alters the yields of the trans → cis photoisomerization reaction compared with solution. Kinetic modeling of the transient spectra for unsubstituted t-Az following nπ* and ππ* excitation reveals steric trapping of excited-state species, as well as an adiabatic excited-state trans → cis isomerization pathway for confined molecules that is not observed in solution. Analysis of the transient spectra following ππ* excitation for a series of 4-alkyl and 4,4'-dialkyl substituted t-Az molecules suggests that additional crowding due to lengthening of the alkyl tails results in deeper trapping of the excited-state species, including distorted trans and cis structures. The variation of the dynamics due to crowding in the confined environment provides new evidence to explain the violation of Kasha's rule for nπ* and ππ* excitation of azobenzenes based on competition between in-plane inversion and out-of-plane rotation channels.

Recent applications of chemometrics in one- and two-dimensional chromatography

The proliferation of increasingly more sophisticated analytical separation systems, often incorporating increasingly more powerful detection techniques, such as high-resolution mass spectrometry, causes an urgent need for highly efficient data-analysis and optimization strategies. This is especially true for comprehensive two-dimensional chromatography applied to the separation of very complex samples. In this contribution, the requirement for chemometric tools is explained and the latest developments in approaches for (pre-)processing and analyzing data arising from one- and two-dimensional chromatography systems are reviewed. The final part of this review focuses on the application of chemometrics for method development and optimization.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

Two-quartet kit* G-quadruplex is formed via double-stranded pre-folded structure

In the promoter of c-KIT proto-oncogene, whose deregulation has been implicated in many cancers, three G-rich regions (kit1, kit* and kit2) are able to fold into G-quadruplexes. While kit1 and kit2 have been studied in depth, little information is available on kit* folding behavior despite its key role in regulation of c-KIT transcription. Notably, kit* contains consensus sites for SP1 and AP2 transcription factors. Herein, a set of complementary spectroscopic and biophysical methods reveals that kit*, d[GGCGAGGAGGGGCGTGGCCGGC], adopts a chair type antiparallel G-quadruplex with two G-quartets at physiological relevant concentrations of KCl. Heterogeneous ensemble of structures is observed in the presence of Na⁺ and NH₄⁺ ions, which however stabilize pre-folded structure. In the presence of K⁺ ions stacking interactions of adenine and thymine residues on the top G-quartet contribute to structural stability together with a G10•C18 base pair and a fold-back motif of the five residues at the 3′-terminal under the bottom G-quartet. The 3′-tail enables formation of a bimolecular pre-folded structure that drives folding of kit* into a single G-quadruplex. Intriguingly, kinetics of kit* G-quadruplex formation matches timescale of transcriptional processes and might demonstrate interplay of kinetic and thermodynamic factors for understanding regulation of c-KIT proto-oncogene expression.

A Heme Propionate Staples the Structure of Cytochrome c for Methionine Ligation to the Heme Iron

Ligand-switch reactions at the heme iron are common in biological systems, but their mechanisms and the features of the polypeptide fold that support dual ligation are not well understood. In cytochrome c (cyt c), two low-stability loops (Ω-loop C and Ω-loop D) are connected by the heme propionate HP6. At alkaline pH, the native Met80 ligand from Ω-loop D switches to a Lys residue from the same loop. Deprotonation of an as yet unknown group triggers the alkaline transition. We have created the two cyt c variants T49V/K79G and T78V/K79G with altered connections of these two loops to HP6. Electronic absorption, NMR, and EPR studies demonstrate that at pH 7.4 ferric forms of these variants are Lys-ligated, whereas ferrous forms maintain the native Met80 ligation. Measurements of protein stability, cyclic voltammetry, pH-jump and gated electron-transfer kinetics have revealed that these Thr to Val substitutions greatly affect the alkaline transition in both ferric and ferrous proteins. The substitutions modify the stability of the Met-ligated species and reduction potentials of the heme iron. The kinetics of ligand-switch processes are also altered, and analyses of these effects implicate redox-dependent differences in metal-ligand interactions and the role of the protein dynamics, including cross-talk between the two Ω-loops. With the two destabilized variants, it is possible to map energy levels for the Met- and Lys-ligated species in both ferric and ferrous proteins and assess the role of the protein scaffold in redox-dependent preferences for these two ligands. The estimated shift in the heme iron reduction potential upon deprotonation of the "trigger" group is consistent with those associated with deprotonation of an HP, suggesting that HP6, on its own or as a part of a hydrogen-bonded cluster, is a likely "trigger" for the Met to Lys ligand switch.

Time-resolved IR spectroscopy reveals mechanistic details of ion transport in the sodium pump Krokinobacter eikastus rhodopsin 2

We report a comparative study on the structural dynamics of the light-driven sodium pump Krokinobacter eikastus rhodopsin 2 wild type under sodium and proton pumping conditions by means of time-resolved IR spectroscopy. The kinetics of KR2 under sodium pumping conditions exhibits a sequential character, whereas the kinetics of KR2 under proton pumping conditions involves several equilibrium states. The sodium translocation itself is characterized by major conformational changes of the protein backbone, such as distortions of the α-helices and probably of the ECL1 domain, indicated by distinct marker bands in the amide I region. Carbonyl stretch modes of specific amino acid residues helped to elucidate structural changes in the retinal Schiff base moiety, including the protonation and deprotonation of D116, which is crucial for a deeper understanding of the mechanistic features in the photocycle of KR2.

Studying NAD(P)H cofactor-binding to alcohol dehydrogenases through global analysis of circular dichroism spectra

The initial step in reactions catalyzed by NAD(P)H-dependent alcohol dehydrogenases (ADHs) is the binding of the cofactor to the active site. To study this process, we measured NAD(P)H concentration-dependent circular dichroism (CD) in the presence of purified enzymes (ADH from horse liver, HLADH; ADH-A from Rhodococcus ruber; YGL157w from Saccharomyces cerevisiae) or enzyme-containing whole cell extract (ADH from Lactobacillus brevis, LbADH). We determined the proportions of binding and non-binding NAD(P)H and the associated dissociation constants (Kd) from matrix least-squares global fitting of law of mass action-derived model. Furthermore, the fitting allowed the back calculation of CD spectra corresponding to the cofactor in its bound conformation. With increasing pH and/or increasing ionic strength, we detected an increase in Kd for the NADH·HLADH complex with the shape of the bound cofactor conformation spectrum remaining unaffected. While the bound cofactor spectrum for the ADH-A·NADH complex was similar to that for HLADH, the corresponding spectra obtained for the NADPH-dependent enzymes YGL157w and LbADH exhibited opposite signs of the most prominent band. In comparison to CD spectra calculated on cofactor geometries from the crystal structures at the sTD-DFT level, we found that the sign of the bound cofactor spectrum correlates with the orientation of the nicotinamide ring of the cofactor in the active site. These results demonstrate the usefulness of the global analysis of cofactor titration CD spectra to study the role of cofactor binding and its geometry in ADH catalysis.

Replacement of the heme axial lysine as a test of conformational adaptability in the truncated hemoglobin THB1

Amino acid replacement is a useful strategy to assess the roles of axial heme ligands in the function of native heme proteins. THB1, the protein product of the Chlamydomonas reinhardtii THB1 gene, is a group 1 truncated hemoglobin that uses a lysine residue in the E helix (Lys53, at position E10 by reference to myoglobin) as an iron ligand at neutral pH. Phylogenetic evidence shows that many homologous proteins have a histidine, methionine or arginine at the same position. In THB1, these amino acids would each be expected to convey distinct reactive properties if replacing the native lysine as an axial ligand. To explore the ability of the group 1 truncated Hb fold to support alternative ligation schemes and distal pocket conformations, the properties of the THB1 variants K53A as a control, K53H, K53M, and K53R were investigated by electronic absorption, EPR, and NMR spectroscopies. We found that His53 is capable of heme ligation in both the Fe(III) and Fe(II) states, that Met53 can coordinate only in the Fe(II) state, and that Arg53 stabilizes a hydroxide ligand in the Fe(III) state. The data illustrate that the group 1 truncated Hb fold can tolerate diverse rearrangement of the heme environment and has a strong tendency to use two protein side chains as iron ligands despite accompanying structural perturbations. Access to various redox pairs and different responses to pH make this protein an excellent test case for energetic and dynamic studies of heme ligation.

Mechanistic studies of a "Declick" reaction

A kinetic analysis of a "declick" reaction is described. Compound 1, previously reported to couple an amine and a thiol (i.e. "click") under mild aqueous conditions to create 2, undergoes release of the unaltered coupling partners upon triggering with dithiothreitol (DTT). In the study reported herein various aniline derivatives possessing para-electron donating and withdrawing groups were used as the amines. UV/vis spectroscopy of the declick reaction shows time-dependent spectra lacking isosbestic points, implying a multi-step mechanism. Global data fitting using numerical integration of rate equations and singular value decomposition afforded the spectra and time-dependence of each species, as well as rate constants for each step. The kinetic analysis reveals a multi-step process with an intermediate where both thiols of DTT have added prior to expulsion of the aniline leaving group, followed by rearrangement to the final product. Hammett plots show a negative rho value on two of the steps, indicating positive charge building (i.e. reduction of a negative charge) in the step leading to the intermediate and its rate-determining breakdown. Overall, the kinetic study reported herein gives a complete mechanistic picture of the declick reaction.

Arguments for an additional long-lived intermediate in the photocycle of the full-length aureochrome 1c receptor: A time-resolved small-angle X-ray scattering study

Aureochromes (AUREO) act as blue-light photoreceptors in algae. They consist of a light-, oxygen-, voltage-sensitive (LOV) domain and a DNA-binding basic region/leucine zipper. Illumination of the flavin cofactor in LOV leads to the formation of an adduct, followed by global structural changes. Here, we first applied UV/vis spectroscopy to characterize the photocycle of full-length aureochrome 1c (PtAUREO1c) from the diatom Phaeodactylum tricornutum. With a time constant of 850 s and a quantum yield of 23%, PtAUREO1c reveals a faster recovery time and a much lower sensitivity toward light than PtAUREO1a, pointing to its role as a high light sensor in vivo. UV/vis spectroscopy offers details on the local recovery of the flavin chromophore. However, kinetic information on the global structural recovery of full-length AUREO or any other multidomain LOV protein is missing. This information is essential not least for the photoreceptors' applications as optogenetic devices. Therefore, we established a procedure to apply small-angle X-ray scattering on PtAUREO1c in a time-resolved manner employing an in-house setup. In combination with UV/vis spectroscopy under similar conditions, we revealed a discrepancy between the recovery of the global protein structure and the adduct lifetime. Accordingly, we propose to supplement the photocycle by an intermediate state (I447), which decays with a time constant of about 800 s and prolongs the lifetime of the signaling state.

Coherent and dissipative quantum process tensor reconstructions in two-dimensional electronic spectroscopy

A major goal of time-resolved spectroscopy is to resolve the dynamical processes that follow photoexcitation. This amounts to identifying all the quantum states involved and the rates of population transfer between them. Unfortunately, such quantum state and kinetic reconstructions are ambiguous using one-dimensional methods such as transient absorption even when all the states of the system are fully resolved. Higher-dimensionality methods like two-dimensional spectroscopy lift some of the ambiguity, but unless the spectral features are well-separated, current inversion methods generally fail. Here, we show that, using both coherence and population signals of the nonlinear response, it is indeed possible to accurately extract both static and dynamic information from the 2D spectrum even when features are highly congested. Coherences report on the positions of the vibronic states of the system, providing a useful constraint for extracting the full kinetic scheme. We model time-resolved 2D photon echo spectra using a sum-over-states approach and show in which regimes the Hamiltonian and kinetic schemes may be recovered. Furthermore, we discuss how such algorithms may be applied to experimental data and where some of the underlying assumptions may fail. The ability to systematically extract the maximal information content of multidimensional spectroscopic data is an important step toward utilizing the full power of these techniques and elucidating the structure and dynamics of increasingly complex molecular systems.

CH2 domain orientation of human immunoglobulin G in solution: Structural comparison of glycosylated and aglycosylated Fc regions using small-angle X-ray scattering

The N-linked glycan in immunoglobulin G is critical for the stability and function of the crystallizable fragment (Fc) region. Alteration of these protein properties upon the removal of the N-linked glycan has often been explained by the alteration of the CH2 domain orientation in the Fc region. To confirm this hypothesis, we examined the small-angle X-ray scattering (SAXS) profile of the glycosylated Fc region (gFc) and aglycosylated Fc region (aFc) in solution. Conformational characteristics of the CH2 domain orientation were validated by comparison with SAXS profiles theoretically calculated from multiple crystal structures of the Fc region with different CH2 domain orientations. The reduced chi-square values from the fitting analyses of gFc and aFc associated with the degree of openness or closure of each crystal structure, as determined from the first principal component that partially governed the variation of the CH2 domain orientation extracted by a singular value decomposition analysis. For both gFc and aFc, the best-fitted SAXS profiles corresponded to ones calculated based on the crystal structure of gFc that formed a "semi-closed" CH2 domain orientation. Collectively, the data indicated that the removal of the N-linked glycan only negligibly affected the CH2 domain orientation in solution. These findings will guide the development of methodology for the production of highly refined functional Fc variants.

Length-Dependent Melting Kinetics of Short DNA Oligonucleotides Using Temperature-Jump IR Spectroscopy

In this work, we utilize Fourier transform infrared and temperature-jump (T-jump) infrared (IR) spectroscopy to investigate the melting thermodynamics and kinetics of a series of five DNA sequences ranging from 6 to 14 base pairs long. IR spectroscopy is well suited for the study of DNA because of its ability to distinguish base-specific information, and the nanosecond time resolution of the T-jump apparatus can access the relevant range of kinetics. Eyring analysis of a two-state model examines both the activation enthalpy and entropy, providing new insights into the energetic driving forces and physical processes behind the association and dissociation while also helping to clarify the commonly observed negative activation energy. Global analysis of the thermodynamic and kinetic data applying a linear dependence of activation barriers on oligo length provides a holistic result by producing reasonable agreement between our data and existing nearest-neighbor (NN) thermodynamic parameters blending the experimental results with established predictive models. By studying the trends in the thermodynamics and kinetics as a function of length, this work demonstrates a direct correlation between the effects additional dinucleotides have on the kinetics and the NN parameters for those dinucleotides. This result further supports the development of a kinetic analogue to the thermodynamic NN parameters.

Vibrational spectroscopic study of pH dependent solvation at a Ge(100)-water interface during an electrode potential triggered surface termination transition

The charge-dependent structure of interfacial water at the n-Ge(100)-aqueous perchlorate interface was studied by controlling the electrode potential. Specifically, a joint attenuated total reflection infrared spectroscopy and electrochemical experiment was used in 0.1M NaClO₄ at pH ≈ 1-10. The germanium surface transformation to an H-terminated surface followed the thermodynamic Nernstian pH dependence and was observed throughout the entire pH range. A singular value decomposition-based spectra deconvolution technique coupled to a sigmoidal transition model for the potential dependence of the main components in the spectra shows the surface transformation to be a two-stage process. The first stage was observed together with the first appearance of Ge-H stretching modes in the spectra and is attributed to the formation of a mixed surface termination. This transition was reversible. The second stage occurs at potentials ≈0.1-0.3 V negative of the first one, shows a hysteresis in potential, and is attributed to the formation of a surface with maximum Ge-H coverage. During the surface transformation, the surface becomes hydrophobic, and an effective desolvation layer, a "hydrophobic gap," developed with a thickness ≈1-3 Å. The largest thickness was observed near neutral pH. Interfacial water IR spectra show a loss of strongly hydrogen-bound water molecules compared to bulk water after the surface transformation, and the appearance of "free," non-hydrogen bound OH groups, throughout the entire pH range. Near neutral pH at negative electrode potentials, large changes at wavenumbers below 1000 cm^-1 were observed. Librational modes of water contribute to the observed changes, indicating large changes in the water structure.

Quenching of chlorophyll triplet states by carotenoids in algal light-harvesting complexes related to fucoxanthin-chlorophyll protein

We have used time-resolved absorption and fluorescence spectroscopy with nanosecond resolution to study triplet energy transfer from chlorophylls to carotenoids in a protective process that prevents the formation of reactive singlet oxygen. The light-harvesting complexes studied were isolated from Chromera velia, belonging to a group Alveolata, and Xanthonema debile and Nannochloropsis oceanica, both from Stramenopiles. All three light-harvesting complexes are related to fucoxanthin-chlorophyll protein, but contain only chlorophyll a and no chlorophyll c. In addition, they differ in the carotenoid content. This composition of the complexes allowed us to study the quenching of chlorophyll a triplet states by different carotenoids in a comparable environment. The triplet states of chlorophylls bound to the light-harvesting complexes were quenched by carotenoids with an efficiency close to 100%. Carotenoid triplet states were observed to rise with a ~5 ns lifetime and were spectrally and kinetically homogeneous. The triplet states were formed predominantly on the red-most chlorophylls and were quenched by carotenoids which were further identified or at least spectrally characterized.

Conformational profiling of a G-rich sequence within the c-KIT promoter

G-quadruplexes (G4) within oncogene promoters are considered to be promising anticancer targets. However, often they undergo complex structural rearrangements that preclude a precise description of the optimal target. Moreover, even when solved structures are available, they refer to the thermodynamically stable forms but little or no information is supplied about their complex multistep folding pathway. To shed light on this issue, we systematically followed the kinetic behavior of a G-rich sequence located within the c-KIT proximal promoter (kit2) in the presence of monovalent cations K⁺ and Na⁺. A very short-lived intermediate was observed to start the G4 folding process in both salt conditions. Subsequently, the two pathways diverge to produce distinct thermodynamically stable species (parallel and antiparallel G-quadruplex in K⁺ and Na⁺, respectively). Remarkably, in K⁺-containing solution a branched pathway is required to drive the wild type sequence to distribute between a monomeric and dimeric G-quadruplex. Our approach has allowed us to identify transient forms whose relative abundance is regulated by the environment; some of them were characterized by a half-life within the timescale of physiological DNA processing events and thus may represent possible unexpected targets for ligands recognition.

Histidine-Lysine Axial Ligand Switching in a Hemoglobin: A Role for Heme Propionates

The hemoglobin of Synechococcus sp. PCC 7002, GlbN, is a monomeric group I truncated protein (TrHb1) that coordinates the heme iron with two histidine ligands at neutral pH. One of these is the distal histidine (His46), a residue that can be displaced by dioxygen and other small molecules. Here, we show with mutagenesis, electronic absorption spectroscopy, and nuclear magnetic resonance (NMR) spectroscopy that at high pH and exclusively in the ferrous state, Lys42 competes with His46 for the iron coordination site. When b heme is originally present, the population of the lysine-bound species remains too small for detailed characterization; however, the population can be increased significantly by using dimethyl-esterified heme. Electronic absorption and NMR spectroscopies showed that the reversible ligand switching process occurs with an apparent pK_a of 9.3 and a Lys-ligated population of ∼60% at the basic pH limit in the modified holoprotein. The switching rate, which is slow on the chemical shift time scale, was estimated to be 20-30 s^-1 by NMR exchange spectroscopy. Lys42-His46 competition and attendant conformational rearrangement appeared to be related to weakened bis-histidine ligation and enhanced backbone dynamics in the ferrous protein. The pH- and redox-dependent ligand exchange process observed in GlbN illustrates the structural plasticity allowed by the TrHb1 fold and demonstrates the importance of electrostatic interactions at the heme periphery for achieving axial ligand selection. An analogy is drawn to the alkaline transition of cytochrome c, in which Lys-Met competition is detected at alkaline pH, but, in contrast to GlbN, in the ferric state only.

Protein dynamics observed by tunable mid-IR quantum cascade lasers across the time range from 10ns to 1s

We have developed a spectrometer based on tunable quantum cascade lasers (QCLs) for recording time-resolved absorption spectra of proteins in the mid-infrared range. We illustrate its performance by recording time-resolved difference spectra of bacteriorhodopsin in the carboxylic range (1800-1700cm^-1) and on the CO rebinding reaction of myoglobin (1960-1840cm^-1), at a spectral resolution of 1cm^-1. The spectrometric setup covers the time range from 4ns to nearly a second with a response time of 10-15ns. Absorption changes as low as 1×10^-4 are detected in single-shot experiments at t>1μs, and of 5×10^-6 in kinetics obtained after averaging 100 shots. While previous time-resolved IR experiments have mostly been conducted on hydrated films of proteins, we demonstrate here that the brilliance of tunable quantum cascade lasers is superior to perform ns time-resolved experiments even in aqueous solution (H₂O).

Coexistence of two main folded G-quadruplexes within a single G-rich domain in the EGFR promoter

EGFR is an oncogene which codifies for a tyrosine kinase receptor that represents an important target for anticancer therapy. Indeed, several human cancers showed an upregulation of the activity of this protein. The promoter of this gene contains some G-rich domains, thus representing a yet unexplored point of intervention to potentially silence this gene. Here, we explore the conformational equilibria of a 30-nt long sequence located at position −272 (EGFR-272). By merging spectroscopic and electrophoretic analysis performed on the wild-type sequence as well as on a wide panel of related mutants, we were able to prove that in potassium ion containing solution this sequence folds into two main G-quadruplex structures, one parallel and one hybrid. They show comparable thermal stabilities and affinities for the metal ion and, indeed, they are always co-present in solution. The folding process is driven by a hairpin occurring in the domain corresponding to the terminal loop which works as an important stabilizing element for both the identified G-quadruplex arrangements.

Disordered Conformation with Low Pii Helix in Phosphoproteins Orchestrates Biomimetic Apatite Formation

The interplay between noncollagenous proteins and biomineralization is widely accepted, yet the contribution of their secondary structure in mineral formation remains to be clarified. This study demonstrates a role for phosvitin, an intrinsically disordered phosphoprotein, in chick embryo skeletal development, and using circular dichroism and matrix least-squares Henderson-Hasselbalch global fitting, unravels three distinct pH-dependent secondary structures in phosvitin. By sequestering phosvitin on a biomimetic 3D insoluble cationic framework at defined pHs, access is gained to phosvitin in various conformational states. Induction of biomimetic mineralization at near physiological conditions reveals that a disordered secondary structure with a low content of P_II helix is remarkably efficient at promoting calcium adsorption, and results in the formation of biomimetic hydroxyapatite through an amorphous calcium phosphate precursor. By extending this finding to phosphorylated full-length human recombinant dentin matrix protein-1 (17-513 AA), this bioinspired approach provides compelling evidence for the role of a disordered secondary structure in phosphoproteins in bone-like apatite formation.

Characterization of G4-G4 Crosstalk in the c-KIT Promoter Region

The proximal promoter of c-KIT contains a peculiar domain that consists of three short G-rich sequences that are close together and can fold into noncanonical DNA secondary structures called G-quadruplexes (G4). Here, we focused on a sequence containing two consecutive G4 (kit2 and kit*). By electrophoretic, surface plasmon resonance, and spectroscopic techniques, we demonstrated that they retain the ability to fold into G4 upon being inserted into the extended sequence. Here, we highlighted the occurrence of crosstalk between the two forming units. This previously unexplored G4-G4 interaction modulates both the conformation and the stability of the overall arrangement of the c-KIT promoter. It is not supported by stacking of single nucleotides but refers to a G4-G4 interaction surface surrounded by a two-nucleotides loop that might represent a reliable unprecedented target for anticancer therapy.

Remote Perturbations in Tertiary Contacts Trigger Ligation of Lysine to the Heme Iron in Cytochrome c

Perturbations in protein structure define the mechanism of allosteric regulation and biological information transfer. In cytochrome c (cyt c), ligation of Met80 to the heme iron is critical for the protein's electron-transfer (ET) function in oxidative phosphorylation and for suppressing its peroxidase activity in apoptosis. The hard base Lys is a better match for the hard ferric iron than the soft base Met is, suggesting the key role of the protein scaffold in favoring Met ligation. To probe the role of the protein structure in the maintenance of Met ligation, mutations T49V and Y67R/M80A were designed to disrupt hydrogen bonding and packing of the heme coordination loop, respectively. Electronic absorption, nuclear magnetic resonance, and electron paramagnetic resonance spectra reveal that ferric forms of both variants are Lys-ligated at neutral pH. A minor change in the tertiary contacts in T49V, away from the heme coordination loop, appears to be sufficient to execute a change in ligation, suggesting a cross-talk between the different regions of the protein structure and a possibility of built-in conformational switches in cyt c. Analyses of thermodynamic stability, kinetics of Lys binding and dissociation, and the pH-dependent changes in ligation provide a detailed characterization of the Lys coordination in these variants and relate these properties to the extent of structural perturbations. The findings emphasize the importance of the hydrogen-bonding network in controlling ligation of the native Met80 to the heme iron.

Model-independent particle species disentanglement by X-ray cross-correlation scattering

Mixtures of different particle species are often investigated using the angular averages of the scattered X-ray intensity. The number of species is deduced by singular value decomposition methods. The full disentanglement of the data into per-species contributions requires additional knowledge about the system under investigation. We propose to exploit higher-order angular X-ray intensity correlations with a new computational protocol, which we apply to synchrotron data from two-species mixtures of two-dimensional static test nanoparticles. Without any other information besides the correlations, we demonstrate the assessment of particle species concentrations in the measured data sets, as well as the full ab initio reconstruction of both particle structures. The concept extends straightforwardly to more species and to the three-dimensional case, whereby the practical application will require the measurements to be performed at an X-ray free electron laser.

Conformational Mobility in Cytochrome P450 3A4 Explored by Pressure-Perturbation EPR Spectroscopy

We used high hydrostatic pressure as a tool for exploring the conformational landscape of human cytochrome P450 3A4 (CYP3A4) by electron paramagnetic resonance and fluorescence spectroscopy. Site-directed incorporation of a luminescence resonance energy transfer donor-acceptor pair allowed us to identify a pressure-dependent equilibrium between two states of the enzyme, where an increase in pressure increased the spatial separation between the two distantly located fluorophores. This transition is characterized by volume change (ΔV°) and P1/2 values of -36.8 ± 5.0 mL/mol and 1.45 ± 0.33 kbar, respectively, which corresponds to a Keq° of 0.13 ± 0.06, so that only 15% of the enzyme adopts the pressure-promoted conformation at ambient pressure. This pressure-promoted displacement of the equilibrium is eliminated by the addition of testosterone, an allosteric activator. Using site-directed spin labeling, we demonstrated that the pressure- and testosterone-sensitive transition is also revealed by pressure-induced changes in the electron paramagnetic resonance spectra of a nitroxide side chain placed at position 85 or 409 of the enzyme. Furthermore, we observed a pressure-induced displacement of the emission maxima of a solvatochromic fluorophore (7-diethylamino-3-((((2-maleimidyl)ethyl)amino)carbonyl) coumarin) placed at the same positions, which suggests a relocation to a more polar environment. Taken together, the results reveal an effector-dependent conformational equilibrium between open and closed states of CYP3A4 that involves a pronounced change at the interface between the region of α-helices A/A' and the meander loop of the enzyme, where residues 85 and 409 are located. Our study demonstrates the high potential of pressure-perturbation strategies for studying protein conformational landscapes.

Observation of Electronic Excitation Transfer Through Light Harvesting Complex II Using Two-Dimensional Electronic-Vibrational Spectroscopy

Light-harvesting complex II (LHCII) serves a central role in light harvesting for oxygenic photosynthesis and is arguably the most important photosynthetic antenna complex. In this work, we present two-dimensional electronic-vibrational (2DEV) spectra of LHCII isolated from spinach, demonstrating the possibility of using this technique to track the transfer of electronic excitation energy between specific pigments within the complex. We assign the spectral bands via comparison with the 2DEV spectra of the isolated chromophores, chlorophyll a and b, and present evidence that excitation energy between the pigments of the complex are observed in these spectra. Finally, we analyze the essential components of the 2DEV spectra using singular value decomposition, which makes it possible to reveal the relaxation pathways within this complex.

Noninvasive Recognition and Biomarkers of Early Allergic Asthma in Cats Using Multivariate Statistical Analysis of NMR Spectra of Exhaled Breath Condensate

Asthma is prevalent in children and cats, and needs means of noninvasive diagnosis. We sought to distinguish noninvasively the differences in 53 cats before and soon after induction of allergic asthma, using NMR spectra of exhaled breath condensate (EBC). Statistical pattern recognition was improved considerably by preprocessing the spectra with probabilistic quotient normalization and glog transformation. Classification of the 106 preprocessed spectra by principal component analysis and partial least squares with discriminant analysis (PLS-DA) appears to be impaired by variances unrelated to eosinophilic asthma. By filtering out confounding variances, orthogonal signal correction (OSC) PLS-DA greatly improved the separation of the healthy and early asthmatic states, attaining 94% specificity and 94% sensitivity in predictions. OSC enhancement of multi-level PLS-DA boosted the specificity of the prediction to 100%. OSC-PLS-DA of the normalized spectra suggest the most promising biomarkers of allergic asthma in cats to include increased acetone, metabolite(s) with overlapped NMR peaks near 5.8 ppm, and a hydroxyphenyl-containing metabolite, as well as decreased phthalate. Acetone is elevated in the EBC of 74% of the cats with early asthma. The noninvasive detection of early experimental asthma, biomarkers in EBC, and metabolic perturbation invite further investigation of the diagnostic potential in humans.

Global analysis of coherence and population dynamics in 2D electronic spectroscopy

2D electronic spectroscopy is a widely exploited tool to study excited state dynamics. A high density of information is enclosed in 2D spectra. A crucial challenge is to objectively disentangle all the features of the third order optical signal. We propose a global analysis method based on the variable projection algorithm, which is able to reproduce simultaneously coherence and population dynamics of rephasing and non-rephasing contributions. Test measures at room temperature on a standard dye are used to validate the procedure and to discuss the advantages of the proposed methodology with respect to the currently employed analysis procedures.

Using Separable Nonnegative Matrix Factorization Techniques for the Analysis of Time-Resolved Raman Spectra

The key challenge of time-resolved Raman spectroscopy is the identification of the constituent species and the analysis of the kinetics of the underlying reaction network. In this work we present an integral approach that allows for determining both the component spectra and the rate constants simultaneously from a series of vibrational spectra. It is based on an algorithm for nonnegative matrix factorization that is applied to the experimental data set following a few pre-processing steps. As a prerequisite for physically unambiguous solutions, each component spectrum must include one vibrational band that does not significantly interfere with the vibrational bands of other species. The approach is applied to synthetic "experimental" spectra derived from model systems comprising a set of species with component spectra differing with respect to their degree of spectral interferences and signal-to-noise ratios. In each case, the species involved are connected via monomolecular reaction pathways. The potential and limitations of the approach for recovering the respective rate constants and component spectra are discussed.

α-Aminoxy Oligopeptides: Synthesis, Secondary Structure, and Cytotoxicity of a New Class of Anticancer Foldamers

α-Aminoxy peptides are peptidomimetic foldamers with high proteolytic and conformational stability. To gain an improved synthetic access to α-aminoxy oligopeptides we used a straightforward combination of solution- and solid-phase-supported methods and obtained oligomers that showed a remarkable anticancer activity against a panel of cancer cell lines. We solved the first X-ray crystal structure of an α-aminoxy peptide with multiple turns around the helical axis. The crystal structure revealed a right-handed 2₈ -helical conformation with precisely two residues per turn and a helical pitch of 5.8 Å. By 2D ROESY experiments, molecular dynamics simulations, and CD spectroscopy we were able to identify the 2₈ -helix as the predominant conformation in organic solvents. In aqueous solution, the α-aminoxy peptides exist in the 2₈ -helical conformation at acidic pH, but exhibit remarkable changes in the secondary structure with increasing pH. The most cytotoxic α-aminoxy peptides have an increased propensity to take up a 2₈ -helical conformation in the presence of a model membrane. This indicates a correlation between the 2₈ -helical conformation and the membranolytic activity observed in mode of action studies, thereby providing novel insights in the folding properties and the biological activity of α-aminoxy peptides.

Excitation Dynamics in Hetero-bichromophoric Calixarene Systems

In this work, the dynamics of electronic energy transfer (EET) in bichromophoric donor-acceptor systems, obtained by functionalizing a calix[4]arene scaffold with two dyes, was experimentally and theoretically characterized. The investigated compounds are highly versatile, due to the possibility of linking the dye molecules to the cone or partial cone structure of the calix[4]arene, which directs the two active units to the same or opposite side of the scaffold, respectively. The dynamics and efficiency of the EET process between the donor and acceptor units was investigated and discussed through a combined experimental and theoretical approach, involving ultrafast pump-probe spectroscopy and density functional theory based characterization of the energetic and spectroscopic properties of the system. Our results suggest that the external medium strongly determines the particular conformation adopted by the bichromophores, with a direct effect on the extent of excitonic coupling between the dyes and hence on the dynamics of the EET process itself.

The oxidative damage to the human telomere: effects of 5-hydroxymethyl-2'-deoxyuridine on telomeric G-quadruplex structures

As part of the genome, human telomeric regions can be damaged by the chemically reactive molecules responsible for oxidative DNA damage. Considering that G-quadruplex structures have been proven to occur in human telomere regions, several studies have been devoted to investigating the effect of oxidation products on the properties of these structures. However only investigations concerning the presence in G-quadruplexes of the main oxidation products of deoxyguanosine and deoxyadenosine have appeared in the literature. Here, we investigated the effects of 5-hydroxymethyl-2'-deoxyuridine (5-hmdU), one of the main oxidation products of T, on the physical-chemical properties of the G-quadruplex structures formed by two human telomeric sequences. Collected calorimetric, circular dichroism and electrophoretic data suggest that, in contrast to most of the results on other damage, the replacement of a T with a 5-hmdU results in only negligible effects on structural stability. Reported results and other data from literature suggest a possible protecting effect of the loop residues on the other parts of the G-quadruplexes.