Search PubMed⌕ Search

Biomedical subjects

So Hirata

Publications and source records attributed to So Hirata.

16 recordsLinked to original sources

Modeling meningioma in vitro in the omics era.

Meningioma biology has been substantially clarified by recent omics-based studies, which have identified recurrent mutations, copy-number alterations, and distinct molecular subgroups. However, although these approaches have provided a valuable framework, they are inherently limited in their ability to establish direct causal relationships. The mechanistic studies are therefore indispensable for translating these molecular observations into biological understanding. Nevertheless, the mechanistic literature has often evolved in a fragmented manner, with individual pathways and model systems studied in relative isolation from the broader multi-omic landscape. In this review, we synthesize these complementary bodies of work into an integrated framework and outline a clear roadmap for future studies. We first review the historical development of established meningioma cell lines, their current molecular characterization, and the recent emergence of 3D models and organoids. Intrinsic challenges in modeling meningioma in vitro are discussed, including the difficulty of establishing immortalized cell lines from predominantly benign tumors, genetic alterations introduced during immortalization, and drift under culture conditions that differ substantially from those of the parental tumors. Next, insights from functional studies centered on these models are integrated within the molecular framework established by large-scale omics analyses. To avoid fragmentation and overemphasis on isolated findings, prior studies are organized into six categories based on major signaling pathways: Hippo, PI3K/Akt/mTOR, MAPK, Wnt/β-catenin, FOXM1, and Notch. Finally, lessons from other cancer models, including experimental approaches to chromosome-scale genomic disturbances, are considered to provide a more integrated view of meningioma biology and to highlight directions for future research.

Meningioma↗

Chromosome-scale genome remodeling in tumor evolution: Copy number alterations and structural variants as two sides of the same coin.

Chromosome-scale genomic rearrangements are a dominant force in tumor evolution. Copy-number alterations (CNAs) and structural variants (SVs) constitute two complementary axes of this process. Although detection technologies now deliver near-comprehensive catalogs, technical resolution has outpaced conceptual integration. In this review, we frame CNAs and SVs as inextricable facets of chromosomal aberrations. They reshape cancer genomes through altered gene dosage and three-dimensional regulatory rewiring. CNAs quantify the gene-dosage imbalance, yet arise through mechanistically distinct routes. Segmental CNAs typically require chromosomal breakage, and therefore often coincide with SV junctions. By contrast, whole-chromosome aneuploidy and whole-genome doubling (WGD) primarily reflect mitotic or cytokinetic failure and can occur without local breakpoints, while nevertheless reshaping the karyotypic landscape and seeding subsequent structural complexity. SVs, in turn, range from unbalanced events that alter copy number to ostensibly balanced exchanges that predominantly rewire regulatory architecture. Despite their diverse and sometimes catastrophic architectures, SVs are ultimately rooted in double-strand break formation and error-prone resolution. By integrating CNAs and SVs within a unified mechanistic and functional framework, we aim to convert catalogs into concepts and distill the organizing principles that govern tumor genome evolution.

Humans↗

High-order electron-correlation methods with scalar relativistic and spin-orbit corrections.

An assortment of computer-generated, parallel-executable programs of ab initio electron-correlation methods has been fitted with the ability to use relativistic reference wave functions. This has been done on the basis of scalar relativistic and spin-orbit effective potentials and by allowing the computer-generated programs to handle complex-valued, spinless orbitals determined by these potentials. The electron-correlation methods that benefit from this extension are high-order coupled-cluster methods (up to quadruple excitation operators) for closed- and open-shell species, coupled-cluster methods for excited and ionized states (up to quadruples), second-order perturbation corrections to coupled-cluster methods (up to triples), high-order perturbation corrections to configuration-interaction singles, and active-space (multireference) coupled-cluster methods for the ground, excited, and ionized states (up to active-space quadruples). A subset of these methods is used jointly such that the dynamical correlation energies and scalar relativistic effects are computed by a lower-order electron-correlation method with more extensive basis sets and all-electron relativistic treatment, whereas the nondynamical correlation energies and spin-orbit effects are treated by a higher-order electron-correlation method with smaller basis sets and relativistic effective potentials. The authors demonstrate the utility and efficiency of this composite scheme in chemical simulation wherein the consideration of spin-orbit effects is essential: ionization energies of rare gases, spectroscopic constants of protonated rare gases, and photoelectron spectra of hydrogen halides.

Journal Article↗

Higher-order equation-of-motion coupled-cluster methods for ionization processes.

Compact algebraic equations defining the equation-of-motion coupled-cluster (EOM-CC) methods for ionization potentials (IP-EOM-CC) have been derived and computer implemented by virtue of a symbolic algebra system largely automating these processes. Models with connected cluster excitation operators truncated after double, triple, or quadruple level and with linear ionization operators truncated after two-hole-one-particle (2h1p), three-hole-two-particle (3h2p), or four-hole-three-particle (4h3p) level (abbreviated as IP-EOM-CCSD, CCSDT, and CCSDTQ, respectively) have been realized into parallel algorithms taking advantage of spin, spatial, and permutation symmetries with optimal size dependence of the computational costs. They are based on spin-orbital formalisms and can describe both alpha and beta ionizations from open-shell (doublet, triplet, etc.) reference states into ionized states with various spin magnetic quantum numbers. The application of these methods to Koopmans and satellite ionizations of N2 and CO (with the ambiguity due to finite basis sets eliminated by extrapolation) has shown that IP-EOM-CCSD frequently accounts for orbital relaxation inadequately and displays errors exceeding a couple of eV. However, these errors can be systematically reduced to tenths or even hundredths of an eV by IP-EOM-CCSDT or CCSDTQ. Comparison of spectroscopic parameters of the FH+ and NH+ radicals between IP-EOM-CC and experiments has also underscored the importance of higher-order IP-EOM-CC treatments. For instance, the harmonic frequencies of the A 2Sigma- state of NH+ are predicted to be 1285, 1723, and 1705 cm(-1) by IP-EOM-CCSD, CCSDT, and CCSDTQ, respectively, as compared to the observed value of 1707 cm(-1). The small adiabatic energy separation (observed 0.04 eV) between the X 2Pi and a 4Sigma- states of NH+ also requires IP-EOM-CCSDTQ for a quantitative prediction (0.06 eV) when the a 4Sigma- state has the low-spin magnetic quantum number (s(z) = 1/2). When the state with s(z) = 3/2 is sought, the energy separations converge much more rapidly with the IP-EOM-CCSD value (0.03 eV) already being close to the observed (0.04 eV).

Journal Article↗

Franck-Condon factors based on anharmonic vibrational wave functions of polyatomic molecules.

Franck-Condon (FC) integrals of polyatomic molecules are computed on the basis of vibrational self-consistent-field (VSCF) or configuration-interaction (VCI) calculations capable of including vibrational anharmonicity to any desired extent (within certain molecular size limits). The anharmonic vibrational wave functions of the initial and final states are expanded unambiguously by harmonic oscillator basis functions of normal coordinates of the respective electronic states. The anharmonic FC integrals are then obtained as linear combinations of harmonic counterparts, which can, in turn, be evaluated by established techniques taking account of the Duschinsky rotations, geometry displacements, and frequency changes. Alternatively, anharmonic wave functions of both states are expanded by basis functions of just one electronic state, permitting the FC integral to be evaluated directly by the Gauss-Hermite quadrature used in the VSCF and VCI steps [Bowman et al., Mol. Phys. 104, 33 (2006)]. These methods in conjunction with the VCI and coupled-cluster with singles, doubles, and perturbative triples [CCSD(T)] method have predicted the peak positions and intensities of the vibrational manifold in the X 2B1 photoelectron band of H2O with quantitative accuracy. It has revealed that two weakly visible peaks are the result of intensity borrowing from nearby states through anharmonic couplings, an effect explained qualitatively by VSCF and quantitatively by VCI, but not by the harmonic approximation. The X 2B2 photoelectron band of H2CO is less accurately reproduced by this method, likely because of the inability of CCSD(T)/cc-pVTZ to describe the potential energy surface of open-shell H2CO+ with the same high accuracy as in H2O+.

Journal Article↗

Advances in methods and algorithms in a modern quantum chemistry program package.

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.

Algorithms↗

Active-space coupled-cluster methods through connected quadruple excitations.

Coupled-cluster methods that include just a subset of all connected triple, quadruple, or both excitation amplitudes, according to the ansatz of and Adamowicz co-workers [Int. Rev. Phys. Chem. 12, 339 (1993); J. Chem. Phys. 99, 1875 (1993); 100, 5792 (1994)] and Piecuch et al. [J. Chem. Phys. 110, 6103 (1999)], have been implemented into parallel execution programs. They are applicable to closed- and open-shell species and they take advantage of real Abelian point-group symmetry. A symbol manipulation program has been invoked to automate the implementation. These methods have been applied to the singlet-triplet separations of five triatomic hydrides (CH2, NH2+, SiH2, PH2+, and AsH2+) with consideration of scalar relativistic effects. They have been shown to be remarkably effective with errors arising from the use of a very small subset of higher-order excitations being no more than a few tenths of 1 kcal/mol.

Journal Article↗

Vibrational and electronic absorption spectroscopy of dibenzo[b,def]chrysene and its ions.

The vibrational and electronic absorption spectra of dibenzo[b,def]chrysene (DBC) and its ions in argon matrixes have been recorded. Assignment of the observed infrared (IR) bands has been made by comparison with the density functional theory (DFT) computations of harmonic vibrational frequencies (with 6-31G(d,p) or 6-311+G(d,p) basis sets). Extensive time-dependent (TD) DFT calculations of vertical excitation energies have aided in the assignment of the experimental electronic absorption transitions. In general, the theoretical predictions are in good agreement with the observed ultraviolet and visible bands. By correlating IR and UV-visible band intensities (after UV photolysis), it has been shown that both DBC cations and anions are formed. The IR band intensity distributions of the DBC ions differ markedly from neutral DBC. A synthetic spectrum composed of neutral, cationic, and anionic DBC contributions compares reasonably well with the interstellar features of the "unidentified infrared" (UIR) bands from the reflection nebula NGC 7023. Finally, it is shown that the electronic absorption bands of the DBC ions lie in close proximity to several of the diffuse interstellar visible absorption bands (DIBs).

Chrysenes↗

Active-space coupled-cluster study of electronic states of Be3.

The active-space coupled-cluster (CC) and equation-of-motion (EOM) CC methods with all single and double excitations with triple excitations defined via active orbitals (CCSDt, EOMCCSDt), as implemented with TENSOR CONTRACTION ENGINE, are applied to the challenging Be3 system, which is characterized by a large number of low-lying excited states dominated by two-electron transitions and significant high-order correlation effects in the ground electronic state. It is demonstrated that the CCSDt and EOMCCSDt methods provide an excellent description of complicated electronic quasidegeneracies present in the Be3 cluster. Different strategies for defining triple excitations within the CCSDtEOMCCSDt approach are discussed.

Journal Article↗

Third- and fourth-order perturbation corrections to excitation energies from configuration interaction singles.

Complete third-order and partial fourth-order Rayleigh-Schrodinger perturbation corrections to excitation energies from configuration interaction singles (CIS) have been derived and termed CIS(3) and CIS(4)(P). They have been implemented by the automated system TENSOR CONTRACTION ENGINE into parallel-execution programs taking advantage of spin, spatial, and index permutation symmetries and applicable to closed- and open-shell molecules. The consistent use of factorization, first introduced by Head-Gordon et al. in the second-order correction to CIS denoted CIS(D), has reduced the computational cost of CIS(3) and CIS(4)(P) from O(n(8)) and O(n(6)) to O(n(6)) and O(n(5)), respectively, with n being the number of orbitals. It has also guaranteed the size extensivity of excited-state energies of these methods, which are in turn the sum of size-intensive excitation energies and the ground-state energies from the standard Møller-Plesset perturbation theory at the respective orders. The series CIS(D), CIS(3), and CIS(4)(P) are usually monotonically convergent at values close to the accurate results predicted by coupled-cluster singles and doubles (CCSD) with a small fraction of computational costs of CCSD for predominantly singly excited states characterized by a 90%-100% overlap between the CIS and CCSD wave functions. When the overlap is smaller, the perturbation theory is incapable of adequately accounting for the mixing of the CIS states through higher-than-singles sectors of the Hamiltonian matrix, resulting in wildly oscillating series with often very large errors in CIS(4)(P). Hence, CIS(3) and CIS(4)(P) have a rather small radius of convergence and a limited range of applicability, but within that range they can be an inexpensive alternative to CCSD.

Journal Article↗

The exchange-correlation potential in ab initio density functional theory.

From coupled-cluster theory and many-body perturbation theory we derive the local exchange-correlation potential of density functional theory in an orbital dependent form. We show the relationship between the coupled-cluster approach and density functional theory, and connections and comparisons with our previous second-order correlation potential [OEP-MBPT(2) (OEP-optimized effective potential)] [I. Grabowski, S. Hirata, S. Ivanov, and R. J. Bartlett, J. Chem. Phys. 116, 4415 (2002)]. Starting from a general theoretical framework based on the density condition in Kohn-Sham theory, we define a rigorous exchange-correlation functional, potential and orbitals. Specifying initially to second-order terms, we show that our ab initio correlation potential provides the correct shape compared to those from reference quantum Monte Carlo calculations, and we demonstrate the superiority of using Fock matrix elements or more general infinite-order semicanonical transformations. This enables us to introduce a method that is guaranteed to converge to the right answer in the correlation and basis set limit, just as does ab initio wave function theory. We also demonstrate that the energies obtained from this generalized second-order method [OEP-MBPT2-f] and [OEP-MBPT2-sc] are often of coupled-cluster accuracy and substantially better than ordinary Hartree-Fock based second-order MBPT=MP2.

Journal Article↗

Combined coupled-cluster and many-body perturbation theories.

Various approximations combining coupled-cluster (CC) and many-body perturbation theories have been derived and implemented into the parallel execution programs that take into account the spin, spatial (real Abelian), and permutation symmetries and that are applicable to closed- and open-shell molecules. The implemented models range from the CCSD(T), CCSD[T], CCSD(2)(T), CCSD(2)(TQ), and CCSDT(2)(Q) methods to the completely renormalized (CR) CCSD(T) and CCSD[T] approaches, where CCSD (CCSDT) stands for the CC method with connected single and double (single, double, and triple) cluster operators, and subscripted or parenthesized 2, T, and Q indicate the perturbation order or the excitation ranks of the cluster operators included in the corrections. The derivation and computer implementation have been automated by the algebraic and symbolic manipulation program TENSOR CONTRACTION ENGINE (TCE). The TCE-synthesized subroutines generate the tensors with the highest excitation rank in a blockwise manner so that they need not be stored in their entirety, while enabling the efficient reuse of other precalculated intermediate tensors defined by prioritizing the memory optimization as well as operation minimization. Consequently, the overall storage requirements for the corrections due to connected triple and quadruple cluster operators scale as O(n(4)) and O(n(6)), respectively (n being a measure of the system size). For systems with modest multireference character of their wave functions, we found that the order of accuracy is CCSD<CR-CCSD(T) approximately CCSD(2)(T) approximately CCSD(T)<CCSDT approximately CCSD(2)(TQ)<CCSDT(2)(Q), whereas CR-CCSD(T) is more effective in cases of larger quasidegeneracy. The operation costs of the TCE-generated CCSD(2)(TQ) and CCSDT(2)(Q) codes scale as rather steep O(n(9)), while the TCE-generated CCSD(T), CCSD(2)(T), and CR-CCSD(T) codes are near operation minimum [a noniterative O(n(7))]. The perturbative correction part of the CCSD(T)/cc-pVDZ calculations for azulene exhibited a 45-fold speedup upon a 64-fold increase in the number of processors from 8 to 512.

Journal Article↗

Higher-order equation-of-motion coupled-cluster methods.

The equation-of-motion coupled-cluster (EOM-CC) methods truncated after double, triple, or quadruple cluster and linear excitation operators (EOM-CCSD, EOM-CCSDT, and EOM-CCSDTQ) have been derived and implemented into parallel execution programs. They compute excitation energies, excited-state dipole moments, and transition moments of closed- and open-shell systems, taking advantage of spin, spatial (real Abelian), and permutation symmetries simultaneously and fully (within the spin-orbital formalisms). The related Lambda equation solvers for coupled-cluster (CC) methods through and up to connected quadruple excitation (CCSD, CCSDT, and CCSDTQ) have also been developed. These developments have been achieved, by virtue of the algebraic and symbolic manipulation program that automated the formula derivation and implementation altogether. The EOM-CC methods and CC Lambda equations introduce a class of second quantized ansatz with a de-excitation operator (Y), a number of excitation operators (X), and a physical (e.g., the Hamiltonian) operator (A), leading to the tensor contraction expressions that can be performed in the order of ((...((yx)x)...)x)a or ((...((ax)x)...)x)y at the minimal peak operation cost, where x, y, and a are basis-set representations (i.e., tensors) of the respective operators X, Y, and A. Any intermediate tensor resulting from either contraction order is shown to have at most six groups of permutable indices, and this knowledge is used to guide the computer-synthesized programs to fully exploit the permutation symmetry of any tensor to minimize the arithmetic and memory costs.

Journal Article↗

Third-order Douglas-Kroll relativistic coupled-cluster theory through connected single, double, triple, and quadruple substitutions: applications to diatomic and triatomic hydrides.

Coupled-cluster methods including through and up to the connected single, double, triple, and quadruple substitutions have been derived and implemented automatically for sequential and parallel executions by an algebraic and symbolic manipulation program TCE (TENSOR CONTRACTION ENGINE) for use in conjunction with a one-component third-order Douglas-Kroll approximation for relativistic corrections. A combination of the converging electron-correlation methods, the accurate relativistic reference wave functions, and the use of systematic basis sets tailored to the relativistic approximation has been shown to predict the experimental singlet-triplet separations within 0.02 eV (0.5 kcal/mol) for five triatomic hydrides (CH2, NH2+, SiH2, PH2+, and AsH2+), the experimental bond lengths (re or r0) within 0.002 angstroms, rotational constants (Be or B0) within 0.02 cm(-1), vibration-rotation constants (alphae) within 0.01 cm(-1), centrifugal distortion constants (De) within 2%, harmonic vibration frequencies (omegae) within 8 cm(-1) (0.4%), anharmonic vibrational constants (xomegae) within 2 cm(-1), and dissociation energies (D0(0)) within 0.02 eV (0.4 kcal/mol) for twenty diatomic hydrides (BH, CH, NH, OH, FH, AlH, SiH, PH, SH, ClH, GaH, GeH, AsH, SeH, BrH, InH, SnH, SbH, TeH, and IH) containing main-group elements across the second through fifth rows of the periodic table. In these calculations, spin-orbit effects on dissociation energies, which were assumed to be additive, were estimated from the measured spin-orbit coupling constants of atoms and diatomic molecules, and an electronic energy in the complete-basis-set, complete-electron-correlation limit has been extrapolated in two ways to verify the robustness of the results: One assuming Gaussian-exponential dependence of total energies on double through quadruple zeta basis sets and the other assuming n(-3) dependence of correlation energies on double through quintuple zeta basis sets.

Journal Article↗

Coupled-cluster singles and doubles for extended systems.

Coupled-cluster theory with connected single and double excitation operators (CCSD) and related approximations, such as linearized CCSD, quadratic configuration interaction with single and double excitation operators, coupled-cluster with connected double excitation operator (CCD), linearized CCD, approximate CCD, and second- and third-order many-body perturbation theories, are formulated and implemented for infinitely extended one-dimensional systems (polymers), on the basis of the periodic boundary conditions and distance-based screening of integrals, density matrix elements, and excitation amplitudes. The variation of correlation energies with the truncation radii of short- and long-range lattice sums and with the number of wave vector sampling points in the first Brillouin zone is examined for polyethylene, polyacetylene, and polyyne, and is shown to be a function of the degree of pi-electron conjugation or the fundamental band gaps. The t2 and t1 amplitudes in the atomic orbital (AO) basis are obtained by first computing the t amplitudes in the Bloch-orbital basis and subsequently back-transforming them into the AO basis. The plot of these AO-based t amplitudes as a function of unit cells also indicates that the t2 amplitudes of polyacetylene and polyyne exhibit appreciably slower decay than those of polyethylene, although the asymptotic decay behavior is invariably 1/r3. The AO-based t1 amplitudes appear to correlate strongly with the electronic structure, and they decay seemingly exponentially for polyethylene whereas they stay at a constant magnitude across the seventh nearest neighbors of polyacetylene and polyyne, which attests to far reaching effects of nondynamical electron correlation mediated by orbital rotation. Nonetheless, the unit cell contributions to the correlation energies taper below 10(-6) hartree after 15 A for all three polymers. The basis set dependence of the decay behavior of t2 amplitudes is also examined for linear hydrogen fluoride polymer (HF)infinity and linear beryllium polymer (Be)infinity employing the STO-3G, 6-31G, and 6-31G* basis sets, and proves to be rather small.

Journal Article↗