Molecular dynamics simulations of N-terminal peptides from a nucleotide binding protein

An investigation into non-covalent functionalization of a single-walled carbon nanotube and a graphene sheet with protein G:A combined experimental and molecular dynamics study

Investigation of non-covalent interaction of hydrophobic surfaces with the protein G (PrG) is necessary due to their frequent utilization in immunosensors and ELISA. It has been confirmed that surfaces, including carbonous-nanostructures (CNS) could orient proteins for a better activation. Herein, PrG interaction with single-walled carbon nanotube (SWCNT) and graphene (Gra) nanostructures was studied by employing experimental and MD simulation techniques. It is confirmed that the PrG could adequately interact with both SWCNT and Gra and therefore fine dispersion for them was achieved in the media. Results indicated that even though SWCNT was loaded with more content of PrG in comparison with the Gra, the adsorption of the PrG on Gra did not induce significant changes in the IgG tendency. Several orientations of the PrG were adopted in the presence of SWCNT or Gra; however, SWCNT could block the PrG-FcR. Moreover, it was confirmed that SWCNT reduced the α-helical structure content in the PrG. Reduction of α-helical structure of the PrG and improper orientation of the PrG-SWCNT could remarkably decrease the PrG tendency to the Fc of the IgG. Importantly, the Gra could appropriately orient the PrG by both exposing the PrG-FcR and also by blocking the fragment of the PrG that had tendency to interact with Fab in IgG.

Molecular simulation-based structural prediction of protein complexes in mass spectrometry: the human insulin dimer

Protein electrospray ionization (ESI) mass spectrometry (MS)-based techniques are widely used to provide insight into structural proteomics under the assumption that non-covalent protein complexes being transferred into the gas phase preserve basically the same intermolecular interactions as in solution. Here we investigate the applicability of this assumption by extending our previous structural prediction protocol for single proteins in ESI-MS to protein complexes. We apply our protocol to the human insulin dimer (hIns2) as a test case. Our calculations reproduce the main charge and the collision cross section (CCS) measured in ESI-MS experiments. Molecular dynamics simulations for 0.075 ms show that the complex maximizes intermolecular non-bonded interactions relative to the structure in water, without affecting the cross section. The overall gas-phase structure of hIns2 does exhibit differences with the one in aqueous solution, not inferable from a comparison with calculated CCS. Hence, care should be exerted when interpreting ESI-MS proteomics data based solely on NMR and/or X-ray structural information.

Sodium hexadecyl sulfate as an interfacial substance adjusting the adsorption of a protein on carbon nanotubes

Carbon nanotubes (CNTs) were functionalized with sodium hexadecyl sulfate (SHS). The lysozyme adsorbed on the SHS-CNTs exhibited a higher activity than that immobilized on the nonfunctionalized CNTs. To explain the experimental results and explore the mechanism of lysozyme adsorption, large-scale molecular dynamics simulations have been performed for a four-component system, including lysozyme, SHS, CNTs in explicit water. It has been found that the assembled SHS molecules form a soft layer on the surface of CNTs. The interactions between lysozyme and SHS induce the rearrangement of SHS molecules, forming a saddle-like structure on the CNT surface. The saddle-like structure fits the shape of the lysozyme, and the active-site cleft of the lysozyme is exposed to the water phase. Whereas, for the lysozyme adsorbed on the nonfunctionalized CNT, due to the hydrophobic interactions, the active-site cleft of the enzyme tends to face the wall of the CNT. The results of this work demonstrate that the SHS molecules as the interfacial substance have a function of adjusting the lysozyme with an appropriate orientation, which is favorable for the lysozyme having a higher activity.

The βαβαβ elementary supersecondary structure of the Rossmann fold from porcine lactate dehydrogenase exhibits characteristics of a molten globule

Protein classifications show that the Rossmann fold, which consists of two betaalphabetaalphabeta motifs (BABAB) related by a rough twofold axis, is the most populated alphabeta fold, and that the betaalphabeta submotif (BAB) is a widespread elementary structural arrangement. Herein, we report MD simulations, circular dichroism and NMR analyses on BAB and BABAB from porcine lactate dehydrogenase to evaluate their intrinsic stability. Our results demonstrate that BAB is not stable in solution and is not a folding nucleus. We also find that BABAB, despite its appearance of a functional and structural unit, is not an independent and thermodynamically stable folding unit. Rather, we show that BABAB retains most native secondary structure but very little tertiary structure, thus displaying characteristics of a molten globule.

GROMACS: Fast, flexible, and free

This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.

Modeling of enzyme-substrate complexes for the metalloproteases MMP-3, ADAM-9 and ADAM-10

The matrix metalloproteases (MMPs) and the ADAMs (A Disintegrin And Metalloprotease domain) are proteolytic enzyme families containing a catalytic zinc ion, that are implicated in a variety of normal and pathological processes involving tissue remodeling and cancer. Synthetic MMP inhibitors have been designed for applications in pathological situations. However, a greater understanding of substrate binding and the catalytic mechanism is required so that more effective and selective inhibitors may be developed for both experimental and clinical purposes. By modeling a natural substrate spanning P4-P4' in complex with the catalytic domains, we aim to compare substrate-specificities between Stromelysin-1 (MMP-3), ADAM-9 and ADAM-10, with the aid of molecular dynamics simulations. Our results show that the substrate retains a favourable antiparallel beta-sheet conformation on the P-side in addition to the well-known orientation of the P'-region of the scissile bond, and that the primary substrate selectivity is dominated by the sidechains in the S1' pocket and the S2/S3 region. ADAM-9 has a hydrophobic residue as the central determinant in the S1' pocket, while ADAM-10 has an amphiphilic residue, which suggests a different primary specificity. The S2/S3 pocket is largely hydrophobic in all three enzymes. Inspired by our molecular dynamics calculations and supported by a large body of literature, we propose a novel, hypothetical, catalytic mechanism where the Zn-ion polarizes the oxygens from the catalytic glutamate to form a nucleophile, leading to a tetrahedral oxyanion anhydride transition state.

Computer simulations aimed at structure prediction of supersecondary motifs in proteins

It is well established that protein structures are more conserved than protein sequences. One-third of all known protein structures can be classified into ten protein folds, which themselves are composed mainly of alpha-helical hairpin, beta hairpin, and betaalphabeta supersecondary structural elements. In this study, we explore the ability of a recent Monte Carlo-based procedure to generate the 3D structures of eight polypeptides that correspond to units of supersecondary structure and three-stranded antiparallel beta sheet. Starting from extended or misfolded compact conformations, all Monte Carlo simulations show significant success in predicting the native topology using a simplified chain representation and an energy model optimized on other structures. Preliminary results on model peptides from nucleotide binding proteins suggest that this simple protein folding model can help clarify the relation between sequence and topology.

Molecular dynamics simulation of the antimicrobial salivary peptide histatin-5 in water and in trifluoroethanol: a microscopic description of the water destructuring effect

The results of 520 ps molecular dynamics simulation of histatin-5, a small peptide present in human saliva and possessing antimicrobial activity, dissolved in water and in 2,2,2-trifluoroethanol, are reported. The simulations indicate that histatin-5 is destabilized in water and begins to unfold after 250 ps, while in organic solvent it maintains a regular secondary structure throughout the trajectory. Analysis of the peptide-solvent hydrogen bonds indicates that 2,2,2-trifluoroethanol is a poorer proton acceptor than water. The fluorine atom of the alcohol is almost never engaged in a hydrogen bond and the organic solvent interacts mainly with the peptide through its hydroxyl group. For some residues analysis of the solvent residence time indicated longer values for 2,2,2-trifluoroethanol than for water. The most striking difference is related to the number of times the solvent enters and leaves the first coordination shell of the peptide. This value was more than one order of magnitude higher for water than for the alcohol, suggesting that this may be the main cause of alpha-helix destabilization perpetrated by water.

Bending of the calmodulin central helix: a theoretical study

The crystal structure of calcium-calmodulin (CaM) reveals a protein with a typical dumbbell structure. Various spectroscopic studies have suggested that the central linker region of CaM, which is alpha-helical in the crystal structure, is flexible in solution. In particular, NMR studies have indicated the presence of a flexible backbone between residues Lys 77 and Asp 80. This flexibility is related directly to the function of the protein because it enables the N- and C-terminal domains of the protein to move toward each other and bind to the CaM-binding domain of a target protein. We have investigated the flexibility of the CaM central helix by a variety of computational techniques: molecular dynamics (MD) simulations, normal mode analysis (NMA), and essential dynamics (ED) analysis. Our MD results reproduce the experimentally determined location of the bend in a simulation of only the CaM central helix, indicating that the bending point is an intrinsic property of the alpha-helix, for which the remainder of the protein is not important. Interestingly, the modes found by the ED analysis of the MD trajectory are very similar to the lowest frequency modes from the NM analysis and to modes found by an ED analysis of different structures in a set of NMR structures. Electrostatic interactions involving residues Arg 74 and Asp 80 seem to be important for these bending motions and unfolding, which is in line with pH-dependent NMR and CD studies.