Harnessing autoimmunity with dominant self-peptide: Modulating the sustainability of tissue-preferential antigen-specific Tregs by governing the binding stability via peptide flanking residues.

Sensitization to self-peptides induces various immunological responses, from autoimmunity to tumor immunity, depending on the peptide sequence; however, the underlying mechanisms remain unclear, and thus, curative therapeutic options considering immunity balance are limited. Herein, two overlapping dominant peptides of myelin proteolipid protein, PLP136-150 and PLP139-151, which induce different forms of experimental autoimmune encephalomyelitis (EAE), monophasic and relapsing EAE, respectively, were investigated. Mice with monophasic EAE exhibited highly resistant to EAE re-induction with any encephalitogenic peptides, whereas mice with relapsing EAE were susceptible, and progressed, to EAE re-induction. This resistance to relapse and re-induction in monophasic EAE mice was associated with the maintenance of potent CD69+CD103+CD4+CD25high regulatory T-cells (Tregs) enriched with antigen specificity, which expanded preferentially in the central nervous system with sustained suppressive activity. This tissue-preferential sustainability of potent antigen-specific Tregs was correlated with the antigenicity of PLP136-150, depending on its flanking residues. That is, the flanking residues of PLP136-150 enable to form pivotally arranged strong hydrogen bonds that secured its binding stability to MHC-class II. These potent Tregs acting tissue-preferentially were induced only by sensitization of PLP136-150, not by its tolerance induction, independent of EAE development. These findings suggest that, for optimal therapy, "benign autoimmunity" can be critically achieved through inverse vaccination with self-peptides by manipulating their flanking residues.

Extended ensemble simulations of a SARS-CoV-2 nsp1-5'-UTR complex.

Nonstructural protein 1 (nsp1) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a 180-residue protein that blocks translation of host mRNAs in SARS-CoV-2-infected cells. Although it is known that SARS-CoV-2's own RNA evades nsp1's host translation shutoff, the molecular mechanism underlying the evasion was poorly understood. We performed an extended ensemble molecular dynamics simulation to investigate the mechanism of the viral RNA evasion. Simulation results suggested that the stem loop structure of the SARS-CoV-2 RNA 5'-untranslated region (SL1) binds to both nsp1's N-terminal globular region and intrinsically disordered region. The consistency of the results was assessed by modeling nsp1-40S ribosome structure based on reported nsp1 experiments, including the X-ray crystallographic structure analysis, the cryo-EM electron density map, and cross-linking experiments. The SL1 binding region predicted from the simulation was open to the solvent, yet the ribosome could interact with SL1. Cluster analysis of the binding mode and detailed analysis of the binding poses suggest residues Arg124, Lys47, Arg43, and Asn126 may be involved in the SL1 recognition mechanism, consistent with the existing mutational analysis.

Critical Roles for Coiled-Coil Dimers of Butyrophilin 3A1 in the Sensing of Prenyl Pyrophosphates by Human Vγ2Vδ2 T Cells.

Vγ2Vδ2 T cells play important roles in human immunity to pathogens and tumors. Their TCRs respond to the sensing of isoprenoid metabolites, such as (E)-4-hydroxy-3-methyl-but-2-enyl pyrophosphate and isopentenyl pyrophosphate, by butyrophilin (BTN) 3A1. BTN3A1 is an Ig superfamily protein with extracellular IgV/IgC domains and intracellular B30.2 domains that bind prenyl pyrophosphates. We have proposed that intracellular α helices form a coiled-coil dimer that functions as a spacer for the B30.2 domains. To test this, five pairs of anchor residues were mutated to glycine to destabilize the coiled-coil dimer. Despite maintaining surface expression, BTN3A1 mutagenesis either abrogated or decreased stimulation by (E)-4-hydroxy-3-methyl-but-2-enyl pyrophosphate. BTN3A2 and BTN3A3 proteins and orthologs in alpacas and dolphins are also predicted to have similar coiled-coil dimers. A second short coiled-coil region dimerizes the B30.2 domains. Molecular dynamics simulations predict that mutation of a conserved tryptophan residue in this region will destabilize the dimer, explaining the loss of stimulation by BTN3A1 proteins with this mutation. The juxtamembrane regions of other BTN/BTN-like proteins with B30.2 domains are similarly predicted to assume α helices, with many predicted to form coiled-coil dimers. An exon at the end of this region and the exon encoding the dimerization region for B30.2 domains are highly conserved. We propose that coiled-coil dimers function as rod-like helical molecular spacers to position B30.2 domains, as interaction sites for other proteins, and as dimerization regions to allow sensing by B30.2 domains. In these ways, the coiled-coil domains of BTN3A1 play critical roles for its function.

Correction: Free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA.

[This corrects the article DOI: 10.1371/journal.pcbi.1006024.].

Large-Scale Parallel Implementation of Hartree-Fock Exchange Energy on Real-Space Grids Using 3D-Parallel Fast Fourier Transform.

Two types of implementation of the Hartree-Fock (HF) exchange energy were developed with the real-space grid approach for the purpose of achieving high efficiency in the parallel execution of the hybrid exchange functional in the density functional theory. First, a parallel implementation of the three-dimensional fast Fourier transform (FFT), referred to as PFFT, was adapted to solve the Poisson equations for the electrostatic potentials of the densities of the orbital pairs. In the other approach, the Poisson equations were solved iteratively through the conjugate gradient (CG) procedures where the operation of Laplacian was parallelized by the domain decomposition scheme. Comparison of the parallel performances for the exchange energy calculation was made between these two approaches, and it was revealed that the calculation with the FFT method is faster than that with CG. The method with FFT is more advantageous than CG because a larger bandwidth can be made available in the collective message passing interface communication associated with the parallel execution of FFT. We also implemented the projection operator to circumvent the laborious calculation of the exchange energy at every self-consistent field step, which made a significant contribution to expedite the convergence. To assess the accuracy of our implementation, the association energies of a hydrated ion were computed, which showed excellent agreement with those given by the Gaussian 09 program employing sophisticated basis sets.

Free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA.

The eukaryotic genome is packaged into a nucleus in the form of chromatin. The fundamental structural unit of chromatin is a protein-DNA complex, the nucleosome, where 146 or 147 base pairs of DNA wrap 1.75 times around a histone core. To function in cellular processes, however, nucleosomal DNA must be unwrapped. Although this unwrapping has been experimentally investigated, details of the process at an atomic level are not yet well understood. Here, we used molecular dynamics simulation with an enhanced sampling method to calculate the free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA. A free energy change of about 11.5 kcal/mol for the unwrapping agrees well with values obtained in single molecule experiments. This simulation revealed a variety of conformational states, indicating there are many potential paths to outer superhelicdal turn unwrapping, but the dominant path is likely asymmetric. At one end of the DNA, the first five bps unwrap, after which a second five bps unwrap at the same end with no increase in free energy. The unwrapping then starts at the other end of the DNA, where 10 bps are unwrapped. During further unwrapping of 15 bps, the unwrapping advances at one of the ends, after which the other end of the DNA unwraps to complete the unwrapping of the outer superhelical turn. These results provide insight into the construction, disruption, and repositioning of nucleosomes, which are continuously ongoing during cellular processes.

Capturing alternative secondary structures of RNA by decomposition of base-pairing probabilities.

BACKGROUND: It is known that functional RNAs often switch their functions by forming different secondary structures. Popular tools for RNA secondary structures prediction, however, predict the single 'best' structures, and do not produce alternative structures. There are bioinformatics tools to predict suboptimal structures, but it is difficult to detect which alternative secondary structures are essential. RESULTS: We proposed a new computational method to detect essential alternative secondary structures from RNA sequences by decomposing the base-pairing probability matrix. The decomposition is calculated by a newly implemented software tool, RintW, which efficiently computes the base-pairing probability distributions over the Hamming distance from arbitrary reference secondary structures. The proposed approach has been demonstrated on ROSE element RNA thermometer sequence and Lysine RNA ribo-switch, showing that the proposed approach captures conformational changes in secondary structures. CONCLUSIONS: We have shown that alternative secondary structures are captured by decomposing base-paring probabilities over Hamming distance. Source code is available from http://www.ncRNA.org/RintW .

Performance evaluation of the zero-multipole summation method in modern molecular dynamics software.

The zero-multiple summation method (ZMM) is a cutoff-based method for calculating electrostatic interactions in molecular dynamics simulations, utilizing an electrostatic neutralization principle as a physical basis. Since the accuracies of the ZMM have been revealed to be sufficient in previous studies, it is highly desirable to clarify its practical performance. In this paper, the performance of the ZMM is compared with that of the smooth particle mesh Ewald method (SPME), where the both methods are implemented in molecular dynamics software package GROMACS. Extensive performance comparisons against a highly optimized, parameter-tuned SPME implementation are performed for various-sized water systems and two protein-water systems. We analyze in detail the dependence of the performance on the potential parameters and the number of CPU cores. Even though the ZMM uses a larger cutoff distance than the SPME does, the performance of the ZMM is comparable to or better than that of the SPME. This is because the ZMM does not require a time-consuming electrostatic convolution and because the ZMM gains short neighbor-list distances due to the smooth damping feature of the pairwise potential function near the cutoff length. We found, in particular, that the ZMM with quadrupole or octupole cancellation and no damping factor is an excellent candidate for the fast calculation of electrostatic interactions. © 2018 Wiley Periodicals, Inc.

Analysis of O2-binding Sites in Proteins Using Gas-Pressure NMR Spectroscopy: Outer Surface Protein A.

Internal cavities in proteins produce conformational fluctuations and enable the binding of small ligands. Here, we report a NMR analysis of O2-binding sites by O2-induced paramagnetic relaxation enhancements (PREs) on amide groups of proteins in solution. Outer surface protein A contains a nonglobular single-layer ß-sheet that connects the N- and C-terminal globular domains. Several cavities have been observed in both domains of the crystallized protein structure. The receptor-binding sites are occluded and line the largest cavity of the C-terminal domain. We observed significant O2-induced PREs for amide protons located around the largest cavity and at the central ß-sheet. We suggested three potential O2-accessible sites in the protein based on the 1/r6 distance dependence of the PRE. Two sites were in or close to the largest cavity and the third site was in the surface crevice of the central ß-sheet. These results provide, to our knowledge, the first evidence of ligand binding to the surface crevice and cavity of the protein in solution. Because O2 generally binds more specifically to hydrophobic rather than hydrophilic cavities within a protein, the results also indicated that the receptor-binding sites lining the largest cavity were in the hydrophobic environment in the ground-state conformation. Molecular dynamics simulations permitted the visualization of the rotational and translational motions of O2 within the largest cavity, egress of O2 from the cavity, and ingress of O2 in the surface crevice of the ß-sheet. These molecular dynamics simulation results qualitatively explained the O2-induced changes in NMR observations. Exploring cavities that are sufficiently dynamic to enable access by small molecules can be a useful strategy for the design of stable proteins and their ligands.

H3 Histone Tail Conformation within the Nucleosome and the Impact of K14 Acetylation Studied Using Enhanced Sampling Simulation.

Acetylation of lysine residues in histone tails is associated with gene transcription. Because histone tails are structurally flexible and intrinsically disordered, it is difficult to experimentally determine the tail conformations and the impact of acetylation. In this work, we performed simulations to sample H3 tail conformations with and without acetylation. The results show that irrespective of the presence or absence of the acetylation, the H3 tail remains in contact with the DNA and assumes an α-helix structure in some regions. Acetylation slightly weakened the interaction between the tail and DNA and enhanced α-helix formation, resulting in a more compact tail conformation. We inferred that this compaction induces unwrapping and exposure of the linker DNA, enabling DNA-binding proteins (e.g., transcription factors) to bind to their target sequences. In addition, our simulation also showed that acetylated lysine was more often exposed to the solvent, which is consistent with the fact that acetylation functions as a post-translational modification recognition site marker.

Comment on "Replica-exchange-with-tunneling for fast exploration of protein landscapes" [J. Chem. Phys. 143, 224102 (2015)].

In "Replica-exchange-with-tunneling for fast exploration of protein landscapes" [F. Yasar et al., J. Chem. Phys. 143, 224102 (2015)], a novel sampling algorithm called "Replica Exchange with Tunneling" was proposed. However, due to its violation of the detailed balance, the algorithm fails to sample from the correct canonical ensemble.

Spotting the difference in molecular dynamics simulations of biomolecules.

Comparing two trajectories from molecular simulations conducted under different conditions is not a trivial task. In this study, we apply a method called Linear Discriminant Analysis with ITERative procedure (LDA-ITER) to compare two molecular simulation results by finding the appropriate projection vectors. Because LDA-ITER attempts to determine a projection such that the projections of the two trajectories do not overlap, the comparison does not suffer from a strong anisotropy, which is an issue in protein dynamics. LDA-ITER is applied to two test cases: the T4 lysozyme protein simulation with or without a point mutation and the allosteric protein PDZ2 domain of hPTP1E with or without a ligand. The projection determined by the method agrees with the experimental data and previous simulations. The proposed procedure, which complements existing methods, is a versatile analytical method that is specialized to find the "difference" between two trajectories.

Interaction-component analysis of the hydration and urea effects on cytochrome c.

Energetics was analyzed for cytochrome c in pure-water solvent and in a urea-water mixed solvent to elucidate the solvation effect in the structural variation of the protein. The solvation free energy was computed through all-atom molecular dynamics simulation combined with the solution theory in the energy representation, and its correlations were examined over sets of protein structures against the electrostatic and van der Waals components in the average interaction energy of the protein with the solvent and the excluded-volume component in the solvation free energy. It was observed in pure-water solvent that the solvation free energy varies in parallel to the electrostatic component with minor roles played by the van der Waals and excluded-volume components. The effect of urea on protein structure was then investigated in terms of the free-energy change upon transfer of the protein solute from pure-water solvent to the urea-water mixed solvent. The decomposition of the transfer free energy into the contributions from urea and water showed that the urea contribution is partially canceled by the water contribution and governs the total free energy of transfer. When correlated against the change in the solute-solvent interaction energy upon transfer and the corresponding changes in the electrostatic, van der Waals, and excluded-volume components, the transfer free energy exhibited strong correlations with the total change in the solute-solvent energy and its van der Waals component. The solute-solvent energy was decomposed into the contributions from the protein backbone and side chain, furthermore, and neither of the contributions was seen to be decisive in the correlation to the transfer free energy.

Ermod: fast and versatile computation software for solvation free energy with approximate theory of solutions.

ERmod is a software package to efficiently and approximately compute the solvation free energy using the method of energy representation. Molecular simulation is to be conducted at two condensed-phase systems of the solution of interest and the reference solvent with test-particle insertion of the solute. The subprogram ermod in ERmod then provides a set of energy distribution functions from the simulation trajectories, and another subprogram slvfe determines the solvation free energy from the distribution functions through an approximate functional. This article describes the design and implementation of ERmod, and illustrates its performance in solvent water for two organic solutes and two protein solutes. Actually, the free-energy computation with ERmod is not restricted to the solvation in homogeneous medium such as fluid and polymer and can treat the binding into weakly ordered system with nano-inhomogeneity such as micelle and lipid membrane. ERmod is available on web at http://sourceforge.net/projects/ermod.

Adaptive lambda square dynamics simulation: an efficient conformational sampling method for biomolecules.

A novel, efficient sampling method for biomolecules is proposed. The partial multicanonical molecular dynamics (McMD) was recently developed as a method that improved generalized ensemble (GE) methods to focus sampling only on a part of a system (GEPS); however, it was not tested well. We found that partial McMD did not work well for polylysine decapeptide and gave significantly worse sampling efficiency than a conventional GE. Herein, we elucidate the fundamental reason for this and propose a novel GEPS, adaptive lambda square dynamics (ALSD), which can resolve the problem faced when using partial McMD. We demonstrate that ALSD greatly increases the sampling efficiency over a conventional GE. We believe that ALSD is an effective method and is applicable to the conformational sampling of larger and more complicated biomolecule systems.

Molecular binding sites are located near the interface of intrinsic dynamics domains (IDDs).

We provide evidence supporting that protein-protein and protein-ligand docking poses are functions of protein shape and intrinsic dynamics. Over sets of 68 protein-protein complexes and 240 nonhomologous enzymes, we recognize common predispositions for binding sites to have minimal vibrations and angular momenta, while two interacting proteins orient so as to maximize the angle between their rotation/bending axes (>65°). The findings are then used to define quantitative criteria to filter out docking decoys less likely to be the near-native poses; hence, the chances to find near-native hits can be doubled. With the novel approach to partition a protein into "domains" of robust but disparate intrinsic dynamics, 90% of catalytic residues in enzymes can be found within the first 50% of the residues closest to the interface of these dynamics domains. The results suggest an anisotropic rather than isotropic distribution of catalytic residues near the mass centers of enzymes.

A theoretical study of the two binding modes between lysozyme and tri-NAG with an explicit solvent model based on the fragment molecular orbital method.

To examine the stabilities and binding characteristics, fragment molecular orbital (FMO) calculations were performed for the two binding modes of hen egg-white lysozyme with tri-N-acetyl-D-glucosamine (tri-NAG). Solvent effects were considered using an explicit solvent model. For comparison with the computational results, we experimentally determined the enthalpic contribution of the binding free-energy. Our calculations showed that the binding mode observed by X-ray analysis was more stable than the other binding mode by -6.2 kcal mol(-1), where it was found that the interaction of protein with solvent molecules was crucial for this stability. The amplitude of this energy difference was of the same order as the experimental enthalpic contribution. Our detailed analysis using the energies divided into each residue was also consistent with a previous mutant study. In addition, the electron density analysis showed that the formal charge of the lysozyme (+8.0 e) was reduced to +5.16 e by charge transfer with solvent molecules.

Interaction of naphthalene derivatives with lipids in membranes studied by the 1H-nuclear Overhauser effect and molecular dynamics simulation.

The location, orientation, and dynamics of hydrophobic small molecules in lipid membranes are studied through combined use of solution-state (1)H-NMR and MD simulation. 1-Naphthol and 1-methylnaphthalene were adopted as the small molecules with or without hydrophilic groups. The nuclear Overhauser effect (NOE) measurement was performed for large unilamellar vesicles (100 nm in diameter) composed of dimyristoylphosphatidylcholine (DMPC) and the naphthalene derivative. The transient NOE-SE (spin-echo) scheme previously reported (J. Phys. Chem. B, 2011, 115, 9106-9115) was employed to quantitatively determine the NOE cross relaxation rate constant between DMPC and the naphthalene derivative. The observed NOE shows that both the naphthalene derivatives distribute over a wide domain across the normal of the essentially planar membrane ranging from the hydrophobic core to the hydrophilic headgroup. The experimental NOE information was further refined in combination with the analysis of time correlation functions in MD simulation. It was found that 1-naphthol exhibits a slight preference for pointing its OH group toward the hydrophilic domain of the membrane and that no definite preference can be concluded for the orientation of 1-methylnaphthalene. When 1-naphthol and 1-methylnaphthalene are compared, the NOE is stronger for 1-naphthol due to the restricted motion of the OH group. The slowdown of the 1-naphthol motion is also evidenced by the (1)H spectral line width.

Evaluation of protein-protein docking model structures using all-atom molecular dynamics simulations combined with the solution theory in the energy representation.

We propose a method to evaluate binding free energy differences among distinct protein-protein complex model structures through all-atom molecular dynamics simulations in explicit water using the solution theory in the energy representation. Complex model structures are generated from a pair of monomeric structures using the rigid-body docking program ZDOCK. After structure refinement by side chain optimization and all-atom molecular dynamics simulations in explicit water, complex models are evaluated based on the sum of their conformational and solvation free energies, the latter calculated from the energy distribution functions obtained from relatively short molecular dynamics simulations of the complex in water and of pure water based on the solution theory in the energy representation. We examined protein-protein complex model structures of two protein-protein complex systems, bovine trypsin/CMTI-1 squash inhibitor (PDB ID: 1PPE) and RNase SA/barstar (PDB ID: 1AY7), for which both complex and monomer structures were determined experimentally. For each system, we calculated the energies for the crystal complex structure and twelve generated model structures including the model most similar to the crystal structure and very different from it. In both systems, the sum of the conformational and solvation free energies tended to be lower for the structure similar to the crystal. We concluded that our energy calculation method is useful for selecting low energy complex models similar to the crystal structure from among a set of generated models.

HNF1A POU Domain Mutations Found in Japanese Liver Cancer Patients Cause Downregulation of HNF4A Promoter Activity with Possible Disruption in Transcription Networks.

Hepatocyte nuclear factor 1A (HNF1A) is the master regulator of liver homeostasis and organogenesis and regulates many aspects of hepatocyte functions. It acts as a tumor suppressor in the liver, evidenced by the increased proliferation in HNF1A knockout (KO) hepatocytes. Hence, we postulated that any loss-of-function variation in the gene structure or composition (mutation) could trigger dysfunction, including disrupted transcriptional networks in liver cells. From the International Cancer Genome Consortium (ICGC) database of cancer genomes, we identified several HNF1A mutations located in the functional Pit-Oct-Unc (POU) domain. In our biochemical analysis, we found that the HNF1A POU-domain mutations Y122C, R229Q and V259F suppressed HNF4A promoter activity and disrupted the binding of HNF1A to its target HNF4A promoter without any effect on the nuclear localization. Our results suggest that the decreased transcriptional activity of HNF1A mutants is due to impaired DNA binding. Through structural simulation analysis, we found that a V259F mutation was likely to affect DNA interaction by inducing large conformational changes in the N-terminal region of HNF1A. The results suggest that POU-domain mutations of HNF1A downregulate HNF4A gene expression. Therefore, to mimic the HNF1A mutation phenotype in transcription networks, we performed siRNA-mediated knockdown (KD) of HNF4A. Through RNA-Seq data analysis for the HNF4A KD, we found 748 differentially expressed genes (DEGs), of which 311 genes were downregulated (e.g., HNF1A, ApoB and SOAT2) and 437 genes were upregulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) mapping revealed that the DEGs were involved in several signaling pathways (e.g., lipid and cholesterol metabolic pathways). Protein-protein network analysis suggested that the downregulated genes were related to lipid and cholesterol metabolism pathways, which are implicated in hepatocellular carcinoma (HCC) development. Our study demonstrates that mutations of HNF1A in the POU domain result in the downregulation of HNF1A target genes, including HNF4A, and this may trigger HCC development through the disruption of HNF4A-HNF1A transcriptional networks.