A perturbative extension to optimized coordinate vibrational self-consistent field (oc-VSCF) is proposed based on the quasi-degenerate perturbation theory (QDPT). A scheme to construct the degenerate space (P space) is developed, which incorporates degenerate configurations and alleviates the divergence of perturbative expansion due to localized coordinates in oc-VSCF (e.g., local O–H stretching modes of water). An efficient configuration selection scheme is also implemented, which screens out the Hamiltonian matrix element between the P space configuration (p) and the complementary Q space configuration (q) based on a difference in their quantum numbers (λpq = ∑s|ps − qs|). It is demonstrated that the second-order vibrational QDPT based on optimized coordinates (oc-VQDPT2) smoothly converges with respect to the order of the mode coupling, and outperforms the conventional one based on normal coordinates. Furthermore, an improved, fast algorithm is developed for optimizing the coordinates. First, the minimization of the VSCF energy is conducted in a restricted parameter space, in which only a portion of pairs of coordinates is selectively transformed. A rational index is devised for this purpose, which identifies the important coordinate pairs to mix from others that may remain unchanged based on the magnitude of harmonic coupling induced by the transformation. Second, a cubic force field (CFF) is employed in place of a quartic force field, which bypasses intensive procedures that arise due to the presence of the fourth-order force constants. It is found that oc-VSCF based on CFF together with the pair selection scheme yields the coordinates similar in character to the conventional ones such that the final vibrational energy is affected very little while gaining an order of magnitude acceleration. The proposed method is applied to ethylene and trans-1,3-butadiene. An accurate, multi-resolution potential, which combines the MP2 and coupled-cluster with singles, doubles, and perturbative triples level of electronic structure theory, is generated and employed in the oc-VQDPT2 calculation to obtain the fundamental tones as well as selected overtones/combination tones coupled to the fundamentals through the Fermi resonance. The calculated frequencies of ethylene and trans-1,3-butadiene are found to be in excellent agreement with the experimental values with a mean absolute error of 8 and 9 cm−1, respectively.

Influenza virus infections are serious public health concerns throughout the world. The development of compounds with novel mechanisms of action is urgently required due to the emergence of viruses with resistance to the currently-approved anti-influenza viral drugs. We performedin silicoscreening using a structure-based drug discovery algorithm called Nagasaki University Docking Engine (NUDE), which is optimised for a GPU-based supercomputer (DEstination for Gpu Intensive MAchine; DEGIMA), by targeting influenza viral PA protein. The compounds selected by NUDE were tested for anti-influenza virus activity using a cell-based assay. The most potent compound, designated as PA-49, is a medium-sized quinolinone derivative bearing a tetrazole moiety, and it inhibited the replication of influenza virus A/WSN/33 at a half maximal inhibitory concentration of 0.47 μM. PA-49 has the ability to bind PA and its anti-influenza activity was promising against various influenza strains, including a clinical isolate of A(H1N1)pdm09 and type B viruses. The docking simulation suggested that PA-49 interrupts the PA–PB1 interface where important amino acids are mostly conserved in the virus strains tested, suggesting the strain independent utility. Because our NUDE/DEGIMA system is rapid and efficient, it may help effective drug discovery against the influenza virus and other emerging viruses.

Objective: While the aetiology of idiopathic multicentric Castleman disease (iMCD) remains unclear, the involvement of autoinflammatory mechanisms has been suggested. Herein we report a Japanese patient with iMCD with a novel heterozygous Ile729Met mutation in exon 10 of the Mediterranean fever (MEFV) gene. Methods: We performed genetic analysis via targeted next-generation sequencing analysis and Sanger sequencing and conducted molecular dynamics simulations to investigate the hydrophobic interactions around the 729th amino acid in human pyrin. Results: In February 2011, a 59-year-old man was diagnosed with IgG4-related disease at our department based on the findings of cervical and hilar lymphadenopathies, typical lung lesions and cervical lymph node biopsy. The patient was followed up without treatment, as he was asymptomatic. However, he had been experiencing prolonged fatigue and fever with high levels of CRP since June 2017. Axillary lymph node biopsy findings led to the diagnosis of iMCD. He was successfully treated with an IL-6 inhibitor and has been in remission for 12 months. Genetic analyses for hereditary autoinflammatory disease–related genes were performed, revealing a novel heterozygous Ile729Met mutation in exon 10 of the MEFV gene. We identified that this novel mutation significantly altered the local interaction of the human pyrin B30.2 domain by molecular dynamics simulation analysis and experimentally had the potential for inflammasome activation with increased inflammatory cytokines. Conclusion: The abnormal function of pyrin due to a mutation in the MEFV gene in this patient may have contributed to the development of MCD by inducing IL-6 production via inflammasome signalling.

Prions are composed solely of the pathological isoform (PrP^Sc) of the normal cellular prion protein (PrP^C). Identification of different PrP^Sc structures is crucially important for understanding prion biology because the pathogenic properties of prions are hypothesized to be encoded in the structures of PrP^Sc. However, these structures remain yet to be identified, because of the incompatibility of PrP^Sc with conventional high-resolution structural analysis methods. Previously, we reported that the region between the first and the second α-helix (H1∼H2) of PrP^C might cooperate with the more C-terminal side region for efficient interactions with PrP^Sc. From this starting point, we created a series of PrP variants with two cysteine substitutions (C;C-PrP) forming a disulfide-crosslink between H1∼H2 and the distal region of the third helix (Ctrm). We then assessed the conversion capabilities of the C;C-PrP variants in N2a cells infected with mouse-adapted scrapie prions (22L-ScN2a). Specifically, Cys substitutions at residues 165, 166, or 168 in H1∼H2 were combined with cysteine scanning along Ctrm residues 220–229. We found that C;C-PrPs are expressed normally with glycosylation patterns and subcellular localization similar to WT PrP, albeit differing in expression levels. Interestingly, some C;C-PrPs converted to protease-resistant isoforms in the 22L-ScN2a cells, but not in Fukuoka1 prion-infected cells. Crosslink patterns of convertible C;C-PrPs indicated a positional change of H1∼H2 toward Ctrm in PrP^Sc–induced conformational conversion. Given the properties of the C;C-PrPs reported here, we propose that these PrP variants may be useful tools for investigating prion strain–specific structures and structure–phenotype relationships of PrP^Sc.