Publishing type：Research paper (scientific journal)   Publisher：Springer Science and Business Media LLC  

DOI： 10.1038/s41467-024-53533-6

researchmap

Other Link： https://www.nature.com/articles/s41467-024-53533-6

Publishing type：Research paper (scientific journal)   Publisher：American Chemical Society (ACS)  

DOI： 10.1021/jacs.4c04648

researchmap

Publishing type：Research paper (scientific journal)   Publisher：American Chemical Society (ACS)  

DOI： 10.1021/acs.jcim.4c00172

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Hydroxycarboxylic acid receptors (HCAR1, HCAR2, and HCAR3) transduce Gi/o signaling upon biding to molecules such as lactic acid, butyric acid and 3-hydroxyoctanoic acid, which are associated with lipolytic and atherogenic activity, and neuroinflammation. Although many reports have elucidated the function of HCAR2 and its potential as a therapeutic target for treating not only dyslipidemia but also neuroimmune disorders such as multiple sclerosis and Parkinson's disease, the structural basis of ligand recognition and ligand-induced Gi-coupling remains unclear. Here we report three cryo-EM structures of the human HCAR2-Gi signaling complex, each bound with different ligands: niacin, acipimox or GSK256073. All three agonists are held in a deep pocket lined by residues that are not conserved in HCAR1 and HCAR3. A distinct hairpin loop at the HCAR2 N-terminus and extra-cellular loop 2 (ECL2) completely enclose the ligand. These structures also reveal the agonist-induced conformational changes propagated to the G-protein-coupling interface during activation. Collectively, the structures presented here are expected to help in the design of ligands specific for HCAR2, leading to new drugs for the treatment of various diseases such as dyslipidemia and inflammation.

DOI： 10.1038/s41467-023-42764-8

PubMed

researchmap

Publishing type：Research paper (scientific journal)   Publisher：MDPI AG  

Tubulin has been recently reported to form a large family consisting of various gene isoforms; however, the differences in the molecular features of tubulin dimers composed of a combination of these isoforms remain unknown. Therefore, we attempted to elucidate the physical differences in the molecular motility of these tubulin dimers using the method of measurable pico-meter-scale molecular motility, diffracted X-ray tracking (DXT) analysis, regarding characteristic tubulin dimers, including neuronal TUBB3 and ubiquitous TUBB5. We first conducted a DXT analysis of neuronal (TUBB3-TUBA1A) and ubiquitous (TUBB5-TUBA1B) tubulin dimers and found that the molecular motility around the vertical axis of the neuronal tubulin dimer was lower than that of the ubiquitous tubulin dimer. The results of molecular dynamics (MD) simulation suggest that the difference in motility between the neuronal and ubiquitous tubulin dimers was probably caused by a change in the major contact of Gln245 in the T7 loop of TUBB from Glu11 in TUBA to Val353 in TUBB. The present study is the first report of a novel phenomenon in which the pico-meter-scale molecular motility between neuronal and ubiquitous tubulin dimers is different.

DOI： 10.3390/ijms242015423

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

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

Publishing type：Research paper (scientific journal)   Publisher：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

Other Link： https://www.nature.com/articles/s41467-023-35844-2

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English  

DOI： 10.1021/jacs.2c07339

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)  

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

Publishing type：Research paper (scientific journal)   Publisher：Wiley  

DOI： 10.1002/asia.202001415

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)   Publisher：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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)   Publisher：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

Language：English   Publishing type：Research paper (scientific journal)  

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

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Springer Verlag  

DOI： 10.1007/s12551-017-0373-4

Scopus

PubMed

researchmap

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

PubMed

researchmap

Language：Japanese   Publisher：THE VITAMIN SOCIETY OF JAPAN  

DOI： 10.20632/vso.92.1_18

CiNii Books

CiNii Research

researchmap

Other Link： https://ndlsearch.ndl.go.jp/books/R000000004-I028798270

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1016/j.bpj.2017.01.029

Web of Science

PubMed

researchmap

Language：Japanese   Publisher：(公社)日本薬学会  

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1021/acs.jmedchem.6b00682

Web of Science

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1063/1.4919636

Web of Science

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1021/ct5008394

Web of Science

PubMed

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1016/j.cplett.2013.05.051

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

DOI： 10.1016/j.jnoncrysol.2007.02.064

Web of Science

researchmap

Language：English   Publishing type：Research paper, summary (international conference)  

Web of Science

researchmap

J-GLOBAL

researchmap

J-GLOBAL

researchmap

Language：English   Publishing type：Research paper, summary (international conference)  

DOI： 10.1016/j.bpj.2015.11.1767

Web of Science

researchmap

J-GLOBAL

researchmap

J-GLOBAL

researchmap

researchmap

Language：English   Publisher：The Biophysical Society of Japan General Incorporated Association  

DOI： 10.2142/biophys.53.S239_1

researchmap

Language：Japanese   Publisher：The Physical Society of Japan (JPS)  

CiNii Books

researchmap

Language：English  

DOI： 10.1063/1.2897786

researchmap

Language：Japanese   Publisher：The Physical Society of Japan (JPS)  

CiNii Books

researchmap