No abstract

The structure of a protein is often not completely accessible by experiments. In silico, replica exchange molecular dynamics (REMD) is the standard sampling method for predicting the secondary and tertiary structures from the amino acid sequence, but it is computationally very expensive. Two recent adaptations from REMD, temperature intervals with global exchange of replicas (TIGER2) and TIGER2A, have been tested here in implicit and explicit solvents. Additionally, explicit, implicit, and hybrid solvent REMD are compared. On the basis of the hybrid REMD (REMDh) method, we present a new hybrid TIGER2h algorithm for faster structural sampling, while retaining good accuracy. The implementations of REMDh, TIGER2, TIGER2A, and TIGER2h are provided for nanoscale molecular dynamics (NAMD). All the methods were tested with two model peptides of known structure, (AAQAA) and HP7, with helix and sheet motifs, respectively. The TIGER2 methods and REMDh were also applied to the unknown structure of the collagen type I telopeptides, which represent bigger proteins with some degree of disorder. We present simulations covering more than 180 μs and analyze the performance and convergence of the distributions of states between the particular methods by dihedral principal component and secondary structure analysis.

In many cases, native states of proteins may be predicted with sufficient accuracy by molecular dynamics simulations (MDSs) with modern force fields. Enhanced sampling methods based on MDS are applied for exploring the phase space of a protein sequence and to overcome barriers on rough conformational energy landscapes. The minimum free energy state is obtained with sampling algorithms providing sufficient convergence and accuracy. A reliable but computationally very expensive method is replica exchange molecular dynamics, with many modifications to this approach presented in the past. Recently, we demonstrated how our temperature intervals with global exchange of replicas hybrid (TIGER2h) solvent sampling algorithm made a good compromise between efficiency and accuracy. There, all states are sampled under full explicit solvent conditions with a freely chosen number of replicas, whereas an implicit solvent is used during the swap decisions. This hybrid method yielded a much better approximation to the agreement with calculations in an explicit solvent than fully implicit solvent simulations. Here, we present an extension of TIGER2h and add a few layers of explicit water molecules around the peptide for the energy calculations, whereas the dynamics in fully explicit water is maintained. We claim that these water layers better reproduce steric effects, the polarization of the solvent, and the resulting reaction field energy than typical implicit solvent models. By investigating the protein−solvent interactions across comprehensive thermodynamic state ensembles, we found a strong conformational dependence of this reaction field energy. All simulations were performed with nanoscale molecular dynamics on two peptides, the α-helical peptide (AAQAA) 3 and the β-hairpin peptide HP7. A production-ready TIGER2hs implementation is supplied, approaching the accuracy of full explicit solvent sampling at a fraction of computational resources.

Synthetic scaffolds containing collagen (Type I) are of increasing interest for bone tissue engineering, especially for highly porous biomaterials in combination with glycosaminoglycans. In experiments the integration of heparin during the fibrillogenesis resulted in different types of collagen fibrils, but models for this aggregation on a molecular scale were only tentative. We conducted molecular dynamic simulations investigating the binding of heparin to collagen and the influence of the telopeptides during collagen aggregation. This aims at explaining experimental findings on a molecular level. Novel structures for N- and C-telopeptides were developed with the TIGER2 replica exchange algorithm and dihedral principle component analysis. We present an extended statistical analysis of the mainly electrostatic interaction between heparin and collagen and identify several binding sites. Finally, we propose a molecular mechanism for the influence of glycosaminoglycans on the morphology of collagen fibrils. Proteins 2017; 85:1119-1130. © 2017 Wiley Periodicals, Inc.

Integrins are transmembrane proteins involved in hemostasis, wound healing, immunity and cancer. In response to intracellular signals and ligand binding, integrins adopt different conformations: the bent (resting) form; the intermediate extended form; and the ligand-occupied active form. An integrin undergoing such conformational dynamics is the heterodimeric platelet receptor αIIbβ3. Although the dramatic rearrangement of the overall structure of αIIbβ3 during the activation process is potentially related to changes in the protein secondary structure, this has not been investigated so far in a membrane environment. Here we examine the Mn 2+ - and drug-induced activation of αIIbβ3 and the impact on the structure of this protein reconstituted into liposomes. By quartz crystal microbalance with dissipation monitoring and activation assays we show that Mn 2+ induces binding of the conformation-specific antibody PAC-1, which only recognizes the extended, active integrin. Circular dichroism spectroscopy reveals, however, that Mn 2+ -treatment does not induce major secondary structural changes of αIIbβ3. Similarly, we found that treatment with clinically relevant drugs (e.g. quinine) led to the activation of αIIbβ3 without significant changes in protein secondary structure. Molecular dynamics simulation studies revealed minor local changes in the beta-sheet probability of several extracellular domains of the integrin. Our experimental setup represents a new approach to study transmembrane proteins, especially integrins, in a membrane environment and opens a new way for testing drug binding to integrins under clinically relevant conditions.