Electrostatic Ratchet for Successive Peptide Synthesis in Nonribosomal Molecular Machine RimK.

A nonribosomal peptide-synthesizing molecular machine, RimK, adds l-glutamic acids to the C-terminus of ribosomal protein S6 (RpsF) in vivo and synthesizes poly-α-glutamates in vitro. However, the mechanism of the successive glutamate addition, which is fueled by ATP, remains unclear. Here, we investigate the successive peptide-synthesizing mechanism of RimK via the molecular dynamics (MD) simulation of glutamate binding.

The Intrinsic Radius as a Key Parameter in the Generalized Born Model to Adjust Protein-Protein Electrostatic Interaction.

The generalized Born (GB) model is an extension of the continuum dielectric theory of Born solvation energy and is a powerful method for accelerating the molecular dynamic (MD) simulations of charged biological molecules in water. While the effective dielectric constant of water that varies as a function of the separation distance between solute molecules is incorporated into the GB model, adjustment of the parameters is indispensable for accurate calculation of the Coulomb (electrostatic) energy. One of the key parameters is the lower limit of the spatial integral of the energy density of the electric field around a charged atom, known as the intrinsic radius ρ.

Torque generation mechanism in F motor of ATP synthase elucidated by free-energy and Coulomb-energy landscapes along the c-ring rotation.

F portion of ATP synthase is a proton-motive rotary motor. The Coulombic attraction between the conserved acidic residues in the c-ring and the arginine in the a-subunit (aR) was early proposed to drive the c-ring rotation relative to the a-subunit, and has been actually observed in our previous molecular dynamics simulation with full atomistic description of F embedded in the membrane. In this study, to quantify the driving force, we conducted the umbrella sampling (US) and obtained the free-energy landscape for the c-ring rotation.

Coulombic Organization in Membrane-Embedded Rotary Motor of ATP Synthase.

The electrochemical potential difference of protons across the membrane is used to synthesize ATP through the proton-motive rotatory motion of the membrane-embedded region of ATP synthase called F. In this study, we illuminate the unsolved proton-motive rotary mechanism of F on the basis of atomistic simulation with full description of protein, lipid, and water molecules, and highlight the underlying Coulombic design. We first show that a water channel is spontaneously formed at the interfacial region between the rotor (-ring) and the stator (-subunit).

Crystal structures of hydroxymethylbilane synthase complexed with a substrate analog: a single substrate-binding site for four consecutive condensation steps.

Hydroxymethylbilane synthase (HMBS), which is involved in the heme biosynthesis pathway, has a dipyrromethane cofactor and combines four porphobilinogen (PBG) molecules to form a linear tetrapyrrole, hydroxymethylbilane. Enzyme kinetic study of human HMBS using a PBG-derivative, 2-iodoporphobilinogen (2-I-PBG), exhibited noncompetitive inhibition with the inhibition constant being 5.4 ± 0.

Biased Brownian Motion of KIF1A and the Role of Tubulin's C-Terminal Tail Studied by Molecular Dynamics Simulation.

KIF1A is a kinesin family protein that moves over a long distance along the microtubule (MT) to transport synaptic vesicle precursors in neurons. A single KIF1A molecule can move toward the plus-end of MT in the monomeric form, exhibiting the characteristics of biased Brownian motion. However, how the bias is generated in the Brownian motion of KIF1A has not yet been firmly established.

Biophysics at Waseda University.

Biophysics in Waseda University was started in 1965 as one of the three key research areas that constitute the Physics Department. In the biophysics group, one theoretical lab and two experimental labs are now working on the cutting-edge themes on biophysics, disseminating the ideas and knowledge of biophysics to undergraduate and graduate students from the viewpoint of physics.

Coupling of Redox and Structural States in Cytochrome P450 Reductase Studied by Molecular Dynamics Simulation.

Cytochrome P450 reductase (CPR) is the key protein that regulates the electron transfer from NADPH to various heme-containing monooxygenases. CPR has two flavin-containing domains: one with flavin adenine dinucleotide (FAD), called FAD domain, and the other with flavin mononucleotide (FMN), called FMN domain. It is considered that the electron transfer occurs via FAD and FMN (NADPH → FAD → FMN → monooxygenase) and is regulated by an interdomain open-close motion.

Hydrophobic Surface Enhances Electrostatic Interaction in Water.

A high dielectric constant is one of the peculiar properties of liquid water, indicating that the electrostatic interaction between charged substances is largely reduced in water. We show by molecular dynamics simulation that the dielectric constant of water is decreased near the hydrophobic surface. We further show that the decrease in the dielectric constant is due to both the decreased water density and the reduced water dipole correlation in the direction perpendicular to the surface.

Structural and dynamical characteristics of tropomyosin epitopes as the major allergens in shrimp.

Ingestion of marine invertebrates often causes food allergy, where the major allergens have been reported to be derived from tropomyosin (TM). Intact or the digestive fragments of food allergens generally show resistance to digestion, which is usually attributable to the structural stability (or rigidity). The difference in the structural and dynamical characteristics between the epitope and the non-epitope regions in TM has not yet been well understood.

Long-range coupling between ATP-binding and lever-arm regions in myosin via dielectric allostery.

A protein molecule is a dielectric substance, so the binding of a ligand is expected to induce dielectric response in the protein molecule, considering that ligands are charged or polar in general. We previously reported that binding of adenosine triphosphate (ATP) to molecular motor myosin actually induces such a dielectric response in myosin due to the net negative charge of ATP. By this dielectric response, referred to as "dielectric allostery," spatially separated two regions in myosin, the ATP-binding region and the actin-binding region, are allosterically coupled.

Electrostatic balance between global repulsion and local attraction in reentrant polymerization of actin.

Actin polymerization depends on the salt concentration, exhibiting a reentrant behavior: the polymerization is promoted by increasing KCl concentration up to 100 mM, and then depressed by further increase above 100 mM. We here investigated the physical mechanism of this reentrant behavior by calculating the polymerization energy, defined by the electrostatic energy change upon binding of an actin subunit to a filament, using an implicit solvent model based on the Poisson-Boltzmann (PB) equation. We found that the polymerization energy as a function of the salt concentration shows a non-monotonic reentrant-like behavior, with the minimum at about 100 mM (1:1 salt).

Over-Destabilization of Protein-Protein Interaction in Generalized Born Model and Utility of Energy Density Integration Cutoff.

The generalize Born (GB) model is frequently used in MD simulations of biomolecular systems in aqueous solution. The GB model is usually based on the so-called Coulomb field approximation (CFA) for the energy density integration. In this study, we report that the GB model with CFA overdestabilizes the long-range electrostatic attraction between oppositely charged molecules (ionic bond forming two-helix system and kinesin-tubulin system) when the energy density integration cutoff, r, which is used to calculate the Born energy, is set to a large value.

Dielectric Allostery of Protein: Response of Myosin to ATP Binding.

Protein uses allostery to execute biological function. The physical mechanism underlying the allostery has long been studied, with the focus on the mechanical response by ligand binding. Here, we highlight the electrostatic response, presenting an idea of "dielectric allostery".

Piezoelectric allostery of protein.

Allostery is indispensable for a protein to work, where a locally applied stimulus is transmitted to a distant part of the molecule. While the allostery due to chemical stimuli such as ligand binding has long been studied, the growing interest in mechanobiology prompts the study of the mechanically stimulated allostery, the physical mechanism of which has not been established. By molecular dynamics simulation of a motor protein myosin, we found that a locally applied mechanical stimulus induces electrostatic potential change at distant regions, just like the piezoelectricity.

Intrinsic disorder accelerates dissociation rather than association.

The intrinsically disordered protein (IDP) has distinct properties both physically and biologically: it often becomes folded when binding to the target and is frequently involved in signal transduction. The physical property seems to be compatible with the biological property where fast association and dissociation between IDP and the target are required. While fast association has been well studied, fueled by the fly-casting mechanism, the dissociation kinetics has received less attention.

Coupling of lever arm swing and biased Brownian motion in actomyosin.

An important unresolved problem associated with actomyosin motors is the role of Brownian motion in the process of force generation. On the basis of structural observations of myosins and actins, the widely held lever-arm hypothesis has been proposed, in which proteins are assumed to show sequential structural changes among observed and hypothesized structures to exert mechanical force. An alternative hypothesis, the Brownian motion hypothesis, has been supported by single-molecule experiments and emphasizes more on the roles of fluctuating protein movement.

Temperature-enhanced association of proteins due to electrostatic interaction: a coarse-grained simulation of actin-myosin binding.

Association of protein molecules constitutes the basis for the interaction network in a cell. Despite its fundamental importance, the thermodynamic aspect of protein-protein binding, particularly the issues relating to the entropy change upon binding, remains elusive. The binding of actin and myosin, which are vital proteins in motility, is a typical example, in which two different binding mechanisms have been argued: the binding affinity increases with increasing temperature and with decreasing salt-concentration, indicating the entropy-driven binding and the enthalpy-driven binding, respectively.

Conformational Ensembles of an Intrinsically Disordered Protein pKID with and without a KIX Domain in Explicit Solvent Investigated by All-Atom Multicanonical Molecular Dynamics.

The phosphorylated kinase-inducible activation domain (pKID) adopts a helix-loop-helix structure upon binding to its partner KIX, although it is unstructured in the unbound state. The N-terminal and C-terminal regions of pKID, which adopt helices in the complex, are called, respectively, αA and αB. We performed all-atom multicanonical molecular dynamics simulations of pKID with and without KIX in explicit solvents to generate conformational ensembles.

Unidirectional Brownian motion observed in an in silico single molecule experiment of an actomyosin motor.

The actomyosin molecular motor, the motor composed of myosin II and actin filament, is responsible for muscle contraction, converting chemical energy into mechanical work. Although recent single molecule and structural studies have shed new light on the energy-converting mechanism, the physical basis of the molecular-level mechanism remains unclear because of the experimental limitations. To provide a clue to resolve the controversy between the lever-arm mechanism and the Brownian ratchet-like mechanism, we here report an in silico single molecule experiment of an actomyosin motor.

Residue network in protein native structure belongs to the universality class of a three-dimensional critical percolation cluster.

Single protein molecules are regarded as contact networks of amino-acid residues. Relationships between the shortest path lengths and the numbers of residues within single molecules in the native structures are examined for various sized proteins. A universal scaling among proteins is obtained, which shows that the residue networks are fractal networks.

Temperature of a Hamiltonian system given as the effective temperature of a nonequilibrium steady-state Langevin thermostat.

In nonequilibrium steady states (NESS) far from equilibrium, it is known that the Einstein relation is violated. Then, the ratio of the diffusion coefficient to the mobility is called an effective temperature, and the physical relevance of this effective temperature has been studied in several works. Although the physical relevance is not yet completely clear, it has been found that the role of an effective temperature in NESS is indeed analogous to that of the temperature in equilibrium systems in a number of respects.

Modeling of a propagation mechanism of infectious prion protein; a hexamer as the minimum infectious unit.

To construct a new model of the propagation mechanism of infectious scrapie-type prion protein (PrP(Sc)), here we conducted a disruption simulation of a PrP(Sc) nonamer using structure-based molecular dynamics simulation method based on a hypothetical PrP(Sc) model structure. The simulation results showed that the nonamer disrupted in cooperative manners into monomers via two significant intermediate states: (1) a nonamer with a partially unfolded surface trimer and (2) a hexamer and three monomers. Dimers and trimers were rarely observed.

Violation of the fluctuation-dissipation theorem in a protein system.

We report the results of molecular dynamics simulations of the protein myosin carried out with an elastic network model. Quenching the system, we observe glassy behavior of a density correlation function and a density response function that are often investigated in structure glasses and spin glasses. In the equilibrium, the fluctuation-response relation, a representative relation of the fluctuation-dissipation theorem, holds that the ratio of the density correlation function to the density response function is equal to the temperature of the environment.