Search

Abstract: We have investigated the effects of surface morphologies of base Nb electrodes on quality parameter, V m, of high critical current density (J c) Nb/AlO x /Nb junctions. The results showed that the V m depended significantly on the surface flatness parameter, root-mean-square (rms), of the sputtered Nb base double layer which consisted of a 200nm-thick Nb layer on a 100nm-thick Nb underlayer. The surface flatness was controlled by changing the Ar gas pressures between 0.1Pa and 1.5Pa during the dc sputter-deposition process. Reducing the Ar gas pressure improved the rms of the base Nb-double layer from 5nm to 1nm. When junctions with underlayers were fabricated at reduced Ar gas pressure, an improvement in V m was realized. The results indicated that smaller rms values were required for obtaining larger V m values. At an rms=1.7nm, a maximum V m of 16.5mV was obtained for junctions with a J c =14kA/cm2 and a V m of 3.0mV was realized for a junction with a J c of 52kA/cm2. [Copyright &y& Elsevier]

In all-atom (AA) molecular dynamics (MD) simulations, the rugged energy profile of the force field makes it challenging to reproduce spontaneous structural changes in biomolecules within a reasonable calculation time. Existing coarse-grained (CG) models, in which the energy profile is set to a global minimum around the initial structure, are unsuitable to explore the structural dynamics between metastable states far away from the initial structure without any bias. In this study, we developed a new hybrid potential composed of an artificial intelligence (AI) potential and minimal CG potential related to the statistical bond length and excluded volume interactions to accelerate the transition dynamics while maintaining the protein character. The AI potential is trained by energy matching using a diverse structural ensemble sampled via multicanonical (Mc) MD simulation and the corresponding AA force field energy, profile of which is smoothed by energy minimization. By applying the new methodology to chignolin and TrpCage, we showed that the AI potential can predict the AA energy with significantly high accuracy, as indicated by a correlation coefficient ( R -value) between the true and predicted energies exceeding 0.89. In addition, we successfully demonstrated that CGMD simulation based on the smoothed hybrid potential can significantly enhance the transition dynamics between various metastable states while preserving protein properties compared to those obtained with conventional CGMD and AAMD.

Background: Mixed lineage leukemia 1-rearranged (MLL1-r) acute leukemia patients respond poorly to currently available treatments and there is a need to develop more effective therapies directly disrupting the Menin‒MLL1 complex. Small-molecule-mediated inhibition of the protein‒protein interaction between Menin and MLL1 fusion proteins is a potential therapeutic strategy for patients with MLL1-r or mutated-nucleophosmin 1 (NPM1c) acute leukemia. In this study, we preclinically evaluated the new compound DS-1594a and its salts., Methods: We evaluated the preclinical efficacy of DS-1594a as well as DS-1594a·HCl (the HCl salt of DS-1594a) and DS-1594a·succinate (the succinic acid salt of DS-1594a, DS-1594b) in vitro and in vivo using acute myeloid leukemia (AML)/acute lymphoblastic leukemia (ALL) models., Results: Our results showed that MLL1-r or NPM1c human leukemic cell lines were selectively and highly sensitive to DS-1594a·HCl, with 50% growth inhibition values < 30 nM. Compared with cytrabine, the standard chemotherapy drug as AML therapy, both DS-1594a·HCl and DS-1594a·succinate mediated the eradication of potential leukemia-initiating cells by enhancing differentiation and reducing serial colony-forming potential in MLL1-r AML cells in vitro. The results were confirmed by flow cytometry, RNA sequencing, RT‒qPCR and chromatin immunoprecipitation sequencing analyses. DS-1594a·HCl and DS-1594a·succinate exhibited significant antitumor efficacy and survival benefit in MOLM-13 cell and patient-derived xenograft models of MLL1-r or NPM1c acute leukemia in vivo., Conclusion: We have generated a novel, potent, orally available small-molecule inhibitor of the Menin-MLL1 interaction, DS-1594a. Our results suggest that DS-1594a has medicinal properties distinct from those of cytarabine and that DS-1594a has the potential to be a new anticancer therapy and support oral dosing regimen for clinical studies (NCT04752163)., (© 2023. The Author(s).)

Histone acetylation is a post-translational modification of histones that is catalyzed by histone acetyltransferases (HATs) and plays an essential role in cellular processes. The HAT domain of EP300/CBP has recently emerged as a potential drug target for cancer therapy. Here, we describe the identification of the novel, highly potent, and selective EP300/CBP HAT inhibitor DS-9300. Our optimization efforts using a structure-based drug design approach based on the cocrystal structures of the EP300 HAT domain in complex with compounds 2 and 3 led to the identification of compounds possessing low-nanomolar EP300 HAT inhibitory potency and the ability to inhibit cellular acetylation of histone H3K27. Optimization of the pharmacokinetic properties in this series resulted in compounds with excellent oral systemic exposure, and once-daily oral administration of 16 (DS-9300) demonstrated potent antitumor effects in a castrated VCaP xenograft mouse model without significant body weight loss.

Changes in intracellular calcium concentrations regulate heart beats. However, the decline in the left ventricular pressure during early diastole is much sharper than that of the Ca 2+ transient, resulting in a rapid supply of blood to the left ventricle during the diastole. At the tissue level, cardiac muscles have a distinct characteristic, known as stretch activation, similar to the function of insect flight muscles. Stretch activation, which is a delayed increase in force following a rapid muscle length increase, has been thought to be related to autonomous control in these muscles. In this numerical simulation study, we introduced a molecular mechanism of stretch activation and investigated the role of this mechanism in the pumping function of the heart, using the previously developed coupling multiple-step active stiffness integration scheme for a Monte Carlo (MC) cross-bridge model and a bi-ventricular finite element model. In the MC cross-bridge model, we introduced a mechanism for trapping the myosin molecule in its post-power stroke state. We then determined the rate constants of transitions for trapping and escaping in a thermodynamically consistent manner. Based on our numerical analysis, we draw the following conclusions regarding the stretch activation mechanism: (i) the delayed force becomes larger than the original isometric force because the population of trapped myosin molecules and their average force increase after stretching; (ii) the delayed force has a duration of more than a few seconds owing to a fairly small rate constant of escape from the trapped state. For the role of stretch activation in heart pumping, we draw the following conclusions: (iii) for the regions in which the contraction force decreases earlier than the neighboring region in the end-systole phase, the trapped myosin molecules prevent further lengthening of the myocytes, which then prevents further shortening of neighboring myocytes; (iv) as a result, the contraction forces are sustained longer, resulting in a larger blood ejection, and their degeneration is synchronized., Competing Interests: TW, JO, SS, and TH are employed by UT-Heart Inc. KY and MW were employed by Fujitsu Japan, Limited. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Yoneda, Kanada, Okada, Watanabe, Sugiura, Hisada and Washio.)

Fragment molecular orbital (FMO) method is a powerful computational tool for structure-based drug design, in which protein-ligand interactions can be described by the inter-fragment interaction energy (IFIE) and its pair interaction energy decomposition analysis (PIEDA). Here, we introduced a dynamically averaged (DA) FMO-based approach in which molecular dynamics simulations were used to generate multiple protein-ligand complex structures for FMO calculations. To assess this approach, we examined the correlation between the experimental binding free energies and DA-IFIEs of six CDK2 inhibitors whose net charges are zero. The correlation between the experimental binding free energies and snapshot IFIEs for X-ray crystal structures was R 2  = 0.75. Using the DA-IFIEs, the correlation significantly improved to 0.99. When an additional CDK2 inhibitor with net charge of -1 was added, the DA FMO-based scheme with the dispersion energies still achieved R 2  = 0.99, whereas R 2 decreased to 0.32 employing all the energy terms of PIEDA., (© 2022 Wiley Periodicals LLC.)

EP300 and its paralog CBP play an important role in post-translational modification as histone acetyltransferases (HATs). EP300/CBP inhibition has been gaining attention as an anticancer treatment target in recent years. Herein, we describe the identification of a novel, highly selective EP300/CBP inhibitor, compound 11 (DS17701585), by scaffold hopping and structure-based optimization of a high-throughput screening hit 1. Compound 11 (DS17701585) shows dose-dependent inhibition of SRY-box transcription factor 2 (SOX2) mRNA expression in a human lung squamous cell carcinoma cell line LK2-xenografted mouse model., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

Compared to all-atom molecular dynamics (AA-MD) simulations, coarse-grained (CG) MD simulations can significantly reduce calculation costs. However, existing CG-MD methods are unsuitable for sampling structures that depart significantly from the initial structure without any biased force. In this study, we developed a new adaptive CG elastic network model (ENM), in which the dynamic cross-correlation coefficient based on short-time AA-MD of at most ns order is considered. By applying Bayesian optimization to search for a suitable parameter among the vast parameter space of adaptive CG-ENM, we succeeded in reducing the searching cost to approximately 10% of those for random sampling and exhaustive sampling. To evaluate the performance of adaptive CG-ENM, we applied the new methodology to adenylate kinase (ADK) and glutamine binding protein (GBP) in the apo state. The results showed that the structural ensembles explored by adaptive CG-ENM could be considerably more diverse than those by conventional ENMs with enhanced sampling such as temperature replica exchange MD and long-time AA-MD of 1 μs. In particular, some of the structures sampled by adaptive ENM are relatively close to the holo-type structures of ADK and GBP. Furthermore, as a challenging task, to demonstrate the advantages of the CG model with lower calculation cost, we applied our new methodology to a larger biomolecule, integrin (αV) in the inactive state. Then, we sampled various structural ensembles, including extended structures that are apparently different from inactive ones.

Histone acetyltransferases (HATs) play a crucial role in post-translational modification. Among them, overexpression, mutation, or hyperfunction of EP300/CBP has been associated with various cancers. In this study, we identified the novel compound 2-chloro-5-[5-[(E)-[1-(3-chlorophenyl)-3-methyl-5-oxo-pyrazol-4-ylidene]methyl]-2-furyl]benzoic acid (1) as an EP300 HAT inhibitor via virtual screening. Further research has been focused on the design, synthesis, and in vitro biological evaluation of virtual hit derivatives. The studies revealed that 4-pyridone-3-carboxylic acid derivatives exhibited bioisosterism of benzoic acid. Replacement proved effective, providing compounds with similar EP300 HAT-inhibitory activity and improved cell growth-inhibitory activity compared to the benzoic acid analogs. Through these studies, we identified a potent and selective EP300/CBP HAT inhibitor., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

In molecular dynamics simulations, the limited time step size has been a barrier to simulating long-time behaviors. Implicit time integration methods allow markedly larger time steps than the standard explicit time method, although they have major drawbacks such as overheads solving linear systems and instability of Newton iterations. To overcome these issues, we propose a semi-implicit time integration scheme, the semi-implicit Hessian correction (SimHec) scheme, for overdamped Langevin dynamics. The method focuses on the Hessian matrices of bonded and nonbonded interactions, where components with large negative Hessian eigenvalues are cut off in the linear approximation of momentum equations to avoid instability. The narrow band Hessian matrix enables an efficient parallelized linear solution with an overlapping approximation. We tested SimHec for the interdomain fluctuations in adenylate kinase and the powerstroke transition of myosin II using a coarse-grained protein model. SimHec reproduced the same dynamics as the explicit method, although the transition dynamics tended to be accelerated and fluctuations in bonded potentials were slightly reduced. These deviations were corrected using a hybrid method, SimHec-H, which adds explicit time steps after the semi-implicit time step. The proposed scheme allowed us to use time steps 50-200 times larger than those in explicit time integration, which resulted in a speedup factor of 7-30 taking the overhead into account.

Approximately 15-30% of patients with lung cancer harbor mutations in the EGFR gene. Major EGFR mutations (>90% of EGFR-mutated lung cancer) are highly sensitive to EGFR tyrosine kinase inhibitors (TKIs). Many uncommon EGFR mutations have been identified, but little is known regarding their characteristics, activation, and sensitivity to various EGFR-TKIs, including allosteric inhibitors. We encountered a case harboring an EGFR-L747P mutation, originally misdiagnosed with EGFR-del19 mutation using a routine diagnostic EGFR mutation test, which was resistant to EGFR-TKI gefitinib. Using this minor mutation and common EGFR-activating mutations, we performed the binding free energy calculations and microsecond-timescale molecular dynamic (MD) simulations, revealing that the L747P mutation considerably stabilizes the active conformation through a salt-bridge formation between K745 and E762. We further revealed why several EGFR inhibitors, including the allosteric inhibitor, were ineffective. Our computational structural analysis strategy would be beneficial for future drug development targeting the EGFR minor mutations.

Here, we have constructed neural network-based models that predict atomic partial charges with high accuracy at low computational cost. The models were trained using high-quality data acquired from quantum mechanics calculations using the fragment molecular orbital method. We have succeeded in obtaining highly accurate atomic partial charges for three representative molecular systems of proteins, including one large biomolecule (approx. 2000 atoms). The novelty of our approach is the ability to take into account the electronic polarization in the system, which is a system-dependent phenomenon, being important in the field of drug design. Our high-precision models are useful for the prediction of atomic partial charges and expected to be widely applicable in structure-based drug designs such as structural optimization, high-speed and high-precision docking, and molecular dynamics calculations.

Biomolecular imaging using X-ray free-electron lasers (XFELs) has been successfully applied to serial femtosecond crystallography. However, the application of single-particle analysis for structure determination using XFELs with 100 nm or smaller biomolecules has two practical problems: the incomplete diffraction data sets for reconstructing 3D assembled structures and the heterogeneous conformational states of samples. A new diffraction template matching method is thus presented here to retrieve a plausible 3D structural model based on single noisy target diffraction patterns, assuming candidate structures. Two concepts are introduced here: prompt candidate diffraction, generated by enhanced sampled coarse-grain (CG) candidate structures, and efficient molecular orientation searching for matching based on Bayesian optimization. A CG model-based diffraction-matching protocol is proposed that achieves a 100-fold speed increase compared to exhaustive diffraction matching using an all-atom model. The conditions that enable multiconformational analysis were also investigated by simulated diffraction data for various conformational states of chromatin and ribosomes. The proposed method can enable multiconformational analysis, with a structural resolution of at least 20 Å for 270-800 Å flexible biomolecules, in experimental single-particle structure analyses that employ XFELs.

Accompanied with an increase of revealed biomolecular structures owing to advancements in structural biology, the molecular dynamics (MD) approach, especially coarse-grained (CG) MD suitable for macromolecules, is becoming increasingly important for elucidating their dynamics and behavior. In fact, CG-MD simulation has succeeded in qualitatively reproducing numerous biological processes for various biomolecules such as conformational changes and protein folding with reasonable calculation costs. However, CG-MD simulations strongly depend on various parameters, and selecting an appropriate parameter set is necessary to reproduce a particular biological process. Because exhaustive examination of all candidate parameters is inefficient, it is important to identify successful parameters. Furthermore, the successful region, in which the desired process is reproducible, is essential for describing the detailed mechanics of functional processes and environmental sensitivity and robustness. We propose an efficient search method for identifying the successful region by using two machine learning techniques, Bayesian optimization and active learning. We evaluated its performance using F1-ATPase, a biological rotary motor, with CG-MD simulations. We successfully identified the successful region with lower computational costs (12.3% in the best case) without sacrificing accuracy compared to exhaustive search. This method can accelerate not only parameter search but also biological discussion of the detailed mechanics of functional processes and environmental sensitivity based on MD simulation studies.

Neurotoxicity caused by nonfibrillar amyloid β (Aβ) oligomers in the brain is suggested to be associated with the onset of Alzheimer's disease (AD). Elucidating the structural features of Aβ oligomers is critical for promoting drug discovery research for AD. One of the Aβ oligomers, known as Aβ*56, is a dodecamer that impairs memory when injected into healthy rats, suggesting that Aβ*56 may contribute to cognitive deficits in AD patients. Another dodecamer structure, formed by 20-residue peptide segments derived from the Aβ peptide (Aβ 17-36 ), has been revealed by X-ray crystallography. The structure of the Aβ 17-36 dodecamer is composed of trimer units and shows the oligomer antibody A11 reactivity, which are characteristic of Aβ*56, indicating that Aβ*56 and the Aβ 17-36 dodecamer share a similar structure. However, the structure of the C-terminal regions (Aβ 37-42 ) remains unclear. The C-terminal region, which is abundant in hydrophobic residues, is thought to play a key role in stabilizing the oligomer structure by forming a hydrophobic core. In this study, we employed dissipative particle dynamics, a coarse-grained simulation method with soft core potentials, utilizing the crystal structure information to unravel Aβ dodecamer structures with C-terminal regions. The simulation results were validated by the reported experimental data. Hence, an analysis of the simulation results can provide structural insights into Aβ oligomers. Our simulations revealed the stabilization mechanism of the dodecamer structure at the molecular level. We showed that C-terminal regions spontaneously form a hydrophobic core in the central cavity, contributing to stabilizing the dodecamer structure. Furthermore, four consecutive hydrophobic residues in the C-terminal region (i.e., Val39-Ala42) are important for core formation.

Dynamics of nuclear proteins in crowded chromatin has only been poorly understood. Here, we address the diffusion, target search, and structural dynamics of three proteins in a model chromatin using coarse-grained molecular simulations run on the K computer. We prepared two structures of chromatin made of 20 nucleosomes with different nucleosome densities and investigated dynamics of two transcription factors, HMGB1 and p53, and one signaling protein, ERK, embedded in the chromatin. We found fast and normal diffusion of the nuclear proteins in the low-density chromatins and slow and subdiffusional movements in the high-density chromatin. The diffusion of the largest transcription factor, p53, is slowed by high-density chromatin most markedly. The on rates and off rates for DNA binding are increased and decreased, respectively, in the high-density chromatin. To our surprise, the DNA sequence search was faster in chromatin with high nucleosome density, though the diffusion is slower. We also found that the three nuclear proteins preferred to bind on the linker DNA and the entry and exit regions of nucleosomal DNA. In addition to these regions, HMGB1 and p53 also bound to the dyad., (Copyright © 2019 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Next generation sequencing (NGS)-based tumor profiling identified an overwhelming number of uncharacterized somatic mutations, also known as variants of unknown significance (VUS). The therapeutic significance of EGFR mutations outside mutational hotspots, consisting of >50 types, in nonsmall cell lung carcinoma (NSCLC) is largely unknown. In fact, our pan-nation screening of NSCLC without hotspot EGFR mutations ( n = 3,779) revealed that the majority (>90%) of cases with rare EGFR mutations, accounting for 5.5% of the cohort subjects, did not receive EGFR-tyrosine kinase inhibitors (TKIs) as a first-line treatment. To tackle this problem, we applied a molecular dynamics simulation-based model to predict the sensitivity of rare EGFR mutants to EGFR-TKIs. The model successfully predicted the diverse in vitro and in vivo sensitivities of exon 20 insertion mutants, including a singleton, to osimertinib, a third-generation EGFR-TKI ( R 2 = 0.72, P = 0.0037). Additionally, our model showed a higher consistency with experimentally obtained sensitivity data than other prediction approaches, indicating its robustness in analyzing complex cancer mutations. Thus, the in silico prediction model will be a powerful tool in precision medicine for NSCLC patients carrying rare EGFR mutations in the clinical setting. Here, we propose an insight to overcome mutation diversity in lung cancer., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

Background: Alectinib has shown a greater efficacy to ALK-rearranged non-small-cell lung cancers in first-line setting; however, most patients relapse due to acquired resistance, such as secondary mutations in ALK including I1171N and G1202R. Although ceritinib or lorlatinib was shown to be effective to these resistant mutants, further resistance often emerges due to ALK-compound mutations in relapse patients following the use of ceritinib or lorlatinib. However, the drug for overcoming resistance has not been established yet., Methods: We established lorlatinib-resistant cells harboring ALK-I1171N or -G1202R compound mutations by performing ENU mutagenesis screening or using an in vivo mouse model. We performed drug screening to overcome the lorlatinib-resistant ALK-compound mutations. To evaluate these resistances in silico, we developed a modified computational molecular dynamic simulation (MP-CAFEE)., Findings: We identified 14 lorlatinib-resistant ALK-compound mutants, including several mutants that were recently discovered in lorlatinib-resistant patients. Some of these compound mutants were found to be sensitive to early generation ALK-TKIs and several BCR-ABL inhibitors. Using our original computational simulation, we succeeded in demonstrating a clear linear correlation between binding free energy and in vitro experimental IC 50 value of several ALK-TKIs to single- or compound-mutated EML4-ALK expressing Ba/F3 cells and in recapitulating the tendency of the binding affinity reduction by double mutations found in this study. Computational simulation revealed that ALK-L1256F single mutant conferred resistance to lorlatinib but increased the sensitivity to alectinib., Interpretation: We discovered lorlatinib-resistant multiple ALK-compound mutations and an L1256F single mutation as well as the potential therapeutic strategies for these ALK mutations. Our original computational simulation to calculate the binding affinity may be applicable for predicting resistant mutations and for overcoming drug resistance in silico. FUND: This work was mainly supported by MEXT/JSPS KAKENHI Grants and AMED Grants., (Copyright © 2019 The Authors. Published by Elsevier B.V. All rights reserved.)

Conformationally restricted analogues of SPD-304, the first small-molecule TNFα inhibitor, in which two heteroaryl groups, indole and chromone, are connected by chiral methyl- or ethyl- cis-cyclopropane, were designed. Synthesis of these molecules was achieved via Suzuki-Miyaura or Stille coupling reactions with chiral bromomethylenecyclopropane or iodovinyl- cis-cyclopropane as the substrate, both of which were prepared from chiral methylenecyclopropane as a common intermediate, constructing the heteroaryl-methyl or -ethyl- cis-cyclopropane structures as key steps. This study presents an efficient synthesis of a series of chiral cis-cyclopropane conjugates with two heteroaryl groups.

High-performance computing approaches that combine molecular-scale and macroscale continuum mechanics have long been anticipated in various fields. Such approaches may enrich our understanding of the links between microscale molecular mechanisms and macroscopic properties in the continuum. However, there have been few successful examples to date owing to various difficulties associated with overcoming the large spatial (from 1 nm to 10 cm) and temporal (from 1 ns to 1 ms) gaps between the two scales. In this paper, we propose an efficient parallel scheme to couple a microscopic model using Langevin dynamics for a protein motor with a finite element continuum model of a beating heart. The proposed scheme allows us to use a macroscale time step that is an order of magnitude longer than the microscale time step of the Langevin model, without loss of stability or accuracy. This reduces the overhead required by the imbalanced loads of the microscale computations and the communication required when switching between scales. An example of the Langevin dynamics model that demonstrates the usefulness of the coupling approach is the molecular mechanism of the actomyosin system, in which the stretch-activation phenomenon can be successfully reproduced. This microscopic Langevin model is coupled with a macroscopic finite element ventricle model. In the numerical simulations, the Langevin dynamics model reveals that a single sarcomere can undergo spontaneous oscillation (15 Hz) accompanied by quick lengthening due to cooperative movements of the myosin molecules pulling on the common Z-line. Also, the coupled simulations using the ventricle model show that the stretch-activation mechanism contributes to the synchronization of the quick lengthening of the sarcomeres at the end of the systolic phase. By comparing the simulation results given by the molecular model with and without the stretch-activation mechanism, we see that this synchronization contributes to maintaining the systolic blood pressure by providing sufficient blood volume without slowing the diastolic process.

5- to 6-member ring enlargements of 3-oxa-2-silacyclopentylmethyl to 4-oxa-3-silacyclohexyl radicals were investigated by EPR spectroscopy and QM computations of model indano-oxasilacyclopentane and oxasilinanyl compounds. Both experimental and computational evidence favored a mechanism via a concerted 1,2-migration of the "tethered" Si-group. Thus, the "forbidden" 1,2-Si-group migration from carbon to carbon becomes allowed when the Si-group is "tethered". The EPR data from 3-oxa-2-silacyclopentylmethyl radicals disclosed ground state conformations having semioccupied p-orbitals close to antiperiplanar with respect to their β-Si-C bonds, but indicated Si-hyperconjugation (β-silicon effect) was insignificant in radicals. Kinetic data was obtained by the steady state EPR method for ring enlargement of indano-3-oxa-2-silacyclopentylmethyl radicals. The scope of the novel rearrangement in terms of other ring types and sizes, as well as the analogous 1,2-migration of "tethered" C-centered groups, was explored computationally.

Lipid chemical mediator resolvins with highly potent anti-inflammatory activity can be leads to develop novel anti-inflammatory drugs; however, they are unstable in oxygen due to their characteristic polyunsaturated structures. To solve the problem, CP-RvE2 has been designed and synthesized in which the cis-olefin of RvE2 was replaced with a cyclopropane. CP-RvE2s were much more stable than RvE2 against autoxidation and equipotent or more potent than RvE2. CP-RvE2s were successfully identified as stable equivalents of RvE2.

Due to hierarchic nature of biomolecular systems, their computational modeling calls for multiscale approaches, in which coarse-grained (CG) simulations are used to address long-time dynamics of large systems. Here, we review recent developments and applications of CG modeling methods, focusing on our methods primarily for proteins, DNA, and their complexes. These methods have been implemented in the CG biomolecular simulator, CafeMol. Our CG model has resolution such that ∼10 non-hydrogen atoms are grouped into one CG particle on average. For proteins, each amino acid is represented by one CG particle. For DNA, one nucleotide is simplified by three CG particles, representing sugar, phosphate, and base. The protein modeling is based on the idea that proteins have a globally funnel-like energy landscape, which is encoded in the structure-based potential energy function. We first describe two representative minimal models of proteins, called the elastic network model and the classic Go̅ model. We then present a more elaborate protein model, which extends the minimal model to incorporate sequence and context dependent local flexibility and nonlocal contacts. For DNA, we describe a model developed by de Pablo's group that was tuned to well reproduce sequence-dependent structural and thermodynamic experimental data for single- and double-stranded DNAs. Protein-DNA interactions are modeled either by the structure-based term for specific cases or by electrostatic and excluded volume terms for nonspecific cases. We also discuss the time scale mapping in CG molecular dynamics simulations. While the apparent single time step of our CGMD is about 10 times larger than that in the fully atomistic molecular dynamics for small-scale dynamics, large-scale motions can be further accelerated by two-orders of magnitude with the use of CG model and a low friction constant in Langevin dynamics. Next, we present four examples of applications. First, the classic Go̅ model was used to emulate one ATP cycle of a molecular motor, kinesin. Second, nonspecific protein-DNA binding was studied by a combination of elaborate protein and DNA models. Third, a transcription factor, p53, that contains highly fluctuating regions was simulated on two perpendicularly arranged DNA segments, addressing intersegmental transfer of p53. Fourth, we simulated structural dynamics of dinucleosomes connected by a linker DNA finding distinct types of internucleosome docking and salt-concentration-dependent compaction. Finally, we discuss many of limitations in the current approaches and future directions. Especially, more accurate electrostatic treatment and a phospholipid model that matches our CG resolutions are of immediate importance.

Molecular motors such as kinesin regulate affinity to a rail protein during the ATP hydrolysis cycle. The regulation mechanism, however, is yet to be determined. To understand this mechanism, we investigated the structural fluctuations of the motor head of the single-headed kinesin called KIF1A in different nucleotide states using molecular dynamics simulations of a Gō-like model. We found that the helix α4 at the microtubule (MT) binding site intermittently exhibits a large structural fluctuation when MT is absent. Frequency of this fluctuation changes systematically according to the nucleotide states and correlates strongly with the experimentally observed binding affinity to MT. We also showed that thermal fluctuation enhances the correlation and the interaction with the nucleotide suppresses the fluctuation of the helix α4. These results suggest that KIF1A regulates affinity to MT by changing the flexibility of the helix α4 during the ATP hydrolysis process: the binding site becomes more flexible in the strong binding state than in the weak binding state., (© 2015 Wiley Periodicals, Inc.)

KIF1A is a single-headed molecular motor that moves processively and unidirectionally along a microtubule by using the chemical energy released by hydrolyzing adenosine triphosphate (ATP) into adenosine diphosphate (ADP) and inorganic phosphate (P(i)). Although the movement of KIF1A seems to have successfully been explained by a simple Brownian motor model of the flashing ratchet type, this model is not suited to discuss the energetics of KIF1A. We introduce an elaborated model of the ratchet type to investigate how the chemical free energy is converted into mechanical work by taking account of the binding and release of reactant (ATP) and product (ADP and P(i)) molecules to and from the motor. The efficiency of energy transduction, the power output, and other quantities are calculated from the analytically obtained steady-state solution of the Fokker-Planck equations. It turns out that the concentrations of the reactant and product molecules that optimize both the efficiency and the power are close to those in the cell.

Background: We investigated the prevalence of depression among patients with COPD treated in long-term in-patient rehabilitation facilities, using the Center for Epidemiologic Studies Depression scale (CES-D). Furthermore, the relationship between the severity of air-flow obstruction (the percent of predicted FEV(1)), BODE (body mass index, degree of air-flow obstruction, dyspnea, exercise capacity) index, health-related quality of life (St George's Respiratory Questionnaire [SGRQ]), and depression were investigated., Methods: We recruited 74 in-patients (64 males, 10 females) with COPD. The mean age of the subjects was 72.7 years (range 52-85 y). Subjects completed the CES-D, and measurements were made of pulmonary function, body mass index, Modified Medical Research Council dyspnea scale, 6-min walk test (6MWT), and SGRQ., Results: Depression was evident in 48.6% (n = 36) of the subjects. A weak correlation was found between the CES-D scores and the percent of predicted FEV(1). The prevalence of depression showed a significant association with BODE stage. Scores for the SGRQ activity and impacts domains, and total SGRQ score were significantly worse in the subjects who were depressed., Conclusions: We found a high prevalence of depression among patients with stable COPD treated in long-term in-patient rehabilitation facilities. Depression among these patients, as measured by the CES-D, was associated with greater impairment in respiratory function and with poorer Modified Medical Research Council dyspnea scale and SGRQ scores. The prevalence of depression increased with BODE stage., (© 2013 Daedalus Enterprises.)

Kinesin is a family of molecular motors that move unidirectionally along microtubules (MT) using ATP hydrolysis free energy. In the family, the conventional two-headed kinesin was experimentally characterized to move unidirectionally through "walking" in a hand-over-hand fashion by coordinated motions of the two heads. Interestingly a single-headed kinesin, a truncated KIF1A, still can generate a biased Brownian movement along MT, as observed by in vitro single molecule experiments. Thus, KIF1A must use a different mechanism from the conventional kinesin to achieve the unidirectional motions. Based on the energy landscape view of proteins, for the first time, we conducted a set of molecular simulations of the truncated KIF1A movements over an ATP hydrolysis cycle and found a mechanism exhibiting and enhancing stochastic forward-biased movements in a similar way to those in experiments. First, simulating stand-alone KIF1A, we did not find any biased movements, while we found that KIF1A with a large friction cargo-analog attached to the C-terminus can generate clearly biased Brownian movements upon an ATP hydrolysis cycle. The linked cargo-analog enhanced the detachment of the KIF1A from MT. Once detached, diffusion of the KIF1A head was restricted around the large cargo which was located in front of the head at the time of detachment, thus generating a forward bias of the diffusion. The cargo plays the role of a diffusional anchor, or cane, in KIF1A "walking."

For simulating proteins at work in millisecond time scale or longer, we develop a coarse-grained (CG) molecular dynamics (MD) method and software, CafeMol. At the resolution of one-particle-per-residue, CafeMol equips four structure-based protein models: (1) the off-lattice Go model, (2) the atomic interaction based CG model for native state and folding dynamics, (3) the multiple-basin model for conformational change dynamics, and (4) the elastic network model for quasiharmonic fluctuations around the native structure. Ligands can be treated either explicitly or implicitly. For mimicking functional motions of proteins driven by some external force, CafeMol has various and flexible means to "switch" the energy functions that induce active motions of the proteins. CafeMol can do parallel computation with modest sized PC clusters. We describe CafeMol methods and illustrate it with several examples, such as rotary motions of F1-ATPase and drug exports from a transporter. The CafeMol source code is available at www.cafemol.org .

Aim: The objective of this study was to utilize commonly applied tools, the Hospital Anxiety and Depression Scale - Depression subscale (HADS-D) and the Center for Epidemiological Studies Depression Scale (CES-D), to screen for depressive symptoms in patients with stable chronic obstructive pulmonary disease (COPD). Furthermore, we sought to identify whether differences existed in the prevalence of depressive symptoms as assessed by CES-D and HADS-D, and predictors of depressive symptoms., Methods: The presence of depressive symptoms in 80 outpatients and 51 inpatients with stable COPD was assessed using the CES-D and HADS-D. Data regarding sex, educational level, body mass index, smoking index and pulmonary function were obtained to evaluate their independent contribution as predictors of depressive symptoms., Results: The prevalence of depressive symptoms was 29.8% based on CES-D and 40.5% based on HADS-D. A MacNemar test of COPD severity and analysis of the results of depressive symptoms based on CES-D and HADS-D revealed significant differences. Logistic regression analysis suggested that 'severity' is a predictor of depressive symptoms as assessed by CES-D, whereas 'body mass index', 'education level' and 'setting' were predictors of depressive symptoms as assessed by HADS-D., Conclusions: The prevalence of depressive symptoms differed when assessed with CES-D and HADS-D. The reasons behind this difference include the fact that HADS-D frequently detected depressive symptoms in patients with mild COPD as well as a tendency for HADS-D to be strongly influenced by education levels. In contrast, the severity of COPD was reflected in CES-D. It is possible that prevalence of depressive symptoms differs in accordance with the applied screening tool., (© 2011 The Authors. Psychiatry and Clinical Neurosciences © 2011 Japanese Society of Psychiatry and Neurology.)

The cytotoxin-associated antigen A (CagA) of Helicobacter pylori prevalent in East Asian countries, where the mortality rate due to gastric cancer is high, has been reported to be structurally different from that in Western countries, where the gastric cancer mortality rate is relatively low. Based on the structural features of the EPIYA motifs located at the carboxyl terminal of the protein, CagA was subdivided into two types: East Asian CagA and Western CagA. A recent study suggested that immunohistochemistry with anti-East Asian-specific antibody (alpha-EAS Ab), which was specifically immunoreactive with East Asian CagA but not with Western CagA, may be useful for diagnosis of the cagA genotype. To further evaluate the value of this diagnostic method in terms of sensitivity, specificity, and accuracy, 143 gastric biopsy specimens with alpha-EAS Ab were analyzed on immunohistochemistry and compared with the sequencing of the cagA gene. It was found that diagnosis of the cagA genotype of H. pylori on immunohistochemistry using the alpha-EAS Ab was highly sensitive (sensitivity 93.2%) and specific (specificity 72.7%), suggesting that immunohistochemical diagnosis of the cagA genotype is useful for diagnosis of the cagA genotype.

Cytotoxin-associated antigen A (CagA) protein produced by Helicobacter pylori is proposed to be associated with the pathogenesis of gastric cancer as well as gastritis and gastroduodenal ulcer. It has been reported that the CagA of H. pylori widespread in East Asian countries, where the mortality rate due to gastric cancer is high, is structurally different from that in Western countries, where the gastric cancer mortality rate is relatively low. In this study, we generated an antibody, East Asian CagA-specific antibody (alpha-EAS Ab), which is specifically immunoreactive with East Asian CagA but not with Western CagA. The CagA was immunohistochemically detected at the surface of the gastric mucosa. Interestingly, positive immunoreactivity was also detected in the nucleus and cytoplasm of the infected gastric epithelium, suggesting that CagA may play some pathogenic role in both the nucleus and cytoplasm. Immunohistochemistry of 47 gastric biopsy specimens detected East Asian CagA-positive H. pylori in 43 cases. In 46 of the 47 cases examined, the data obtained by immunohistochemistry were completely consistent with those obtained by sequencing of the cagA gene of the isolated strain, suggesting that our immunohistochemical method is reliable and useful for diagnosis of the cagA genotype.

The processive motion of two-headed molecular motors is studied theoretically by introducing a model that takes into account the coordinated motion of the constituent heads and the detachment process of heads from linear molecular tracks. The mean velocity, the mean run length, and the mean run time of the motor along the track are calculated numerically based on the Langevin equation. It turns out that the model, with appropriate choice of model parameters, can explain qualitatively the dependence of these quantities on the external load and adenosin triphosphate concentration observed experimentally for kinesin motors. Furthermore, we discuss how the motility and processivity of the motor are affected by various model parameters, which may be tested by experiments.