Site Of Cytochrome Research Articles

QSAR modelling, molecular docking, molecular dynamic and ADMET prediction of pyrrolopyrimidine derivatives as novel Bruton’s tyrosine kinase (BTK) inhibitors

In recent years, there has been a focus on developing and discovering novel Bruton's tyrosine kinase (BTK) inhibitors, as they offer an effective treatment strategy for B-cell malignancies. BTK plays a crucial role in B cell receptor (BCR)-mediated activation and proliferation by regulating downstream factors such as the NF-κB and MAP kinase pathways. To address this challenge and propose potential therapeutic options for B-cell lymphomas, researchers conducted 2D-QSAR and ADMET studies on pyrrolopyrimidine derivatives that act as inhibitors of the BCR site in cytochrome b. These studies aim to improve and identify new compounds that could serve as more potent potential BTK inhibitors, which would lead to the identification of new drug candidates in this field.In our study, we used 2D-QSAR (multiple linear regression, multiple nonlinear regression, and artificial neural networks), molecular docking, molecular dynamics, and ADMET properties to investigate the potential of 35 pyrrolopyrimidine derivatives as BTK inhibitors. A molecular docking study and molecular dynamics simulations of molecule 13 over 10 ns revealed that it establishes multiple hydrogen bonds with several residues and exhibits frequent stability throughout the simulation period. Based on the results obtained by molecular modeling, we proposed six new compounds (Pred1, Pred2, Pred3, Pred4, Pred5, and Pred6) with highly significant predicted activity by MLR models. A study based on the in silico evaluation of the predicted ADMET properties of the new candidate molecules is strongly recommended to classify these molecules as promising candidates for new anticancer agents specifically designed to target Bruton's tyrosine kinase (BTK) inhibition.

Probing the Proton-Loading Site of Cytochrome C Oxidase Using Time-Resolved Fourier Transform Infrared Spectroscopy.

Crystal structure analyses at atomic resolution and FTIR spectroscopic studies of cytochrome c oxidase have yet not revealed protonation or deprotonation of key sites of proton transfer in a time-resolved mode. Here, a sensitive technique to detect protolytic transitions is employed. In this work, probing a proton-loading site of cytochrome c oxidase from Paracoccus denitrificans with time-resolved Fourier transform infrared spectroscopy is presented for the first time. For this purpose, variants with single-site mutations of N131V, D124N, and E278Q, the key residues in the D-channel, were studied. The reaction of mutated CcO enzymes with oxygen was monitored and analyzed. Seven infrared bands in the “fast” kinetic spectra were found based on the following three requirements: (1) they are present in the “fast” phases of N131V and D124N mutants, (2) they have reciprocal counterparts in the “slow” kinetic spectra in these mutants, and (3) they are absent in “fast” kinetic spectra of the E278Q mutant. Moreover, the double-difference spectra between the first two mutants and E278Q revealed more IR bands that may belong to the proton-loading site protolytic transitions. From these results, it is assumed that several polar residues and/or water molecule cluster(s) share a proton as a proton-loading site. This site can be propionate itself (holding only a fraction of H+), His403, and/or water cluster(s).

Arabidopsis thaliana Hcc1 is a Sco-like metallochaperone for CuA assembly in Cytochrome c Oxidase.

The assembly of the CuA site in Cytochrome c Oxidase (COX) is a critical step for aerobic respiration in COX-dependent organisms. Several gene products have been associated with the assembly of this copper site, the most conserved of them belonging to the Sco family of proteins, which have been shown to perform different roles in different organisms. Plants express two orthologs of Sco proteins: Hcc1 and Hcc2. Hcc1 is known to be essential for plant development and for COX maturation, but its precise function has not been addressed until now. Here, we report the biochemical, structural and functional characterization of Arabidopsis thaliana Hcc1 protein (here renamed Sco1). We solved the crystal structure of the Cu+1 -bound soluble domain of this protein, revealing a tri coordinated environment involving a CxxxCxn H motif. We show that AtSco1 is able to work as a copper metallochaperone, inserting two Cu+1 ions into the CuA site in a model of CoxII. We also show that AtSco1 does not act as a thiol-disulfide oxido-reductase. Overall, this information sheds new light on the biochemistry of Sco proteins, highlighting the diversity of functions among them despite their high structural similarities. DATABASE: PDB entry 6N5U (Crystal structure of Arabidopsis thaliana ScoI with copper bound).

A Protonated Water Cluster as a Transient Proton‐Loading Site in Cytochrome c Oxidase

AbstractCytochrome c oxidase (CcO) is a redox‐driven proton pump that powers aerobic respiratory chains. We show here by multi‐scale molecular simulations that a protonated water cluster near the active site is likely to serve as the transient proton‐loading site (PLS) that stores a proton during the pumping process. The pKa of this water cluster is sensitive to the redox states of the enzyme, showing distinct similarities to other energy converting proton pumps.

A Protonated Water Cluster as a Transient Proton-Loading Site in Cytochrome c Oxidase.

Cytochrome c oxidase (CcO) is a redox-driven proton pump that powers aerobic respiratory chains. We show here by multi-scale molecular simulations that a protonated water cluster near the active site is likely to serve as the transient proton-loading site (PLS) that stores a proton during the pumping process. The pKa of this water cluster is sensitive to the redox states of the enzyme, showing distinct similarities to other energy converting proton pumps.

Distinct properties of semiquinone species detected at the ubiquinol oxidation Qo site of cytochrome bc1 and their mechanistic implications.

The two-electron ubiquinol oxidation or ubiquinone reduction typically involves semiquinone (SQ) intermediates. Natural engineering of ubiquinone binding sites of bioenergetic enzymes secures that SQ is sufficiently stabilized, so that it does not leave the site to membranous environment before full oxidation/reduction is completed. The ubiquinol oxidation Qo site of cytochrome bc1 (mitochondrial complex III, cytochrome b6f in plants) has been considered an exception with catalytic reactions assumed to involve highly unstable SQ or not to involve any SQ intermediate. This view seemed consistent with long-standing difficulty in detecting any reaction intermediates at the Qo site. New perspective on this issue is now offered by recent, independent reports on detection of SQ in this site. Each of the described SQs seems to have different spectroscopic properties leaving space for various interpretations and mechanistic considerations. Here, we comparatively reflect on those properties and their consequences on the SQ stabilization, the involvement of SQ in catalytic reactions, including proton transfers, and the reactivity of SQ with oxygen associated with superoxide generation activity of the Qo site.

Enhancing Mn(II)-Binding and Manganese Peroxidase Activity in a Designed Cytochrome c Peroxidase through Fine-Tuning Secondary-Sphere Interactions.

Noncovalent second-shell interactions are important in controlling metal-binding affinity and activity in metalloenzymes, but fine-tuning these interactions in designed metalloenzymes has not been fully explored. As a result, most designed metalloenzymes have low metal-binding affinity and activity. Here we identified three mutations in the second coordination shell of an engineered Mn(II)-binding site in cytochrome c peroxidase (called MnCcP.1, containing Glu45, Glu37, and Glu181 ligands) that mimics the native manganese peroxidase (MnP), and explored their effects on both Mn(II)-binding affinity and MnP activity. First, removing a hydrogen bond to Glu45 through Tyr36Phe mutation enhanced Mn(II)-binding affinity, as evidenced by a 2.8-fold decrease in the KM of Mn(II) oxidation. Second, introducing a salt bridge through Lys179Arg mutation improved Glu35 and Glu181 coordination to Mn(II), decreasing KM 2.6-fold. Third, eliminating a steric clash that prevented Glu37 from orienting toward Mn(II) resulted in an 8.6-fold increase in kcat/KM, arising primarily from a 3.6-fold decrease in KM, with a KM value comparable to that of the native enzyme (0.28 mM vs 0.19 mM for Pleurotus eryngii MnP PS3). We further demonstrated that while the effects of Tyr36Phe and Lys179Arg mutations are additive, because involved in secondary-shell interactions to different ligands, other combinations of mutations were antagonistic because they act on different aspects of the Mn(II) coordination at the same residues. Finally, we showed that these MnCcP variants are functional models of MnP that mimic its activity in both Mn(II) oxidation and degradation of a phenolic lignin model compound and kraft lignin. In addition to achieving KM in a designed protein that is similar to the that of native enzyme, our results offer molecular insight into the role of noncovalent interactions around metal-binding sites for improving metal binding and overall activity; such insight can be applied to rationally enhance these properties in other metalloenzymes and their models.

Mitochondrial Disease-related Mutation G167P in Cytochrome b of Rhodobacter capsulatus Cytochrome bc1 (S151P in Human) Affects the Equilibrium Distribution of [2Fe-2S] Cluster and Generation of Superoxide

Cytochrome bc1 is one of the key enzymes of many bioenergetic systems. Its operation involves a large scale movement of a head domain of iron-sulfur protein (ISP-HD), which functionally connects the catalytic quinol oxidation Qo site in cytochrome b with cytochrome c1. The Qo site under certain conditions can generate reactive oxygen species in the reaction scheme depending on the actual position of ISP-HD in respect to the Qo site. Here, using a bacterial system, we show that mutation G167P in cytochrome b shifts the equilibrium distribution of ISP-HD toward positions remote from the Qo site. This renders cytochrome bc1 non-functional in vivo. This effect is remediated by addition of alanine insertions (1Ala and 2Ala) in the neck region of the ISP subunit. These insertions, which on their own shift the equilibrium distribution of ISP-HD in the opposite direction (i.e. toward the Qo site), also act in this manner in the presence of G167P. Changes in the equilibrium distribution of ISP-HD in G167P lead to an increased propensity of cytochrome bc1 to generate superoxide, which becomes evident when the concentration of quinone increases. This result corroborates the recently proposed model in which "semireverse" electron transfer back to the Qo site, occurring when ISP-HD is remote from the site, favors reactive oxygen species production. G167P suggests possible molecular effects of S151P (corresponding in sequence to G167P) identified as a mitochondrial disease-related mutation in human cytochrome b. These effects may be valid for other human mutations that change the equilibrium distribution of ISP-HD in a manner similar to G167P.

Reversible Switching of Redox-Active Molecular Orbitals and Electron Transfer Pathways in Cu(A) Sites of Cytochrome c Oxidase.

The Cu(A) site of cytochrome c oxidase is a redox hub that participates in rapid electron transfer at low driving forces with two redox cofactors in nearly perpendicular orientations. Spectroscopic and electrochemical characterizations performed on first and second-sphere mutants have allowed us to experimentally detect the reversible switching between two alternative electronic states that confer different directionalities to the redox reaction. Specifically, the M160H variant of a native Cu(A) shows a reversible pH transition that allows to functionally probe both states in the same protein species. Alternation between states exerts a dramatic impact on the kinetic redox parameters, thereby suggesting this effect as the mechanism underlying the efficiency and directionality of Cu(A) electron transfer in vivo. These findings may also prove useful for the development of molecular electronics.

Protonation of the binuclear active site in cytochrome c oxidase decreases the reduction potential of CuB

One of the remaining mysteries regarding the respiratory enzyme cytochrome c oxidase is how proton pumping can occur in all reduction steps in spite of the low reduction potentials observed in equilibrium titration experiments for two of the active site cofactors, CuB(II) and Fea3(III). It has been speculated that, at least the copper cofactor can acquire two different states, one metastable activated state occurring during enzyme turnover, and one relaxed state with lower energy, reached only when the supply of electrons stops. The activated state should have a transiently increased CuB(II) reduction potential, allowing proton pumping. The relaxed state should have a lower reduction potential, as measured in the titration experiments. However, the structures of these two states are not known. Quantum mechanical calculations show that the proton coupled reduction potential for CuB is inherently high in the active site as it appears after reaction with oxygen, which explains the observed proton pumping. It is suggested here that, when the flow of electrons ceases, a relaxed resting state is formed by the uptake of one extra proton, on top of the charge compensating protons delivered in each reduction step. The extra proton in the active site decreases the proton coupled reduction potential for CuB by almost half a volt, leading to agreement with titration experiments. Furthermore, the structure for the resting state with an extra proton is found to have a hydroxo-bridge between CuB(II) and Fea3(III), yielding a magnetic coupling that can explain the experimentally observed EPR silence.

CO-dynamics in the active site of cytochrome c oxidase.

The transfer of CO from heme a3 to the Cu(B) site in Cytochrome c oxidase (CcO) after photolysis is studied using molecular dynamics simulations using an explicitly reactive, parametrized potential energy surface based on density functional theory calculations. After photodissociation from the heme-Fe, the CO ligand rebinds to the Cu(B) site on the sub-picosecond time scale. Depending on the simulation protocol the characteristic time ranges from 260 fs to 380 fs which compares with an estimated 450 fs from experiment based on the analysis of the spectral changes as a function of time delay after the photodissociating pulse. Following photoexcitation ≈90% of the ligands are found to rebind to either the Cu(B) (major component, 85%) or the heme-Fe (minor component, 2%) whereas about 10% remain in an unbound state. The infrared spectra of unbound CO in the active site is broad and featureless and no appreciable shift relative to gas-phase CO is found, which is in contrast to the situation in myoglobin. These observations explain why experimentally, unbound CO in the binuclear site of CcO has not been found as yet.

Predicting Charged-Ligand Binding from Molecular Simulations

We seek to use molecular simulations to predict and understand protein-ligand interactions. Our approach uses so-called alchemical free energy calculations to quantitatively predict binding affinities. Here, we will discuss recent work on blind predictions of protein-ligand binding in a model binding site in cytochrome C peroxidase, as well as related work in other binding sites. We discuss the outcome of our predictions, what we learned, and the challenges facing computational methods in predicting biomolecular interactions. We also discuss ongoing related work in other binding sites, and implications of our work for binding studies more broadly.

Computational Prediction and in Vitro Analysis of Potential Physiological Ligands of the Bile Acid Binding Site in Cytochrome c Oxidase

A conserved bile acid site has been crystallographically defined in the membrane domain of mammalian and Rhodobacter sphaeroides cytochrome c oxidase (RsCcO). Diverse amphipathic ligands were shown previously to bind to this site and affect the electron transfer equilibrium between heme a and a3 cofactors by blocking the K proton uptake path. Current studies identify physiologically relevant ligands for the bile acid site using a novel three-pronged computational approach: ROCS comparison of ligand shape and electrostatics, SimSite3D comparison of ligand binding site features, and SLIDE screening of potential ligands by docking. Identified candidate ligands include steroids, nicotinamides, flavins, nucleotides, retinoic acid, and thyroid hormones, which are predicted to make key protein contacts with the residues involved in bile acid binding. In vitro oxygen consumption and ligand competition assays on RsCcO wildtype and its Glu101Ala mutant support regulatory activity and specificity of some of these ligands. An ATP analog and GDP inhibit RsCcO under low substrate conditions, while fusidic acid, cholesteryl hemisuccinate, retinoic acid, and T3 thyroid hormone are more potent inhibitors under both high and low substrate conditions. The sigmoidal kinetics of RsCcO inhibition in the presence of certain nucleotides is reminiscent of previously reported ATP inhibition of mammalian CcO, suggesting regulation involving the conserved core subunits of both mammalian and bacterial oxidases. Ligand binding to the bile acid site is noncompetitive with respect to cytochrome c and appears to arrest CcO in a semioxidized state with some resemblance to the "resting" state of the enzyme.

Blind Prediction of Charged Ligand Binding Affinities in a Model Binding Site

Predicting absolute protein–ligand binding affinities remains a frontier challenge in ligand discovery and design. This becomes more difficult when ionic interactions are involved because of the large opposing solvation and electrostatic attraction energies. In a blind test, we examined whether alchemical free-energy calculations could predict binding affinities of 14 charged and 5 neutral compounds previously untested as ligands for a cavity binding site in cytochrome c peroxidase. In this simplified site, polar and cationic ligands compete with solvent to interact with a buried aspartate. Predictions were tested by calorimetry, spectroscopy, and crystallography. Of the 15 compounds predicted to bind, 13 were experimentally confirmed, while 4 compounds were false negative predictions. Predictions had a root-mean-square error of 1.95kcal/mol to the experimental affinities, and predicted poses had an average RMSD of 1.7Å to the crystallographic poses. This test serves as a benchmark for these thermodynamically rigorous calculations at predicting binding affinities for charged compounds and gives insights into the existing sources of error, which are primarily electrostatic interactions inside proteins. Our experiments also provide a useful set of ionic binding affinities in a simplified system for testing new affinity prediction methods.

Separated topologies—A method for relative binding free energy calculations using orientational restraints

Orientational restraints can improve the efficiency of alchemical free energy calculations, but they are not typically applied in relative binding calculations, which compute the affinity difference been two ligands. Here, we describe a new "separated topologies" method, which computes relative binding free energies using orientational restraints and which has several advantages over existing methods. While standard approaches maintain the initial and final ligand in a shared orientation, the separated topologies approach allows the initial and final ligands to have distinct orientations. This avoids a slowly converging reorientation step in the calculation. The separated topologies approach can also be applied to determine the relative free energies of multiple orientations of the same ligand. We illustrate the approach by calculating the relative binding free energies of two compounds to an engineered site in Cytochrome C Peroxidase.

An investigation into a cardiolipin acyl chain insertion site in cytochrome c

Mitochondrial cytochrome c associates with the phosphoplipid cardiolipin (CL) through a combination of electrostatic and hydrophobic interactions. The latter occurs by insertion into cytochrome c of an acyl chain, resulting in the dissociation of the axial Met-80 heme-iron ligand. The resulting five coordinate cytochrome c/CL complex has peroxidatic properties leading to peroxidation of CL and dissociation of the complex. These events are considered to be pre-apoptotic and culminate with release of cytochrome c from the mitochondria into the cytoplasm. Two distinct surface regions on cytochrome c have been suggested to mediate CL acyl chain insertion and this study has probed one of these regions. We have constructed a series of alanine mutants aimed at disrupting a surface cleft formed between residues 67–71 and 82–85. The physicochemical properties, peroxidase activity, CL binding, and kinetics of carbon monoxide (CO) binding to the ferrous cytochrome c/CL complex have been assessed for the individual mutants. Our findings reveal that the majority of mutants are capable of binding CL in the same apparent stoichiometry as the wild-type protein, with the extent to which the Met-80 ligand is bound in the ferrous cytochrome c/CL complex being mutant specific at neutral pH. Mutation of the species conserved Arg-91 residue, that anchors the cleft, results in the greatest changes to physicochemical properties of the protein leading to a change in the CL binding ratio required to effect structural changes and to the ligand-exchange properties of the ferrous cytochrome c/CL complex.

Electronic structures of the Cu2S2 core of the CuA site in cytochrome c oxidase and nitrous oxide reductase

AbstractThe CuA site functions as an electron transfer intermediate in cytochrome c oxidase (CcO) and nitrous oxide reductase. X‐ray crystal structures have shown that the CuA site is a binuclear copper site, bridged by two cysteinyl thiolate groups. Each copper ion is coordinated by a histidine side chain and a weaker axial ligand: a methionine side chain or a backbone carbonyl group. The ground state of the oxidized CuA site is strongly related to the singly occupied orbital, as the oxidized state of the CuA site is doublet. Spectroscopic studies have suggested the σu* ground state of the CuA site due to the direct CuCu interaction, facilitating a rapid electron transfer over long distances with low driving forces. In this study, to elucidate what factors are responsible for the stabilization of the σu* ground state of the CuA site, the electronic structures of the Cu2S2 core of the CuA have been studied, using the density functional theory (M06). Our computational results show that the ground state of the Cu2S2 core is the πu state even with a shorter CuCu distance, and that an addition of electrostatic and orbital interactions by ligand coordination stabilizes the σu* state rather than the πu state, because electrostatic repulsion and antibonding orbital interaction between the dσ* orbital of the CuCu dimer and the lone pairs of the coordinating His residues rise the orbital energy of the σu* orbital of the CuA site. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

Thermodynamics of copper and zinc distribution in the cyanobacterium Synechocystis PCC 6803.

Copper is supplied to plastocyanin for photosynthesis and cytochrome c oxidase for respiration in the thylakoids of Synechocystis PCC 6803 by the membrane-bound P-type ATPases CtaA and PacS, and the metallochaperone Atx1. We have determined the Cu(I) affinities of all of the soluble proteins and domains in this pathway. The Cu(I) affinities of the trafficking proteins range from 5 × 10(16) to 5 × 10(17) M(-1) at pH 7.0, consistent with values for homologues. Unusually, Atx1 binds Cu(I) significantly tighter than the metal-binding domains (MBDs) of CtaA and PacS (CtaA(N) and PacS(N)), and equilibrium copper exchange constants of approximately 0.2 are obtained for transfer to the MBDs. Dimerization of Atx1 increases the affinity for Cu(I), but the loop 5 His61 residue has little influence. The MBD of the zinc exporter ZiaA (ZiaA(N)) exhibits an almost identical Cu(I) affinity, and Cu(I) exchange with Atx1, as CtaA(N) and PacS(N), and the relative stabilities of the complexes must enable the metallochaperone to distinguish between the MBDs. The binding of potentially competing zinc to the trafficking proteins has been studied. ZiaA(N) has the highest Zn(II) affinity and thermodynamics could be important for zinc removal from the cell. Plastocyanin has a Cu(I) affinity of 2.6 × 10(17) M(-1), 15-fold tighter than that of the Cu(A) site of cytochrome c oxidase, highlighting the need for specific mechanisms to ensure copper delivery to both of these targets. The narrow range of Cu(I) affinities for the cytoplasmic copper proteins in Synechocystis will facilitate relocation when copper is limiting.

A theoretical investigation of the functional role of the axial methionine ligand of the Cu A site in cytochrome c oxidase

The functional roles of the amino acid residues of the Cu A site in bovine cytochrome c oxidase (C cO) were investigated by utilizing hybrid quantum mechanics (QM)/molecular mechanics (MM) calculations. The energy levels of the molecular orbitals (MOs) involving Cu d zx orbitals unexpectedly increased, as compared with those found previously with a simplified model system lacking the axial Met residue (i.e., Cu 2S 2N 2). This elevation of MO energies stemmed from the formation of the anti-bonding orbitals, which are generated by hybridization between the d zx orbitals of Cu ions and the p-orbitals of the S and O atoms of the axial ligands. To clarify the roles of the axial Met ligand, the inner-sphere reorganization energies of the Cu A site were computed, with the Met residue assigned to either the QM or MM region. The reorganization energy slightly increased when the Met residue was excluded from the QM region. The existing experimental data and the present structural modeling study also suggested that the axial Met residue moderately increased the redox potential of the Cu A site. Thus, the role of the Met may be to regulate the electron transfer rate through the fine modulation of the electronic structure of the Cu A “platform”, created by two Cys/His residues coordinated to the Cu ions. This regulation would provide the optimum redox potential/reorganization energy of the Cu A site, and thereby facilitate the subsequent cooperative reactions, such as the proton pump and the enzymatic activity, of C cO. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

2I1512 A computational study of functional role of axial methionine ligand of Cu_A site in cytochrome c oxidase(Heme proteins 2,The 48th Annual Meeting of the Biophysical Society of Japan)

2I1512 A computational study of functional role of axial methionine ligand of Cu_A site in cytochrome c oxidase(Heme proteins 2,The 48th Annual Meeting of the Biophysical Society of Japan)