Assessing protein–ligand binding modes with computational tools: the case of PDE4B

In a first step in the discovery of novel potent inhibitor structures for the PDE4B family with limited side effects, we present a protocol to rank newly designed molecules through the estimation of their IC\(_{50}\) values. Our protocol is based on reproducing the linear relationship between the logarithm of experimental IC\(_{50}\) values [\(\log\)(IC\(_{50}\))] and their calculated binding free energies (\(\Delta G_\mathrm{binding}\)). From 13 known PDE4B inhibitors, we show here that (1) binding free energies obtained after a docking process by AutoDock are not accurate enough to reproduce this linear relationship; (2) MM-GB/SA post-processing of molecular dynamics (MD) trajectories of the top ranked AutoDock pose improves the linear relationship; (3) by taking into account all representative structures obtained by AutoDock and by averaging MM-GB/SA computations on a series of 40 independent MD trajectories, a linear relationship between \(\log\)(IC\(_{50}\)) and the lowest \(\Delta G_\mathrm{binding}\) is achieved with \(R^2=0.944\).     

Uncertainty from sampling: workshop to launch a Nordtest handbook on sampling uncertainty estimation and control

A workshop on uncertainty in sampling was held in Hillerød, Denmark, on 12–13 April 2007 to launch a new handbook on sampling quality assurance and uncertainty estimation. The participants of the workshop were approximately 60 delegates from 15 European countries, representing institutions performing sampling, users of the data, research institutions, as well as accreditation bodies. Materials from the workshop, including examples, tools, and calculation aids for the work can be found at http://www.samplersguide.com. The Nordtest handbook Uncertainty from sampling will be made available on the Nordtest web site at http://www.nordicinnovation.net/nordtest.cfm under NT technical reports, report number NT tec 604. Until the final report is available on the Nordtest web site, an advance draft of the Nordtest handbook is available from http://www.samplersguide.com.     

CDOCKER and $$\lambda$$-dynamics for prospective prediction in D3R Grand Challenge 2

The opportunity to prospectively predict ligand bound poses and free energies of binding to the Farnesoid X Receptor in the D3R Grand Challenge 2 provided a useful exercise to evaluate CHARMM based docking (CDOCKER) and \(\lambda\)-dynamics methodologies for use in “real-world” applications in computer aided drug design. In addition to measuring their current performance, several recent methodological developments have been analyzed retrospectively to highlight best procedural practices in future applications. For pose prediction with CDOCKER, when the protein structure used for rigid receptor docking was close to the crystallographic holo structure, reliable poses were obtained. Benzimidazoles, with a known holo receptor structure, were successfully docked with an average RMSD of 0.97 \(\AA\). Other non-benzimidazole ligands displayed less accuracy largely because the receptor structures we chose for docking were too different from the experimental holo structures. However, retrospective analysis has shown that when these ligands were re-docked into their holo structures, the average RMSD dropped to 1.18 \(\AA\) for all ligands. When sulfonamides and spiros were docked with the apo structure, which agrees more with their holo structure than the structures we chose, five out of six ligands were correctly docked. These docking results emphasize the need for flexible receptor docking approaches. For \(\lambda\)-dynamics techniques, including multisite \(\lambda\)-dynamics (MS\(\lambda\)D), reasonable agreement with experiment was observed for the 33 ligands investigated; root mean square errors of 2.08 and 1.67 kcal/mol were obtained for free energy sets 1 and 2, respectively. Retrospectively, soft-core potentials, adaptive landscape flattening, and biasing potential replica exchange (BP-REX) algorithms were critical to model large substituent perturbations with sufficient precision and within restrictive timeframes, such as was required with participation in Grand Challenge 2. These developments, their associated benefits, and proposed procedures for their use in future applications are discussed.     

Identification of Microorganisms by High Resolution Tandem Mass Spectrometry with Accurate Statistical Significance

Correct and rapid identification of microorganisms is the key to the success of many important applications in health and safety, including, but not limited to, infection treatment, food safety, and biodefense. With the advance of mass spectrometry (MS) technology, the speed of identification can be greatly improved. However, the increasing number of microbes sequenced is challenging correct microbial identification because of the large number of choices present. To properly disentangle candidate microbes, one needs to go beyond apparent morphology or simple ‘fingerprinting’; to correctly prioritize the candidate microbes, one needs to have accurate statistical significance in microbial identification. We meet these challenges by using peptidome profiles of microbes to better separate them and by designing an analysis method that yields accurate statistical significance. Here, we present an analysis pipeline that uses tandem MS (MS/MS) spectra for microbial identification or classification. We have demonstrated, using MS/MS data of 81 samples, each composed of a single known microorganism, that the proposed pipeline can correctly identify microorganisms at least at the genus and species levels. We have also shown that the proposed pipeline computes accurate statistical significances, i.e., E-values for identified peptides and unified E-values for identified microorganisms. The proposed analysis pipeline has been implemented in MiCId, a freely available software for Microorganism Classification and Identification. MiCId is available for download at http://www.ncbi.nlm.nih.gov/CBBresearch/Yu/downloads.html.

Open image in new window

Graphical Abstract ?

     

In vitro and in silico inhibition of angiotensin-converting enzyme by carbohydrates and cyclitols

Fifteen carbohydrates (d-mannose, d-glucose, d-galactose, methyl-α-d-glucose, l-rhamnose, d-xylose, d-fructose, d-arabinose, dulcitol, mannitol, β-maltose, α-lactose, melibiose, sucrose, and raffinose) and four cyclitols [l-(+)-bornesitol, myo-inositol, per-O-acetyl-1-l-(+)-bornesitol, and quinic acid] were assayed for in vitro ACE inhibition. Of these molecules, per-O-Acetyl-1-l-(+)-bornesitol, quinic acid, methyl-α-d-glucose, d-rhamnose, raffinose, and the disaccharides were determined to be either inactive or weak ACE inhibitors, whereas l-(+)-bornesitol, d-galactose, d-glucose, and myo-inositol exhibited significant ACE inhibition. Molecular docking studies were performed to investigate interactions between active compounds and human ACE (Protein Data Bank, PDB 1O83). The results of various calculations showed that all active sugars bind to the same enzyme region, which is a tunnel directed towards the active site. With the exception of myo-inositol (K _i = 13.95 μM, IC₅₀ = 449.2 μM), the active compounds presented similar K _i and IC₅₀ values. d-Galactose (K _i = 19.6 μM, IC₅₀ = 35.7 μM) and l-(+)-bornesitol (K _i = 25.3 μM, IC₅₀ = 41.4 μM) were the most active compounds, followed by d-glucose (K _i = 32.9 μM, IC₅₀ = 85.7 μM). Our docking calculations are in agreement with the experimental data and show a new binding region for sugar-like molecules, which may be explored for the development of new ACE inhibitors.     

A novel view of modelling interactions between synthetic and biological polymers via docking

Multipoint interactions between synthetic and natural polymers provide a promising platform for many topical applications, including therapeutic blockage of virus-specific targets. Docking may become a useful tool for modelling of such interactions. However, the rigid docking cannot be correctly applied to synthetic polymers with flexible chains. The application of flexible docking to these polymers as whole macromolecule ligands is also limited by too many possible conformations. We propose to solve this problem via stepwise flexible docking. Step 1 is docking of separate polymer components: (1) backbone units (BU), multi-repeated along the chain, and (2) side groups (SG) consisting of functionally active elements (SG _F ) and bridges (SG _B ) linking SG _F with BU. At this step, probable binding sites locations and binding energies for the components are scored. Step 2 is docking of component-integrating models: [BU] _m , SG = SG _F –SG _B , BU–SG, BU–BU(SG)–BU, BU(SG)–[BU] _m –BU(SG), and [BU _var (SG _var )] _m . Every modelling level yields new information, including how the linkage of various components influences on the ligand—target contacts positioning, orientation, and binding energy in step-by-step approximation to polymeric ligand motifs. Step 3 extrapolates the docking results to real-scale macromolecules. This approach has been demonstrated by studying the interactions between hetero-SG modified anionic polymers and the N-heptad repeat region tri-helix core of the human immunodeficiency virus type 1 (HIV-1) envelope glycoprotein gp41, the key mediator of HIV-1 fusion during virus entry. The docking results are compared to real polymeric compounds, acting as HIV-1 entry inhibitors in vitro. This study clarifies the optimal macromolecular design for the viral fusion inhibition and drug resistance prevention.     

GalaxyDock BP2 score: a hybrid scoring function for accurate protein–ligand docking

Protein–ligand docking is a useful tool for providing atomic-level understanding of protein functions in nature and design principles for artificial ligands or proteins with desired properties. The ability to identify the true binding pose of a ligand to a target protein among numerous possible candidate poses is an essential requirement for successful protein–ligand docking. Many previously developed docking scoring functions were trained to reproduce experimental binding affinities and were also used for scoring binding poses. However, in this study, we developed a new docking scoring function, called GalaxyDock BP2 Score, by directly training the scoring power of binding poses. This function is a hybrid of physics-based, empirical, and knowledge-based score terms that are balanced to strengthen the advantages of each component. The performance of the new scoring function exhibits significant improvement over existing scoring functions in decoy pose discrimination tests. In addition, when the score is used with the GalaxyDock2 protein–ligand docking program, it outperformed other state-of-the-art docking programs in docking tests on the Astex diverse set, the Cross2009 benchmark set, and the Astex non-native set. GalaxyDock BP2 Score and GalaxyDock2 with this score are freely available at http://galaxy.seoklab.org/softwares/galaxydock.html.     

Azahar: a PyMOL plugin for construction,visualization and analysis of glycan molecules

Glycans are key molecules in many physiological and pathological processes. As with other molecules, like proteins, visualization of the 3D structures of glycans adds valuable information for understanding their biological function. Hence, here we introduce Azahar, a computing environment for the creation, visualization and analysis of glycan molecules. Azahar is implemented in Python and works as a plugin for the well known PyMOL package (Schrodinger in The PyMOL molecular graphics system, version 1.3r1, 2010). Besides the already available visualization and analysis options provided by PyMOL, Azahar includes 3 cartoon-like representations and tools for 3D structure caracterization such as a comformational search using a Monte Carlo with minimization routine and also tools to analyse single glycans or trajectories/ensembles including the calculation of radius of gyration, Ramachandran plots and hydrogen bonds. Azahar is freely available to download from http://www.pymolwiki.org/index.php/Azahar and the source code is available at https://github.com/agustinaarroyuelo/Azahar.     

The UPS in diabetes and obesity

Abstract

Type 2 diabetes is caused by defects in both insulin signaling and insulin secretion. Though the role of the ubiquitin proteasome system (UPS) in the pathogenesis of type 2 diabetes remains largely unexplored, the few examples present in the literature are interesting and suggest targets for drug development. Studies indicate that insulin resistance can be induced by stimulating the degradation of important molecules in the insulin signaling pathway, in particular the insulin receptor substrate proteins IRS1, IRS2 and the kinase AKT1 (Akt). In addition, a defect in insulin secretion could occur due to UPS-mediated degradation of IRS2 in the β-cells of the pancreas. The UPS also appears to be involved in regulating lipid synthesis in adipocytes and lipid production by the liver and could influence the development of obesity. Other possible mechanisms for inducing defects in insulin signaling and secretion remain to be explored, including the role of ubiquitylation in insulin receptor internalization and trafficking.

Publication history

Republished from Current BioData's Targeted Proteins database (TPdb; http://www.targetedproteinsdb.com).

     

Theoretical mapping of new L–(N)–L family of species with donor–acceptor bonding between N and ligand L

We perform a computational mapping study of a family of new inorganic species, based on idea of donor–acceptor type bonding between N⁺ and a ligand L with a terminal electron lone pair. The nitrogen ion is seen as being in an atomic ¹D state, with empty 2p acceptor orbitals [I.S.K. Kerkines., A. Papakondylis, A. Mavridis, J. Phys. Chem. A, 2002, 106, 4435]. We consider a series of small ligands, such as PN, CCH⁻, CCCN⁻, , and others. Chemical bonding analysis confirms the suggested bonding picture as characteristic for experimentally known and as well as for most of the predicted species. A number of these new compounds is found to be thermodynamically stable with respect to the existing or . They are candidates for new synthetic targets.     

Biaryl-based anion receptors bearing thiourea groups: fluoride anion receptor

The synthesis, characterization and binding studies with anions for biaryl-based anion receptors bearing thiourea groups have been described. The results revealed that receptors (1 and 2) showed good selectivity and binding affinity for F^?, and among them binaphthyl-based receptor (1a) showed the best binding affinity for F^? in comparison to other tested anions (Cl^?, Br^?, I^?, $ {\text{NO}}_{3}^{ - } ,\;{\text{HSO}}_{4}^{ - } , $ AcO^? and $ {\text{H}}_{2} {\text{PO}}_{4}^{ - } $ ). This is probably due to the fact that the moderate rigidity of binaphthyl skeleton in 1a is able to provide the better geometry of two thiourea groups for incorporating F^? into the binding pocket. The higher basicity of F^? also participated in this selectivity.     

Synthesis of Some N1-(3-Hydroxy-2-pyridyl)benzamides

Benzoxazoles are of special interest because of having various types of biological properties such as antihistaminic, antihelmintic antifungal and antibacterial activities. ¹ Hisano, T., Ichikawa, M., Tsumoto, K. and Tasaki, M. 1982. Chem. Pharm. Bull, 30: 2996 [CSA][Crossref], [Web of Science ®] , [Google Scholar], ² Prudhomme, M., Guyot, J. and Jeminet, G. 1986. J. Antibiotics, 39: 934 [CSA][Crossref], [PubMed], [Web of Science ®] , [Google Scholar], ³ Ersan, S., Nacak, S., Berkem, R. and Özden, T. 1997. Arzneim. Forsch, 47: 963 [CSA][PubMed], [Web of Science ®] , [Google Scholar], ⁴ Sener, E., Yalc?n, ?., Temeiz, Ö., Ören, ?., Ak?n, A. and Ucartürk, N. 1997. Farmaco, 52: 99 [CSA][PubMed] , [Google Scholar], ⁵ Ören, ?., Temiz, Ö., Yalc?n, ?., Sener, E., Ak?n and Ucartürk, N. 1997. Forsch, 47: 1393 [CSA] [Google Scholar], ⁶ Temiz, Ö., Ören, ?., Sener, E., Yalc?n, ?. and Ucartürk, N. 1998. Farmaco, 53: 337 [CSA][CROSSREF][Crossref], [PubMed] , [Google Scholar], ⁷ Yalc?n, ?., Ören, ?., Sener, E., Ak?n, A. and Ucartürk, N. 1992. J. Med. Chem, 27: 401 [CSA][CROSSREF][Crossref], [Web of Science ®] , [Google Scholar], ⁸ Temiz, Ö., Ören, ?. and Altanlar, N. 2002. Farmaco, 57: 175 [CSA][CROSSREF][Crossref], [PubMed] , [Google Scholar], ⁹ Sener, E., Temeiz, Ö., Yalc?n, ?. and Altanlar, N. 2000. Farmaco, 55: 397 [CSA][CROSSREF][Crossref], [PubMed] , [Google Scholar], ¹⁰ Ören, ?., Temiz, Ö., Yalc?n, ?., Sener, E. and Altanlar, N. 1998. Eur. J. Pharm. Sci., 7: 153 [CSA][Crossref], [Web of Science ®] , [Google Scholar] Benzamide derivatives, as the possible metabolites of benzoxazoles, show various types of biological activities. ¹¹ Sener, E., Bingöl, K., Ören, ?., Arpac?, Ö., Yalc?n, K. ?. and Altanlar, N. 2000. Farmaco, 55: 469 [CSA][CROSSREF][Crossref], [PubMed] , [Google Scholar], ¹² Sener, E., Yalc?n, K. ?. and Altanlar, N. 2002. Farmaco, 57: 451 [CSA][CROSSREF][Crossref], [PubMed] , [Google Scholar] Some N-(2-hydroxyphenyl)- benzamides, which showed significant activity compared to phenylacetamides and furamides,11 have been synthesized by treating 2-aminophenol with carboxylic acid chlorides under weak basic solution. ¹¹ Sener, E., Bingöl, K., Ören, ?., Arpac?, Ö., Yalc?n, K. ?. and Altanlar, N. 2000. Farmaco, 55: 469 [CSA][CROSSREF][Crossref], [PubMed] , [Google Scholar], ¹² Sener, E., Yalc?n, K. ?. and Altanlar, N. 2002. Farmaco, 57: 451 [CSA][CROSSREF][Crossref], [PubMed] , [Google Scholar] However, synthesis and especially biological activity of N-(3-hydroxy-2-pyridyl)benzamides were not studied well. Owing to the versatility of benzamides we have extended the reaction of 2-amino-3-pyridinol with different carboxylic acid chlorides in order to preparation of some new benzamide derivatives containing a hydroxypyridyl ring.     

Calcium-induced tetramerization and zinc chelation shield human calprotectin from degradation by host and bacterial extracellular proteases

Calprotectin (CP, S100A8/S100A9 oligomer, MRP-8/14 oligomer, calgranulins A and B) is a protein component of the innate immune system that contributes to the metal-withholding response by sequestering bioavailable transition metal ions at sites of infection. Human CP employs Ca(ii) ions to modulate its quaternary structure, transition metal binding properties, and antimicrobial activity. In this work, we report the discovery that Ca(ii)-induced self-association of human CP to afford heterotetramers protects the protein scaffold from degradation by host serine proteases. We present the design and characterization of two new human CP-Ser variants, S100A8(C42S)(I60E)/S100A9(C3S) and S100A8(C42S)(I60K)/S100A9(C3S), that exhibit defective tetramerization properties. Analytical size exclusion chromatography and analytical ultracentrifugation reveal that both variants, hereafter I60E and I60K, persist as heterodimers in the presence of Ca(ii) only, and form heterotetramers in the presence of Mn(ii) only and both Ca(ii) and Mn(ii). Coordination to Fe(ii) also causes I60E and I60K to form heterotetramers in both the absence and presence of Ca(ii). The Ca(ii)-bound I60E and I60K heterodimers are readily degraded by trypsin, chymotrypsin, and human neutrophil elastase, whereas the Ca(ii)-bound CP-Ser heterotetramers and the Ca(ii)- and Mn(ii)-bound I60E and I60K heterotetramers are resistant to degradation by these host proteases. The staphylococcal extracellular protease GluC cuts the S100A8 subunit of CP-Ser at the C-terminal end of residue 89 to afford a ΔSHKE variant. The GluC cleavage site is in close proximity to the His₃Asp metal-binding site, which coordinates Zn(ii) with high affinity, and Zn(ii) chelation protects the S100A8 subunit from GluC cleavage. Taken together, these results provide new insight into how Ca(ii) ions and transition metals modulate the chemistry and biology of CP, and indicate that coordination to divalent cations transforms human CP into a protease-resistant form and enables innate immune function in the hostile conditions of an infection site.     

New four-stages symmetric six-step method with improved phase properties for second order problems with periodical and/or oscillating solutions

In this paper, we build, for the first time in the literature, a new four-stages symmetric six-step finite difference pair with optimized properties. The method:

1. 1.
   is a symmetric non-linear six-step method,
    
2. 2.
   is of four stages
    
3. 3.
   is of fourteenth algebraic order,
    
4. 4.
   has eliminated the phase-lag,
    
5. 5.
   has eliminated the first and second derivatives of the phase-lag.
    

An analysis of the new proposed method is given in details in this paper. More specifically, we present:

1. 1.
   the building of the new four-stages symmetric six-step method,
    
2. 2.
   the computation of the local truncation error of the new proposed method,
    
3. 3.
   the comparative local truncation error analysis of the new proposed method with other finite difference pairs of the same family.
    
4. 4.
   the stability and the interval of periodicity analysis and
    
5. 5.
   finally, the investigation and evaluation of the computational efficiency of the new proposed scheme for the approximate solution of the Schrödinger equation.
    

The theoretical, computational and numerical results for the new proposed method show its effectiveness compared with other known or recently obtained finite difference pairs in the literature.

     

Role of the UPS in Liddle syndrome

Abstract

Hypertension is a serious medical problem affecting a large population worldwide. Liddle syndrome is a hereditary form of early onset hypertension caused by mutations in the epithelial Na⁺ channel (ENaC). The mutated region, called the PY (Pro-Pro-x-Tyr) motif, serves as a binding site for Nedd4-2, an E3 ubiquitin ligase from the HECT family. Nedd4-2 binds the ENaC PY motif via its WW domains, normally leading to ENaC ubiquitylation and endocytosis, reducing the number of active channels at the plasma membrane. In Liddle syndrome, this endocytosis is impaired due to the inability of the mutated PY motif in ENaC to properly bind Nedd4-2. This leads to accumulation of active channels at the cell surface and increased Na⁺ (and fluid) absorption in the distal nephron, resulting in elevated blood volume and blood pressure. Small molecules/compounds that destabilize cell surface ENaC, or enhance Nedd4-2 activity in the kidney, could potentially serve to alleviate hypertension.

Publication history

Republished from Current BioData's Targeted Proteins database (TPdb; http://www.targetedproteinsdb.com).

     

Assessment of the ��B97 family for excited-state calculations

We benchmark three recently proposed range-separated hybrids, namely ??B97, ??B97X and ??B97XD in the framework of time-dependent density functional theory simulations of electronic absorption spectra. Comparisons are made with both theoretical estimates obtained by highly correlated approaches and experimental wavelengths of maximal absorption measured for important classes of $\pi \rightarrow \pi^\ast$ and $n \rightarrow \pi^\ast$ chromogens. The amplitude of the errors induced by the lack of vibronic coupling in our computational model is also evaluated for five dyes. The performances of the ??B97 group are systematically compared to the results of other global and range-separated hybrids. It turns out that ??B97XD provides, in general, more accurate estimates than ??B97X and ??B97.     

A local-MP2 approach to the ab initio study of electron correlation in crystals and to the simulation of vibrational spectra: the case of Ice XI

Wave-function-based computational techniques for describing electron correlation effects in periodic systems have a long history. Among early attempts, the method proposed by Colle and Salvetti (Theor Chim Acta 37:329, 1975) more than 30 years ago is noteworthy for its simplicity, power, and far-reaching consequences. The renewed interest in this topic is due to the well-known failure of techniques based on Density Functional Theory when it comes to obtaining very accurate estimates of some important quantities. Here we present the essential features of an ab initio code, ( CRYSCOR ) recently implemented (Phys Rev B 76: 075101, 075102, 2007), which solves the MP2 equations for crystals by adopting a local-correlation approach and using as a reference the Hartree-Fock solution provided by the CRYSTAL program. As an example of application, we discuss the MP2 corrections to the frequency of some vibrational modes in a proton-ordered structure of water ice (Ice XI).     

Targeting the UPS as therapy in multiple myeloma

Abstract

The coordinated regulation of cellular protein synthesis and degradation is essential for normal cellular functioning. The ubiquitin proteasome system mediates the intracellular protein degradation that is required for normal cellular homeostasis. The 26S proteasome is a multi-enzyme protease that degrades redundant proteins; conversely, inhibition of proteasomal degradation results in intracellular aggregation of unwanted proteins and cell death. This observation led to the development of proteasome inhibitors as therapeutics for use in cancer. The clinical applicability of targeting proteasomes is exemplified by the recent FDA approval of the first proteasome inhibitor, bortezomib, for the treatment of relapsed/refractory multiple myeloma. Although bortezomib represents a major advance in the treatment of this disease, it can be associated with toxicity and the development of drug resistance. Importantly, extensive preclinical studies suggest that combination therapies can both circumvent drug resistance and reduce toxicity. In addition, promising novel proteasome inhibitors, which are distinct from bortezomib, and exhibit equipotent anti-multiple myeloma activities, are undergoing clinical evaluation in order to improve patient outcome in multiple myeloma.

Publication history

Republished from Current BioData's Targeted Proteins database (TPdb; http://www.targetedproteinsdb.com).

     

Structure,reactivity and luminescence studies of triphenylsiloxide complexes of tetravalent lanthanides

Among the 14 lanthanide elements (Ce–Lu), until recently, the tetravalent oxidation state was readily accessible in solution only for cerium while Pr(iv), Nd(iv), Dy(iv) and Tb(iv) had only been detected in the solid state. The triphenylsiloxide ligand recently allowed the isolation of molecular complexes of Tb(iv) and Pr(iv) providing an unique opportunity of investigating the luminescent properties of Ln(iv) ions. Here we have expanded the coordination studies of the triphenylsiloxide ligand with Ln(iii) and Ln(iv) ions and we report the first observed luminescence emission spectra of Pr(iv) complexes which are assigned to a ligand-based emission on the basis of the measured lifetime and computational studies. Binding of the ligand to the Pr(iv) ion leads to an unprecedented large shift of the ligand triplet state which is relevant for future applications in materials science.

The first observed luminescence emission spectra of Pr(iv) complexes are assigned to a ligand-based emission. Binding of the triphenylsiloxide ligand to the Pr(iv) ion leads to an unprecedented large red shift of its triplet state.