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-04-07	DFTB+ 25.1 is now available
This is the latest version of DFTB+ currently available. Both MPI and non-MPI versions are available. Sample files can be found in /apl/dftb+/25.1/mpi/samples/ and /apl/dftb+/25.1/non-mpi/samples/. Read more about DFTB+ 25.1 is now available

2026-03-31	OpenMolcas v26.02 installed
OpenMolcas v26.02 was installed in /apl/openmolcas/v26.02 directory. Read more about OpenMolcas v26.02 installed

2026-03-25	OpenMM 8.5.0 installed
OpenMM 8.5.0 was installed. Job samples are available in /apl/openmm/8.5.0/samples. Read more about OpenMM 8.5.0 installed

2026-03-16	NWChem 7.3.1 installed.
NWChem 7.3.1 was installed. GCC-15 and Intel Classic Compiler 2023 versions are available. Job samples are available in /apl/nwchem/7.3.1/samples and /apl/nwchem/7.3.1-intel/samples. For installation details, please check following links. Read more about NWChem 7.3.1 installed.

2026-03-06	GENESIS 2.1.6.1 installed.
GENESIS 2.1.6.1 was installed. CPU and CPU+GPU versions are available. Read more about GENESIS 2.1.6.1 installed.

2026-02-18	Molpro 2025.4.1 is Now Available
Molpro 2025.4.1, the latest release of Molpro, is now available. Read more about Molpro 2025.4.1 is Now Available

2026-02-12	Siesta 5.4.2 installed
Siesta 5.4.2 was installed. As with the previous build (5.4.1), this build also employs gcc and openmpi. Job sample is available in /apl/siesta/5.4.2/samples. Read more about Siesta 5.4.2 installed

2026-02-05	LAMMPS 22Jul2025 Update 3 installed
LAMMPS 22Jul2025 Update 3 was installed. Read more about LAMMPS 22Jul2025 Update 3 installed

2026-02-05	Boltz 2.2.1, BoltzGen 0.2.0, BindCraft 1.5.2 installed
Boltz 2.2.1, BoltzGen 0.2.0, and BindCraft 1.5.2 were installed. Read more about Boltz 2.2.1, BoltzGen 0.2.0, BindCraft 1.5.2 installed

2026-01-14	CENSO 3.0.5 installed
CENSO 3.0.5 was installed. Job samples are available in /apl/censo/3.0.5/samples. Read more about CENSO 3.0.5 installed