Picometer-scale conformational heterogeneity separates functional from nonfunctional states of a photoreceptor protein

Spectroscopic Validation of Crystallographic Structures of a Protein Active Site by Chiroptical Spectroscopy

Out-of-plane distortions of a cofactor molecule in a protein active site are functionally important, and in photoreceptors, it has been proposed that they are crucial for spectral tuning and energy storage in photocycle intermediates. However, these subtle structural features are often beyond the grasp of structural biology. This issue is strikingly exemplified by photoactive yellow protein: its 14 independently determined crystal structures exhibit considerable differences in the dihedral angles defining the chromophore geometry, even though most of these are at excellent resolution. Here we developed a strategy to verify cofactor distortions in crystal structures by using quantum chemical calculations and chiroptical spectroscopy, particularly Raman optical activity and electronic circular dichroism spectroscopies. Based on this approach, we identify seven crystal structures with the chromophore geometries inconsistent with the experimentally observed data. The strategy implemented here promises to be widely applicable to uncovering cofactor distortions at active sites and to studies of reaction intermediates.

Highly scalable maximum likelihood and conjugate Bayesian inference for ERGMs on graph sets with equivalent vertices

The exponential family random graph modeling (ERGM) framework provides a highly flexible approach for the statistical analysis of networks (i.e., graphs). As ERGMs with dyadic dependence involve normalizing factors that are extremely costly to compute, practical strategies for ERGMs inference generally employ a variety of approximations or other workarounds. Markov Chain Monte Carlo maximum likelihood (MCMC MLE) provides a powerful tool to approximate the maximum likelihood estimator (MLE) of ERGM parameters, and is generally feasible for typical models on single networks with as many as a few thousand nodes. MCMC-based algorithms for Bayesian analysis are more expensive, and high-quality answers are challenging to obtain on large graphs. For both strategies, extension to the pooled case—in which we observe multiple networks from a common generative process—adds further computational cost, with both time and memory scaling linearly in the number of graphs. This becomes prohibitive for large networks, or cases in which large numbers of graph observations are available. Here, we exploit some basic properties of the discrete exponential families to develop an approach for ERGM inference in the pooled case that (where applicable) allows an arbitrarily large number of graph observations to be fit at no additional computational cost beyond preprocessing the data itself. Moreover, a variant of our approach can also be used to perform Bayesian inference under conjugate priors, again with no additional computational cost in the estimation phase. The latter can be employed either for single graph observations, or for observations from graph sets. As we show, the conjugate prior is easily specified, and is well-suited to applications such as regularization. Simulation studies show that the pooled method leads to estimates with good frequentist properties, and posterior estimates under the conjugate prior are well-behaved. We demonstrate the usefulness of our approach with applications to pooled analysis of brain functional connectivity networks and to replicated x-ray crystal structures of hen egg-white lysozyme.

Structural Basis of Design and Engineering for Advanced Plant Optogenetics

In optogenetics, light-sensitive proteins are specifically expressed in target cells and light is used to precisely control the activity of these proteins at high spatiotemporal resolution. Optogenetics initially used naturally occurring photoreceptors to control neural circuits, but has expanded to include carefully designed and engineered photoreceptors. Several optogenetic constructs are based on plant photoreceptors, but their application to plant systems has been limited. Here, we present perspectives on the development of plant optogenetics, considering different levels of design complexity. We discuss how general principles of light-driven signal transduction can be coupled with approaches for engineering protein folding to develop novel optogenetic tools. Finally, we explore how the use of computation, networks, circular permutation, and directed evolution could enrich optogenetics.

Spectroscopic ruler for measuring active-site distortions based on Raman optical activity of a hydrogen out-of-plane vibration

Photoactive yellow protein (PYP), from the phototrophic bacterium Halorhodospira halophila, is a small water-soluble photoreceptor protein and contains p-coumaric acid (pCA) as a chromophore. PYP has been an attractive model for studying the physical chemistry of protein active sites. Here, we explore how Raman optical activity (ROA) can be used to extract quantitative information on distortions of the pCA chromophore at the active site in PYP. We use ¹³C8-pCA to assign an intense signal at 826 cm^-1 in the ROA spectrum of PYP to a hydrogen out-of-plane vibration of the ethylenic moiety of the chromophore. Quantum-chemical calculations based on density functional theory demonstrate that the sign of this ROA band reports the direction of the distortion in the dihedral angle about the ethylenic C=C bond, while its amplitude is proportional to the dihedral angle. These results document the ability of ROA to quantify structural deformations of a cofactor molecule embedded in a protein moiety.

Correlating structural and photochemical heterogeneity in cyanobacteriochrome NpR6012g4

Phytochrome photoreceptors control plant growth, development, and the shade avoidance response that limits crop yield in high-density agricultural plantings. Cyanobacteriochromes (CBCRs) are distantly related photosensory proteins that control cyanobacterial metabolism and behavior in response to light. Photoreceptors in both families reversibly photoconvert between two photostates via photoisomerization of linear tetrapyrrole (bilin) chromophores. Spectroscopic and biochemical studies have demonstrated heterogeneity in both photostates, but the structural basis for such heterogeneity remains unclear. We report solution NMR structures for both photostates of the red/green CBCR NpR6012g4 from Nostoc punctiforme In addition to identifying structural changes accompanying photoconversion, these structures reveal structural heterogeneity for residues Trp655 and Asp657 in the red-absorbing NpR6012g4 dark state, yielding two distinct environments for the phycocyanobilin chromophore. We use site-directed mutagenesis and fluorescence and absorbance spectroscopy to assign an orange-absorbing population in the NpR6012g4 dark state to the minority configuration for Asp657. This population does not undergo full, productive photoconversion, as shown by time-resolved spectroscopy and absorption spectroscopy at cryogenic temperature. Our studies thus elucidate the spectral and photochemical consequences of structural heterogeneity in a member of the phytochrome superfamily, insights that should inform efforts to improve photochemical or fluorescence quantum yields in the phytochrome superfamily.

Excitation-Wavelength-Dependent Photocycle Initiation Dynamics Resolve Heterogeneity in the Photoactive Yellow Protein from Halorhodospira halophila

Photoactive yellow proteins (PYPs) make up a diverse class of blue-light-absorbing bacterial photoreceptors. Electronic excitation of the p-coumaric acid chromophore covalently bound within PYP results in triphasic quenching kinetics; however, the molecular basis of this behavior remains unresolved. Here we explore this question by examining the excitation-wavelength dependence of the photodynamics of the PYP from Halorhodospira halophila via a combined experimental and computational approach. The fluorescence quantum yield, steady-state fluorescence emission maximum, and cryotrapping spectra are demonstrated to depend on excitation wavelength. We also compare the femtosecond photodynamics in PYP at two excitation wavelengths (435 and 475 nm) with a dual-excitation-wavelength-interleaved pump-probe technique. Multicompartment global analysis of these data demonstrates that the excited-state photochemistry of PYP depends subtly, but convincingly, on excitation wavelength with similar kinetics with distinctly different spectral features, including a shifted ground-state beach and altered stimulated emission oscillator strengths and peak positions. Three models involving multiple excited states, vibrationally enhanced barrier crossing, and inhomogeneity are proposed to interpret the observed excitation-wavelength dependence of the data. Conformational heterogeneity was identified as the most probable model, which was supported with molecular mechanics simulations that identified two levels of inhomogeneity involving the orientation of the R52 residue and different hydrogen bonding networks with the p-coumaric acid chromophore. Quantum calculations were used to confirm that these inhomogeneities track to altered spectral properties consistent with the experimental results.

The Low Barrier Hydrogen Bond in the Photoactive Yellow Protein: A Vacuum Artifact Absent in the Crystal and Solution

There has been considerable debate on the existence of a low-barrier hydrogen bond (LBHB) in the photoactive yellow protein (PYP). The debate was initially triggered by the neutron diffraction study of Yamaguchi et al. ( Proc. Natl. Acad. Sci., U. S. A. , 2009 , 106 , 440 - 444 ) who suggested a model in which a neutral Arg52 residue triggers the formation of the LBHB in PYP. Here, we present an alternative model that is consistent within the error margins of the Yamaguchi structure factors. The model explains an increased hydrogen bond length without nuclear quantum effects and for a protonated Arg52. We tested both models by calculations under crystal, solution, and vacuum conditions. Contrary to the common assumption in the field, we found that a single PYP in vacuum does not provide an accurate description of the crystal conditions but instead introduces strong artifacts, which favor a LBHB and a large ¹H NMR chemical shift. Our model of the crystal environment was found to stabilize the two Arg52 hydrogen bonds and crystal water positions for the protonated Arg52 residue in free MD simulations and predicted an Arg52 pK_a upshift with respect to PYP in solution. The crystal and solution environments resulted in almost identical ¹H chemical shifts that agree with NMR solution data. We also calculated the effect of the Arg52 protonation state on the LBHB in 3D nuclear equilibrium density calculations. Only the charged crystal structure in vacuum supports a LBHB if Arg52 is neutral in PYP at the previously reported level of theory ( J. Am. Chem. Soc. , 2014 , 136 , 3542 - 3552 ). We attribute the anomalies in the interpretation of the neutron data to a shift of the potential minimum, which does not involve nuclear quantum effects and is transferable beyond the Yamaguchi structure.

Accessing protein conformational ensembles using room-temperature X-ray crystallography

Modern protein crystal structures are based nearly exclusively on X-ray data collected at cryogenic temperatures (generally 100 K). The cooling process is thought to introduce little bias in the functional interpretation of structural results, because cryogenic temperatures minimally perturb the overall protein backbone fold. In contrast, here we show that flash cooling biases previously hidden structural ensembles in protein crystals. By analyzing available data for 30 different proteins using new computational tools for electron-density sampling, model refinement, and molecular packing analysis, we found that crystal cryocooling remodels the conformational distributions of more than 35% of side chains and eliminates packing defects necessary for functional motions. In the signaling switch protein, H-Ras, an allosteric network consistent with fluctuations detected in solution by NMR was uncovered in the room-temperature, but not the cryogenic, electron-density maps. These results expose a bias in structural databases toward smaller, overpacked, and unrealistically unique models. Monitoring room-temperature conformational ensembles by X-ray crystallography can reveal motions crucial for catalysis, ligand binding, and allosteric regulation.

Proline 68 enhances photoisomerization yield in photoactive yellow protein

In proteins and enzymes, the local environment of an active cofactor plays an important role in controlling the outcome of a functional reaction. In photoactive yellow protein (PYP), it ensures photoisomerization of the chromophore, a prerequisite for formation of a signaling state. PYP is the prototype of a PAS domain, and the preferred model system for the studies of molecular mechanisms of biological light sensing. We investigated the effect of replacing proline-68, positioned near but not in direct contact with the chromophore, with other neutral amino acids (alanine, glycine, and valine), using ultrafast spectroscopy probing the visible and the mid-IR spectral regions, and molecular simulation to understand the interactions tuning the efficiency of light signaling. Transient absorption measurements indicate that the quantum yield of isomerization in the mutants is lower than the yield observed for the wild type. Subpicosecond mid-IR spectra and molecular dynamics simulations of the four proteins reveal that the hydrogen bond interactions around the chromophore and the access of water molecules in the active site of the protein determine the efficiency of photoisomerization. The mutants provide additional hydrogen bonds to the chromophore, directly and by allowing more water molecules access to its binding pocket. We conclude that proline-68 in the wild type protein optimizes the yield of photochemistry by maintaining a weak hydrogen bond with the chromophore, at the same time restraining the entrance of water molecules close to the alkylic part of pCa. This study provides a molecular basis for the structural optimization of biological light sensing.

Characterizing rare-event property distributions via replicate molecular dynamics simulations of proteins

As computational resources increase, molecular dynamics simulations of biomolecules are becoming an increasingly informative complement to experimental studies. In particular, it has now become feasible to use multiple initial molecular configurations to generate an ensemble of replicate production-run simulations that allows for more complete characterization of rare events such as ligand-receptor unbinding. However, there are currently no explicit guidelines for selecting an ensemble of initial configurations for replicate simulations. Here, we use clustering analysis and steered molecular dynamics simulations to demonstrate that the configurational changes accessible in molecular dynamics simulations of biomolecules do not necessarily correlate with observed rare-event properties. This informs selection of a representative set of initial configurations. We also employ statistical analysis to identify the minimum number of replicate simulations required to sufficiently sample a given biomolecular property distribution. Together, these results suggest a general procedure for generating an ensemble of replicate simulations that will maximize accurate characterization of rare-event property distributions in biomolecules.