Design and implementation of a hybrid cloud system for large-scale human genomic research.

In the field of genomic medical research, the amount of large-scale information continues to increase due to advances in measurement technologies, such as high-performance sequencing and spatial omics, as well as the progress made in genomic cohort studies involving more than one million individuals. Therefore, researchers require more computational resources to analyze this information. Here, we introduce a hybrid cloud system consisting of an on-premise supercomputer, science cloud, and public cloud at the Kyoto University Center for Genomic Medicine in Japan as a solution. This system can flexibly handle various heterogeneous computational resource-demanding bioinformatics tools while scaling the computational capacity. In the hybrid cloud system, we demonstrate the way to properly perform joint genotyping of whole-genome sequencing data for a large population of 11,238, which can be a bottleneck in sequencing data analysis. This system can be one of the reference implementations when dealing with large amounts of genomic medical data in research centers and organizations.

Parallel Fock matrix construction with distributed shared memory model for the FMO-MO method.

A parallel Fock matrix construction program for FMO-MO method has been developed with the distributed shared memory model. To construct a large-sized Fock matrix during FMO-MO calculations, a distributed parallel algorithm was designed to make full use of local memory to reduce communication, and was implemented on the Global Array toolkit. A benchmark calculation for a small system indicates that the parallelization efficiency of the matrix construction portion is as high as 93% at 1,024 processors. A large FMO-MO application on the epidermal growth factor receptor (EGFR) protein (17,246 atoms and 96,234 basis functions) was also carried out at the HF/6-31G level of theory, with the frontier orbitals being extracted by a Sakurai-Sugiura eigensolver. It takes 11.3 h for the FMO calculation, 49.1 h for the Fock matrix construction, and 10 min to extract 94 eigen-components on a PC cluster system using 256 processors.

Fragment molecular orbital study of the electronic excitations in the photosynthetic reaction center of Blastochloris viridis.

All electron calculations were performed on the photosynthetic reaction center of Blastochloris viridis, using the fragment molecular orbital (FMO) method. The protein complex of 20,581 atoms and 77,754 electrons was divided into 1398 fragments, and the two-body expansion of FMO/6-31G* was applied to calculate the ground state. The excited electronic states of the embedded electron transfer system were separately calculated by the configuration interaction singles approach with the multilayer FMO method. Despite the structural symmetry of the system, asymmetric excitation energies were observed, especially on the bacteriopheophytin molecules. The asymmetry was attributed to electrostatic interaction with the surrounding proteins, in which the cytoplasmic side plays a major role.

Parallel Fock matrix construction program for molecular orbital calculation--specific computer with a hierarchical network.

A parallel Fock matrix construction program for a hierarchical network has been developed on the molecular orbital calculation-specific EHPC system. To obtain high parallelization efficiency on the hierarchical network system, a multilevel dynamic load-balancing scheme was adopted, which provides equal load balance and localization of communications on a tree-structured hierarchical network. The parallelized Fock matrix construction routine was implemented into a GAMESS program on the EHPC system, which has a tree-structured hierarchical network. Benchmark results on a 63-processor system showed high parallelization efficiency even on the tree-structured hierarchical network.

DNA and estrogen receptor interaction revealed by fragment molecular orbital calculations.

Molecular orbital calculations of the complex between DNA-ERE (estrogen response element) and ER (estrogen receptor)-DBD (DNA-binding domain) were performed using the fragment molecular orbital (FMO) method, which enables large-scale MO (molecular orbital) calculations by reducing the computational cost and by significantly increasing efficiency for parallel computation. Such a large system, which contains 3354 atoms, is impractical via conventional MO methods due to the immense computational cost. Details of the interaction between DNA-ERE and ER-DBD were revealed in this study as follows by using the FMO calculations to analyze the interfragment interaction energies (IFIEs) and the electrostatic potentials (ESPs). An area with a high positive ESP is identified on the DNA-binding side of ER-DBD and is the main driving force behind access to the DNA. The position of the ER-DBD monomer can be fixed on a phosphate group of DNA-ERE by the strong electrostatic interactions, whereas the rotation cannot be fixed. In contrast, both the position and rotation of the ER-DBD dimer can be fixed and can therefore form the stable (ER-DBD)2...DNA-ERE complex. Dimerization of the ER-DBD monomers, each of which have a charge of +5 , is mainly due to large attractive interaction energies of the second Zn fragments. The base pairs in the consensus sequence of DNA-ERE interact only with the recognition helix located in the major groove due to the large shielding effect of the phosphate groups of DNA. The recognition helix has weaker interactions with the base pairs than the electrostatic interactions with the phosphate groups. Thus, the DNA-binding machinery of the ER-DBD dimer, which can secure the recognition helix in the major groove of DNA, is crucial for interactions between the recognition helix and base pairs.

PEACH 4 with ABINIT-MP: a general platform for classical and quantum simulations of biological molecules.

A program package for molecular simulations of biological molecules was developed. The package, "PEACH version 4 with ABINIT-MP version 20021029," was constructed by incorporating ABINIT-MP, a program for the fragment molecular orbital (FMO) method [Chem.Phys. Lett. 313 (1999) 701], into PEACH, a program package for classical molecular dynamics simulations (MD). A few capabilities of the package were demonstrated. First, high parallel efficiency of FMO was demonstrated in a single point calculation of a protein. Second,FMO-MD simulations [Chem. Phys. Lett. 372 (2003) 342] of a peptide were performed with and without explicit solvent, and the simulations showed the influence of the solvent on the electronic state of the peptide.