Comparing Properties of Common Bioinorganic Ligands with Switchable Variants of Cytochrome c

Quantum chemical package Jaguar: A survey of recent developments and unique features

This paper is dedicated to the quantum chemical package Jaguar, which is commercial software developed and distributed by Schrödinger, Inc. We discuss Jaguar's scientific features that are relevant to chemical research as well as describe those aspects of the program that are pertinent to the user interface, the organization of the computer code, and its maintenance and testing. Among the scientific topics that feature prominently in this paper are the quantum chemical methods grounded in the pseudospectral approach. A number of multistep workflows dependent on Jaguar are covered: prediction of protonation equilibria in aqueous solutions (particularly calculations of tautomeric stability and pKa), reactivity predictions based on automated transition state search, assembly of Boltzmann-averaged spectra such as vibrational and electronic circular dichroism, as well as nuclear magnetic resonance. Discussed also are quantum chemical calculations that are oriented toward materials science applications, in particular, prediction of properties of optoelectronic materials and organic semiconductors, and molecular catalyst design. The topic of treatment of conformations inevitably comes up in real world research projects and is considered as part of all the workflows mentioned above. In addition, we examine the role of machine learning methods in quantum chemical calculations performed by Jaguar, from auxiliary functions that return the approximate calculation runtime in a user interface, to prediction of actual molecular properties. The current work is second in a series of reviews of Jaguar, the first having been published more than ten years ago. Thus, this paper serves as a rare milestone on the path that is being traversed by Jaguar's development in more than thirty years of its existence.

Applications of the Newly Developed Force-Field Parameters Uncover a Dynamic Nature of Ω-Loop C in the Lys-Ligated Alkaline Form of Cytochrome c

Lys-ligated cytochromes make up an emerging family of heme proteins. Density functional theory calculations on the amine/imidazole-ligated c-type ferric heme were employed to develop force-field parameters for molecular dynamics (MD) simulations of structural and dynamic features of these proteins. The new force-field parameters were applied to the alkaline form of yeast iso-1 cytochrome c to rationalize discrepancies resulting from distinct experimental conditions in prior structural studies and to provide insights into the mechanisms of the alkaline transition. Our simulations have revealed the dynamic nature of Ω-loop C in the Lys-ligated protein and its unfolding in the Lys-ligated conformer having this loop in the same position as in the native Met-ligated protein. The proximity of Tyr67 or Tyr74 to the Lys ligand of ferric heme iron suggests a possible mechanism of the backward alkaline transition where a proton donor Tyr assists in Lys dissociation. The developed force-field parameters will be useful in structural and dynamic characterization of other native or engineered Lys-ligated heme proteins.

Influence of the Interdomain Interface on Structural and Redox Properties of Multiheme Proteins

Multiheme proteins are important in energy conversion and biogeochemical cycles of nitrogen and sulfur. A diheme cytochrome c4 (c4) was used as a model to elucidate roles of the interdomain interface on properties of iron centers in its hemes A and B. Isolated monoheme domains c4-A and c4-B, together with the full-length diheme c4 and its Met-to-His ligand variants, were characterized by a variety of spectroscopic and stability measurements. In both isolated domains, the heme iron is Met/His-ligated at pH 5.0, as in the full-length c4, but becomes His/His-ligated in c4-B at higher pH. Intradomain contacts in c4-A are minimally affected by the separation of c4-A and c4-B domains, and isolated c4-A is folded. In contrast, the isolated c4-B is partially unfolded, and the interface with c4-A guides folding of this domain. The c4-A and c4-B domains have the propensity to interact even without the polypeptide linker. Thermodynamic cycles have revealed properties of monomeric folded isolated domains, suggesting that ferrous (FeII), but not ferric (FeIII) c4-A and c4-B, is stabilized by the interface. This study illustrates the effects of the interface on tuning structural and redox properties of multiheme proteins and enriches our understanding of redox-dependent complexation.

Naturally Occurring I81N Mutation in Human Cytochrome c Regulates Both Inherent Peroxidase Activity and Interactions with Neuroglobin

Human cytochrome c (hCyt c) is a crucial heme protein and plays an indispensable role in energy conversion and intrinsic apoptosis pathways. The sequence and structure of Cyt c were evolutionarily conserved and only a few naturally occurring mutants were detected in humans. Among those variable sites, position 81 was proposed to act as a peroxidase switch in the initiation stages of apoptosis. In this study, we show that Ile81 not only suppresses the intrinsic peroxidase activity but also is essential for Cyt c to interact with neuroglobin (Ngb), a potential protein partner. The kinetic assays showed that the peroxidase activity of the naturally occurring variant I81N was enhanced up to threefold under pH 5. The local stability of the Ω-loop D (residues 70-85) in the I81N variant was decreased. Moreover, the Alphafold2 program predicted that Ile81 forms stable contact with human Ngb. Meanwhile, the Ile81 to Asn81 missense mutation abolishes the interaction interface, resulting in a ∼40-fold decrease in binding affinity. These observations provide an insight into the structure-function relationship of the conserved Ile81 in vertebrate Cyt c.