Distribution and Structure Analysis of Fibril-Forming Peptides Focusing on Concentration Dependency.

We focus on the concentration dependency of fibril-forming peptides, which have the potential of aggregation by themselves. In this study, we performed replica-exchange molecular dynamics simulations of Lys-Phe-Phe-Glu (KFFE) fragments, which are known to form fibrils in experiments under different concentration environments. The analysis by static structure factors suggested that the density fluctuation of the KFFE fragments becomes large as the concentration increases.

Absolute Protein Binding Free Energy Simulations for Ligands with Multiple Poses, a Thermodynamic Path That Avoids Exhaustive Enumeration of the Poses.

We propose a free energy calculation method for receptor-ligand binding, which have multiple binding poses that avoids exhaustive enumeration of the poses. For systems with multiple binding poses, the standard procedure is to enumerate orientations of the binding poses, restrain the ligand to each orientation, and then, calculate the binding free energies for each binding pose. In this study, we modify a part of the thermodynamic cycle in order to sample a broader conformational space of the ligand in the binding site.

Structural Analysis of a Trimer of β-Microgloblin Fragment by Molecular Dynamics Simulations.

A peptide β-m, which is a fragment from residue 21 to residue 31 of β-microgloblin, is experimentally known to self-assemble and form amyloid fibrils. In order to understand the mechanism of amyloid fibril formations, we applied the replica-exchange molecular dynamics method to the system consisting of three fragments of β-m. From the analyses on the temperature dependence, we found that there is a clear phase transition temperature in which the peptides aggregate with each other.

N-C-H···O Hydrogen bonds in protein-ligand complexes.

In the context of drug design, C-H···O hydrogen bonds have received little attention so far, mostly because they are considered weak relative to other noncovalent interactions such as O-H···O hydrogen bonds, π/π interactions, and van der Waals interactions. Herein, we demonstrate the significance of hydrogen bonds between C-H groups adjacent to an ammonium cation and an oxygen atom (N-C-H···O hydrogen bonds) in protein-ligand complexes. Quantum chemical calculations revealed details on the strength and geometrical requirements of these N-C-H···O hydrogen bonds, and a subsequent survey of the Protein Data Bank (PDB) based on these criteria suggested that numerous protein-ligand complexes contain such N-C-H···O hydrogen bonds.

Enhanced sampling method in molecular simulations using genetic algorithm for biomolecular systems.

We propose a molecular simulation method using genetic algorithm (GA) for biomolecular systems to obtain ensemble averages efficiently. In this method, we incorporate the genetic crossover, which is one of the operations of GA, to any simulation method such as conventional molecular dynamics (MD), Monte Carlo, and other simulation methods. The genetic crossover proposes candidate conformations by exchanging parts of conformations of a target molecule between a pair of conformations during the simulation.

Computational analysis for selectivity of histone deacetylase inhibitor by replica-exchange umbrella sampling molecular dynamics simulations.

We performed protein-ligand docking simulations with a ligand T247, which has been reported as a selective inhibitor of a histone deacetylase HDAC3, by the replica-exchange umbrella sampling method in order to estimate the free energy profiles along ligand docking pathways of HDAC3-T247 and HDAC2-T247 systems. The simulation results showed that the docked state of the HDAC3-T247 system is more stable than that of the HDAC2-T247 system although the amino-acid sequences and structures of HDAC3 and HDAC2 are very similar. By comparing structures obtained from the simulations of both systems, we found the difference between structures of hydrophobic residues at the entrance of the catalytic site.

Structure-function insights into direct lipid transfer between membranes by Mmm1-Mdm12 of ERMES.

The endoplasmic reticulum (ER)-mitochondrial encounter structure (ERMES) physically links the membranes of the ER and mitochondria in yeast. Although the ER and mitochondria cooperate to synthesize glycerophospholipids, whether ERMES directly facilitates the lipid exchange between the two organelles remains controversial. Here, we compared the x-ray structures of an ERMES subunit Mdm12 from with that of Mdm12 from and found that both Mdm12 proteins possess a hydrophobic pocket for phospholipid binding.

Conformational effects of N-glycan core fucosylation of immunoglobulin G Fc region on its interaction with Fcγ receptor IIIa.

Antibody-dependent cellular cytotoxicity (ADCC) is promoted through interaction between the Fc region of immunoglobulin G1 (IgG1) and Fcγ receptor IIIa (FcγRIIIa), depending on N-glycosylation of these glycoproteins. In particular, core fucosylation of IgG1-Fc N-glycans negatively affects this interaction and thereby compromises ADCC activity. To address the mechanisms of this effect, we performed replica-exchange molecular dynamics simulations based on crystallographic analysis of a soluble form of FcγRIIIa (sFcγRIIIa) in complex with IgG1-Fc.

[Molecular Simulations of Protein Systems toward Drug Discovery].

  Molecular simulations have been widely used in biomolecular systems such as proteins, DNA, etc. The search for stable conformations of proteins by molecular simulations is important to understand the function and stability of proteins. However, finding the stable state by conformational search is difficult, because the energy landscape of the system is characterized by many local minima separated by high energy barriers.

Exploration of conformational spaces of high-mannose-type oligosaccharides by an NMR-validated simulation.

Exploration of the conformational spaces of flexible biomacromolecules is essential for quantitatively understanding the energetics of their molecular recognition processes. We employed stable isotope- and lanthanide-assisted NMR approaches in conjunction with replica-exchange molecular dynamics (REMD) simulations to obtain atomic descriptions of the conformational dynamics of high-mannose-type oligosaccharides, which harbor intracellular glycoprotein-fate determinants in their triantennary structures. The experimentally validated REMD simulation provided quantitative views of the dynamic conformational ensembles of the complicated, branched oligosaccharides, and indicated significant expansion of the conformational space upon removal of a terminal mannose residue during the functional glycan-processing pathway.

Amino-acid-dependent main-chain torsion-energy terms for protein systems.

Many commonly used force fields for protein systems such as AMBER, CHARMM, GROMACS, OPLS, and ECEPP have amino-acid-independent force-field parameters for main-chain torsion-energy terms. Here, we propose a new type of amino-acid-dependent torsion-energy terms in the force fields. As an example, we applied this approach to AMBER ff03 force field and determined new amino-acid-dependent parameters for ψ (N-C(α)-C-N) and ζ (C(β)-C(α)-C-N) angles for each amino acid by using our optimization method, which is one of the knowledge-based approach.

Protein structure predictions by parallel simulated annealing molecular dynamics using genetic crossover.

We propose a conformational search method to find a global minimum energy structure for protein systems. The simulated annealing is a powerful method for local conformational search. On the other hand, the genetic crossover can search the global conformational space.

Kinetic investigations of the process of encapsulation of small hydrocarbons into a cavitand-porphyrin.

Exchange of guest molecules into capsule shaped host molecules is the most fundamental process in host-guest chemistry. Several examples of quantitative measurements of guest exchange rates have been reported. However, there have been no reports on the activation energies of these processes.