Computational Analysis of Activation of Dimerized Epidermal Growth Factor Receptor Kinase Using the String Method and Markov State Model.

Epidermal growth factor receptor (EGFR) activation is accompanied by dimerization. During the activation of the intracellular kinase domain, two EGFR kinases form an asymmetric dimer, and one side of the dimer (receiver) is activated. Using the string method and Markov state model (MSM), we performed a computational analysis of the structural changes in the activation of the EGFR dimer in this study. The string method reveals the minimum free-energy pathway (MFEP) from the inactive to active structure. The MSM was constructed from numerous trajectories of molecular dynamics simulations around the MFEP, which revealed the free-energy map of structural changes. In the activation of the receiver kinase, the unfolding of the activation loop (A-loop) is followed by the rearrangement of the C-helix, as observed in other kinases. However, unlike other kinases, the free-energy map of EGFR at the asymmetric dimer showed that the active state yielded the highest stability and revealed how interactions at the dimer interface induced receiver activation. As the H-helix of the activator approaches the C-helix of the receiver during activation, the A-loop unfolds. Subsequently, L782 of the receiver enters the pocket between the G- and H-helices of the activator, leading to a rearrangement of the hydrophobic residues around L782 of the receiver, which constitutes a structural rearrangement of the C-helix of the receiver from an outward to an inner position. The MSM analysis revealed long-time scale trajectories via kinetic Monte Carlo.

Comparison of the Molecular Motility of Tubulin Dimeric Isoforms: Molecular Dynamics Simulations and Diffracted X-ray Tracking Study.

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.

Dynamic Solution Structures of Whole Human NAP1 Dimer Bound to One and Two Histone H2A-H2B Heterodimers Obtained by Integrative Methods.

Nucleosome assembly protein 1 (NAP1) binds to histone H2A-H2B heterodimers, mediating their deposition on and eviction from the nucleosome. Human NAP1 (hNAP1) consists of a dimerization core domain and intrinsically disordered C-terminal acidic domain (CTAD), both of which are essential for H2A-H2B binding. Several structures of NAP1 proteins bound to H2A-H2B exhibit binding polymorphisms of the core domain, but the distinct structural roles of the core and CTAD domains remain elusive. Here, we have examined dynamic structures of the full-length hNAP1 dimer bound to one and two H2A-H2B heterodimers by integrative methods. Nuclear magnetic resonance (NMR) spectroscopy of full-length hNAP1 showed CTAD binding to H2A-H2B. Atomic force microscopy revealed that hNAP1 forms oligomers of tandem repeated dimers; therefore, we generated a stable dimeric hNAP1 mutant exhibiting the same H2A-H2B binding affinity as wild-type hNAP1. Size exclusion chromatography (SEC), multi-angle light scattering (MALS) and small angle X-ray scattering (SAXS), followed by modelling and molecular dynamics simulations, have been used to reveal the stepwise dynamic complex structures of hNAP1 binding to one and two H2A-H2B heterodimers. The first H2A-H2B dimer binds mainly to the core domain of hNAP1, while the second H2A-H2B binds dynamically to both CTADs. Based on our findings, we present a model of the eviction of H2A-H2B from nucleosomes by NAP1.

Development of the force field for cyclosporine A.

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.

Hybrid in vitro/in silico analysis of low-affinity protein-protein interactions that regulate signal transduction by Sema6D.

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.

Identification of changes in the microflora composition of Japanese horse mackerel (Trachurus japonicus) during storage to identify specific spoilageorganisms.

Japanese horse mackerel (Trachurus japonicus) is an important marine resource, and its loss and waste should be reduced. This study aimed to identify the changes in the microflora composition during storage and specific spoilage organisms (SSOs) in Japanese horse mackerel, for spoilage prevention. They were stored at either 20 °C or 4 °C aerobically, and the bacterial viable counts, concentration of total volatile basic nitrogen (TVB-N), and microflora composition for each group were analyzed. Samples stored at 20 °C for 48 h showed similar viable counts to those stored at 4 °C for 168 h; however, the TVB-N concentrations increased at 20 °C, but not at 4 °C. 16S rRNA metagenome analysis showed that Shewanella became dominant genus in the microflora regardless of the storage temperature. However, dominant amplicon sequence variants (ASVs), which are a more detailed classification level than the genus, differed depending on the storage temperatures; therefore, dominant ASVs at 20 °C were assumed to be potential SSOs. Shewanella sp. Strain NFH-SH190041, which was genetically closely related to the dominant ASVs at 20 °C, was isolated, and its spoilage ability was verified. The strain NFH-SH190041 may be considered a novel SSO of Japanese horse mackerel because its 16S rRNA sequence is clearly different from those of known species.

3D-RISM-AI: A Machine Learning Approach to Predict Protein-Ligand Binding Affinity Using 3D-RISM.

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.

Supramolecular Mechanosensitive Potassium Channel Formed by Fluorinated Amphiphilic Cyclophane.

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.

Mechanism of Vitamin D Receptor Ligand-Binding Domain Regulation Studied by gREST Simulations.

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.

Elimination of Finite-Size Effects on Binding Free Energies via the Warp-Drive Method.

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.

Ionic scattering factors of atoms that compose biological molecules.

Ionic scattering factors of atoms that compose biological molecules have been computed by the multi-configuration Dirac-Fock method. These ions are chemically unstable and their scattering factors had not been reported except for O-. Yet these factors are required for the estimation of partial charges in protein molecules and nucleic acids. The electron scattering factors of these ions are particularly important as the electron scattering curves vary considerably between neutral and charged atoms in the spatial-resolution range explored in structural biology. The calculated X-ray and electron scattering factors have then been parameterized for the major scattering curve models used in X-ray and electron protein crystallography and single-particle cryo-EM. The X-ray and electron scattering factors and the fitting parameters are presented for future reference.

Energetics and conformational pathways of functional rotation in the multidrug transporter AcrB.

The multidrug transporter AcrB transports a broad range of drugs out of the cell by means of the proton-motive force. The asymmetric crystal structure of trimeric AcrB suggests a functionally rotating mechanism for drug transport. Despite various supportive forms of evidence from biochemical and simulation studies for this mechanism, the link between the functional rotation and proton translocation across the membrane remains elusive. Here, calculating the minimum free energy pathway of the functional rotation for the complete AcrB trimer, we describe the structural and energetic basis behind the coupling between the functional rotation and the proton translocation at atomic resolution. Free energy calculations show that protonation of Asp408 in the transmembrane portion of the drug-bound protomer drives the functional rotation. The conformational pathway identifies vertical shear motions among several transmembrane helices, which regulate alternate access of water in the transmembrane as well as peristaltic motions that pump drugs in the periplasm.

Solution structure of the isolated histone H2A-H2B heterodimer.

During chromatin-regulated processes, the histone H2A-H2B heterodimer functions dynamically in and out of the nucleosome. Although detailed crystal structures of nucleosomes have been established, that of the isolated full-length H2A-H2B heterodimer has remained elusive. Here, we have determined the solution structure of human H2A-H2B by NMR coupled with CS-Rosetta. H2A and H2B each contain a histone fold, comprising four α-helices and two ß-strands (α1-ß1-α2-ß2-α3-αC), together with the long disordered N- and C-terminal H2A tails and the long N-terminal H2B tail. The N-terminal αN helix, C-terminal ß3 strand, and 310 helix of H2A observed in the H2A-H2B nucleosome structure are disordered in isolated H2A-H2B. In addition, the H2A α1 and H2B αC helices are not well fixed in the heterodimer, and the H2A and H2B tails are not completely random coils. Comparison of hydrogen-deuterium exchange, fast hydrogen exchange, and {(1)H}-(15)N hetero-nuclear NOE data with the CS-Rosetta structure indicates that there is some conformation in the H2A 310 helical and H2B Lys11 regions, while the repression domain of H2B (residues 27-34) exhibits an extended string-like structure. This first structure of the isolated H2A-H2B heterodimer provides insight into its dynamic functions in chromatin.

Functional rotation induced by alternating protonation states in the multidrug transporter AcrB: all-atom molecular dynamics simulations.

The multidrug transporter AcrB actively exports a wide variety of noxious compounds using proton-motive force as an energy source in Gram-negative bacteria. AcrB adopts an asymmetric structure comprising three protomers with different conformations that are sequentially converted during drug export; these cyclic conformational changes during drug export are referred to as functional rotation. To investigate functional rotation driven by proton-motive force, all-atom molecular dynamics simulations were performed. Using different protonation states for the titratable residues in the middle of the transmembrane domain, our simulations revealed the correlation between the specific protonation states and the side-chain configurations. Changing the protonation state for Asp408 induced a spontaneous structural transition, which suggests that the proton translocation stoichiometry may be one proton per functional rotation cycle. Furthermore, our simulations demonstrate that alternating the protonation states in the transmembrane domain induces functional rotation in the porter domain, which is primarily responsible for drug transport.

Side-chain conformational changes of transcription factor PhoB upon DNA binding: a population-shift mechanism.

Using molecular dynamics (MD) simulations and analyses of NMR relaxation order parameters, we investigated conformational changes of side chains in hydrophobic cores upon DNA binding for the DNA binding/transactivation domain of the transcription factor PhoB, in which backbone conformational changes upon DNA binding are small. The simulation results correlated well with experimental order parameters for the backbone and side-chain methyl groups, showing that the order parameters generally represent positional fluctuations of the backbone and side-chain methyl groups. However, topological effects of the side chains on the order parameters were also found and could be eliminated using normalized order parameters for each amino acid type. Consistent with the NMR experiments, the normalized order parameters from the MD simulations showed that the side chains in one of the two hydrophobic cores (the soft core) were highly flexible in comparison with those in the other hydrophobic core (the hard core) before DNA binding and that the flexibility of the hydrophobic cores, particularly of the soft core, was reduced upon DNA binding. Principal component analysis of methyl group configurations revealed strikingly different side-chain dynamics for the soft and hard cores. In the hard core, side-chain configurations were simply distributed around one or two average configurations. In contrast, the side chains in the soft core dynamically varied their configurations in an equilibrium ensemble that included binding configurations as minor components before DNA binding. DNA binding led to a restriction of the side-chain dynamics and a shift in the equilibrium toward binding configurations, in clear correspondence with a population-shift model.

Structure and functional characterization of Vibrio parahaemolyticus thermostable direct hemolysin.

Thermostable direct hemolysin (TDH) is a major virulence factor of Vibrio parahaemolyticus that causes pandemic foodborne enterocolitis mediated by seafood. TDH exists as a tetramer in solution, and it possesses extreme hemolytic activity. Here, we present the crystal structure of the TDH tetramer at 1.5 A resolution. The TDH tetramer forms a central pore with dimensions of 23 A in diameter and approximately 50 A in depth. Pi-cation interactions between protomers comprising the tetramer were indispensable for hemolytic activity of TDH. The N-terminal region was intrinsically disordered outside of the pore. Molecular dynamic simulations suggested that water molecules permeate freely through the central and side channel pores. Electron micrographs showed that tetrameric TDH attached to liposomes, and some of the tetramer associated with liposome via one protomer. These findings imply a novel membrane attachment mechanism by a soluble tetrameric pore-forming toxin.

Structure of the small ubiquitin-like modifier (SUMO)-interacting motif of MBD1-containing chromatin-associated factor 1 bound to SUMO-3.

Post-translational modification by small ubiquitin-like modifier (SUMO) proteins has been implicated in the regulation of a variety of cellular events. The functions of sumoylation are often mediated by downstream effector proteins harboring SUMO-interacting motifs (SIMs) that are composed of a hydrophobic core and a stretch of acidic residues. MBD1-containing chromatin-associated factor 1 (MCAF1), a transcription repressor, interacts with SUMO-2/3 and SUMO-1, with a preference for SUMO-2/3. We used NMR spectroscopy to solve the solution structure of the SIM of MCAF1 bound to SUMO-3. The hydrophobic core of the SIM forms a parallel beta-sheet pairing with strand beta2 of SUMO-3, whereas its C-terminal acidic stretch seems to mediate electrostatic interactions with a surface area formed by basic residues of SUMO-3. The significance of these electrostatic interactions was shown by mutations of both SUMO-3 and MCAF1. The present structural and biochemical data suggest that the acidic stretch of the SIM of MCAF1 plays an important role in the binding to SUMO-3.

Water-mediated interactions between DNA and PhoB DNA-binding/transactivation domain: NMR-restrained molecular dynamics in explicit water environment.

The solution structure of the complex between the transcription factor PhoB DNA-binding/transactivation domain and DNA was determined by NMR spectroscopy and simulated annealing in a periodic boundary box of explicit water with the particle mesh Ewald method. The refined structures provided better convergence and better local geometry compared with the structures determined in vacuum. The hydrogen bond interactions between the PhoB domain and DNA in the aqueous environment were fully formed. The complex structure was found to be very similar to the crystal structure, particularly at the PhoB-DNA interface, much more so than expected from the vacuum structure. These results indicate the importance of the proper treatment of electrostatic and hydration influences in describing protein-DNA interactions. The hydration structures observed for the refined structures contained most of the crystal waters as a subset. We observed that various water-mediated PhoB-DNA interactions contributed to the molecular recognition between PhoB and DNA.

Uniaxial alignment of the smallest diamagnetic susceptibility axis using time-dependent magnetic fields.

A diamagnetic particle with magnetic susceptibilities chi3 < chi2 = chi1 < 0 was subjected to a rotating magnetic field to obtain an alignment of the chi3 axis (the smallest susceptibility axis) in the direction perpendicular to the plane of the rotating magnetic field. A polymer short fiber, whose fiber axis coincides with the chi3 axis, was suspended in a fluid with the same density, and then a rotating magnetic field generated by a rotation of a pair of permanent magnets was applied. The fiber axis, rotating following the applied field, finally ended up with an alignment perpendicular to the plane of the rotating magnetic field. The experimental data on the time course of the alignment was in good agreement with the numerical calculation based on the equation of rotation.