Autoencoder-based detection of the residues involved in G protein-coupled receptor signaling.

Regulator binding and mutations alter protein dynamics. The transmission of the signal of these alterations to distant sites through protein motion results in changes in protein expression and cell function. The detection of residues involved in signal transmission contributes to an elucidation of the mechanisms underlying processes as vast as cellular function and disease pathogenesis.

Mutual information analysis of the dynamic correlation between side chains in proteins.

Protein dynamics play an essential role in function regulation. In recent years, many experimental and theoretical studies have shown that changes in protein fluctuations in the backbone and side chains fulfill a pivotal role associated with amino acid mutations, chemical modifications, and ligand binding. The dynamic correlations between protein side chains have not been sufficiently studied, and no reliable analysis method has been available so far.

Path Ensembles for Pin1-Catalyzed Cis-Trans Isomerization of a Substrate Calculated by Weighted Ensemble Simulations.

Pin1 enzyme protein recognizes specifically phosphorylated serine/threonine (pSer/pThr) and catalyzes the slow interconversion of the peptidyl-prolyl bond between cis and trans forms. Structural dynamics between the cis and trans forms are essential to reveal the underlying molecular mechanism of the catalysis. In this study, we apply the weighted ensemble (WE) simulation method to obtain comprehensive path ensembles for the Pin1-catalyzed isomerization process.

Rebuilding Ring-Type Assembly of Peroxiredoxin by Chemical Modification.

Direct control of the protein quaternary structure (QS) is challenging owing to the complexity of the protein structure. As a protein with a characteristic QS, peroxiredoxin from K1 (ApPrx) forms a decamer, wherein five dimers associate to form a ring. Here, we disrupted and reconstituted ApPrx QS via amino acid mutations and chemical modifications targeting hot spots for protein assembly.

Atomistic detailed free-energy landscape of intrinsically disordered protein studied by multi-scale divide-and-conquer molecular dynamics simulation.

Calcineurin (CaN) is a eukaryotic serine/threonine protein phosphatase activated by both Ca and calmodulin (CaM), including intrinsically disordered region (IDR). The region undergoes folding into an α-helix form in the presence Ca -loaded CaM. To sample the ordered structure of the IDR by conventional all atom model (AAM) molecular dynamics (MD) simulation, the IDR and Ca -loaded CaM must be simultaneously treated.

Epigenetic Protection of Vertebrate Lymphoid Progenitor Cells by Dnmt1.

DNA methylation is a universal epigenetic mechanism involved in regulation of gene expression and genome stability. The DNA maintenance methylase DNMT1 ensures that DNA methylation patterns are faithfully transmitted to daughter cells during cell division. Because loss of DNMT1 is lethal, a pan-organismic analysis of DNMT1 function is lacking.

Mutational effects of Cys113 on structural dynamics of Pin1.

Pin1 is a peptidyl-prolyl isomerase (PPIase) which catalyzes isomerization of pS/pT-P bond. Its activity is related to various cellular functions including suppression of Alzheimer's disease. A cysteine residue C113 is known to be important for its PPIase activity; a mutation C113A reduced the activity by 130-fold.

Loosening of Side-Chain Packing Associated with Perturbations in Peripheral Dynamics Induced by the D76N Mutation of β-Microglobulin Revealed by Pressure-NMR and Molecular Dynamic Simulations.

β-Microglobulin (βm) is the causative protein of dialysis-related amyloidosis, and its D76N variant is less stable and more prone to aggregation. Since their crystal structures are indistinguishable from each other, enhanced amyloidogenicity induced by the mutation may be attributed to changes in the structural dynamics of the molecule. We examined pressure and mutation effects on the βm molecule by NMR and MD simulations, and found that the mutation induced the loosening of the inter-sheet packing of βm, which is relevant to destabilization and subsequent amyloidogenicity.

Autoencoder-Based Detection of Dynamic Allostery Triggered by Ligand Binding Based on Molecular Dynamics.

Dynamic allostery on proteins, triggered by regulator binding or chemical modifications, transmits information from the binding site to distant regions, dramatically altering protein function. It is accompanied by subtle changes in side-chain conformations of the protein, indicating that the changes in dynamics, and not rigid or large conformational changes, are essential to understand regulation of protein function. Although a lot of experimental and theoretical studies have been dedicated to investigate this issue, the regulation mechanism of protein function is still being debated.

Allosteric activation of cytochrome P450 3A4 by efavirenz facilitates midazolam binding.

1. The purpose of this study is to investigate the heteroactivation mechanism of CYP3A4 by efavirenz, which enhances metabolism of midazolam in vivo, in terms of its binding to CYP3A4 with in vitro spectroscopic methods. 2.

On-the-fly analysis of molecular dynamics simulation trajectories of proteins using the Bayesian inference method.

Robust and reliable analyses of long trajectories from molecular dynamics simulations are important for investigations of functions and mechanisms of proteins. Structural fitting is necessary for various analyses of protein dynamics, thus removing time-dependent translational and rotational movements. However, the fitting is often difficult for highly flexible molecules.

Free Energy Reconstruction from Logarithmic Mean-Force Dynamics Using Multiple Nonequilibrium Trajectories.

Efficient and reliable estimation of the mean force (MF), the derivatives of the free energy with respect to a set of collective variables (CVs), has been a challenging problem because free energy differences are often computed by integrating the MF. Among various methods for computing free energy differences, logarithmic mean-force dynamics (LogMFD) [ Morishita et al., Phys.

Substrate recognition of N,N'-diacetylchitobiose deacetylase from Pyrococcus horikoshii.

Enzymes of carbohydrate esterase (CE) family 14 catalyze hydrolysis of N-acetyl groups at the non-reducing end of the N-acetylglucosamine (GlcNAc) residue of chitooligosaccharides or related compounds. N,N'-diacetylchitobiose deacetylase (Dac) belongs to the CE-14 family and plays a role in the chitinolytic pathway in archaea by deacetylating N,N'-diacetylchitobiose (GlcNAc2), which is the end product of chitinase. In this study, we revealed the structural basis of reaction specificity in CE-14 deacetylases by solving a crystal structure of Dac from Pyrococcus horikoshii (Ph-Dac) in complex with a novel reaction intermediate analog.

Development of massive multilevel molecular dynamics simulation program, Platypus (PLATform for dYnamic Protein Unified Simulation), for the elucidation of protein functions.

A massively parallel program for quantum mechanical-molecular mechanical (QM/MM) molecular dynamics simulation, called Platypus (PLATform for dYnamic Protein Unified Simulation), was developed to elucidate protein functions. The speedup and the parallelization ratio of Platypus in the QM and QM/MM calculations were assessed for a bacteriochlorophyll dimer in the photosynthetic reaction center (DIMER) on the K computer, a massively parallel computer achieving 10 PetaFLOPs with 705,024 cores. Platypus exhibited the increase in speedup up to 20,000 core processors at the HF/cc-pVDZ and B3LYP/cc-pVDZ, and up to 10,000 core processors by the CASCI(16,16)/6-31G** calculations.

A method for predicting protein conformational pathways by using molecular dynamics simulations guided by difference distance matrices.

Here, an efficient method that predicts natural transition pathways between two endpoint states of an allosteric protein has been proposed. This method helps create structures that bridge these endpoints through multiple iterative and unbiased molecular dynamics simulations with explicit water. Difference distance matrices provide an approach for identifying states involving concerted slow motion.

Direct Interaction between the Voltage Sensors Produces Cooperative Sustained Deactivation in Voltage-gated H+ Channel Dimers.

The voltage-gated H(+) channel (Hv) is a voltage sensor domain-like protein consisting of four transmembrane segments (S1-S4). The native Hv structure is a homodimer, with the two channel subunits functioning cooperatively. Here we show that the two voltage sensor S4 helices within the dimer directly cooperate via a π-stacking interaction between Trp residues at the middle of each segment.

Molecular dynamics study of the phosphorylation effect on the conformational states of the C-terminal domain of RNA polymerase II.

The carboxyl-terminal domain (CTD) of RNA polymerase II in eukaryotes regulates mRNA processing processes by recruiting various regulation factors. A main function of the CTD relies on the heptad consensus sequence (YSPTSPS). The CTD dynamically changes its conformational state to recognize and bind different regulation factors.

Simple and accurate scheme to compute electrostatic interaction: zero-dipole summation technique for molecular system and application to bulk water.

The zero-dipole summation method was extended to general molecular systems, and then applied to molecular dynamics simulations of an isotropic water system. In our previous paper [I. Fukuda, Y.

A long-range electrostatic potential based on the Wolf method charge-neutral condition.

Molecular simulations rely heavily on a long range electrostatic Coulomb interaction. The Coulomb potential decays inversely with distance, indicating infinite effective range. In practice, molecular simulations do not directly take into account such an infinite interaction.

Free Energy Landscapes of Alanine Dipeptide in Explicit Water Reproduced by the Force-Switching Wolf Method.

Precise and rapid calculation of long-range interactions is of crucial importance for molecular dynamics (MD) and Monte Carlo simulations. Instead of the Ewald method or its high speed variant, PME, we applied our novel method, called the force-switching Wolf method, to computation of the free energy landscapes of a short peptide in explicit water. Wolf and co-workers showed that long-range electrostatic energy under a periodic boundary condition can be well reproduced even by truncating the contribution from the distant charges, when the charge neutrality is taken into account.

Molecular dynamics scheme for precise estimation of electrostatic interaction via zero-dipole summation principle.

We propose a novel idea, zero-dipole summation, for evaluating the electrostatic energy of a classical particle system, and have composed an algorithm for effectively utilizing the idea for molecular dynamics. It conceptually prevents the nonzero-charge and nonzero-dipole states artificially generated by a simple cutoff truncation. The resulting energy formula is nevertheless represented by a simple pairwise function sum, which enables facile application to high-performance computation.

Theory for trivial trajectory parallelization of multicanonical molecular dynamics and application to a polypeptide in water.

Trivial trajectory parallelization of multicanonical molecular dynamics (TTP-McMD) explores the conformational space of a biological system with multiple short runs of McMD starting from various initial structures. This method simply connects (i.e.

Simple and effective application of the Wang-Landau method for multicanonical molecular dynamics simulation.

We propose a novel application of the Wang-Landau method (WLM) for multicanonical molecular dynamics (McMD) simulations. Originally, WLM was developed for Monte Carlo (MC) simulations. Fundamentally, WLM remarkably reduces simulation efforts because it estimates the optimal multicanonical energy function automatically.

Enhanced free-energy calculation using multiscale simulation.

We propose a multiscale simulation method combining the efficiency of a coarse-grained model (CGM) and the accuracy of an all-atom model (AAM) for free-energy landscape calculation of protein systems. A protein's conformation space is quickly searched first using CGM. Then the obtained information is incorporated into AAM simulations.

Membrane attachment facilitates ligand access to the active site in monoamine oxidase A.

Monoamine oxidase membrane enzymes are responsible for the catalytic breakdown of extra- and intracellular neurotransmitters and are targets for the development of central nervous system drugs. We analyzed the dynamics of rat MAOA by performing multiple independent molecular dynamics simulations of membrane-bound and membrane-free forms to clarify the relationship between the mechanics of the enzyme and its function, with particular emphasis on the significance of membrane attachment. Principal component analysis of the simulation trajectories as well as correlations in the fluctuations of the residues pointed to the existence of three domains that define the global dynamics of the protein.