Dynamic Asymmetry and the Role of the Conserved Active-Site Thiol in Rabbit Muscle Creatine Kinase

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

Vibrational Lifetime of the SCN Protein Label in H2O and D2O Reports Site-Specific Solvation and Structure Changes During PYP's Photocycle

The application of vibrational labels such as thiocyanate  (-S-C≡N) for studying protein structure and dynamics is thriving. Absorption spectroscopy is usually employed to obtain wavenumber and line shape of the label. An observable of great significance might be the vibrational lifetime, which can be obtained by pump probe or 2D-IR spectroscopy. Due to the insulating effect of the heavy sulfur atom in the case of the SCN label, the lifetime of the C≡N oscillator is expected to be particularly sensitive to its surrounding as it is not dominated by through-bond relaxation. We therefore investigate the vibrational lifetime of the SCN label at various positions in the blue light sensor protein Photoactive Yellow Protein (PYP) in the ground state and signaling state of the photoreceptor. We find that the vibrational lifetime of the C≡N stretching mode is strongly affected both by its protein environment and by the degree of exposure to the solvent. Even for label positions where the line shape and wavenumber observed by FTIR are barely changing upon activation of the photoreceptor, we find that the lifetime can change considerably. To obtain an unambiguous measure for the solvent exposure of the labeled site, we show that it is imperative to compare the lifetimes in H₂O and D₂O. Importantly, the lifetimes shorten in H₂O as compared to D₂O for water exposed labels, while they stay largely the same for buried labels. We quantify this effect by defining a solvent exclusion coefficient (SEC). The response of the label's vibrational lifetime to its solvent exposure renders it a suitable universal probe for protein investigations. This applies even to systems that are otherwise hard to address, such as transient or short-lived states, which could be created during a protein's working cycle (as here in PYP) or during protein folding. It is also applicable to flexible systems (intrinsically disordered proteins), protein-protein and protein-membrane interactions.

Following local light-induced structure changes and dynamics of the photoreceptor PYP with the thiocyanate IR label

Photoactive Yellow Protein (PYP) is a bacterial blue light receptor that enters a photocycle after excitation. The intermediate states are formed on time scales ranging from femtoseconds up to hundreds of milliseconds, after which the signaling state with a lifetime of about 1 s is reached. To investigate structural changes and dynamics, we incorporated the SCN IR label at distinct positions of the photoreceptor via cysteine mutation and cyanylation. FT-IR measurements of the SCN label at different sites of the well-established dark state structure of PYP characterized the spectral response of the label to differences in the environment. Under constant blue light irradiation, we observed the formation of the signaling state with significant changes of wavenumber and lineshape of the SCN bands. Thereby we deduced light-induced structural changes in the local environment of the labels. These results were supported by molecular dynamics simulations on PYP providing the solvent accessible surface area (SASA) at the different positions. To follow protein dynamics via the SCN label during the photocycle, we performed step-scan FT-IR measurements with a time resolution of 10 μs. Global analysis yielded similar time constants of τ1 = 70 μs, τ2 = 640 μs, and τ3 > 20 ms for the wild type and τ1 = 36 μs, τ2 = 530 μs, and τ3 > 20 ms for the SCN-labeled mutant PYP-A44C*, a mutant which provided a sufficiently large SCN difference signal to measure step-scan FT-IR spectra. In comparison to the protein (amide, E46) and chromophore bands the dynamics of the SCN label show a different behavior. This result indicates that the local kinetics sensed by the label are different from the global protein kinetics.

Mitochondrial Proteolipid Complexes of Creatine Kinase

Isoforms of creatine kinase (CK) generate and use phosphocreatine, a concentrated and highly diffusible cellular "high energy" intermediate, for the main purpose of energy buffering and transfer in order to maintain cellular energy homeostasis. The mitochondrial CK isoform (mtCK) localizes to the mitochondrial intermembrane and cristae space, where it assembles into peripherally membrane-bound, large cuboidal homooctamers. These are part of proteolipid complexes wherein mtCK directly interacts with cardiolipin and other anionic phospholipids, as well as with the VDAC channel in the outer membrane. This leads to a stabilization and cross-linking of inner and outer mitochondrial membrane, forming so-called contact sites. Also the adenine nucleotide translocator of the inner membrane can be recruited into these proteolipid complexes, probably mediated by cardiolipin. The complexes have functions mainly in energy transfer to the cytosol and stimulation of oxidative phosphorylation, but also in restraining formation of reactive oxygen species and apoptosis. In vitro evidence indicates a putative role of mtCK in mitochondrial phospholipid distribution, and most recently a role in thermogenesis has been proposed. This review summarizes the essential structural and functional data of these mtCK complexes and describes in more detail the more recent advances in phospholipid interaction, thermogenesis, cancer and evolution of mtCK.

Cyanylated Cysteine Reports Site-Specific Changes at Protein-Protein-Binding Interfaces Without Perturbation

To investigate the cyanylated cysteine vibrational probe group’s ability to report on binding-induced changes along a protein–protein interface, the probe group was incorporated at several sites in a peptide of the calmodulin (CaM)-binding domain of skeletal muscle myosin light chain kinase. Isothermal titration calorimetry was used to determine the binding thermodynamics between calmodulin and each peptide. For all probe positions, the binding affinity was nearly identical to that of the unlabeled peptide. The CN stretching infrared band was collected for each peptide free in solution and bound to calmodulin. Binding-induced shifts in the IR spectral frequencies were correlated with estimated solvent accessibility based on molecular dynamics simulations. This work generally suggests (1) that site-specific incorporation of this vibrational probe group does not cause major perturbations to its local structural environment and (2) that this small probe group might be used quite broadly to map dynamic protein-binding interfaces. However, site-specific perturbations due to artificial labeling groups can be somewhat unpredictable and should be evaluated on a site-by-site basis through complementary measurements. A fully quantitative, simulation-based interpretation of the rich probe IR spectra is still needed but appears to be possible given recent advances in simulation techniques.

A Direct, Quantitative Connection between Molecular Dynamics Simulations and Vibrational Probe Line Shapes

A quantitative connection between molecular dynamics simulations and vibrational spectroscopy of probe-labeled systems would enable direct translation of experimental data into structural and dynamical information. To constitute this connection, all-atom molecular dynamics (MD) simulations were performed for two SCN probe sites (solvent-exposed and buried) in a calmodulin-target peptide complex. Two frequency calculation approaches with substantial nonelectrostatic components, a quantum mechanics/molecular mechanics (QM/MM)-based technique and a solvatochromic fragment potential (SolEFP) approach, were used to simulate the infrared probe line shapes. While QM/MM results disagreed with experiment, SolEFP results matched experimental frequencies and line shapes and revealed the physical and dynamic bases for the observed spectroscopic behavior. The main determinant of the CN probe frequency is the exchange repulsion between the probe and its local structural neighbors, and there is a clear dynamic explanation for the relatively broad probe line shape observed at the "buried" probe site. This methodology should be widely applicable to vibrational probes in many environments.

Dynamic traction of lattice-confined platinum atoms into mesoporous carbon matrix for hydrogen evolution reaction

Constructing atomically dispersed platinum (Pt) electrocatalysts is essential to build high-performance and cost-effective electrochemical water-splitting systems. We present a novel strategy to realize the traction and stabilization of isolated Pt atoms in the nitrogen-containing porous carbon matrix (Pt@PCM). In comparison with the commercial Pt/C catalyst (20 weight %), the as-prepared Pt@PCM catalyst exhibits significantly boosted mass activity (up to 25 times) for hydrogen evolution reaction. Results of extended x-ray absorption fine structure investigation and density functional theory calculation suggest that the active sites are associated with the lattice-confined Pt centers and the activated carbon (C)/nitrogen (N) atoms at the adjacency of the isolated Pt centers. This strategy may provide insights into constructing highly efficient single-atom catalysts for different energy-related applications.

Vibrational Probes: From Small Molecule Solvatochromism Theory and Experiments to Applications in Complex Systems

The vibrational frequency of a chosen normal mode is one of the most accurately measurable spectroscopic properties of molecules in condensed phases. Accordingly, infrared absorption and Raman scattering spectroscopy have provided valuable information on both distributions and ensemble-average values of molecular vibrational frequencies, and these frequencies are now routinely used to investigate structure, conformation, and even absolute configuration of chemical and biological molecules of interest. Recent advancements in coherent time-domain nonlinear vibrational spectroscopy have allowed the study of heterogeneous distributions of local structures and thermally driven ultrafast fluctuations of vibrational frequencies. To fully utilize IR probe functional groups for quantitative bioassays, a variety of biological and chemical techniques have been developed to site-specifically introduce vibrational probe groups into proteins and nucleic acids. These IR-probe-labeled biomolecules and chemically reactive systems are subject to linear and nonlinear vibrational spectroscopic investigations and provide information on the local electric field, conformational changes, site-site protein contacts, and/or function-defining features of biomolecules. A rapidly expanding library of data from such experiments requires an interpretive method with atom-level chemical accuracy. However, despite prolonged efforts to develop an all-encompassing theory for describing vibrational solvatochromism and electrochromism as well as dynamic fluctuations of instantaneous vibrational frequencies, purely empirical and highly approximate theoretical models have often been used to interpret experimental results. They are, in many cases, based on the simple assumption that the vibrational frequency of an IR reporter is solely dictated by electric potential or field distribution around the vibrational chromophore. Such simplified description of vibrational solvatochromism generally referred to as vibrational Stark effect theory has been considered to be quite appealing and, even in some cases, e.g., carbonyl stretch modes in amide, ester, ketone, and carbonate compounds or proteins, it works quantitatively well, which makes it highly useful in determining the strength of local electric field around the IR chromophore. However, noting that the vibrational frequency shift results from changes of solute-solvent intermolecular interaction potential along its normal coordinate, Pauli exclusion repulsion, polarization, charge transfer, and dispersion interactions, in addition to the electrostatic interaction between distributed charges of both vibrational chromophore and solvent molecules, are to be properly included in the theoretical description of vibrational solvatochromism. Since the electrostatic and nonelectrostatic intermolecular interaction components have distinctively different distance and orientation dependences, they affect the solvatochromic vibrational properties in a completely different manner. Over the past few years, we have developed a systematic approach to simulating vibrational solvatochromic data based on the effective fragment potential approach, one of the most accurate and rigorous theories on intermolecular interactions. We have further elucidated the interplay of local electric field with the general vibrational solvatochromism of small IR probes in either solvents or complicated biological systems, with emphasis on contributions from non-Coulombic intermolecular interactions to vibrational frequency shifts and fluctuations. With its rigorous foundation and close relation to quantitative interpretation of experimental data, this and related theoretical approaches and experiments will be of use in studying and quantifying the structure and dynamics of biomolecules with unprecedented time and spatial resolution when combined with time-resolved vibrational spectroscopy and chemically sensitive vibrational imaging techniques.

Probing the effectiveness of spectroscopic reporter unnatural amino acids: a structural study

The X-ray crystal structures of superfolder green fluorescent protein (sfGFP) containing the spectroscopic reporter unnatural amino acids (UAAs) 4-cyano-L-phenylalanine (pCNF) or 4-ethynyl-L-phenylalanine (pCCF) at two unique sites in the protein have been determined. These UAAs were genetically incorporated into sfGFP in a solvent-exposed loop region and/or a partially buried site on the β-barrel of the protein. The crystal structures containing the UAAs at these two sites permit the structural implications of UAA incorporation for the native protein structure to be assessed with high resolution and permit a direct correlation between the structure and spectroscopic data to be made. The structural implications were quantified by comparing the root-mean-square deviation (r.m.s.d.) between the crystal structure of wild-type sfGFP and the protein constructs containing either pCNF or pCCF in the local environment around the UAAs and in the overall protein structure. The results suggest that the selective placement of these spectroscopic reporter UAAs permits local protein environments to be studied in a relatively nonperturbative fashion with site-specificity.

Azidoethoxyphenylalanine as a Vibrational Reporter and Click Chemistry Partner in Proteins

An unnatural amino acid, 4-(2-azidoethoxy)-L-phenylalanine (AePhe, 1), was designed and synthesized in three steps from known compounds in 54% overall yield. The sensitivity of the IR absorption of the azide of AePhe was established by comparison of the frequency of the azide asymmetric stretch vibration in water and dimethyl sulfoxide. AePhe was successfully incorporated into superfolder green fluorescent protein (sfGFP) at the 133 and 149 sites by using the amber codon suppression method. The IR spectra of these sfGFP constructs indicated that the azide group at the 149 site was not fully solvated despite the location in sfGFP and the three-atom linker between the azido group and the aromatic ring of AePhe. An X-ray crystal structure of sfGFP-149-AePhe was solved at 1.45 Å resolution and provides an explanation for the IR data as the flexible linker adopts a conformation which partially buries the azide on the protein surface. Both sfGFP-AePhe constructs efficiently undergo a bioorthogonal strain-promoted click cycloaddition with a dibenzocyclooctyne derivative.

The substrate-free and -bound crystal structures of the duplicated taurocyamine kinase from the human parasite Schistosoma mansoni

The taurocyamine kinase from the blood fluke Schistosoma mansoni (SmTK) belongs to the phosphagen kinase (PK) family and catalyzes the reversible Mg(2+)-dependent transfer of a phosphoryl group between ATP and taurocyamine. SmTK is derived from gene duplication, as are all known trematode TKs. Our crystallographic study of SmTK reveals the first atomic structure of both a TK and a PK with a bilobal structure. The two unliganded lobes present a canonical open conformation and interact via their respective C- and N-terminal domains at a helix-mediated interface. This spatial arrangement differs from that observed in true dimeric PKs, in which both N-terminal domains make contact. Our structures of SmTK complexed with taurocyamine or l-arginine compounds explain the mechanism by which an arginine residue of the phosphagen specificity loop is crucial for substrate specificity. An SmTK crystal was soaked with the dead end transition state analog (TSA) components taurocyamine-NO3 (2-)-MgADP. One SmTK monomer was observed with two bound TSAs and an asymmetric conformation, with the first lobe semiclosed and the second closed. However, isothermal titration calorimetry and enzyme kinetics experiments showed that the two lobes function independently. A small angle x-ray scattering model of SmTK-TSA in solution with two closed active sites was generated.