Phytochromes are biological red/far-red light sensors found in many organisms. Photoisomerization of the linear methine-bridged tetrapyrrole triggers transient proton translocation events in the chromophore binding pocket (CBP) leading to major conformational changes of the protein matrix that are in turn associated with signaling. By combining pH-dependent resonance Raman and UV–visible absorption spectroscopy, we analyzed protonation-dependent equilibria in the CBP of Cph1 involving the proposed Pr-I and Pr-II substates that prevail below and above pH 7.5, respectively. The protonation pattern and vibrational properties of these states were further characterized by means of hybrid quantum mechanics/molecular mechanics calculations. From this combined experimental–theoretical study, we were able to identify His260 as the key residue controlling pH-dependent equilibria. This residue is not only responsible for the conformational heterogeneity of CBP in the Pr state of prokaryotic phytochromes, discussed extensively in the past, but it constitutes the sink and source of protons in the proton release/uptake mechanism involving the tetrapyrrole chromophore which finally leads to the formation of the Pfr state. Thus, this work provides valuable information that may guide further experiments toward the understanding of the specific role of protons in controlling structure and function of phytochromes in general.

Carbon monoxide dehydrogenase (CODH) is a key enzyme for reversible CO interconversion. To elucidate structural and mechanistic details of CO binding at the CODH active site (C-cluster), cyanide is frequently used as an iso-electronic substitute and inhibitor. However, previous studies revealed conflicting results on the structure of the cyanide-bound complex and the mechanism of cyanide-inhibition. To address this issue in this work, we have employed IR spectroscopy, crystallography, site directed mutagenesis, and theoretical methods to analyse the cyanide complex of the CODH from Carboxydothermus hydrogenoformans (CODHIICh). IR spectroscopy demonstrates that a single cyanide binds to the Ni ion. Whereas the inhibitor could be partially removed at elevated temperature, irreversible degradation of the C-cluster occurred in the presence of an excess of cyanide on the long-minute time scale, eventually leading to the formation of [Fe(CN)6]4− and [Ni(CN)4]2− complexes. Theoretical calculations based on a new high-resolution structure of the cyanide-bound CODHIICh indicated that cyanide binding to the Ni ion occurs upon dissociation of the hydroxyl ligand from the Fe1 subsite of the C-cluster. The hydroxyl group is presumably protonated by Lys563 which, unlike to His93, does not form a hydrogen bond with the cyanide ligand. A stable deprotonated ε-amino group of Lys563 in the cyanide complex is consistent with the nearly unchanged CN stretching in the Lys563Ala variant of CODHIICh. These findings support the view that the proton channel connecting the solution phase with the active site displays a strict directionality, controlled by the oxidation state of the C-cluster.

Understanding the interaction and immobilization of [NiFe] hydrogenases on functionalized surfaces is important in the field of biotechnology and, in particular, for the development of biofuel cells. In this study, we investigated the adsorption behavior of the standard [NiFe] hydrogenase of Desulfovibrio gigas on amino-terminated alkanethiol self-assembled monolayers (SAMs) with different levels of protonation. Classical all-atom molecular dynamics (MD) simulations revealed a strong correlation between the adsorption behavior and the level of ionization of the chemically modified electrode surface. While the hydrogenase undergoes a weak but stable initial adsorption process on SAMs with a low degree of protonation, a stronger immobilization is observable on highly ionized SAMs, affecting protein reorientation and conformation. These results were validated by complementary surface-enhanced infrared absorption (SEIRA) measurements on the comparable [NiFe] standard hydrogenases from Desulfovibrio vulgaris Miyazaki F and allowed in this way for a detailed insight into the adsorption mechanism at the atomic level.

Sulfite oxidase (SO) is an enzyme catalyzing the terminal step of the metabolism of sulfur-containing amino acids that is essential for almost all living organisms. The catalytic activity of SO in vertebrates strongly depends on the efficiency of the intramolecular electron transfer (IET) between the catalytic Moco domain and the cytochrome b5 (cyt b5) domain. The IET process is assumed to be mediated by large domain motions of the cyt b5 domains within the enzyme. Thus, the interaction of SO with charged surfaces may affect the mobility of the cyt b5 domain required for IET and consequently hinder SO activation. In this study, we present a molecular dynamics approach to investigating the ionic strength dependence of the initial surface adsorption of SO in two different conformations—the crystallographic structure and the model structure for an activated SO—onto mixed amino- and hydroxyl-terminated SAMs. The results show for both conformations at low ionic strengths a strong adsorption of the cyt b5 units onto the SAM, which inhibits the domain motion event required for IET. Under higher ion concentrations, however, the interaction with the surface is weakened by the negatively charged ions acting as a buffer and competing in adsorption with the cathodic cyt b5 domains. This competition prevents the immobilization of the cytochrome b5 units onto the surface, allowing the intramolecular domain motions favoring IET. Our predictions support the interpretation of recent experimental spectroelectrochemical studies on SO.

A homology structural model was generated for plant phytochrome phyA utilizing the crystal structure of the sensory module of cyanobacterial phytochrome Cph1 (Cph1Δ2). As chromophores, either the native phytochromobilin cofactor (PΦB) or phycocyanobilin (PCB), the natural cofactor in Cph1, was incorporated. These homology models were further optimized by molecular dynamics (MD) simulations revealing a satisfying overall agreement with the crystal structure of Cph1Δ2. Notable differences in the PΦB adduct of phyA result from a restructuring of the small helical segment α7 that leads to displacements of a few amino acids away from the cofactor. This repositioning of residues also include aspartate 218 such that, instead of its carbonyl function as in Cph1Δ2, an additional water molecule forms hydrogen bonds with the ring B and C NH groups. To validate the phyA structural model in the chromophore binding pocket, Raman spectra of the cofactor were calculated by means of the quantum mechanics/molecular mechanics (QM/MM) hybrid methodology and compared with the experimental resonance Raman (RR) spectra. The satisfactory overall agreement between calculated and experimental spectra is taken as an indication for the good quality of the structural model. Moreover, the methine bridge stretching modes and the effects of isotopic labeling at selected positions of the chromophore are very well reproduced to allow confirming even details of the methine bridge geometry as predicted by the homology model. Specifically, it is demonstrated that the experimental RR spectra are consistent with a torsional angle of ring D with respect to ring C that is distinctly higher for phyA−PCB (45°) and phyA−PΦB (42°) than for Cph1Δ2 (30°). Raman spectra calculated from different points of the MD trajectory display variations of the mode frequencies and intensities reflecting the structural fluctuations from snapshot to snapshot. The snapshot spectrum of the lowest energy structure and the sum of all snapshot spectra afford an equally good description of the experimental data. Particularly large variations between the snapshots are noted for the N−H in-plane bending mode of the pyrrole rings B and C, which reflect alterations of the hydrogen bond interactions brought about by fluctuations of water molecules in the cofactor cavity. This overestimation of the water molecule mobility is a consequence of the deficiency of the current QM/MM methodology that, due to the lack of appropriate protein force fields, cannot adequately account for the electrostatics in the cofactor pocket.

Bone morphogenetic protein-2 (BMP-2) plays a crucial role in osteoblast differentiation and proliferation. Its effective therapeutic use for ectopic bone and cartilage regeneration depends, among other factors, on the interaction with the carrier at the implant site. In this study, we used classical molecular dynamics (MD) and a hybrid approach of steered molecular dynamics (SMD) combined with MD simulations to investigate the initial stages of the adsorption of BMP-2 when approaching two implant surfaces, hydrophobic graphite and hydrophilic titanium dioxide rutile. Surface adsorption was evaluated for six different orientations of the protein, two end-on and four side-on, in explicit water environment. On graphite, we observed a weak but stable adsorption. Depending on the initial orientation, hydrophobic patches as well as flexible loops of the protein were involved in the interaction with graphite. On the contrary, BMP-2 adsorbed only loosely to hydrophilic titanium dioxide. Despite a favorable interaction energy between protein and the TiO2 surface, the rapid formation of a two-layer water structure prevented the direct interaction between protein and titanium dioxide. The first water adlayer had a strong repulsive effect on the protein, while the second attracted the protein toward the surface. For both surfaces, hydrophobic graphite and hydrophilic titanium dioxide, denaturation of BMP-2 induced by adsorption was not observed on the nanosecond time scale.

The Self-Consistent Charge Density Functional Tight Binding (SCC-DFTB) method has been extended for the calculation of vibrational Raman spectra employing the Fourier Transform of Time-Correlation Function (FTTCF) formalism. As Witek and co-workers have already shown for a set of various organic molecules, the minimal basis SCC-DFTB approach performs surprisingly good in terms of polarizability calculations. Therefore, we were encouraged to use this electronic structure method for the purpose of Raman spectra calculations via FTTCF. The molecular polarizability was accessed via second order numeric derivatives of the SCC-DFTB energy with respect to the components of an external electric field “on-the-fly” during a molecular dynamics (MD) simulation. The finite electric field approach delivers Raman spectra that are in overall good agreement for most of 10 small organic model compounds examined in the gas phase compared to a standard Normal Mode Analysis (NMA) approach at the same (SCC-DFTB) and at a higher level of theory (BLYP aug-cc-pVTZ). With the use of reparametrized SCC-DFTB repulsive potentials, a distinct improvement of the Raman spectra from the SCC-DFTB/FTTCF protocol of conjugated hydrocarbons has been observed. Further QM/MM test calculations of l-phenylalanine in aqueous solution revealed larger deviations concerning vibrational frequencies and relative intensities for several stretching and bending modes in the benzene ring as compared to experimental results. Our SCC-DFTB/FTTCF approach was also tested against a hybrid method, in which polarizability calculations at the B3YLP 6-31G(d) level were performed on a trajectory at the SCC-DFTB level. We found that our SCC-DFTB/FTTCF protocol is not only much more efficient but in terms of the resulting Raman spectra also of similar accuracy compared to the hybrid approach. In our opinion, the more accurately calculated polarizabilities at the B3YLP 6-31G(d) level cannot compensate for the usually insufficient sampling of phase space when employing high level QM methods in a FTTCF framework.

A quantum mechanics (QM)/molecular mechanics (MM) hybrid method was applied to the Pr state of the cyanobacterial phytochrome Cph1 to calculate the Raman spectra of the bound PCB cofactor. Two QM/MM models were derived from the atomic coordinates of the crystal structure. The models differed in the protonation site of His260 in the chromophore-binding pocket such that either the δ-nitrogen (M-HSD) or the ɛ-nitrogen (M-HSE) carried a hydrogen. The optimized structures of the two models display small differences specifically in the orientation of His260 with respect to the PCB cofactor and the hydrogen bond network at the cofactor-binding site. For both models, the calculated Raman spectra of the cofactor reveal a good overall agreement with the experimental resonance Raman (RR) spectra obtained from Cph1 in the crystalline state and in solution, including Cph1 adducts with isotopically labeled PCB. However, a distinctly better reproduction of important details in the experimental spectra is provided by the M-HSD model, which therefore may represent an improved structure of the cofactor site. Thus, QM/MM calculations of chromoproteins may allow for refining crystal structure models in the chromophore-binding pocket guided by the comparison with experimental RR spectra. Analysis of the calculated and experimental spectra also allowed us to identify and assign the modes that sensitively respond to chromophore-protein interactions. The most pronounced effect was noted for the stretching mode of the methine bridge A-B adjacent to the covalent attachment site of PCB. Due a distinct narrowing of the A-B methine bridge bond angle, this mode undergoes a large frequency upshift as compared with the spectrum obtained by QM calculations for the chromophore in vacuo. This protein-induced distortion of the PCB geometry is the main origin of a previous erroneous interpretation of the RR spectra based on QM calculations of the isolated cofactor.