X-ray Crystallography Research Articles

Structure-based mapping of the histone-binding pocket of KDM4D using functionalized tetrazole and pyridine core compounds

KDM4 histone demethylases became an exciting target for inhibitor development as the evidence linking them directly to tumorigenesis mounts. In this study, we set out to better understand the binding cavity using an X-ray crystallographic approach to provide a detailed landscape of possible interactions within the under-investigated region of KDM4. Our design strategy was based on utilizing known KDM binding motifs, such as nicotinic acid and tetrazolylhydrazides, as core motifs that we decided to enrich with flexible tails to map the distal histone binding site. The resulting X-ray structures of the novel compounds bound to KDM4D, a representative of the KDM4 family, revealed the interaction pattern with distal residues in the histone-binding site. The most prominent protein rearrangement detected upon ligand binding is the loop movement that blocks the accessibility to the histone binding site. Apart from providing new sites that potential inhibitors can target, the novel compounds may prove helpful in exploring the capacity of ligands to bind in sites distal to the cofactor-binding site of other KDMs or 2-oxoglutarate (2OG)-dependent oxygenases. The case study proves that combining a strong small binding motif with flexible tails to probe the binding pocket will facilitate lead discovery in classical drug-discovery campaigns, given the ease of accessing X-ray quality crystals.

Conformational disorder in the crystal structure of methyl 2-acetamido-2-deoxy-β-D-glucopyranosyl-(1→4)-2-acetamido-2-deoxy-β-D-glucopyranoside (methyl β-chitobioside) methanol monosolvate.

Methyl 2-acetamido-2-deoxy-β-D-glucopyranosyl-(1→4)-2-acetamido-2-deoxy-β-D-glucopyranoside (methyl β-chitobioside), (IV), crystallizes from aqueous methanol at room temperature to give a structure (C17H30N2O22·CH3OH) containing conformational disorder in the exocyclic hydroxymethyl group of one of its βGlcNAc residues. As observed in other X-ray structures of disaccharides containing β-(1→4) O-glycosidic linkages, inter-residue hydrogen bonding between O3H of the βGlcNAc bearing the OCH3 aglycone and O5 of the adjacent βGlcNAc is observed based on the 2.79 Å internuclear distance between the O atoms. The structure of (IV) was compared to that determined previously for 2-acetamido-2-deoxy-β-D-glucopyranosyl-(1→4)-2-acetamido-2-deoxy-β-D-glucopyranose (β-chitobiose), (III). The O-glycosidic linkage torsion angles, phi (φ) and psi (ψ), in (III) and (IV) differ by 6-8°. The N-acetyl side chain conformation in (III) and (IV) shows some context dependence, with the C1-C2-N-Ccar torsion angle 10-15° smaller for the βGlcNAc residue involved in the internal O-glycosidic linkage. In (IV), conformational disorder is observed in the exocyclic hydroxymethyl (-CH2OH) group in the βGlcNAc residue bearing the OCH3 aglycone, and a fitting of the electron density indicates an approximate 50:50 distribution of the gauche-gauche (gg) and gauche-trans (gt) conformers in the lattice. Similar behavior is not observed in (III), presumably due to the different packing structure in the vicinity of the -CH2OH substituent that affects its ability to hydrogen bond to proximal donors/acceptors. Unlike (IV), a re-examination of the previously reported electron density of (III) revealed conformational disorder in the N-acetyl side chain attached to the reducing-end βGlcNAc residue caused by rotation about the C2-N bond.

Intersite communication in dimeric enzymes highlighted by structural and thermodynamic analysis of didansyltyrosine binding to thymidylate synthases

Intersite communication in dimeric enzymes, triggered by ligand binding, represents both a challenge and an opportunity in enzyme inhibition strategy. Though often understestimated, it can impact on the in vivo biological mechansim of an inhibitor and on its pharmacokinetics. Thymidylate synthase (TS) is a homodimeric enzyme present in almost all living organisms that plays a crucial role in DNA synthesis and cell replication. While its inhibition is a valid strategy in the therapy of several human cancers, designing specific inhibitors of bacterial TSs poses a challenge to the development of new anti-infective agents. N,O-didansyl-l-tyrosine (DDT) inhibits both Escherichia coli TS (EcTS) and Lactobacillus casei TS (LcTS). The available X-ray structure of the DDT:dUMP:EcTS ternary complex indicated an unexpected binding mode for DDT to EcTS, involving a rearrangement of the protein and addressing the matter of communication between the two active sites of an enzyme dimer. Combining molecular-level information on DDT binding to EcTS and LcTS extracted from structural and FRET-based fluorometric evidence with a thermodynamic characterization of these events obtained by fluorometric and calorimetric titrations, this study unveiled a negative cooperativity between the DDT bindings to the two monomers of each enzyme dimer. This result, complemented by the species-specific thermodynamic signatures of the binding events, implied that communication across the protein dimer was triggered by the first DDT binding. These findings could challenge the conventional understanding of TS inhibition and open the way for the development of novel TS inhibitors with a different mechanism of action and enhanced efficacy and specificity.

DNA binding/cleavage study of mixed ligand complexes of di-/triorganotin 4-fluroanthranilates with 1,10-phenanthroline/benzimidazole: Synthesis, characterization, X-ray structures of [Bu2Sn(FA)2.phen], [Me3Sn(FA).bz], and DFT calculations

The mixed ligand complexes of di-/tri-organotin 4-fluoroanthranilates with 1,10-phenanthroline (phen), viz. [R2Sn(FA)2.phen] {R=Bu (1), Oct (2), Me (3)} and benzimidazole (bz), viz. [R3Sn(FA).bz] {R=Me (4), Ph (5)} have been synthesized and satisfactorily characterized by analytical and spectroscopic techniques, especially FT-IR, NMR (1H, 13C, and 119Sn NMR) spectroscopy and ESI-MS spectrometry. Single crystal-Xray diffraction measurements (Sc-XRD) confirmed a distorted pentagonal-bipyramidal geometry for complex (1), where one carboxylate group is monodentate, and the second (–COO) group behaves as a bidentate, while complex (4) affirms trigonal–bipyramidal geometry around the tin centre. The bond angles, especially ∠C–Sn–C obtained by optimizing complexes’ gaseous phase geometry showed significant conformity with the data obtained from Sc-XRD and NMR. The Kb (binding constant) and Ksv (Stern-Volmer quenching constant) values, ranging from 103 to 107, have been obtained through UV–visible absorption and fluorescence titrations with CT-DNA. These results strongly indicate partial minor groove binding/electrostatic interactions between the complex and CT-DNA, which has also been evidenced by circular dichroism studies. The studied complexes cleaved efficiently pBR322 plasmid DNA in NC (nicked circular form), SC (supercoiled form), and linear form

Synthesis, Crystal Structure, Hirshfeld Surface Analysis, and Computational Approach of a New Pyrazolo[3,4-g]isoquinoline Derivative as Potent against Leucine-Rich Repeat Kinase 2 (LRRK2).

Ethyl-2-((8-cyano-3,5,9a-trimethyl-1-(4-oxo-4,5-dihydrothiazol-2-yl)-4-phenyl-3a,4,9,9a-tetrahydro-1H-pyrazolo[3,4-g]isoquinolin-7-yl)thio)acetate (5) was synthesized, and its structure was characterized by IR, MS, and NMR (1H and 13C) and verified by a single-crystal X-ray structure determination. Compound 5 adopts a "pincer" conformation. In the crystal, the hydrogen bonds of -H···O, C-H···O, and O-H···S form thick layers of molecules that are parallel to (101). The layers are linked by C-H···π(ring) interactions. The Hirshfeld surface analysis shows that intermolecular hydrogen bonding plays a more important role than both intramolecular hydrogen bonding and π···π stacking in the crystal. The intramolecular noncovalent interactions in 5 were studied by QTAIM, NCI, and DFT-NBO calculations. Based on structural activity relationship studies, leucine-rich repeat kinase 2 (LRRK2) was found to bind 5 and was further subjected to molecular docking studies, molecular dynamics, and ADMET analysis to probe potential drug candidacy.

Structure-activity relationships of middle-size cyclic peptides, KRAS inhibitors derived from an mRNA display

Cyclic peptides are attracting attention as therapeutic agents due to their potential for oral absorption and easy access to tough intracellular targets. LUNA18, a clinical KRAS inhibitor, was transformed—without scaffold hopping—from the initial hit by using an mRNA display library that met our criteria for drug-likeness. In drug discovery using mRNA display libraries, hit compounds always possess a site linked to an mRNA tag. Here, we describe our examination of the Structure-Activity Relationship (SAR) using X-ray structures for chemical optimization near the site linked to the mRNA tag, equivalent to the C-terminus. Structural modifications near the C-terminus demonstrated a relatively wide range of tolerance for side chains. Furthermore, we show that a single atom modification is enough to change the pharmacokinetic (PK) profile. Since there are four positions where side chain modification is permissible in terms of activity, it is possible to flexibly adjust the pharmacokinetic profile by structurally optimizing the side chain. The side chain transformation findings demonstrated here may be generally applicable to hits obtained from mRNA display libraries.

Cobalt(II)-Substituted Cysteamine Dioxygenase Oxygenation Proceeds through a Cobalt(III)-Superoxo Complex.

We investigated the metal-substituted catalytic activity of human cysteamine dioxygenase (ADO), an enzyme pivotal in regulating thiol metabolism and contributing to oxygen homeostasis. Our findings demonstrate the catalytic competence of cobalt(II)- and nickel(II)-substituted ADO in cysteamine oxygenation. Notably, Co(II)-ADO exhibited superiority over Ni(II)-ADO despite remaining significantly less active than the natural enzyme. Structural analyses through X-ray crystallography and cobalt K-edge excitation confirmed successful metal substitution with minimal structural perturbations. This provided a robust structural basis, supporting a conserved catalytic mechanism tailored to distinct metal centers. This finding challenges the proposed high-valent ferryl-based mechanism for thiol dioxygenases, suggesting a non-high-valent catalytic pathway in the native enzyme. Further investigation of the cysteamine-bound or a peptide mimic of N-terminus RGS5 bound Co(II)-ADO binary complex revealed the metal center's high-spin (S = 3/2) state. Upon reaction with O2, a kinetically and spectroscopically detectable intermediate emerged with a ground spin state of S = 1/2. This intermediate exhibits a characteristic 59Co hyperfine splitting (A = 67 MHz) structure in the EPR spectrum alongside UV-vis features, consistent with known low-spin Co(III)-superoxo complexes. This observation, unique for protein-bound thiolate-ligated cobalt centers in a protein, unveils the capacities for O2 activation in such metal environments. These findings provide valuable insights into the non-heme iron-dependent thiol dioxygenase mechanistic landscape, furthering our understanding of thiol metabolism regulation. The exploration of metal-substituted ADO sheds light on the intricate interplay between metal and catalytic activity in this essential enzyme.

Structural analysis of the FERM domain of human protein tyrosine phosphatase non-receptor type 21.

Protein tyrosine phosphatase non-receptor type 21 (PTPN21) is a cytosolic protein tyrosine phosphatase that regulates cell growth and invasion. Due to its oncogenic properties, PTPN21 has recently emerged as a potential therapeutic target for cancer. In this study, the three-dimensional structure of the PTPN21 FERM domain was determined at 2.1 Å resolution by X-ray crystallography. The crystal structure showed that this domain harbors canonical FERM folding and consists of three subdomains that are tightly packed via highly conserved intramolecular hydrophobic interactions. Consistent with this, the PTPN21 FERM domain shares high structural homology with several other FERM domains. Moreover, structural superimposition demonstrated two putative protein-binding sites of the PTPN21 FERM domain, which are presumed to be associated with interaction with its binding partner, kinesin family member 1C. Thus, these data suggest that the FERM domain of PTPN21 serves as a module that mediates protein-protein interaction, like other FERM domains.

Exploring crystal structure-physical property relationships with pressure.

From its conception, X-ray crystallography has provided a unique understanding of the structure, bonding and electronic state of materials, which, in turn, unlocks a means of examining the properties and function of crystalline systems. Using state-of-the-art single-crystal X-ray diffraction, along with UV-Vis spectroscopy and DFT calculations, Zwolenik et al. [(2024). IUCrJ, 11, 519-527] have provided a comprehensive study of the structure-optical property relationship of 1,3-diacetylpyrene with methodologies that are increasingly accessible to non-specialist laboratories.

Shy is a proteobacterial steroid hydratase which catalyzes steroid side chain degradation without requiring a catalytically inert partner domain

Shy (side chain hydratase) and Sal (side chain aldolase), are involved in successive reactions in the pathway of bile acid side chain catabolism in Proteobacteria. Untagged Shy co-purified with His-tagged Sal indicating that the two enzymes form a complex. Shy contains a MaoC and a DUF35 domain. When co-expressed with Sal, the DUF35 domain but not the MaoC domain of Shy was observed to co-purify with Sal, indicating Sal interacts with Shy through its DUF35 domain. The MaoC domain of Shy (ShyMaoC) remained catalytically viable and could hydrate cholyl-enoyl-CoA with similar catalytic efficiency as in the Shy-Sal complex. Sal expressed with the DUF35 domain of Shy (Sal-ShyDUF35) was similarly competent for the retroaldol cleavage of cholyl-3-OH-CoA. ShyMaoC showed a preference for C5 side chain bile acid substrates, exhibiting low activity towards C3 side chain substrates. The ShyMaoC structure was determined by X-ray crystallography, showing a hot dog fold with a short central helix surrounded by a twisted anti-parallel β-sheet. Modeling and mutagenesis studies suggest that the bile acid substrate occupies the large open cleft formed by the truncated central helix and repositioning of the active site housing. ShyMaoC therefore contains two substrate binding sites per homodimer, making it distinct from previously characterized MaoC steroid hydratases that are (pseudo)-heterodimers with one substrate binding site per dimer. The characterization of Shy provides insight into how MaoC family hydratases have adapted to accommodate large polycyclic substrates that can facilitate future engineering of these enzymes to produce novel steroid pharmaceuticals.

Binding of steroid substrates reveals the key to the productive transition of the cytochrome P450 OleP.

OleP is a bacterial cytochrome P450 involved in oleandomycin biosynthesis as it catalyzes regioselective epoxidation on macrolide intermediates. OleP has recently been reported to convert lithocholic acid (LCA) into murideoxycholic acid through a highly regioselective reaction and to unspecifically hydroxylate testosterone (TES). Since LCA and TES mainly differ by the substituent group at the C17, here we used X-ray crystallography, equilibrium binding assays, and molecular dynamics simulations to investigate the molecular basis of the diverse reactivity observed with the two steroids. We found that the differences in the structure of TES and LCA affect the capability of these molecules to directly form hydrogen bonds with N-terminal residues of OleP internal helix I. The establishment of these contacts, by promoting the bending of helix I, fosters an efficient trigger of the open-to-closed structural transition that occurs upon substrate binding to OleP and contributes to the selectivity of the subsequent monooxygenation reaction.

Unusually Short H⋅⋅⋅H Contacts in Intramolecularly Cyclized Helically Fused Anthracenes.

The intramolecular coupling of dichloro-substituted helically fused anthracenes using the Yamamoto coupling yielded cyclized products with sterically congested molecular structures. The X-ray analysis and DFT calculations showed that the aromatic framework adopted a nonplanar structure with a twisted conformation about the newly formed single bond, which acts as a chiral axis. Interestingly, the X-ray structure obtained through the Hirshfeld atom refinement revealed short interatomic distances between the inner hydrogen atoms (1.648-1.692 Å), much shorter than the sum of their van der Waals radii. Owing to these unusually short contacts, the 1H NMR spectrum exhibited a significant deshielding (12.5 ppm) and a large nuclear Overhauser effect (44 %). Additionally, the IR spectrum displayed a high-frequency shift of the C-H stretching vibration. These observations, along with the noncovalent interaction plot indicative of a characteristic steric environment, strongly support the presence of steric hindrance. Moreover, dynamic NMR measurement of the mesityl-substituted derivative yielded a barrier to helical inversion of 84 kJ mol-1. The optical properties and crystal packing of the cyclized products are also reported.

Comprehensive encoding of conformational and compositional protein structural ensembles through the mmCIF data structure.

In the folded state, biomolecules exchange between multiple conformational states crucial for their function. However, most structural models derived from experiments and computational predictions only encode a single state. To represent biomolecules accurately, we must move towards modeling and predicting structural ensembles. Information about structural ensembles exists within experimental data from X-ray crystallography and cryo-electron microscopy. Although new tools are available to detect conformational and compositional heterogeneity within these ensembles, the legacy PDB data structure does not robustly encapsulate this complexity. We propose modifications to the macromolecular crystallographic information file (mmCIF) to improve the representation and interrelation of conformational and compositional heterogeneity. These modifications will enable the capture of macromolecular ensembles in a human and machine-interpretable way, potentially catalyzing breakthroughs for ensemble-function predictions, analogous to the achievements of AlphaFold with single-structure prediction.

Deciphering Structural Traits for Thermal and Kinetic Stability across Protein Family Evolution through Ancestral Sequence Reconstruction.

Natural proteins are frequently marginally stable, and an increase in environmental temperature can easily lead to unfolding. As a result, protein engineering to improve protein stability is an area of intensive research. Nonetheless, since there is usually a high degree of structural homology between proteins from thermophilic organisms and their mesophilic counterparts, the identification of structural determinants for thermoadaptation is challenging. Moreover, in many cases, it has become clear that the success of stabilization strategies is often dependent on the evolutionary history of a protein family. In the last few years, the use of ancestral sequence reconstruction (ASR) as a tool for elucidation of the evolutionary history of functional traits of a protein family has gained strength. Here, we used ASR to trace the evolutionary pathways between mesophilic and thermophilic kinases that participate in the biosynthetic pathway of vitamin B1 in bacteria. By combining biophysics approaches, X-ray crystallography, and molecular dynamics simulations, we found that the thermal stability of these enzymes correlates with their kinetic stability, where the highest thermal/kinetic stability is given by an increase in small hydrophobic amino acids that allow a higher number of interatomic hydrophobic contacts, making this type of interaction the main support for stability in this protein architecture. The results highlight the potential benefits of using ASR to explore the evolutionary history of protein sequence and structure to identify traits responsible for the kinetic and thermal stability of any protein architecture.

Palladium(II) complexes of amino acid derived ONO-tridentate enamine ligand for Suzuki-Miyaura cross-coupling reaction

Three palladium(II) complexes with an ONO-tridentate enamine ligand were synthesized by reacting the glycine derivative ligand with palladium(II) acetate in the presence of selected ancillary ligands. The structures of the Pd(II) complexes were determined through melting point, C, H, N elemental analysis, ATR-FTIR, and NMR. Additionally, the molecular structures of two complexes were confirmed using single-crystal X-ray crystallography studies, revealing that the ligand coordinates to palladium(II) via ONO-chelation, resulting in a square planar geometry around the metal center. Catalytic activities of these compounds were assessed in the Suzuki-Miyaura cross-coupling reaction of phenylboronic acid with the different aryl halides.

Protecting Group-Free Total Synthesis of (-)-Boscartin H.

Herein, we report the first protecting group-free total synthesis of (-)-boscartin H, which features a 5-12-5-fused tricyclic structure. The key steps, which include a diastereoselective THF-ring-forming/aldol reaction sequence and ring-closing metathesis, afforded high stereoselectivity with (-)-boscartin H obtained in 3.6% overall yield using a 11-step long linear sequence. In addition, X-ray crystallography clearly confirmed the stereochemistry of boscartin H.

Automated multiconformer model building for X-ray crystallography and cryo-EM.

In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift toward modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior Rfree and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g., Coot) and fit can be further improved by refinement using standard pipelines (e.g., Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.

Automated multiconformer model building for X-ray crystallography and cryo-EM

In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift toward modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior Rfree and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g., Coot) and fit can be further improved by refinement using standard pipelines (e.g., Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.

Targeted discovery of unusual diterpenoids with anti-fungal activity from the root of Euphorbia lathyris

Lathyrisone A (1), a diterpene with an undescribed tricyclic 6/6/6 fused carbon skeleton, along with spirolathyrisins B-D (3–5), three diterpenes with a rare [4.5.0] spirocyclic carbon skeleton, and one known compound (2) were isolated from the roots of Euphorbia lathyris. Their chemical structures were characterized by extensive spectroscopic analysis, X-ray crystallography, ECD and quantum chemistry calculation. A plausible biosynthetic pathway for compounds 1–5 was proposed, which suggested it is a competitive pathway for ingenol biosynthesis in the plant. The anti-fungal activities of these compounds were tested, especially, compound 2 showed stronger anti-fungal activities against Fusarium oxysporum and Alternaria alternata than the positive control fungicide thiophanate-methyl. The preliminary structure–activity relationship of compounds 1–5 was also discussed. These results not only expanded the chemical diversities of E. lathyris, but also provided a lead compound for the control of plant pathogens.

SIRT5 mutants reveal the role of conserved asparagine and glutamine residues in the NAD+-binding pocket.

SIRT5, one of the mammalian sirtuins, specifically recognizes succinyl-lysine residues on proteins and catalyzes the desuccinylation reaction. In this study, we characterized SIRT5 mutants with hydrophobic amino acid substitutions at Q140 and N141, in addition to the catalytic residue H158, known as an active site residue, by the Michaelis-Menten analysis and X-ray crystallography. Kinetic analysis showed that the catalytic efficiency (kcat/Km) of the Q140L and N141V mutants decreased to 0.02 times and 0.0038 times that of the wild-type SIRT5, respectively, with the activity of the N141V mutant becoming comparable to that of the H158M mutant. Our findings indicate that N141 contributes significantly to the desuccinylation reaction.