Conformational study of Met-enkephalin based on the ECEPP force fields

Role of single-point mutations and deletions on transition temperatures in ideal proteinogenic heteropolymer chains in the gas phase

A coarse-grained statistical mechanics-based model for ideal heteropolymer proteinogenic chains of non-interacting residues is presented in terms of the size K of the chain and the set of helical propensities [Formula: see text] associated with each residue j along the chain. For this model, we provide an algorithm to compute the degeneracy tensor [Formula: see text] associated with energy level [Formula: see text] where [Formula: see text] is the number of residues with a native contact in a given conformation. From these results, we calculate the equilibrium partition function [Formula: see text] and characteristic temperature [Formula: see text] at which a transition from a low to a high entropy states is observed. The formalism is applied to analyze the effect on characteristic temperatures [Formula: see text] of single-point mutations and deletions of specific amino acids [Formula: see text] along the chain. Two probe systems are considered. First, we address the case of a random heteropolymer of size K and given helical propensities [Formula: see text] on a conformational phase space. Second, we focus our attention to a particular set of neuropentapeptides, [Met-5] and [Leu-5] enkephalins whose thermodynamic stability is a key feature on their coupling to [Formula: see text] and [Formula: see text] receptors and the triggering of biochemical responses.

Transition pathway and its free-energy profile: a protocol for protein folding simulations

We propose a protocol that provides a systematic definition of reaction coordinate and related free-energy profile as the function of temperature for the protein-folding simulation. First, using action-derived molecular dynamics (ADMD), we investigate the dynamic folding pathway model of a protein between a fixed extended conformation and a compact conformation. We choose the pathway model to be the reaction coordinate, and the folding and unfolding processes are characterized by the ADMD step index, in contrast to the common a priori reaction coordinate as used in conventional studies. Second, we calculate free-energy profile as the function of temperature, by employing the replica-exchange molecular dynamics (REMD) method. The current method provides efficient exploration of conformational space and proper characterization of protein folding/unfolding dynamics from/to an arbitrary extended conformation. We demonstrate that combination of the two simulation methods, ADMD and REMD, provides understanding on molecular conformational changes in proteins. The protocol is tested on a small protein, penta-peptide of met-enkephalin. For the neuropeptide met-enkephalin system, folded, extended, and intermediate sates are well-defined through the free-energy profile over the reaction coordinate. Results are consistent with those in the literature.

Exploring Multidimensional Free Energy Landscapes Using Time-Dependent Biases on Collective Variables

A new implementation of the adaptive biasing force (ABF) method is described. This implementation supports a wide range of collective variables and can be applied to the computation of multidimensional energy profiles. It is provided to the community as part of a code that implements several analogous methods, including metadynamics. ABF and metadynamics have not previously been tested side by side on identical systems. Here, numerical tests are carried out on processes including conformational changes in model peptides and translocation of a halide ion across a lipid membrane through a peptide nanotube. On the basis of these examples, we discuss similarities and differences between the ABF and metadynamics schemes. Both approaches provide enhanced sampling and free energy profiles in quantitative agreement with each other in different applications. The method of choice depends on the dimension of the reaction coordinate space, the height of the barriers, and the relaxation times of degrees of freedom in the orthogonal space, which are not explicitly described by the chosen collective variables.

Computational study of the Trp-cage miniprotein based on the ECEPP/3 force field

Using a newly developed Monte Carlo global optimization method called basin paving, we have performed an ab initio computation for the structure of Trp-cage based on the ECEPP/3 force field in vacuo. The lowest energy minimum has been located. Its corresponding configuration is comparable to the native structure of Trp-cage (PDB code 1L2Y) with a backbone root mean square deviation of 2.24 A.