Hiroaki Tokiwa

The active form of vitamin D₃, 1α,25-dihydroxyvitamin D₃ [1,25(OH)₂D₃], is a major regulator of calcium homeostasis through activation of the vitamin D receptor (VDR). We have previously synthesized vitamin D derivatives with large adamantane (AD) rings at position 24, 25, or 26 of the side chain to study VDR agonist and/or antagonist properties. One of them—ADTK1, with an AD ring and 23,24-triple bond—shows a high VDR affinity and cell-selective VDR activity. In this study, we synthesized novel vitamin D derivatives (ADKM1-6) with an alkyl group substituted at position 25 of ADTK1 to develop more cell-selective VDR ligands. ADKM2, ADKM4, and ADKM6 had VDR transcriptional activity comparable to 1,25(OH)₂D₃ and ADTK1, although their VDR affinities were weaker. Interestingly, ADKM2 has selective VDR activity in kidney- and skin-derived cells—a unique phenotype that differs from ADTK1. Furthermore, ADKM2, ADKM4, and ADKM6 induced osteoblast differentiation in human dedifferentiated fat cells more effectively than ADTK1. The development of vitamin D derivatives with bulky modifications such as AD at position 24, 25, or 26 of the side chain is useful for increased stability and tissue selectivity in VDR-targeting therapy.

In recent years, London dispersion interactions, which are the attractive component of the van der Waals potential, have been found to play an important role in controlling the regio- and/or stereoselectivity of various reactions. Particularly, the dispersion interactions between substrates and catalysts (or ligands) are dominant in various selective catalyzes. In contrast, repulsive steric interactions, rather than the attractive dispersion interactions, between bulky substituents are predominant in most of the noncatalytic reactions. Herein, we demonstrate the first example of London dispersion-controlled noncatalytic (2 + 2) cyclodimerization of substituted benzynes to selectively afford proximal biphenylenes in high yields and regioselectivities, depending on the extent of dispersion interactions in the substituents. This method can be applied for the synthesis of novel helical biphenylenes, which would be fascinating for chemists as these compounds are potential skeletons for ligands, catalysts, and medicines.

Infection of hosts by morbilliviruses is facilitated by the interaction between viral hemagglutinin (H-protein) and the signaling lymphocytic activation molecule (SLAM). Recently, the functional importance of the n-terminal region of human SLAM as a measles virus receptor was demonstrated. However, the functional roles of this region in the infection process by other morbilliviruses and host range determination remain unknown, partly because this region is highly flexible, which has hampered accurate structure determination of this region by X-ray crystallography. In this study, we analyzed the interaction between the H-protein from canine distemper virus (CDV-H) and SLAMs by a computational chemistry approach. Molecular dynamics simulations and fragment molecular orbital analysis demonstrated that the unique His28 in the N-terminal region of SLAM from Macaca is a key determinant that enables the formation of a stable interaction with CDV-H, providing a basis for CDV infection in Macaca. The computational chemistry approach presented should enable the determination of molecular interactions involving regions of proteins that are difficult to predict from crystal structures because of their high flexibility.

The selective PPARα modulator (SPPARMα) is expected to medicate dyslipidemia with minimizing adverse effects. Recently, pemafibrate was screened from the ligand library as an SPPARMα bearing strong potency. Several clinical pieces of evidence have proved the usefulness of pemafibrate as a medication; however, how pemafibrate works as a SPPARMα at the molecular level is not fully known. In this study, we investigate the molecular mechanism behind its novel SPPARMα character through a combination of approaches of X-ray crystallography, isothermal titration calorimetry (ITC), and fragment molecular orbital (FMO) analysis. ITC measurements have indicated that pemafibrate binds more strongly to PPARα than to PPARγ. The crystal structure of PPARα-ligand binding domain (LBD)/pemafibrate/steroid receptor coactivator-1 peptide (SRC1) determined at 3.2 Å resolution indicates that pemafibrate binds to the ligand binding pocket (LBP) of PPARα in a Y-shaped form. The structure also reveals that the conformation of the phenoxyalkyl group in pemafibrate is flexible in the absence of SRC1 coactivator peptide bound to PPARα; this gives a freedom for the phenoxyalkyl group to adopt structural changes induced by the binding of coactivators. FMO calculations have indicated that the accumulation of hydrophobic interactions provided by the residues at the LBP improve the interaction between pemafibrate and PPARα compared with the interaction between fenofibrate and PPARα.

Ligands for retinoid X receptors (RXRs), “rexinoids”, are attracting interest as candidates for therapy of type 2 diabetes, Alzheimer’s and Parkinson’s diseases. However, current screening methods for rexinoids are slow and require special apparatus or facilities. Here, we created 7-hydroxy-2-oxo-6-(3,5,5,8,8-pentamethyl-5,6,7,8-tetrahydronaphthalen-2-yl)-2H-chromene-3-carboxylic acid (10, CU-6PMN) as a new fluorescent RXR agonist and developed a screening system of rexinoids using 10. Compound 10 was designed based on the fact that umbelliferone emits strong fluorescence in a hydrophilic environment, but the fluorescence intensity decreases in hydrophobic environments such as the interior of proteins. The developed assay using 10 enabled screening of rexinoids to be performed easily within a few hours by monitoring changes of fluorescence intensity with widely available fluorescence microplate readers, without the need for processes such as filtration.

A total synthesis of an anticancer xanthone natural product termicalcicolanone A utilizing multiple nucleophilic aromatic substitutions and pericyclic reactions has been developed. The pyrano[3,2- b]xanthen-6-one scaffold was constructed via NHC-catalyzed aroylation to produce the benzophenone intermediate, Claisen cyclization to form the pyran ring, and intramolecular 1,4-addition to construct the xanthone framework. The prenyl group was introduced in the final stages of the synthesis through regioselective Claisen rearrangement. The synthesis has been achieved in 19 steps.

In modern praxis, a knowledge-driven design of pharmaceutical compounds relies heavily on protein structure data. Nonetheless, quantification of the interaction between protein and ligand is of great importance in the theoretical evaluation of the ability of a pharmaceutical compound to comply with certain expectations. The FMO (fragment molecular orbital) method is handy in this regard. However, the physical complexity and the number of the interactions within a protein-ligand complex render analysis of the results somewhat complicated. This situation prompted us to develop the toolkit AnalysisFMO, which should enable a more efficient and convenient workflow with FMO data generated by quantum-chemical packages such as GAMESS, PAICS, and ABINIT-MP. AnalysisFMO consists of two separate units, RbAnalysisFMO, and the PyMOL plugins. The former can extract inter-fragment interaction energies (IFIEs) or pair interaction energies (PIEs) from the FMO output files generated by the aforementioned quantum-chemical packages. The PyMOL plugins enable visualization of IFIEs or PIEs in the protein structure in PyMOL. We demonstrate the use of this tool on a lectin protein from Burkholderia cenocepacia in which FMO analysis revealed the existence of a new interaction between Gly84 and fucose. Moreover, we found that second-shell interactions are crucial in forming the sugar binding site. In the case of bilirubin oxidase from Myrothecium verrucaria (MvBO), we predict that interactions between Asp105 and three His residues (His401, His403, and His136) are essential for optimally positioning the His residues to coordinate Cu atoms to form a Type 2 and two Type 3 Cu ions.

Protein residue networks (PRN) from energetic and geometric data are probably not identical. PRNs constructed from ab initio pair interaction energies are analyzed for the first time and compared to PRN based on center of mass separation. We use modern, previously unused algorithms such as global and local efficiencies to quantitatively confirm that both types of PRNs do exhibit small-world character. The main novelty finding is that interaction energy-based PRNs preserve small-world character even when clustered. A node hierarchy independent of the cutoff energy used for the edge creation is characteristic for them. Efficiency centrality identifies hubs responsible for such behavior. The interaction energy-based PRNs seem to comply with the scale-free network model with respect to efficiency centrality distribution as opposed to distance based PRNs. Community detection is introduced into protein network research as an extension beyond cluster analysis to study tertiary and quaternary structures.

1‐[(3,5,5,8,8‐pentamethyl‐5,6,7,8‐tetrahydronaphthalen‐2‐yl)amino]benzotriazole‐5‐carboxylic acid (CBt‐PMN), a partial agonist of retinoid X receptor (RXR), has attracted attention due to its potential to treat type 2 diabetes and central nervous system diseases with reduced adverse effects of existing full agonists. Herein, we report the crystal structure of CBt‐PMN‐bound ligand‐binding domain of human RXRα (hRXRα) and its biochemical characterization. Interestingly, the structure is a tetramer in nature, in which CBt‐PMNs are clearly found binding in two different conformations. The dynamics of the hRXRα/CBt‐PMN complex examined using molecular dynamics simulations suggest that the flexibility of the AF‐2 interface depends on the conformation of the ligand. These facts reveal that the dual conformation of CBt‐PMN in the complex is probably the reason behind its partial agonistic activity.