Home | Research Center for Computational Science

Electrochemical cascade access to hetero[8]circulenes as potent organophotocatalysts for diverse C–X bond formations

The chemistry of hetero[8]circulenes has been limited to three types due to challenges in making unsymmetrical variants. We introduce an electrochemical approach to a hetero[8]circulene bearing five hexagons and three pentagons. Exploiting electrochemistry’s selectivity, differential oxidation enables selective inter- and intramolecular couplings, giving dioxaza[8]circulenes in up to 83% yield with good functional-group tolerance. A one-pot protocol from commercially available substrates forges six new bonds. Structural, optical, and electrochemical characterizations (X-ray crystallography, spectrophotometric analysis) with DFT calculations elucidate properties. Motivated by these features, we demonstrate transition-metal-free organophotocatalysis under 365-nm LED irradiation, effecting diverse C–X (X = C, B, P) bond formations in up to 97% yield.

A. S. Gabr, M. S. H. Salem*, M. I. Khalid, R. Takahashi, Y. Nishimoto, M. Yasuda, S. Takizawa*, Nat. Commun., 16, 5682 (2025).

Electrochemical cascade access to hetero[8]circulenes as potent organophotocatalysts for diverse C–X bond formations

The chemistry of hetero[8]circulenes has been limited to three types due to challenges in making unsymmetrical variants. We introduce an electrochemical approach to a hetero[8]circulene bearing five hexagons and three pentagons. Exploiting electrochemistry’s selectivity, differential oxidation enables selective inter- and intramolecular couplings, giving dioxaza[8]circulenes in up to 83% yield with good functional-group tolerance. A one-pot protocol from commercially available substrates forges six new bonds. Structural, optical, and electrochemical characterizations (X-ray crystallography, spectrophotometric analysis) with DFT calculations elucidate properties. Motivated by these features, we demonstrate transition-metal-free organophotocatalysis under 365-nm LED irradiation, effecting diverse C–X (X = C, B, P) bond formations in up to 97% yield.

A. S. Gabr, M. S. H. Salem*, M. I. Khalid, R. Takahashi, Y. Nishimoto, M. Yasuda, S. Takizawa*, Nat. Commun., 16, 5682 (2025).

Molecular Simulation × AI Reveals Structural Changes of Transporter Proteins

Transporter proteins exist in the plasma membrane and transport various substrates into and out of the cell. It is believed that substrates are transported by a mechanism involving conformational changes called alternating access. Therefore, it is necessary to elucidate the structural changes in order to understand the mechanism. We have elucidated the conformational changes of the oxalate transporter, which transports oxalic acid, a substance contained in vegetables and nuts that causes urinary tract stones, by combining molecular dynamics (MD) simulations and the structure prediction AI AlphaFold. MD simulations predicted the inward-open conformation, which is experimentally difficult to analyze, and AlphaFold successfully predicted mutations that stabilize the inward-open conformation.

Accelerated Molecular Dynamics and AlphaFold Uncover a Missing Conformational State of Transporter Protein OxlT
J. Ohnuki, T. Jaunet-Lahary, A. Yamashita, and K. Okazaki, J. Phys. Chem. Lett. 15, 725–732 (2024).

Electronic structure informatics to search for hightly efficient photovoltaic materials using singlet fission

The transition from fossil fuels to environment-friendly and renewable sources of energy has become one of the central challenges for today’s society and across scientific disciplines. Though global energy studies report that there has been record-breaking growth in sustainable energy production, further improvements are desirable. Solar energy is a promising solution for the near future, owing to its availability and potentially high output.

In this study, an "electronic structure informatics" approach combining the semiempirical molecular orbital method, time-dependent density functional theory, and machine learning has been applied to rapidly search for materials that can be used as photovoltaic cell materials with high efficiency by multiple exciton generation using singlet fission from over 4 million derivatives based on indigo, a famous dye of natural origin. The obtained promising candidate molecules and their possible candidates were identified and compared with those not. Random forest classification of the chemical structures of the obtained promising and non-promising candidate molecules led to the discovery of chemical structure rules that effectively show singlet fission and the identification of candidate materials with symmetric chemical structures and sufficient synthesizability in the order of several hundred.

Weber F., Mori H., npj Comput. Mater., 8, 176 (2022).

Electronic structure informatics to search for hightly efficient photovoltaic materials using singlet fission

The transition from fossil fuels to environment-friendly and renewable sources of energy has become one of the central challenges for today’s society and across scientific disciplines. Though global energy studies report that there has been record-breaking growth in sustainable energy production, further improvements are desirable. Solar energy is a promising solution for the near future, owing to its availability and potentially high output.

In this study, an "electronic structure informatics" approach combining the semiempirical molecular orbital method, time-dependent density functional theory, and machine learning has been applied to rapidly search for materials that can be used as photovoltaic cell materials with high efficiency by multiple exciton generation using singlet fission from over 4 million derivatives based on indigo, a famous dye of natural origin. The obtained promising candidate molecules and their possible candidates were identified and compared with those not. Random forest classification of the chemical structures of the obtained promising and non-promising candidate molecules led to the discovery of chemical structure rules that effectively show singlet fission and the identification of candidate materials with symmetric chemical structures and sufficient synthesizability in the order of several hundred.

Weber F., Mori H., npj Comput. Mater., 8, 176 (2022).

The mechanism of divalent cation block of ion channel activity

Nervous activity is subjected to complex regulation by ion channels. In this study, we reproduced the divalent cations block in molecular dynamics simulation and elucidated the blocking mechanism. It was observed that introducing a hydrophilic environment into the channel lumen allows the stable existence of divalent cations. This retention of divalent cations makes it difficult for ions to permeate the ion channel. This result also provides new knowledge about the reproduction of polyvalent ions in molecular dynamics simulations.

K. Irie, Y. Oda, T. Sumikama, A. Oshima, Y. Fujiyoshi, Nat. Commun., 14, 4236 (2023).

The mechanism of divalent cation block of ion channel activity

Nervous activity is subjected to complex regulation by ion channels. In this study, we reproduced the divalent cations block in molecular dynamics simulation and elucidated the blocking mechanism. It was observed that introducing a hydrophilic environment into the channel lumen allows the stable existence of divalent cations. This retention of divalent cations makes it difficult for ions to permeate the ion channel. This result also provides new knowledge about the reproduction of polyvalent ions in molecular dynamics simulations.

K. Irie, Y. Oda, T. Sumikama, A. Oshima, Y. Fujiyoshi, Nat. Commun., 14, 4236 (2023).

Au-Ag nanoclusters showing intense phosphorescence in cells

Luminescent metal nanoclusters are expected to exhibit unique physical properties in the cluster structure depending on the ligand structure, metal type, number of nuclei and arrangement. In this study, carbon-centered gold-silver (CAu₆Ag₂) clusters with N-heterocyclic carbene (NHC) ligands were designed and synthesized, and it was found that these clusters emit strong phosphorescence in solution, and the contribution of NHC ligands to phosphorescence emission was revealed by theoretical calculation. The luminescence rate constant was calculated by an analysis including spin-orbit interactions, and the quantum yield was discussed in terms of the energy barrier to the minimum energy crossing point. Furthermore, the phosphorescent gold-silver clusters with long luminescence lifetime were used for cellular imaging, which revealed the pathway of uptake into the cell and selective localization to specific organelles, confirming their superior functionality, which is different from the non-selective uptake of conventional phosphine ligands.
Z. Lei, M. Endo, H. Ube, T. Shiraogawa, P. Zhao, K. Nagata, X.-L. Pei, T. Eguchi, T. Kamachi, M. Ehara, T. Ozawa, M. Shionoya Nature Commun. 13, 4288-1-9 (2022).

Au-Ag nanoclusters showing intense phosphorescence in cells

Luminescent metal nanoclusters are expected to exhibit unique physical properties in the cluster structure depending on the ligand structure, metal type, number of nuclei and arrangement. In this study, carbon-centered gold-silver (CAu₆Ag₂) clusters with N-heterocyclic carbene (NHC) ligands were designed and synthesized, and it was found that these clusters emit strong phosphorescence in solution, and the contribution of NHC ligands to phosphorescence emission was revealed by theoretical calculation. The luminescence rate constant was calculated by an analysis including spin-orbit interactions, and the quantum yield was discussed in terms of the energy barrier to the minimum energy crossing point. Furthermore, the phosphorescent gold-silver clusters with long luminescence lifetime were used for cellular imaging, which revealed the pathway of uptake into the cell and selective localization to specific organelles, confirming their superior functionality, which is different from the non-selective uptake of conventional phosphine ligands.
Z. Lei, M. Endo, H. Ube, T. Shiraogawa, P. Zhao, K. Nagata, X.-L. Pei, T. Eguchi, T. Kamachi, M. Ehara, T. Ozawa, M. Shionoya Nature Commun. 13, 4288-1-9 (2022).

“Bucket brigade” of remdesivir in RNA-dependent RNA polymerase of SARS-CoV-2

　The process of remdesivir uptake into the protein (RNA-dependent RNA polymerase) that replicates the gene of SARS-CoV-2 was revealed by molecular dynamics simulations. Remdesivir has a phosphate group with a negative charge, and the binding site of RNA polymerase has Mg²⁺ ions. In addition, RNA polymerase has lysine residues, which are positively charged, in a line toward the binding site. It was found that the lysine residue attracts the phosphate group of remdesivir and passes it to the next lysine residue sequentially, transporting the drug to the binding site like a bucket brigade.

S. Tanimoto, S. G. Itoh, and H. Okumura: Biophys. J. (2021), DOI: 10.1016/j.bpj.2021.07.026, “Bucket brigade” using lysine residues in RNA-dependent RNA polymerase of SARS-CoV-2

“Bucket brigade” of remdesivir in RNA-dependent RNA polymerase of SARS-CoV-2

　The process of remdesivir uptake into the protein (RNA-dependent RNA polymerase) that replicates the gene of SARS-CoV-2 was revealed by molecular dynamics simulations. Remdesivir has a phosphate group with a negative charge, and the binding site of RNA polymerase has Mg²⁺ ions. In addition, RNA polymerase has lysine residues, which are positively charged, in a line toward the binding site. It was found that the lysine residue attracts the phosphate group of remdesivir and passes it to the next lysine residue sequentially, transporting the drug to the binding site like a bucket brigade.

S. Tanimoto, S. G. Itoh, and H. Okumura: Biophys. J. (2021), DOI: 10.1016/j.bpj.2021.07.026, “Bucket brigade” using lysine residues in RNA-dependent RNA polymerase of SARS-CoV-2

Scientific results using RCCS

Update Information

2026-07-08	Foundry 0.2.0 available
Foundry 0.2.0 (RFdiffusion3, RFdiffusion3NA, RosettaFold3, ProteinMPNN) is now available. Job sample can be found in /apl/foundry/0.2.0/samples. Read more about Foundry 0.2.0 available

2026-06-30	Amber26 (PMEMD26, AmberTools26) installed
Amber26 (PMEMD26, AmberTools26) was installed. PMEMD and AmberTools are installed as separated packages. Read more about Amber26 (PMEMD26, AmberTools26) installed

2026-06-26	CENSO 3.0.8 installed
CENSO 3.0.8 was installed. Job sample is available in /apl/censo/3.0.8/samples. Please check https://ccportal.ims.ac.jp/en/node/4006 for installation procedure. Read more about CENSO 3.0.8 installed

2026-06-24	Molpro 2026.1.0 installed
Molpro 2026.1.0 was installed. Job sample is available in /apl/molpro/2026.1.0/samples. Read more about Molpro 2026.1.0 installed

2026-06-24	AlphaFold 3.0.3 installed
AlphaFold 3.0.3 was installed. Job sample can be found in /apl/alphafold/samples/3.0.3. Users are required to request and download the AlphaFold3 model parameters by themselves. Read more about AlphaFold 3.0.3 installed

2026-06-10	Gromacs 2026.2 installed
Gromacs 2026.2 was installed. CPU and GPU versions are available. Read more about Gromacs 2026.2 installed

2026-05-12	LAMMPS 22Jul2025 Update 4 installed
LAMMPS 22Jul2025 Update 4 was installed. CPU versions (GCC, Intel) and GPU version (CUDA) are available. Read more about LAMMPS 22Jul2025 Update 4 installed

2026-04-27	NOUS scheduled maintenance: May 14 (Thu) 9:00-17:00
Due to maintenance, NOUS will be unavailable from 9:00 to 17:00 on May 14 (Thu). New applications and usage reports cannot be submitted during the maintenance. Read more about NOUS scheduled maintenance: May 14 (Thu) 9:00-17:00

2026-04-24	AlphaFold 3.0.2 installed
AlphaFold 3.0.2 was installed. Sample job is available in /apl/alphafold/samples/3.0.2. Read more about AlphaFold 3.0.2 installed

2026-04-22	CP2K 2026.1 installed
CP2K 2026.1 was installed under /apl/cp2k/2026.1. Sample files are available at /apl/cp2k/2026.1/samples. Installation details can be found at this page. Read more about CP2K 2026.1 installed