Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Binding of pentagalloyl glucose to two globular proteins occurs via multiple surface sites

The interaction between pentagalloyl glucose (PGG) and two globular proteins, bovine serum albumin (BSA) and ribulose-1,5-bisphosphate carboxylase oxygenase (rubisco), was investigated by isothermal titration calorimetry (ITC). ITC data fit to a binding model consisting of two sets of multiple binding sites, which reveal similarities in the mode of binding of PGG to BSA and rubisco. In both cases, the interaction is characterized by a high number of binding sites, which suggests that binding occurs by a surface adsorption mechanism that leads to coating of the protein surface, which promotes aggregation and precipitation of the PGG-protein complex. This model was confirmed by turbidimetry analysis of the PGG-BSA interaction. Analysis of tryptophan fluorescence quenching during the interaction of PGG with BSA suggests that binding of PGG leads to some conformational changes that are energetically closer to the unfolded state of the BSA structure, because small red shifts in the resulting emission spectra were observed.

c-Jun N-terminal kinase (JNK)-mediated modulation of brain mitochondria function: new target proteins for JNK signalling in mitochondrion-dependent apoptosis

The molecular mechanisms underlying the initiation and control of the release of cytochrome c during mitochondrion-dependent apoptosis are thought to involve the phosphorylation of mitochondrial Bcl-2 and Bcl-x(L). Although the c-Jun N-terminal kinase (JNK) has been proposed to mediate the phosphorylation of Bcl-2/Bcl-x(L) the mechanisms linking the modification of these proteins and the release of cytochrome c remain to be elucidated. This study was aimed at establishing interdependency between JNK signalling and mitochondrial apoptosis. Using an experimental model consisting of isolated, bioenergetically competent rat brain mitochondria, these studies show that (i) JNK catalysed the phosphorylation of Bcl-2 and Bcl-x(L) as well as other mitochondrial proteins, as shown by two-dimensional isoelectric focusing/SDS/PAGE; (ii) JNK-induced cytochrome c release, in a process independent of the permeability transition of the inner mitochondrial membrane (imPT) and insensitive to cyclosporin A; (iii) JNK mediated a partial collapse of the mitochondrial inner-membrane potential (Deltapsim) in an imPT- and cyclosporin A-independent manner; and (iv) JNK was unable to induce imPT/swelling and did not act as a co-inducer, but as an inhibitor of Ca-induced imPT. The results are discussed with regard to the functional link between the Deltapsim and factors influencing the permeability transition of the inner and outer mitochondrial membranes. Taken together, JNK-dependent phosphorylation of mitochondrial proteins including, but not limited to, Bcl-2/Bcl-x(L) may represent a potential of the modulation of mitochondrial function during apoptosis.

Lipid binding interactions of antimicrobial plant seed defence proteins: puroindoline-a and β-purothionin

The indolines and thionins are basic, amphiphilic and cysteine-rich proteins found in cereals; puroindoline-a (Pin-a) and β-purothionin (β-Pth) are members of these families in wheat (Triticum aestivum). Pin-a and β-Pth have been suggested to play a significant role in seed defence against microbial pathogens, making the interaction of these proteins with model bacterial membranes an area of potential interest. We have examined the binding of these proteins to lipid monolayers composed of 1,2-dipalmitoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DPPG) using a combination of neutron reflectometry, Brewster angle microscopy, and infrared spectroscopy. Results showed that both Pin-a and β-Pth interact strongly with condensed phase DPPG monolayers, but the degree of penetration was different. β-Pth was shown to penetrate the lipid acyl chain region of the monolayer and remove lipids from the air/liquid interface during the adsorption process, suggesting this protein may be able to both form membrane spanning ion channels and remove membrane phospholipids in its lytic activity. Conversely, Pin-a was shown to interact mainly with the head-group region of the condensed phase DPPG monolayer and form a 33 Å thick layer below the lipid film. The differences between the interfacial structures formed by these two proteins may be related to the differing composition of the Pin-a and β-Pth hydrophobic regions.

Dynamics of interacting edge defects in copolymer lamellae

It is known that terraces at the air-polymer interface of lamella forming diblock copolymers do not make discontinuous jumps in height. Despite the underlying discretized structure, the height profiles are smoothly varying. The width of a transition region of a terrace edge in isolation is typically several hundreds of nanometres, resulting from a balance between surface tension, chain stretching penalties, and the enthalpy of mixing. What is less well known in these systems is what happens when two transition regions interact with one another. In this study, we investigate the dynamics of the interactions between copolymer lamellar edges. We find that the data can be well described by a model that assumes a repulsion between adjacent edges. While the model is simplistic, and does not include molecular level details, its agreement with the data suggest that some of the the underlying assumptions provide insight into the complex interplay between defects.