記述言語：英語   掲載種別：研究論文（学術雑誌）  

Mutations of G protein-coupled receptors (GPCRs) cause various human diseases, but the mechanistic details are limited. Here, we establish p.E303K in the gene encoding the endothelin receptor type A (ETAR/EDNRA) as a recurrent mutation causing mandibulofacial dysostosis with alopecia (MFDA), with craniofacial changes similar to those caused by p.Y129F. Mouse models carrying either of these missense mutations exhibited a partial maxillary-to-mandibular transformation, which was rescued by deleting the ligand endothelin 3 (ET3/EDN3). Pharmacological experiments confirmed the causative ETAR mutations as gain of function, dependent on ET3. To elucidate how an amino acid substitution far from the ligand binding site can increase ligand affinity, we used molecular dynamics (MD) simulations. E303 is located at the intracellular end of transmembrane domain 6, and its replacement by a lysine increased flexibility of this portion of the helix, thus favoring G protein binding and leading to G protein-mediated enhancement of agonist affinity. The Y129F mutation located under the ligand binding pocket reduced the sodium-water network, thereby affecting the extracellular portion of helices in favor of ET3 binding. These findings provide insight into the pathogenesis of MFDA and into allosteric mechanisms regulating GPCR function, which may provide the basis for drug design targeting GPCRs.

DOI： 10.1172/JCI151536

PubMed

researchmap

掲載種別：研究論文（学術雑誌）   出版者・発行元：Springer Science and Business Media LLC  

Abstract

Toll-like receptor 3 (TLR3) is a member of the TLR family, which plays an important role in the innate immune system and is responsible for recognizing viral double-stranded RNA (dsRNA). Previous biochemical and structural studies have revealed that a minimum length of approximately 40–50 base pairs of dsRNA is necessary for TLR3 binding and dimerization. However, efficient TLR3 activation requires longer dsRNA and the molecular mechanism underlying its dsRNA length-dependent activation remains unknown. Here, we report cryo-electron microscopy analyses of TLR3 complexed with longer dsRNA. TLR3 dimers laterally form a higher multimeric complex along dsRNA, providing the basis for cooperative binding and efficient signal transduction.

DOI： 10.1038/s41467-023-35844-2

researchmap

その他リンク： https://www.nature.com/articles/s41467-023-35844-2

記述言語：英語   掲載種別：研究論文（学術雑誌）  

Based on many crystal structures of ligand complexes, much study has been devoted to understanding the molecular recognition of SARS-CoV-2 3C-like protease (3CLpro), a potent drug target for COVID-19. In this research, to extend this present static view, we examined the kinetic process of binding/unbinding of an eight-residue substrate peptide to/from 3CLpro by evaluating the path ensemble with the weighted ensemble simulation. The path ensemble showed the mechanism of how a highly flexible peptide folded into the bound form. At the early stage, the dominant motion was the diffusion on the protein surface showing a broad distribution, whose center was led into the cleft of the chymotrypsin fold. We observed a definite sequential formation of the hydrogen bonds at the later stage occurring in the cleft, initiated between Glu166 (3CLpro) and P3_Val (peptide), followed by binding to the oxyanion hole and completed by the sequence-specific recognition at P1_Gln.

DOI： 10.1021/acs.jcim.2c00946

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

UNLABELLED: SARS-CoV-2 3C-like protease (3CLpro), a potential therapeutic target for COVID-19, consists of a chymotrypsin fold and a C-terminal α-helical domain (domain III), the latter of which mediates dimerization required for catalytic activation. To gain further understanding of the functional dynamics of SARS-CoV-2 3CLpro, this review extends the scope to the comparative study of many crystal structures of proteases having the chymotrypsin fold (clan PA of the MEROPS database). First, the close correspondence between the zymogen-enzyme transformation in chymotrypsin and the allosteric dimerization activation in SARS-CoV-2 3CLpro is illustrated. Then, it is shown that the 3C-like proteases of family Coronaviridae (the protease family C30), which are closely related to SARS-CoV-2 3CLpro, have the same homodimeric structure and common activation mechanism via domain III mediated dimerization. The survey extended to order Nidovirales reveals that all 3C-like proteases belonging to Nidovirales have domain III, but with various chain lengths, and 3CLpro of family Mesoniviridae (family C107) has the same homodimeric structure as that of C30, even though they have no sequence similarity. As a reference, monomeric 3C proteases belonging to the more distant family Picornaviridae (family C3) lacking domain III are compared with C30, and it is shown that the 3C proteases are rigid enough to maintain their structures in the active state. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12551-022-01020-x.

DOI： 10.1007/s12551-022-01020-x

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

Semaphorins constitute a large family of secreted and membrane-bound proteins that signal through cell-surface receptors, plexins. Semaphorins generally use low-affinity protein-protein interactions to bind with their specific plexin(s) and regulate distinct cellular processes such as neurogenesis, immune response, and organogenesis. Sema6D is a membrane-bound semaphorin that interacts with class A plexins. Sema6D exhibited differential binding affinities to class A plexins in prior cell-based assays, but the molecular mechanism underlying this selectivity is not well understood. Therefore, we performed hybrid in vitro/in silico analysis to examine the binding mode of Sema6D to class A plexins and to identify residues that give rise to the differential affinities and thus contribute to the selectivity within the same class of semaphorins. Our biophysical binding analysis indeed confirmed that Sema6D has a higher affinity for Plexin-A1 than for other class A plexins, consistent with the binding selectivity observed in the previous cell-based assays. Unexpectedly, our present crystallographic analysis of the Sema6D-Plexin-A1 complex showed that the pattern of polar interactions is not interaction-specific because it matches the pattern in the prior structure of the Sema6A-Plexin-A2 complex. Thus, we performed in silico alanine scanning analysis and discovered hotspot residues that selectively stabilized the Sema6D-Plexin-A1 pair via Van der Waals interactions. We then validated the contribution of these hotspot residues to the variation in binding affinity with biophysical binding analysis and molecular dynamics simulations on the mutants. Ultimately, our present results suggest that shape complementarity in the binding interfaces is a determinant for binding selectivity.

DOI： 10.1002/pro.4452

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

Hydration free energy (HFE) is a key factor in improving protein-ligand binding free energy (BFE) prediction accuracy. The HFE itself can be calculated using the three-dimensional reference interaction model (3D-RISM); however, the BFE predictions solely evaluated using 3D-RISM are not correlated to the experimental BFE for abundant protein-ligand pairs. In this study, to predict the BFE for multiple sets of protein-ligand pairs, we propose a machine learning approach incorporating the HFEs obtained using 3D-RISM, termed 3D-RISM-AI. In the learning process, structural metrics, intra-/intermolecular energies, and HFEs obtained via 3D-RISM of ∼4000 complexes in the PDBbind database (ver. 2018) were used. The BFEs predicted using 3D-RISM-AI were well correlated to the experimental data (Pearson's correlation coefficient of 0.80 and root-mean-square error of 1.91 kcal/mol). As important factors for the prediction, the difference in the solvent accessible surface area between the bound and unbound structures and the hydration properties of the ligands were detected during the learning process.

DOI： 10.1021/acs.jpcb.2c03384

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

The coronavirus membrane protein (M) is the most abundant viral structural protein and plays a central role in virus assembly and morphogenesis. However, the process of M protein-driven virus assembly are largely unknown. Here, we report the cryo-electron microscopy structure of the SARS-CoV-2 M protein in two different conformations. M protein forms a mushroom-shaped dimer, composed of two transmembrane domain-swapped three-helix bundles and two intravirion domains. M protein further assembles into higher-order oligomers. A highly conserved hinge region is key for conformational changes. The M protein dimer is unexpectedly similar to SARS-CoV-2 ORF3a, a viral ion channel. Moreover, the interaction analyses of M protein with nucleocapsid protein (N) and RNA suggest that the M protein mediates the concerted recruitment of these components through the positively charged intravirion domain. Our data shed light on the M protein-driven virus assembly mechanism and provide a structural basis for therapeutic intervention targeting M protein.

DOI： 10.1038/s41467-022-32019-3

PubMed

researchmap

記述言語：英語  

DOI： 10.1021/jacs.2c07339

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

Inspired by mechanosensitive potassium channels found in nature, we developed a fluorinated amphiphilic cyclophane composed of fluorinated rigid aromatic units connected via flexible hydrophilic octa(ethylene glycol) chains. Microscopic and emission spectroscopic studies revealed that the cyclophane could be incorporated into the hydrophobic layer of the lipid bilayer membranes and self-assembled to form a supramolecular transmembrane ion channel. Current recording measurements using cyclophane-containing planer lipid bilayer membranes successfully demonstrated an efficient transmembrane ion transport. We also demonstrated that the ion transport property was sensitive to the mechanical forces applied to the membranes. In addition, ion transport assays using pH-sensitive fluorescence dye revealed that the supramolecular channel possesses potassium ion selectivity. We also performed all-atom hybrid quantum-mechanical/molecular mechanical simulations to assess the channel structures at atomic resolution and the mechanism of selective potassium ion transport. This research demonstrated the first example of a synthetic mechanosensitive potassium channel, which would open a new door to sensing and manipulating biologically important processes and purification of key materials in industries.

DOI： 10.1021/jacs.2c04118

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

Membrane permeability of cyclic peptides is an important factor in drug design. To investigate the membrane permeability of cyclic peptides using molecular dynamics (MD) simulations, the accurate force fields for unnatural amino acids present in the cyclic peptides are required. Therefore, we developed the CHARMM force fields of the unnatural amino acids present in cyclosporin A (CsA), a cyclic peptide used as an immune suppressor. Especially for N-methyl amino acids, which contribute to the membrane permeability of cyclic peptides, we developed a grid correction map (CMAP) of the energy surface using the φ and ψ dihedral angles in the main chain of CsA. To validate the developed force field, we performed MD simulations, including the generalized replica exchange with solute tempering method, of CsA in water and chloroform solvents. The conformations of CsA in water and chloroform sampled using the developed force field were consistent with those of the experimental results of the solution nuclear magnetic resonance spectroscopy.

DOI： 10.2142/biophysico.bppb-v19.0045

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

The accumulation of epigenetic alterations is one of the major causes of tumorigenesis. Aberrant DNA methylation patterns cause genome instability and silencing of tumor suppressor genes in various types of tumors. Therefore, drugs that target DNA methylation-regulating factors have great potential for cancer therapy. Ubiquitin-like containing PHD and RING finger domain 1 (UHRF1) is an essential factor for DNA methylation maintenance. UHRF1 is overexpressed in various cancer cells and down-regulation of UHRF1 in these cells reactivates the expression of tumor suppressor genes, thus UHRF1 is a promising target for cancer therapy. We have previously shown that interaction between the tandem Tudor domain (TTD) of UHRF1 and DNA ligase 1 (LIG1) di/trimethylated on Lys126 plays a key role in the recruitment of UHRF1 to replication sites and replication-coupled DNA methylation maintenance. An arginine binding cavity (Arg-binding cavity) of the TTD is essential for LIG1 interaction, thus the development of inhibitors that target the Arg-binding cavity could potentially repress UHRF1 function in cancer cells. To develop such an inhibitor, we performed in silico screening using not only static but also dynamic metrics based on all-atom molecular dynamics simulations, resulting in efficient identification of 5-amino-2,4-dimethylpyridine (5A-DMP) as a novel TTD-binding compound. Crystal structure of the TTD in complex with 5A-DMP revealed that the compound stably bound to the Arg-binding cavity of the TTD. Furthermore, 5A-DMP inhibits the full-length UHRF1:LIG1 interaction in Xenopus egg extracts. Our study uncovers a UHRF1 inhibitor which can be the basis of future experiments for cancer therapy.

DOI： 10.1016/j.bmc.2021.116500

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

The 3C-like protease (3CLpro) of SARS-CoV-2 is a potential therapeutic target for COVID-19. Importantly, it has an abundance of structural information solved as a complex with various drug candidate compounds. Collecting these crystal structures (83 Protein Data Bank (PDB) entries) together with those of the highly homologous 3CLpro of SARS-CoV (101 PDB entries), we constructed the crystal structure ensemble of 3CLpro to analyze the dynamic regulation of its catalytic function. The structural dynamics of the 3CLpro dimer observed in the ensemble were characterized by the motions of four separate loops (the C-loop, E-loop, H-loop, and Linker) and the C-terminal domain III on the rigid core of the chymotrypsin fold. Among the four moving loops, the C-loop (also known as the oxyanion binding loop) causes the order (active)-disorder (collapsed) transition, which is regulated cooperatively by five hydrogen bonds made with the surrounding residues. The C-loop, E-loop, and Linker constitute the major ligand binding sites, which consist of a limited variety of binding residues including the substrate binding subsites. Ligand binding causes a ligand size dependent conformational change to the E-loop and Linker, which further stabilize the C-loop via the hydrogen bond between the C-loop and E-loop. The T285A mutation from SARS-CoV 3CLpro to SARS-CoV-2 3CLpro significantly closes the interface of the domain III dimer and allosterically stabilizes the active conformation of the C-loop via hydrogen bonds with Ser1 and Gly2; thus, SARS-CoV-2 3CLpro seems to have increased activity relative to that of SARS-CoV 3CLpro.

DOI： 10.1016/j.jmb.2021.167324

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

The mineralocorticoid receptor (MR) is a nuclear receptor whose endogenous ligands are mineralocorticoids, a type of steroid hormone. The activating S810L mutation is known to cause severe early-onset and pregnancy-related hypertension. Progesterone binds to the wild-type (WT) MR as a passive antagonist with fast dissociation; however, it binds to the S810L mutant as a full agonist with slow dissociation. The switch in the biological activity of progesterone is considered to be one of the causes of the disease. First, we used steered molecular dynamics simulations to analyze the dissociation process of progesterone for the WT and the S810L mutant. Progesterone in the WT dissociated from the ligand-binding pocket with a weak force in comparison with progesterone in the S810L mutant due to the large inflow of water molecules into the pocket. Therefore, we used conventional molecular dynamics simulations for the ligand-free structures of the WT and the S810L mutant to investigate the effect of the mutation on the inflow of water. In the WT, water molecules enter the ligand-binding pocket in two ways: in the vicinity of (i) Arg817 and (ii) Ser810. In contrast, few water molecules enter the pocket in the S810L mutant because of the large size and hydrophobic nature of the Leu810 side chain. Fast dissociation is a common feature among passive antagonists of MR; therefore, we inferred that the water inflow could be responsible for the dissociation kinetics of progesterone in the WT and the S810L mutant.

DOI： 10.1021/acs.jcim.1c00389

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

The vitamin D receptor ligand-binding domain (VDR-LBD) undergoes conformational changes upon ligand binding. In this nuclear receptor family, agonistic or antagonistic activities are controlled by the conformation of the helix (H)12. However, all crystal structures of VDR-LBD reported to date correspond to the active H12 conformation, regardless of agonist/antagonist binding. To understand the mechanism of VDR-LBD regulation structurally, conformational samplings of agonist- and antagonist-bound rat VDR-LBD were performed using the generalized replica exchange with solute tempering (gREST) method. The gREST simulations demonstrated different structural responses of rat VDR-LBD to agonist or antagonist binding, whereas in conventional molecular dynamics simulations, the conformation was the same as that of the crystal structures, regardless of agonist/antagonist binding. In the gREST simulations, a spontaneous conformational change of H12 was observed only for the antagonist complex. The different responses to agonist/antagonist binding were attributed to hydrophobic core formation at the ligand-binding pocket and cooperative rearrangements of H11. The gREST method can be applied to the examination of structure-activity relationships for multiple VDR-LBD ligands.

DOI： 10.1021/acs.jcim.1c00534

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

Antibody labeling has been conducted extensively for structure determination using both X-ray crystallography and electron microscopy (EM). However, establishing target-specific antibodies is a prerequisite for applying antibody-assisted structural analysis. To expand the applicability of this strategy, an alternative method has been developed to prepare an antibody complex by inserting an exogenous epitope into the target. It has already been demonstrated that the Fab of the NZ-1 monoclonal antibody can form a stable complex with a target containing a PA12 tag as an inserted epitope. Nevertheless, it was also found that complex formation through the inserted PA12 tag inevitably caused structural changes around the insertion site on the target. Here, an attempt was made to improve the tag-insertion method, and it was consequently discovered that an alternate tag (PA14) could replace various loops on the target without inducing large structural changes. Crystallographic analysis demonstrated that the inserted PA14 tag adopts a loop-like conformation with closed ends in the antigen-binding pocket of the NZ-1 Fab. Due to proximity of the termini in the bound conformation, the more optimal PA14 tag had only a minor impact on the target structure. In fact, the PA14 tag could also be inserted into a sterically hindered loop for labeling. Molecular-dynamics simulations also showed a rigid structure for the target regardless of PA14 insertion and complex formation with the NZ-1 Fab. Using this improved labeling technique, negative-stain EM was performed on a bacterial site-2 protease, which enabled an approximation of the domain arrangement based on the docking mode of the NZ-1 Fab.

DOI： 10.1107/S2059798321002527

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

TFIIH is a crucial transcription and DNA repair factor consisting of the seven-subunit core. The core subunit p62 contains a pleckstrin homology domain (PH-D), which is essential for locating TFIIH at transcription initiation and DNA damage sites, and two BSD (BTF2-like transcription factors, synapse-associated proteins and DOS2-like proteins) domains. A recent cryo-electron microscopy (cryo-EM) structure of human TFIIH visualized most parts of core, except for the PH-D. Here, by nuclear magnetic resonance spectroscopy we have established the solution structure of human p62 PH-D connected to the BSD1 domain by a highly flexible linker, suggesting the flexibility of PH-D in TFIIH. Based on this dynamic character, the PH-D was modeled in the cryo-EM structure to obtain the whole human TFIIH core structure, which indicates that the PH-D moves around the surface of core with a specific but limited spatial distribution; these dynamic structures were refined by molecular dynamics (MD) simulations. Furthermore, we built models, also refined by MD simulations, of TFIIH in complex with five p62-binding partners, including transcription factors TFIIEα, p53 and DP1, and nucleotide excision repair factors XPC and UVSSA. The models explain why the PH-D is crucially targeted by these factors, which use their intrinsically disordered acidic regions for TFIIH recruitment.

DOI： 10.1093/nar/gkaa1045

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

Transmembrane anion transport is an important biological process in maintaining cellular functions. Thus, synthetic anion transporters are widely developed for their biological applications. Imidazolinium was introduced as anion recognition site to a multiblock amphiphilic structure that consists of octa(ethylene glycol) and aromatic units. Ion transport assay using halide-sensitive lucigenin and pH-sensitive 8-hydroxypyrene-1,3,6-trisulfonate (HPTS) revealed that imidazolinium-based multiblock amphiphile (IMA) transports anions and showed high selectivity for nitrate, which plays crucial roles in many biological events. Temperature-dependent ion transport assay using 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) indicated that IMA works as a mobile carrier. 1 H NMR titration experiments indicated that the C2 proton of the imidazolinium ring recognizes anions via a (C-H)+ ⋅⋅⋅X- hydrogen bond. Furthermore, all-atom molecular dynamics simulations revealed a dynamic feature of IMA within the membranes during ion transportation.

DOI： 10.1002/asia.202001106

PubMed

researchmap

掲載種別：研究論文（学術雑誌）   出版者・発行元：Wiley  

DOI： 10.1002/asia.202001415

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Elsevier BV  

Ubiquitin-like with PHD and RING finger domains 1 (UHRF1) is an essential factor for the maintenance of mammalian DNA methylation and harbors several reader modules for recognizing epigenetic marks. The tandem Tudor domain (TTD) of UHRF1 has a peptide-binding groove that functions as a binding platform for intra- or intermolecular interactions. Besides the groove interacting with unphosphorylated linker 2 and spacer of UHRF1, it also interacts with di/tri-methylated histone H3 at Lys9 and DNA ligase 1 (LIG1) at Lys126. Here we focus on the phosphorylation of Ser298 in linker 2, which was implied to regulate the ligand-binding property of the TTD. Although the protein expression level of UHRF1 is unchanged throughout the cell cycle, Ser298 phosphorylated form of UHRF1 is notably increased in the G2/M phase, which is revealed by immunoprecipitation followed by Western blotting. Molecularly, while unphosphorylated linker 2 covers the peptide-binding groove to prevent access of other interactors, small-angle X-ray scattering, thermal stability assay and molecular dynamics simulation revealed that the phosphate group of Ser298 dissociates linker 2 from the peptide-binding groove of the TTD to permit the other interactors to access to the groove. Our data reveal a mechanism in which Ser298 phosphorylation in linker 2 triggers a change of the TTD's structure and may affect multiple functions of UHRF1 by facilitating associations with LIG1 at DNA replication sites and histone H3K9me2/me3 at heterochromatic regions.

DOI： 10.1016/j.jmb.2020.05.006

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

Biological membranes play pivotal roles in the cellular activities. Transmembrane proteins are the central molecules that conduct membrane-mediated biochemical functions such as signal transduction and substance transportation. Not only the molecular functions but also the supramolecular properties of the transmembrane proteins such as self-assembly, delocalization, orientation and signal response are essential for controlling cellular activities. Here we report anisotropic ligand responses of a synthetic multipass transmembrane ion channel. An unsymmetrical molecular structure allows for oriented insertion of the synthetic amphiphile to a bilayer by addition to a pre-formed membrane. Complexation with a ligand prompts ion transportation by forming a supramolecular channel, and removal of the ligand deactivates the transportation function. Biomimetic regulation of the synthetic channel by agonistic and antagonistic ligands is also demonstrated not only in an artificial membrane but also in a biological membrane of a living cell.

DOI： 10.1038/s41467-020-16770-z

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：AIP Publishing  

A new method is developed for calculating hydration free energies (HFEs) of polyatomic solutes. The solute insertion is decomposed into the creation of a cavity in water matching the geometric characteristics of the solute at the atomic level (process 1) and the incorporation of solute-water van der Waals and electrostatic interactions (process 2). The angle-dependent integral equation theory combined with our morphometric approach and the three-dimensional interaction site model theory are applied to processes 1 and 2, respectively. Neither a stage of training nor parameterization is necessitated. For solutes with various sizes including proteins, the HFEs calculated by the new method are compared to those obtained using a molecular dynamics simulation based on solution theory in energy representation (the ER method developed by Matubayasi and co-workers), currently the most reliable tool. The agreement is very good especially for proteins. The new method is characterized by the following: The calculation can rapidly be finished; a solute possessing a significantly large total charge can be handled without difficulty; and since it yields not only the HFE but also its many physically insightful energetic and entropic components, it is best suited to the elucidation of mechanisms of diverse phenomena such as the receptor-ligand binding, different types of molecular recognition, and protein folding, denaturation, and association.

DOI： 10.1063/1.5093110

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

The combination of molecular dynamics (MD) simulations and small-angle X-ray scattering (SAXS), called the MD-SAXS method, is efficient for investigating protein dynamics. To overcome the time-scale limitation of all-atom MD simulations, coarse-grained (CG) representations are often utilized for biomolecular simulations. In this study, we propose a method to combine CG MD simulations with SAXS, termed the CG-MD-SAXS method. In the CG-MD-SAXS method, the scattering factors of CG particles for proteins and nucleic acids are evaluated using high-resolution structural data in the Protein Data Bank, and the excluded volume and the hydration shell are modeled using two adjustable parameters to incorporate solvent effects. To avoid overfitting, only the two parameters are adjusted for an entire structure ensemble. To verify the developed method, theoretical SAXS profiles for various proteins, DNA/RNA, and a protein-RNA complex are compared with both experimental profiles and theoretical profiles obtained by the all-atom representation. In the present study, we applied the CG-MD-SAXS method to the Swi5-Sfr1 complex and three types of nucleosomes to obtain reliable ensemble models consistent with the experimental SAXS data.

DOI： 10.2142/biophysico.16.0_377

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Frontiers Media SA  

V1-ATPase exemplifies the ubiquitous rotary motor, in which a central shaft DF complex rotates inside a hexagonally arranged catalytic A3B3 complex, powered by the energy from ATP hydrolysis. We have recently reported a number of crystal structures of the Enterococcus hirae A3B3DF (V1) complex corresponding to its nucleotide-bound intermediate states, namely the forms waiting for ATP hydrolysis (denoted as catalytic dwell), ATP binding (ATP-binding dwell), and ADP release (ADP-release dwell) along the rotatory catalytic cycle of ATPase. Furthermore, we have performed microsecond-scale molecular dynamics simulations and free-energy calculations to investigate the conformational transitions between these intermediate states and to probe the long-time dynamics of the molecular motor. In this article, the molecular structure and dynamics of the V1-ATPase are reviewed to bring forth a unified model of the motor's remarkable rotational mechanism.

DOI： 10.3389/fphys.2019.00046

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）  

The accurate calculation of protein-ligand binding free energies is necessary for computer-aided drug design. The alchemical perturbation method frequently used for binding free energy calculations under periodic boundary conditions suffers from finite-size effects related to the cell-size dependence of the charging free energy at different cell sizes. The finite-size effect on the binding free energy of charged ligands is not negligible in comparison to the binding free energy itself. In this study, we propose an effective perturbation protocol for calculating the binding free energy termed the "warp-drive" method for eliminating the finite-size effect. When the warp-drive method is applied, a solution system consisting of a protein-ligand complex and an unbound ligand located at a distant position is used. Diminished partial charges of the bound ligand simultaneously emerge in the other unbound ligand, and in turn, the total charge of the system does not change at all intermediate states. To assess the performance of the warp-drive method, charging free energies for systematically varied cell sizes are examined and compared to those calculated via alchemical perturbation. In contrast to that of alchemical perturbation, the charging free energy obtained via the warp-drive method does not exhibit finite-size effects, even for typical cell sizes without any corrections, and this result is in good agreement with that calculated on the basis of alchemical perturbation levels measured from large cells with full corrections of the finite-size effect. This finding reveals an advantage of the warp-drive method, as alchemical perturbation is computationally costly due to the large cell sizes and specificities involved in correction schemes depending on the total charge of proteins and components of solvent molecules.

DOI： 10.1021/acs.jctc.8b00280

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Springer Verlag  

Protein functions require specific structures frequently coupled with conformational changes. The scale of the structural dynamics of proteins spans from the atomic to the molecular level. Theoretically, all-atom molecular dynamics (MD) simulation is a powerful tool to investigate protein dynamics because the MD simulation is capable of capturing conformational changes obeying the intrinsically structural features. However, to study long-timescale dynamics, efficient sampling techniques and coarse-grained (CG) approaches coupled with all-atom MD simulations, termed multiscale MD simulations, are required to overcome the timescale limitation in all-atom MD simulations. Here, we review two examples of rotary motor proteins examined using free energy landscape (FEL) analysis and CG-MD simulations. In the FEL analysis, FEL is calculated as a function of reaction coordinates, and the long-timescale dynamics corresponding to conformational changes is described as transitions on the FEL surface. Another approach is the utilization of the CG model, in which the CG parameters are tuned using the fluctuation matching methodology with all-atom MD simulations. The long-timespan dynamics is then elucidated straightforwardly by using CG-MD simulations.

DOI： 10.1007/s12551-017-0373-4

Scopus

PubMed

researchmap

DOI： 10.1007/978-981-13-2200-6_15

PubMed

researchmap

記述言語：日本語   出版者・発行元：公益社団法人 日本ビタミン学会  

DOI： 10.20632/vso.92.1_18

CiNii Books

CiNii Research

researchmap

その他リンク： https://ndlsearch.ndl.go.jp/books/R000000004-I028798270

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：CELL PRESS  

Enterococcus hirae V1-ATPase is a molecular motor composed of the A3B3 hexamer ring and the central stalk. In association with ATP hydrolysis, three catalytic AB pairs in the A3B3 ring undergo conformational changes, which lead to a 120° rotation of the central stalk. To understand how the conformational changes of three catalytic pairs induce the 120° rotation of the central stalk, we performed multiscale molecular dynamics (MD) simulations in which coarse-grained and all-atom MD simulations were combined using a fluctuation matching methodology. During the rotation, a catalytic AB pair spontaneously adopted an intermediate conformation, which was not included in the initial inputs of the simulations and was essentially close to the "bindable-like" structure observed in a recently solved crystal structure. Furthermore, the creation of a space between the bindable-like and tight pairs was required for the central stalk to rotate without steric hindrance. These cooperative rearrangements of the three catalytic pairs are crucial for the rotation of the central stalk.

DOI： 10.1016/j.bpj.2017.01.029

Web of Science

PubMed

researchmap

記述言語：日本語   出版者・発行元：(公社)日本薬学会  

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：AMER CHEMICAL SOC  

Vitamin D receptor (VDR) controls the expression of numerous genes through the conformational change caused by binding 1α,25-dihydroxyvitamin D3. Helix 12 in the ligand-binding domain (LBD) is key to regulating VDR activation. The structures of apo VDR-LBD and the VDR-LBD/antagonist complex are unclear. Here, we reveal their unprecedented structures in solution using a hybrid method combining small-angle X-ray scattering and molecular dynamics simulations. In apo rat VDR-LBD, helix 12 is partially unraveled, and it is positioned around the canonical active position and fluctuates. Helix 11 greatly bends toward the outside at Q396, creating a kink. In the rat VDR-LBD/antagonist complex, helix 12 does not generate the activation function 2 surface, and loop 11-12 is remarkably flexible compared to that in the apo rat VDR-LBD. On the basis of these structural insights, we propose a "folding-door model" to describe the mechanism of agonism/antagonism of VDR-LBD.

DOI： 10.1021/acs.jmedchem.6b00682

Web of Science

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：AMER INST PHYSICS  

The hydration free energy (HFE) is a crucially important physical quantity to discuss various chemical processes in aqueous solutions. Although an explicit-solvent computation with molecular dynamics (MD) simulations is a preferable treatment of the HFE, huge computational load has been inevitable for large, complex solutes like proteins. In the present paper, we propose an efficient computation method for the HFE. In our method, the HFE is computed as a sum of 〈UUV〉/2 (〈UUV〉 is the ensemble average of the sum of pair interaction energy between solute and water molecule) and the water reorganization term mainly reflecting the excluded volume effect. Since 〈UUV〉 can readily be computed through a MD of the system composed of solute and water, an efficient computation of the latter term leads to a reduction of computational load. We demonstrate that the water reorganization term can quantitatively be calculated using the morphometric approach (MA) which expresses the term as the linear combinations of the four geometric measures of a solute and the corresponding coefficients determined with the energy representation (ER) method. Since the MA enables us to finish the computation of the solvent reorganization term in less than 0.1 s once the coefficients are determined, the use of the MA enables us to provide an efficient computation of the HFE even for large, complex solutes. Through the applications, we find that our method has almost the same quantitative performance as the ER method with substantial reduction of the computational load.

DOI： 10.1063/1.4919636

Web of Science

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：AMER CHEMICAL SOC  

The finite-size effect in periodic system is examined for the charging free energy of protein in explicit solvent over a variety of charged states. The key to the finite-size correction is the self-energy, which is defined as the interaction energy of the solute with its own periodic images and the neutralizing background. By employing the thermodynamic-integration method with systematically varied sizes of the unit cell of molecular dynamics (MD) simulations, we show for ubiquitin that the self-energy corrects the finite-size effect on the charging free energy within 1 kcal/mol at total charges of -5e, -1e, neutral, and +1e and within 5 kcal/mol even for a highly charged state with +8e. We then sought the additional correction from the solvation effect using the numerical solution to the Poisson equation of the protein with implicit solvent. This correction reduces the cell-size dependence of the charging free energy at +8e to 3 kcal/mol and is well expressed as the self-energy divided by the dielectric constant of solvent water.

DOI： 10.1021/ct5008394

Web of Science

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：ELSEVIER SCIENCE BV  

On the basis of the free energy landscape theory, we develop a framework to calculate the structural relaxation time in supercooled liquids and glasses. By the framework, the relaxation time is obtained by an escaping time from a basin in a given free energy surface. In order to demonstrate its usefulness, we apply the framework to monodisperse hard-sphere glass systems. Then we show that the relaxation time increases drastically with the density. Additionally, we discuss an explicit picture of the cooperatively rearranging region by analyzing the spatial distribution of an activation free energy of one particle. (C) 2013 Elsevier B.V. All rights reserved.

DOI： 10.1016/j.cplett.2013.05.051

Web of Science

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：ELSEVIER SCIENCE BV  

We present a phenomenological understanding of the fast and slow relaxations in super cooled liquids on the basis of the free energy landscape (FEL). The dynamics of a representative point in the FEL is assumed to obey the time-dependent Ginzburg-Landau type (TDGL) equation with the FEL as the driving potential. We analyze the dynamical structure factor for a toy model and demonstrate that the fast and slow relaxations are determined by jump motion among basins and stochastic motion within a basin, respectively and that the relaxation time of the fast process is determined by the curvature of the FEL around a local minimum. (C) 2007 Elsevier B.V. All rights reserved.

DOI： 10.1016/j.jnoncrysol.2007.02.064

Web of Science

researchmap

J-GLOBAL

researchmap

J-GLOBAL

researchmap

記述言語：英語   掲載種別：研究発表ペーパー・要旨（国際会議）   出版者・発行元：CELL PRESS  

DOI： 10.1016/j.bpj.2015.11.1767

Web of Science

researchmap

J-GLOBAL

researchmap

J-GLOBAL

researchmap

researchmap

記述言語：英語   出版者・発行元：一般社団法人 日本生物物理学会  

DOI： 10.2142/biophys.53.S239_1

researchmap

記述言語：日本語   出版者・発行元：一般社団法人日本物理学会  

CiNii Books

researchmap

記述言語：英語  

DOI： 10.1063/1.2897786

researchmap

記述言語：日本語   出版者・発行元：一般社団法人日本物理学会  

CiNii Books

researchmap