Lowest Free Energy Conformation Research Articles

SWOTein: a structure-based approach to predict stability Strengths and Weaknesses of prOTEINs.

Although structured proteins adopt their lowest free energy conformation in physiological conditions, the individual residues are generally not in their lowest free energy conformation. Residues that are stability weaknesses are often involved in functional regions, whereas stability strengths ensure local structural stability. The detection of strengths and weaknesses provides key information to guide protein engineering experiments aiming to modulate folding and various functional processes. We developed the SWOTein predictor which identifies strong and weak residues in proteins on the basis of three types of statistical energy functions describing local interactions along the chain, hydrophobic forces and tertiary interactions. The large-scale analysis of the different types of strengths and weaknesses demonstrated their complementarity and the enhancement of the information they provide. Moreover, a good average correlation was observed between predicted and experimental strengths and weaknesses obtained from native hydrogen exchange data. SWOTein application to three test cases further showed its suitability to predict and interpret strong and weak residues in the context of folding, conformational changes and protein-protein binding. In summary, SWOTein is both fast and accurate and can be applied at small and large scale to analyze and modulate folding and molecular recognition processes. The SWOTein webserver provides the list of predicted strengths and weaknesses and a protein structure visualization tool that facilitates the interpretation of the predictions. It is freely available for academic use at http://babylone.ulb.ac.be/SWOTein/.

Specific Binding of the Nucleocapsid Protein Transforms the Folding Landscape of the HIV-1 TAR RNA Hairpin

Retroviral nucleocapsid (NC) proteins are nucleic acid chaperones that play a key role in the viral life cycle. During reverse transcription, HIV-1 NC destabilizes nucleic acids, including the transactivation response (TAR) hairpin to facilitate their structural re-arrangement into the lowest free energy conformations. By combining single molecule optical tweezers measurements with mfold-based free energy calculations of the unfolding landscape, we determine the equilibrium TAR stability and characterize the unfolding transition state in both the presence and absence of NC. Our results show that protein binding occurs specifically at guanines that border structural defects in the upper part of TAR hairpin. This results in preferential destabilization by NC of this part of TAR, thereby shifting the transition state closer to the bottom of the TAR stem, leading to a shorter critical unfolding length. Thus, NC facilitates TAR RNA annealing to its complementary DNA hairpin while having little effect on the stability of the final annealed duplexes. These results provide the first direct measurement of alterations in an RNA folding landscape that are induced by protein-RNA interactions, and are a novel case study of the effect of a protein that specifically destabilizes particular elements of the nucleic acid secondary structure.

Prediction of Protein Tertiary Structure using Genetic Algorithm

Proteins are essential for the biological processes in the human body. They can only perform their functions when they fold into their tertiary structure .Protein structure can be determined experimentally and computationally. Experimental methods are time consuming and high-priced and it is not always feasible to identify the protein structure experimentally. In order to predict the protein structure using computational methods, the problem is formulated as an optimization problem and the goal is to find the lowest free energy conformation. In this paper, Genetic Algorithm (GA) based optimization is used. This algorithm is adapted to search the protein conformational search space to find the lowest free energy conformation. Interestingly, the algorithm was able to find the lowest free energy conformation for a test protein (i.e. Met enkephalin) using ECEPP force fields.

Non-arrhenius behavior in the unfolding of a short, hydrophobic alpha-helix. Complementarity of molecular dynamics and lattice model simulations.

The unfolding of the last, C-terminal residue of AcNH(2)-(l-Leu)(11)-NHMe in its alpha-helical form has been investigated by measuring the variation of free energy involved in the alpha(R) to beta conformational transition. These calculations were performed using large-scale molecular dynamics simulations in conjunction with the umbrella sampling method. For different temperatures ranging from 280 to 370 K, the free energy of activation was estimated. Concurrently, unfolding simulations of a homopolypeptide formed by twelve hydrophobic residues were carried out, employing a three-dimensional lattice model description of the peptide, with a temperature-dependent interaction potential. Using a Monte Carlo approach, the lowest free energy conformation, an analogue of a right-handed alpha-helix, was determined in the region where the peptide chain is well ordered. The free energy barrier separating this state from a distinct, compact conformation, analogue to a beta-strand, was determined over a large enough range of temperatures. The results of these molecular dynamics and lattice model simulations are consistent and indicate that the kinetics of the unfolding of a hydrophobic peptide exhibits a non-Arrhenius behavior closely related to the temperature dependence of the hydrophobic effect. These results further illuminate the necessity to include a temperature dependence in potential energy functions designed for coarse-grained models of proteins.

The Dynamics, Structure, and Conformational Free Energy of Proline-Containing Antifreeze Glycoprotein

Recent NMR studies of the solution structure of the 14-amino acid antifreeze glycoprotein AFGP-8 have concluded that the molecule lacks long-range order. The implication that an apparently unstructured molecule can still have a very precise function as a freezing inhibitor seems startling at first consideration. To gain insight into the nature of conformations and motions in AFGP-8, we have undertaken molecular dynamics simulations augmented with free energy calculations using a continuum solvation model. Starting from 10 different NMR structures, 20 ns of dynamics of AFGP were explored. The dynamics show that AFGP structure is composed of four segments, joined by very flexible pivots positioned at alanine 5, 8, and 11. The dynamics also show that the presence of prolines in this small AFGP structure facilitates the adoption of the poly-proline II structure as its overall conformation, although AFGP does adopt other conformations during the course of dynamics as well. The free energies calculated using a continuum solvation model show that the lowest free energy conformations, while being energetically equal, are drastically different in conformations. In other words, this AFGP molecule has many structurally distinct and energetically equal minima in its energy landscape. In addition, conformational, energetic, and hydrogen bond analyses suggest that the intramolecular hydrogen bonds between the N-acetyl group and the protein backbone are an important integral part of the overall stability of the AFGP molecule. The relevance of these findings to the mechanism of freezing inhibition is discussed.

Conformational study of angiotensin II.

Conformational free energy calculations using an empirical potential ECEPP/3 (Empirical Conformational Energy Program for Peptides, Version 3) were carried out on angiotensin II (AII) of sequence Asp-Arg-Val-Tyr-Ile-His-Pro-Phe to find the stable conformations of the free state in the unhydrated and the hydrated states. A conformational analysis of the unhydrated state was carried out using the buildup procedure. The free energy calculation using the hydration shell model was also carried out to obtain the stable conformation of the hydrated state. The calculated stable conformations of AII in both states have a partially right-handed alpha-helical structure stabilized by short- and medium-range interactions. The similarity between the lowest free energy conformations of the unhydrated and hydrated states suggests that the hydration might not be important to stabilize the overall conformation of AII in a free state. The absence of any intramolecular interaction of the Tyr side chain suggests the possible interaction of this residue with the receptor. In this study, we found that the low free energy conformations contain both the parallel-plate and the perpendicular-plate geometries of the His and Phe rings, suggesting the coexistence of both conformations.