Ab initio O(N) elongation-counterpoise method for BSSE-corrected interaction energy analyses in biosystems

Local electronic structure analysis by ab initio elongation method: A benchmark using DNA block polymers

The ab initio elongation (ELG) method based on a polymerization concept is a feasible way to perform linear-scaling electronic structure calculations for huge aperiodic molecules while maintaining computational accuracy. In the method, the electronic structures are sequentially elongated by repeating (1) the conversion of canonical molecular orbitals (CMOs) to region-localized MOs (RLMOs), that is, active RLMOs localized onto a region close to an attacking monomer or frozen RLMOs localized onto the remaining region, and the subsequent (2) partial self-consistent-field calculations for an interaction space composed of the active RLMOs and the attacking monomer. For each ELG process, one can obtain local CMOs for the interaction space and the corresponding local orbital energies. Local site information, such as the local highest-occupied/lowest-unoccupied MOs, can be acquired with linear-scaling efficiency by correctly including electronic effects from the frozen region. In this study, we performed a local electronic structure analysis using the ELG method for various DNA block polymers with different sequential patterns. This benchmark aimed to confirm the effectiveness of the method toward the efficient detection of a singular local electronic structure in unknown systems as a future practical application. We discussed the high-throughput efficiency of our method and proposed a strategy to detect singular electronic structures by combining with a machine learning technique.

On the Stabilization of Carbynes Encapsulated in Penta-Graphene Nanotubes: a DFT Study

We carried out density functional theory calculations to investigate the electronic and structural properties of linear carbon chains (carbynes) encapsulated in armchair and zigzag penta-graphene (PGNT) nanotubes. Results showed that PGNTs-wrapped carbyne can present negative formation energies that tend to stabilize that encapsulated carbon chains. These chains were stabilized in their cumulene and polyyne forms, with slight dependence on tube diameter. As a general trend, the PGNT band structures are kept almost unchanged upon carbyne encapsulation. This finding indicates weak orbital interactions between the PGNT and the carbyne. No net charge was found in chains encapsulated on zigzag PGNTs. Schematic representation of a carbyne encapsulated in a pentagraphene nanotube.

Ab initio multi-level layered elongation method and its application to local interaction analysis between DNA bulge and ligand molecules

A multi-level layered elongation method was developed for efficiently analyzing the electronic states of local structures in large bio/nano-systems at the full ab initio level of theory. The original elongation method developed during the last three decades in our group has focused on the system in one direction from one terminal to the other terminal to sequentially construct the electronic states of a polymer, called a theoretical synthesis of polymers. In this study, an important region termed the central (C) part is targeted in a large polymer and the remainder are terminal (T) parts. The electronic structures along with polymer elongation are calculated repeatedly from both end T parts to the C central part at the same time. The important C part is treated with large basis sets (high level) and the other regions are treated with small basis sets (low level) in the ab initio theoretical framework. The electronic structures besides the C part can be reused for other systems with different structures at the C part, which renders the method computationally efficient. This multi-level layered elongation method was applied to the investigation on DNA single bulge recognition of small molecules (ligands). The reliability and validity of our approach were examined in comparison with the results obtained by direct calculations using a conventional quantum chemical method for the entire system. Furthermore, stabilization energies by the formation of the complex of bulge DNA and a ligand were estimated with basis set superposition error corrections incorporated into the elongation method.

Absorption spectra of Müller cell intermediate filaments: Experimental results and theoretical models

Experimental spectra of Müller cell (MC) intermediate filaments (IFs) isolated from porcine retina are reported in this work. The absorption spectra recorded at different MC IF concentrations were used to estimate their absorption cross-sections at different wavelengths. The average absorption cross-section of a single MC IF was ca. (0.97 … 2.01) × 10^-10 cm² in the 650-445 nm spectral range. To interpret these experimental absorption spectra, we made ab initio calculations of the optical spectra of α-helix polypeptides, and also used a simplified theoretical approach that modeled an IF by a conductive wire. The energy spectra of the refractive index, extinction coefficient (absorption cross-section), energy loss and reflectivity functions for different photon polarizations, with strong anisotropy with respect to the system axis, were calculated ab initio for polyglycine α-helix molecule containing 1000 glycine residues. Strong anisotropy of these parameters was explained by photons interacting with different electronic transitions. Note that similarly strong anisotropy was also obtained for the optical absorption cross-sections in the simplified model. Both modeling approaches were used for calculating the absorption cross-sections of interest. As a result, the absorption cross-section for photons propagating axially along MC IFs was much larger than their geometrical cross-section. The latter result was explained taking into account the density of electronic states, with numerous electrons contributing to the transition intensity at a given energy. We found that the simple conductive wire model describes the MC IF absorption spectrum better than the ab initio spectra. The latter conclusion was explained by the limitations of ab initio analysis, which only took into account one α-helix with 1000 aminoacids, whereas each porcine Müller cell IF is assembled of thousands of protein molecules, reaching the total length of ca. 100 μm. The presently reported results contribute to the understanding of the quantum mechanism of high-contrast vision of vertebrate eyes.

Subsystem Analysis for the Fragment Molecular Orbital Method and Its Application to Protein-Ligand Binding in Solution

A subsystem analysis is derived incorporating interfragment interactions into the fragment properties, such as energies or charges. The relative stabilities of three alanine isomers, the α-helix, the β-turn, and the extended form are studied and the differences in fragment properties are elucidated. The analysis is further elaborated for studies of binding energies. The binding of the Trp-cage protein (PDB: 1L2Y ) to two ligands is studied in detail. Binding energies defined for each fragment can be used as a convenient descriptor for analyzing contributions to binding in solution.