Thermal conductivity and conductance of protein in aqueous solution: Effects of geometrical shape.

Considering the importance of elucidating the heat transfer in living cells, we evaluated the thermal conductivity κ and conductance G of hydrated protein through all-atom non-equilibrium molecular dynamics simulation. Extending the computational scheme developed in earlier studies for spherical protein to cylindrical one under the periodic boundary condition, we enabled the theoretical analysis of anisotropic thermal conduction and also discussed the effects of protein size correction on the calculated results. While the present results for myoglobin and green fluorescent protein (GFP) by the spherical model were in fair agreement with previous computational and experimental results, we found that the evaluations for κ and G by the cylindrical model, in particular, those for the longitudinal direction of GFP, were enhanced substantially, but still keeping a consistency with experimental data. We also studied the influence by salt addition of physiological concentration, finding insignificant alteration of thermal conduction of protein in the present case.

Remarked suppression of Aß42 protomer-protomer dissociation reaction elucidated by molecular dynamics simulation.

Multimeric protein complexes are molecular apparatuses to regulate biological systems and often determine their fate. Among proteins forming such molecular assemblies, amyloid proteins have drawn attention over a half-century since amyloid fibril formation of these proteins is supposed to be a common pathogenic cause for neurodegenerative diseases. This process is triggered by the accumulation of fibril-like aggregates, while the microscopic mechanisms are mostly elusive due to technical limitation of experimental methodologies in individually observing each of diverse aggregate species in the aqueous solution. We then addressed this problem by employing atomistic molecular dynamics simulations for the paradigmatic amyloid protein, amyloid-ß (Aß42 ). Seven different dimeric forms of oligomeric Aß42 fibril-like aggregate in aqueous solution, ranging from tetramer to decamer, were considered. We found additive effects of the size of these fibril-like aggregates on their thermodynamic stability and have clarified kinetic suppression of protomer-protomer dissociation reactions at and beyond the point of pentamer dimer formation. This observation was obtained from the specific combination of the Aß42 protomer structure and the physicochemical condition that we here examined, while it is worthwhile to recall that several amyloid fibrils take dimeric forms of their protomers. We could thus conclude that the stable formation of fibril-like protomer dimer should be involved in a turning point where rapid growth of amyloid fibrils is triggered.

Computational prediction of heteromeric protein complex disassembly order using hybrid Monte Carlo/molecular dynamics simulation.

The physicochemical entities comprising the biological phenomena in the cell form a network of biochemical reactions and the activity of such a network is regulated by multimeric protein complexes. Mass spectroscopy (MS) experiments and multimeric protein docking simulations based on structural bioinformatics techniques have revealed the molecular-level stoichiometry and static configuration of subcomplexes in their bound forms, thus revealing the subcomplex population and formation orders. Meanwhile, these methodologies are not designed to straightforwardly examine the temporal dynamics of multimeric protein assembly and disassembly, essential physicochemical properties to understand the functional expression mechanisms of proteins in the biological environment. To address this problem, we have developed an atomistic simulation in the framework of the hybrid Monte Carlo/molecular dynamics (hMC/MD) method and succeeded in observing the disassembly of a homomeric pentamer of the serum amyloid P component protein in an experimentally consistent order. In this study, we improved the hMC/MD method to examine the disassembly processes of the tryptophan synthase tetramer, a paradigmatic heteromeric protein complex in MS studies. We employed the likelihood-based selection scheme to determine a dissociation-prone subunit pair at every hMC/MD simulation cycle and achieved highly reliable predictions of the disassembly orders without a priori knowledge of the MS experiments and structural bioinformatics simulations. The success rate for the experimentally-observed disassembly order is over 0.9. We similarly succeeded in reliable predictions for three other tetrameric protein complexes. These achievements indicate the potential applicability of our hMC/MD approach as a general-purpose methodology to obtain microscopic and physicochemical insights into multimeric protein complex formation.

Elucidating microscopic events driven by GTP hydrolysis reaction in the Ras-GAP system with semi-reactive molecular dynamics simulations: the alternative role of a phosphate binding loop for mechanical energy storage.

ATPase and GTPase have been widely found as chemical energy-mechanical work transducers, whereas the physicochemical mechanisms are not satisfactorily understood. We addressed the problem by examining John Ross' conjecture that repulsive Coulomb interaction between ADP/GDP and inorganic phosphate (Pi) does the mechanical work upon the system. We effectively simulated the consequence of a GTP hydrolysis reaction in a complex system of Rat sarcoma (Ras) and GTPase activation protein (GAP) in the framework of classical molecular dynamics by switching force field parameters between the reactant and product systems. We then observed a ca. 5 kcal mol-1 increase of potential energy about the phosphate-binding loop (P-loop) in the Ras protein, indicating that the mechanical work generated via the GTP hydrolysis is converted into the local interaction energy and stored in the P-loop. Interestingly, this local energy storage in the P-loop depends on neither impulsive nor consecutive collisions of GDP and Pi with the P-loop. Instead, GTP-GDP conversion itself does work on the Ras system, elevating the potential energy. These observations encourage us to challenge a conjecture previously given by Ross. We assert that triphosphate nucleotide hydrolyses do mechanical work by producing emergent steric interaction accompanied by relaxation, namely, a shift of the biomolecular system to the non-equilibrium state on the reshaped potential energy landscape. Recalling the universality of the P-loop motif among GTPases and ATPases, the observations that we obtained through this study would progress the physicochemical understanding of the operating principles of GTP/ATP hydrolysis-driven biological nano-machines.

Combined mechanism of conformational selection and induced fit in U1A-RNA molecular recognition.

In this study, we demonstrate that U1A-RNA molecular recognition is mediated by a combined mechanism of conformational selection and induced fit. The binding of U1A to RNA has been discussed in the context of induced fit that involves the reorientation of the α-helix in the C-terminal region (Helix-C) of U1A to permit RNA access only when U1A correctly recognizes RNA. However, according to our molecular dynamics simulations, even in the absence of RNA, Helix-C spontaneously reoriented to permit RNA access. Nonetheless, such a conformational change was still incomplete. Helix-C was often partially or even fully unfolded and in an infrequent RNA-accessible conformation, which can be detected using state-of-the-art nuclear magnetic resonance methodology. These results suggest that the formation of an energetically stabilized complex is promoted by specific interactions between U1A and RNA. In conclusion, in the recognition of RNA by U1A protein, we propose a combined mechanism that requires the reorientation of Helix-C and the subsequent contact with RNA through conformational selection, although the stabilization of the U1A-RNA complex is caused by induced fit. We further propose a modification to the conventional assumption regarding the mechanism of U1A-RNA molecular recognition.

Simulation toolkits at the molecular scale for trans-scale thermal signaling.

Thermogenesis is a physiological activity of releasing heat that originates from intracellular biochemical reactions. Recent experimental studies discovered that externally applied heat changes intracellular signaling locally, resulting in global changes in cell morphology and signaling. Therefore, we hypothesize an inevitable contribution of thermogenesis in modulating biological system functions throughout the spatial scales from molecules to individual organisms. One key issue examining the hypothesis, namely, the "trans-scale thermal signaling," resides at the molecular scale on the amount of heat released via individual reactions and by which mechanism the heat is employed for cellular function operations. This review introduces atomistic simulation tool kits for studying the mechanisms of thermal signaling processes at the molecular scale that even state-of-the-art experimental methodologies of today are hardly accessible. We consider biological processes and biomolecules as potential heat sources in cells, such as ATP/GTP hydrolysis and multiple biopolymer complex formation and disassembly. Microscopic heat release could be related to mesoscopic processes via thermal conductivity and thermal conductance. Additionally, theoretical simulations to estimate these thermal properties in biological membranes and proteins are introduced. Finally, we envisage the future direction of this research field.

Inhibition of Melittin Activity Using a Small Molecule with an Indole Ring.

We investigated d-amino acids as potential inhibitors targeting l-peptide toxins. Among the l- and d-amino acids tested, we found that d-tryptophan (d-Trp) acted as an inhibitor of melittin-induced hemolysis. We then evaluated various Trp derivatives and found that 5-chlorotryptamine (5CT) had the largest inhibitory effect on melittin. The indole ring, amino group, and steric hindrance of an inhibitor played important roles in the inhibition of melittin activity. Despite the small size and simple molecular structure of 5CT, its IC50 was approximately 13 µg/mL. Fluorescence quenching, circular dichroism measurements, and size-exclusion chromatography revealed that 5CT interacted with Trp19 in melittin and affected the formation of the melittin tetramer involved in hemolysis. Molecular dynamics simulation of melittin also indicated that the interaction of 5CT with Trp19 in melittin affected the formation of the tetramer.

Reaction Pathway Sampling and Free-Energy Analyses for Multimeric Protein Complex Disassembly by Employing Hybrid Configuration Bias Monte Carlo/Molecular Dynamics Simulation.

Physicochemical characterization of multimeric biomacromolecule assembly and disassembly processes is a milestone to understand the mechanisms for biological phenomena at the molecular level. Mass spectroscopy (MS) and structural bioinformatics (SB) approaches have become feasible to identify subcomplexes involved in assembly and disassembly, while they cannot provide atomic information sufficient for free-energy calculation to characterize transition mechanism between two different sets of subcomplexes. To combine observations derived from MS and SB approaches with conventional free-energy calculation protocols, we here designed a new reaction pathway sampling method by employing hybrid configuration bias Monte Carlo/molecular dynamics (hcbMC/MD) scheme and applied it to simulate the disassembly process of serum amyloid P component (SAP) pentamer. The results we obtained are consistent with those of the earlier MS and SB studies with respect to SAP subcomplex species and the initial stage of SAP disassembly processes. Furthermore, we observed a novel dissociation event, ring-opening reaction of SAP pentamer. Employing free-energy calculation combined with the hcbMC/MD reaction pathway trajectories, we moreover obtained experimentally testable observations on (1) reaction time of the ring-opening reaction and (2) importance of Asp42 and Lys117 for stable formation of SAP oligomer.

Chloride Ions Stabilize Human Adult Hemoglobin in the T-State, Competing with Allosteric Interaction of Oxygen Molecules.

In the context of a molecular-level understanding of the allostery mechanisms, human adult hemoglobin (HbA) has been extensively studied for over half a century. Chloride ions (Cl-) have been known as one of HbA allosteric effectors, which stabilizes the T-state preferable to release oxygen molecules. The functional mechanisms were individually proposed by Ueno and Perutz several decades ago. Ueno considered that the site-specific Cl- binding is essential, while Perutz proposed the non-site-specific interaction between HbA and Cl-. Each speculation explains the mechanism plausibly since each was tightly associated with its reasonable experimental observation. However, both mechanisms themselves still seem to make their speculations controversial. In the present study, we have theoretically reconsidered these apart from their approaches. Our atomistic molecular dynamics simulations then showed that the increase of Cl- concentration suppresses the conformational conversion from the T-state. Interestingly, chloride ions loosely interact with the amino acid residues inside the HbA central cavity, suggesting that both Perutz's and Ueno's speculations are involved in understanding the microscopic roles of Cl-. In conclusion, we theoretically certified that the effect of Cl- competes against that of solvated O2, i.e., the destabilization of T-state through the non-site-specific interaction, implying the concerted regulation of HbA under physiological conditions.

ATP Converts Aß42 Oligomer into Off-Pathway Species by Making Contact with Its Backbone Atoms Using Hydrophobic Adenosine.

Adenosine triphosphate (ATP) is newly expected to be involved in the clearance of amyloid ß 1-42 (Aß42) fibril and its precursors, Aß42 oligomer. Meanwhile, the microscopic mechanism of the role in dissolving the protein aggregate still remains elusive. Aiming to elucidate the mechanism, we examined effects of ATP on the conformational change and thermodynamic stability of the protomer dimer of Aß42 pentamer and tetramer, Aß42(9), by employing all-atom molecular dynamics simulations. We observed interprotomer twisting and intraprotomer peeling of Aß42(9). These conformational changes remarkably accelerate dissociation of the protomer dimer. However, the presence of ATP itself has no positive effect on dissociation processes of the protomer dimer and a monomer from the dimer, indicating its irrelevance to decomposition of the Aß42 oligomer. Rather, it could be supposed that ATP prevents additional binding and rebinding of Aß42 monomers to the Aß42 oligomer and it then converts Aß42 oligomer into an off-pathway species which is excluded from Aß42 fibril growth processes. Interestingly, hydrophobic adenosine in ATP makes contact with Aß42(9) on its backbone atoms, with respect to both Aß42 monomers on the edge of Aß42(9) and dissociated Aß42 monomers in Aß42(9). These roles of ATP would be applied without regard to the structural polymorphism of the Aß42 fibril.

Visualization analysis of inter-fragment interaction energies of CRP-cAMP-DNA complex based on the fragment molecular orbital method.

A visualization method for inter-fragment interaction energies (IFIEs) of biopolymers is presented on the basis of the fragment molecular orbital (FMO) method. The IFIEs appropriately illustrate the information about the interaction energies between the fragments consisting of amino acids, nucleotides and other molecules. The IFIEs are usually analyzed in a matrix form called an IFIE matrix. Analyzing the IFIE matrix, we detect important fragments for the function of biomolecular systems and quantify the strength of interaction energies based on quantum chemistry, including the effects of charge transfer, electronic polarization and dispersion force. In this study, by analyzing a protein-DNA complex, we report a visual representation of the IFIE matrix, a so-called IFIE map. We comprehensively examine what information the IFIE map contains concerning structures and stabilities of the protein-DNA complex.

Na+ Binding Is Ineffective in Forming a Primary Substrate Pocket of Thrombin.

Thrombin is a serine protease involved in the blood coagulation reaction, and it shows maximum enzymatic activity in the presence of Na+. It has been supposed that Na+ binding promotes conversion from the inactive form, with a collapsed primary substrate pocket (S1 pocket), to the active form, with a properly formed S1 pocket. However, the evidence supporting this activation mechanism was derived from the X-ray crystallographic structures solved under nonphysiological conditions and using thrombin mutants; thus, it still remains elusive whether the activation mechanism is actually attributed to Na+ binding. To address the problem, we employed all-atom molecular dynamics simulations for both active and inactive forms of thrombin in the presence and absence of Na+ binding and examined the effect of Na+ binding on S1-pocket formation. In contrast to the conventional supposition, we revealed that Na+ binding does not prevent S1-pocket collapse virtually, but rather, the bound Na+ can move to the S1 pocket, thus blocking substrate access directly. Additionally, it was clarified that Na+ binding does not promote S1-pocket formation. According to these insights, we concluded that Na+ binding is irrelevant to the interconversion between the inactive and active forms of thrombin.

Bound Na(+) is a Negative Effecter for Thrombin-Substrate Stereospecific Complex Formation.

Thrombin has been studied as a paradigmatic protein of Na(+)-activated allosteric enzymes. Earlier structural studies suggest that Na(+)-binding promotes the thrombin-substrate association reaction. However, it is still elusive because (1) the structural change, driven by Na(+)-binding, is as small as the thermal fluctuation, and (2) the bound Na(+) is close to Asp189 in the primary substrate binding pocket (S1-pocket), possibly preventing substrate access via repulsive interaction. It still remains a matter of debate whether Na(+)-binding actually promotes the reaction. To solve this problem, we examined the effect of Na(+) on the reaction by employing molecular dynamics (MD) simulations. By executing independent 210 MD simulations of apo and holo systems, we obtained 80 and 26 trajectories undergoing substrate access to S1-pocket, respectively. Interestingly, Na(+)-binding results in a 3-fold reduction of the substrate access. Furthermore, we examined works for the substrate access and release, and found that Na(+)-binding is disadvantageous for the presence of the substrate in the S1-pocket. These observations provide the insight that the bound Na(+) is essentially a negative effecter in thrombin-substrate stereospecific complex formation. The insight rationalizes an enigmatic feature of thrombin, relatively low Na(+)-binding affinity. This is essential to reduce the disadvantage of Na(+)-binding in the substrate-binding.

Toward understanding allosteric activation of thrombin: a conjecture for important roles of unbound Na(+) molecules around thrombin.

We shed light on important roles of unbound Na(+) molecules in enzymatic activation of thrombin. Molecular mechanism of Na(+)-activation of thrombin has been discussed in the context of allostery. However, the recent challenge to redesign K(+)-activated thrombin revealed that the allosteric interaction is insufficient to explain the mechanism. Under these circumstances, we have examined the roles of unbound Na(+) molecule in maximization of thrombin-substrate association reaction rate. We performed all-atomic molecular dynamics (MD) simulations of thrombin in the presence of three different cations; Li(+), Na(+), and Cs(+). Although these cations are commonly observed in the vicinity of the S1-pocket of thrombin, smaller cations are distributed more densely and extensively than larger ones. This suggests the two observation rules: (i) thrombin surrounded by Na(+) is at an advantage in the initial step of association reaction, namely, the formation of an encounter complex ensemble, and (ii) the presence of Na(+) molecules does not necessarily have an advantage in the final step of association reaction, namely, the formation of the stereospecific complex. In conclusion, we propose a conjecture that unbound Na(+) molecules also affect the maximization of rate constant of thrombin-substrate association reaction through optimally forming an encounter complex ensemble.

Dewetting of S1-Pocket via Water Channel upon Thrombin-Substrate Association Reaction.

Upon protein-substrate association reaction, dewetting of the substrate-binding pocket is one of the rate-limiting processes. However, understanding the microscopic mechanism still remains challenging because of practical limitations of experimental methodologies. We have addressed the problem here by using molecular dynamics (MD) simulation of the thrombin-substrate association reaction. During the MD simulation, ArgP1 in a substrate accessed thrombin's substrate-binding pocket and formed specific hydrogen bonds (H-bonds) with Asp189 in thrombin, while the catalytic serine of thrombin was still away from the substrate's active site. It is assumed that the thrombin-substrate association reaction is regulated by a stepwise mechanism. Furthermore, in the earlier stage of ArgP1 access to the pocket, we observed that ArgP1 was spatially separated from Asp189 by two water molecules in the pocket. These water molecules transferred from the pocket, followed by the specific H-bond formation between thrombin and the substrate. Interestingly, they were not evacuated directly from the pocket to the bulk solvent, but moved to the water channel of thrombin. This observation indicates that the channel plays functional roles in dewetting upon the association reaction.

Non-site-specific allosteric effect of oxygen on human hemoglobin under high oxygen partial pressure.

Protein allostery is essential for vital activities. Allosteric regulation of human hemoglobin (HbA) with two quaternary states T and R has been a paradigm of allosteric structural regulation of proteins. It is widely accepted that oxygen molecules (O2) act as a "site-specific" homotropic effector, or the successive O2 binding to the heme brings about the quaternary regulation. However, here we show that the site-specific allosteric effect is not necessarily only a unique mechanism of O2 allostery. Our simulation results revealed that the solution environment of high O2 partial pressure enhances the quaternary change from T to R without binding to the heme, suggesting an additional "non-site-specific" allosteric effect of O2. The latter effect should play a complementary role in the quaternary change by affecting the intersubunit contacts. This analysis must become a milestone in comprehensive understanding of the allosteric regulation of HbA from the molecular point of view.

Oxygen entry through multiple pathways in T-state human hemoglobin.

The heme oxygen (O2) binding site of human hemoglobin (HbA) is buried in the interior of the protein, and there is a debate over the O2 entry pathways from solvent to the binding site. As a first step to understand HbA O2 binding process at the atomic level, we detected all significant multiple O2 entry pathways from solvent to the binding site in the α and ß subunits of the T-state tetramer HbA by utilizing ensemble molecular dynamics (MD) simulation. By executing 128 independent 8 ns MD trajectories in O2-rich aqueous solvent, we simulated the O2 entry processes and obtained 141 and 425 O2 entry events in the α and ß subunits of HbA, respectively. We developed the intrinsic pathway identification by clustering method to achieve a persuasive visualization of the multiple entry pathways including both the shapes and relative importance of each pathway. The rate constants of O2 entry estimated from the MD simulations correspond to the experimentally observed values, suggesting that O2 ligands enter the binding site through multiple pathways. The obtained multiple pathway map can be utilized for future detailed analysis of HbA O2 binding process.

Spontaneous adjustment mechanism in an RNA-binding protein: cooperation between energetic stabilization and target search enhancement.

We propose a novel concept associated with the relationship between structure and function in biomolecular systems. We performed a 75 nanoseconds molecular dynamics (MD) simulation for an RNA-binding protein, neuro-oncological ventral antigen (NOVA), and examined its physico-chemical properties. NOVA dissociated from the NOVA-RNA complex showed a large conformational change: formation of intra-molecular hydrogen bonds between the C-terminal region and the loop structure located at the middle of amino acid sequence. The free energy analysis suggests that the deformed structure is more stabilized in macromolecular crowding environment where the dielectric constant is smaller than 5. The solvent accessible surface area (SASA) analysis indicates that NOVA enhances the efficiency of association with RNA by changing the relative SASA for the target sequence in RNA molecules. Based on the obtained results, we propose a novel concept of spontaneous adjustment mechanism to explain the structural and energetic changes observed for NOVA in the free state.

Identification and expression analysis of rainbow trout pumilio-1 and pumilio-2.

Pumilio is a sequence-specific RNA-binding protein that regulates translation from the relevant mRNA. The PUF-domain, the RNA-binding motif of Pumilio, is highly conserved across species. In the present study, we have identified two pumilio genes (pumilio-1 and pumilio-2) in rainbow trout and analyzed their expression patterns in its tissues. Pumilio-1 mRNA and pumilio-2A mRNA code for typical full length Pumilio proteins that contain a PUF-domain, whereas pumilio-2B mRNA is a splice variant of pumilio-2 and encodes a protein that lacks the PUF-domain. We have also identified a novel 72-bp exon that has not been reported in other animal species but is conserved in fish species. The insertion of this novel exon leads to the expression of an isoform of the Pumilio-2 protein with a slightly altered conformation of the PUF-domain. Pumilio-1 mRNA and pumilio-2A mRNA (irrespective of the presence of the 72-bp exon) are expressed in both the brain and ovaries at high levels, whereas pumilio-2B mRNA is expressed at low levels in all the rainbow trout tissues examined. Western blot analysis also indicates that the full length Pumilio proteins are expressed predominantly in the brain and ovaries. These data suggest that the Pumilio proteins have physiological roles and are involved in regulatory mechanisms in rainbow trout.