We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use 13C NMR spin-lattice relaxation times T1 of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4-hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius-type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force-field (termed as AMBER99SB-ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.

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.

The bidirectional force transmission process of integrin through the cell membrane is still not well understood. Several possible mechanisms have been discussed in literature on the basis of experimental data, and in this study, we investigate these mechanisms by free and steered molecular dynamics simulations. For the first time, constant velocity pulling on the complete integrin molecule inside a dipalmitoyl‐phosphatidylcholine membrane is conducted. From the results, the most likely mechanism for inside‐out and outside‐in signaling is the switchblade model with further separation of the transmembrane helices.

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.