Comparison of the United- and All-Atom Representations of (Halo)alkanes Based on Two Condensed-Phase Force Fields Optimized against the Same Experimental Data Set

The level of accuracy that can be achieved by a force field is influenced by choices made in the interaction-function representation and in the relevant simulation parameters. These choices, referred to here as functional-form variants (FFVs), include for example the model resolution, the charge-derivation procedure, the van der Waals combination rules, the cutoff distance, and the treatment of the long-range interactions. Ideally, assessing the effect of a given FFV on the intrinsic accuracy of the force-field representation requires that only the specific FFV is changed and that this change is performed at an optimal level of parametrization, a requirement that may prove extremely challenging to achieve in practice. Here, we present a first attempt at such a comparison for one specific FFV, namely the choice of a united-atom (UA) versus an all-atom (AA) resolution in a force field for saturated acyclic (halo)alkanes. Two force-field versions (UA vs AA) are optimized in an automated way using the CombiFF approach against 961 experimental values for the pure-liquid densities ρ_liq and vaporization enthalpies ΔH_vap of 591 compounds. For the AA force field, the torsional and third-neighbor Lennard-Jones parameters are also refined based on quantum-mechanical rotational-energy profiles. The comparison between the UA and AA resolutions is also extended to properties that have not been included as parameterization targets, namely the surface-tension coefficient γ, the isothermal compressibility κ_T, the isobaric thermal-expansion coefficient α_P, the isobaric heat capacity c_P, the static relative dielectric permittivity ϵ, the self-diffusion coefficient D, the shear viscosity η, the hydration free energy ΔG_wat, and the free energy of solvation ΔG_che in cyclohexane. For the target properties ρ_liq and ΔH_vap, the UA and AA resolutions reach very similar levels of accuracy after optimization. For the nine other properties, the AA representation leads to more accurate results in terms of η; comparably accurate results in terms of γ, κ_T, α_P, ϵ, D, and ΔG_che; and less accurate results in terms of c_P and ΔG_wat. This work also represents a first step toward the calibration of a GROMOS-compatible force field at the AA resolution.

Application of a Distance-Dependent Sigmoidal Dielectric Constant to the REMC/SAAP3D Simulations of Chignolin, Trp-Cage, and the G10q Mutant

The replica-exchange Monte Carlo method based on the single amino acid potential (SAAP) force field, i.e., REMC/SAAP3D, was recently developed by our group for the molecular simulation of short peptides. In this study, the method has been improved by applying a distance-dependent dielectric (DDD) constant and extended to the peptides containing D-amino acid (AA) residues. For chignolin (10 AAs), a sigmoidal DDD model reasonably allocated the native-like β-hairpin structure with all-atom root mean square deviation (RMSD) = 2.0 Å as a global energy minimum. The optimal DDD condition was subsequently applied for Trp-cage (20 AAs) and its G10q mutant. The native-like α-rich folded structures with main-chain RMSD = 3.7 and 3.8 Å were obtained as global energy minima for Trp-cage and G10q, respectively. The results suggested that the REMC/SAAP3D method with the sigmoidal DDD model is useful for structural prediction for the short peptides comprised of up to 20 AAs. In addition, the relative contributions of SAAP to the total energy (%SAAP) were evaluated by energetic component analysis. The ratios of %SAAP were about 40 and 20% for chignolin and Trp-cage (or G10q), respectively. It was proposed that SAAP is more important for the secondary structure formation than for the assembly to a higher-order folded structure, in which the attractive van der Waals interaction may play a more important role.

Temperature evolution of Trp-cage folding pathways: An analysis by dividing the probability flux field into stream tubes

Owing to its small size and very fast folding rate, the Trp-cage miniprotein has become a benchmark system to study protein folding. Two folding pathways were found to be characteristic of this protein: pathway I, in which the hydrophobic collapse precedes the formation of α-helix, and pathway II, in which the events occur in the reverse order. At the same time, the relative contribution of these pathways at different temperatures as well as the nature of transition from one pathway to the other remain unclear. To gain insight into this issue, we employ a recently proposed hydrodynamic description of protein folding, in which the process of folding is considered as a motion of a "folding fluid" (Chekmarev et al., Phys. Rev. Lett. 100(1), 018107 2008). Using molecular dynamics simulations, we determine the field of probability fluxes of transitions in a space of collective variables and divide it into stream tubes. Each tube contains a definite fraction of the total folding flow and can be associated with a certain pathway. Specifically, three temperatures were considered, T = 285K, T = 315K, and T = 325K. We have found that as the temperature increases, the contribution of pathway I, which is approximately 90% of the total folding flow at T = 285K, decreases to approximately 10% at T = 325K, i.e., pathway II becomes dominant. At T = 315K, both pathways contribute approximately equally. All these temperatures are found below the calculated melting point, which suggests that the Trp-cage folding mechanism is determined by kinetic factors rather than thermodynamics.

A hydrodynamic view of the first-passage folding of Trp-cage miniprotein

We study folding of Trp-cage miniprotein in the conditions when the native state of the protein is stable and unfolding events are improbable, which corresponds to physiological conditions. Using molecular dynamics simulations with an implicit solvent model, an ensemble of folding trajectories from unfolded (practically extended) states of the protein to the native state was generated. To get insight into the folding kinetics, the free energy surface and kinetic network projected on this surface were constructed. This, "conventional" analysis of the folding reaction was followed by a recently proposed hydrodynamic description of protein folding (Chekmarev et al. in Phys Rev Lett 100(1):018107, 2008), in which the process of the first-passage folding is viewed as a stationary flow of a folding "fluid" from the unfolded to native state. This approach is conceptually different from the previously used approaches and thus allows an alternative view of the folding dynamics and kinetics of Trp-cage, the conclusions about which are very diverse. In agreement with most previous studies, we observed two characteristic folding pathways: in one pathway (I), the collapse of the hydrophobic core precedes the formation of the [Formula: see text]-helix, and in the other pathway (II), these events occur in the reverse order. We found that although pathway II is complicated by a repeated partial protein unfolding, it contributes to the total folding flow as little as ≈10%, so that the folding kinetics remain essentially single-exponential.

Tuning the Attempt Frequency of Protein Folding Dynamics via Transition-State Rigidification: Application to Trp-Cage

The attempt frequency or prefactor (k0) of the transition-state rate equation of protein folding kinetics has been estimated to be on the order of 10(6) s(-1), which is many orders of magnitude smaller than that of chemical reactions. Herein we use the mini-protein Trp-cage to show that it is possible to significantly increase the value of k0 for a protein folding reaction by rigidifying the transition state. This is achieved by reducing the conformational flexibility of a key structural element (i.e., an α-helix) formed in the transition state via photoisomerization of an azobenzene cross-linker. We find that this strategy not only decreases the folding time of the Trp-cage peptide by more than an order of magnitude (to ∼100 ns at 25°C) but also exposes parallel folding pathways, allowing us to provide, to the best of our knowledge, the first quantitative assessment of the curvature of the transition-state free-energy surface of a protein.