Solvent quality and nonbiological oligomer folding: revisiting conventional paradigms.

We report on extensive molecular dynamics atomistic simulations of a -substituted poly-phenylacetylene (pPA) foldamer dispersed in three solvents, water HO, cyclohexane cCH, and -hexane CH, and for three oligomer lengths 12mer, 16mer and 20mer. At room temperature, we find a tendency of the pPA foldamer to collapse into a helical structure in all three solvents but with rather different stability character, stable in water, marginally stable in -hexane, unstable in cyclohexane. In the case of water, the initial and final number of hydrogen bonds of the foldamer with water molecules is found to be unchanged, with no formation of intrachain hydrogen bonding, thus indicating that hydrogen bonding plays no role in the folding process.

Amino-Acid Characteristics in Protein Native State Structures.

The molecular machines of life, proteins, are made up of twenty kinds of amino acids, each with distinctive side chains. We present a geometrical analysis of the protrusion statistics of side chains in more than 4000 high-resolution protein structures. We employ a coarse-grained representation of the protein backbone viewed as a linear chain of C atoms and consider just the heavy atoms of the side chains.

A Tale of Two Chains: Geometries of a Chain Model and Protein Native State Structures.

Linear chain molecules play a central role in polymer physics with innumerable industrial applications. They are also ubiquitous constituents of living cells. Here, we highlight the similarities and differences between two distinct ways of viewing a linear chain.

III. Geometrical framework for thinking about globular proteins: Turns in proteins.

We have shown recently that the notion of poking pairwise interactions along a chain provides a unifying framework for understanding the formation of both secondary and the tertiary protein structure based on symmetry and geometry. α-helices and β-sheets are found to be special geometries that have systematic poking contacts in a repetitive manner with the contacts being local along the α-helix and non-local along a pair of adjacent strands within a β-sheet. Pairwise poking interactions also govern tertiary structure formation, but they are weaker and there are no special geometrical constraints as in secondary structure formation.

II. Geometrical framework for thinking about globular proteins: The power of poking.

Recently, we presented a framework for understanding protein structure based on the idea that simple constructs of holding hands or touching of objects can be used to rationalize the common characteristics of globular proteins. We developed a consistent approach for understanding the formation of the two key common building blocks of helices and sheets as well as the compatible assembly of secondary structures into the tertiary structure through the notion of poking pairwise interactions. Here we benchmark our predictions with a detailed analysis of structural data of over 4000 proteins from the Protein Data Bank.

A geometrical framework for thinking about proteins.

We present a model, based on symmetry and geometry, for proteins. Using elementary ideas from mathematics and physics, we derive the geometries of discrete helices and sheets. We postulate a compatible solvent-mediated emergent pairwise attraction that assembles these building blocks, while respecting their individual symmetries.

Spontaneous dimensional reduction and ground state degeneracy in a simple chain model.

Chain molecules play a key role in the polymer field and in living cells. Our focus is on a new homopolymer model of a linear chain molecule subject to an attractive self-interaction promoting compactness. We analyze the model using simple analytic arguments complemented by extensive computer simulations.

Building blocks of protein structures: Physics meets biology.

The native state structures of globular proteins are stable and well packed indicating that self-interactions are favored over protein-solvent interactions under folding conditions. We use this as a guiding principle to derive the geometry of the building blocks of protein structures-α helices and strands assembled into β sheets-with no adjustable parameters, no amino acid sequence information, and no chemistry. There is an almost perfect fit between the dictates of mathematics and physics and the rules of quantum chemistry.

Marginally compact phase and ordered ground states in a model polymer with side spheres.

We present the results of a quantitative study of the phase behavior of a model polymer chain with side spheres using two independent computer simulation techniques. We find that the mere addition of side spheres results in key modifications of standard polymer behavior. One obtains a marginally compact phase at low temperatures; the structures in this phase are reduced in dimensionality and are ordered, they include strands assembled into sheets and a variety of helices, and at least one of the transitions on lowering the temperature to access these ordered states is found to be first order.

Local sequence-structure relationships in proteins.

We seek to understand the interplay between amino acid sequence and local structure in proteins. Are some amino acids unique in their ability to fit harmoniously into certain local structures? What is the role of sequence in sculpting the putative native state folds from myriad possible conformations? In order to address these questions, we represent the local structure of each C atom of a protein by just two angles, θ and μ, and we analyze a set of more than 4,000 protein structures from the PDB. We use a hierarchical clustering scheme to divide the 20 amino acids into six distinct groups based on their similarity to each other in fitting local structural space.

Can the roles of polar and non-polar moieties be reversed in non-polar solvents?

Using thermodynamic integration, we study the solvation free energy of 18 amino acid side chain equivalents in solvents with different polarities, ranging from the most polar water to the most non-polar cyclohexane. The amino acid side chain equivalents are obtained from the 20 natural amino acids by replacing the backbone part with a hydrogen atom, and discarding proline and glycine that have special properties. A detailed analysis of the relative solvation free energies suggests how it is possible to achieve a robust and unambiguous hydrophobic scale for the amino acids.

Chain stiffness bridges conventional polymer and bio-molecular phases.

Chain molecules play important roles in industry and in living cells. Our focus here is on distinct ways of modeling the stiffness inherent in a chain molecule. We consider three types of stiffnesses-one yielding an energy penalty for local bends (energetic stiffness) and the other two forbidding certain classes of chain conformations (entropic stiffness).

Local symmetry determines the phases of linear chains: a simple model for the self-assembly of peptides.

We discuss the relation between the emergence of new phases with broken symmetry within the framework of simple models of biopolymers. We start with a classic model for a chain molecule of spherical beads tethered together, with the steric constraint that non-consecutive beads cannot overlap, and with a pairwise attractive square well potential accounting for the hydrophobic effect and promoting compaction. We then discuss the consequences of the successive breaking of spurious symmetries.

The elixir phase of chain molecules.

A phase of matter is a familiar notion for inanimate physical matter. The nature of a phase of matter transcends the microscopic material properties. For example, materials in the liquid phase have certain common properties independent of the chemistry of the constituents: liquids take the shape of the container; they flow; and they can be poured-alcohol, oil, and water as well as a Lennard-Jones computer model exhibit similar behavior when poised in the liquid phase.

Universal effects of solvent species on the stabilized structure of a protein.

We investigate the effects of solvent specificities on the stability of the native structure (NS) of a protein on the basis of our free-energy function (FEF). We use CPB-bromodomain (CBP-BD) and apoplastocyanin (apoPC) as representatives of the protein universe and water, methanol, ethanol, and cyclohexane as solvents. The NSs of CBP-BD and apoPC consist of 66% α-helices and of 35% β-sheets and 4% α-helices, respectively.

Vibrational entropy estimation can improve binding affinity prediction for non-obligatory protein complexes.

Predicting the binding affinity between protein monomers is of paramount importance for the understanding of thermodynamical and structural factors that guide the formation of a complex. Several numerical techniques have been developed for the calculation of binding affinities with different levels of accuracy. Approaches such as thermodynamic integration and Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) methodologies which account for well defined physical interactions offer good accuracy but are computationally demanding.

Unraveling protein folding mechanism by analyzing the hierarchy of models with increasing level of detail.

Taking protein G with 56 residues for a case study, we investigate the mechanism of protein folding. In addition to its native structure possessing α-helix and β-sheet contents of 27% and 39%, respectively, we construct a number of misfolded decoys with a wide variety of α-helix and β-sheet contents. We then consider a hierarchy of 8 different models with increasing level of detail in terms of the number of entropic and energetic physical factors incorporated.

Effective stiffness and formation of secondary structures in a protein-like model.

We use Wang-Landau and replica exchange techniques to study the effect of an increasing stiffness on the formation of secondary structures in protein-like systems. Two possible models are considered. In both models, a polymer chain is formed by tethered beads where non-consecutive backbone beads attract each other via a square-well potential representing the tendency of the chain to fold.

From polymers to proteins: the effect of side chains and broken symmetry on the formation of secondary structures within a Wang-Landau approach.

We use a micro-canonical Wang-Landau technique to study the equilibrium properties of a single flexible homopolymer where consecutive monomers are represented by impenetrable hard spherical beads tangential to each other, and non-consecutive monomers interact via a square-well potential. To mimic the characteristics of a protein-like system, the model is then refined in two different directions. Firstly, by allowing partial overlap between consecutive beads, we break the spherical symmetry and thus provide a severe constraint on the possible conformations of the chain.

Folding pathways of a knotted protein with a realistic atomistic force field.

We report on atomistic simulation of the folding of a natively-knotted protein, MJ0366, based on a realistic force field. To the best of our knowledge this is the first reported effort where a realistic force field is used to investigate the folding pathways of a protein with complex native topology. By using the dominant-reaction pathway scheme we collected about 30 successful folding trajectories for the 82-amino acid long trefoil-knotted protein.

The role of non-native interactions in the folding of knotted proteins: insights from molecular dynamics simulations.

For several decades, the presence of knots in naturally-occurring proteins was largely ruled out a priori for its supposed incompatibility with the efficiency and robustness of folding processes. For this very same reason, the later discovery of several unrelated families of knotted proteins motivated researchers to look into the physico-chemical mechanisms governing the concerted sequence of folding steps leading to the consistent formation of the same knot type in the same protein location. Besides experiments, computational studies are providing considerable insight into these mechanisms.

Unfolding thermodynamics of cysteine-rich proteins and molecular thermal-adaptation of marine ciliates.

Euplotes nobilii and Euplotes raikovi are phylogenetically closely allied species of marine ciliates, living in polar and temperate waters, respectively. Their evolutional relation and the sharply different temperatures of their natural environments make them ideal organisms to investigate thermal-adaptation. We perform a comparative study of the thermal unfolding of disulfide-rich protein pheromones produced by these ciliates.

The role of non-native interactions in the folding of knotted proteins.

Stochastic simulations of coarse-grained protein models are used to investigate the propensity to form knots in early stages of protein folding. The study is carried out comparatively for two homologous carbamoyltransferases, a natively-knotted N-acetylornithine carbamoyltransferase (AOTCase) and an unknotted ornithine carbamoyltransferase (OTCase). In addition, two different sets of pairwise amino acid interactions are considered: one promoting exclusively native interactions, and the other additionally including non-native quasi-chemical and electrostatic interactions.

Dominant folding pathways of a WW domain.

We investigate the folding mechanism of the WW domain Fip35 using a realistic atomistic force field by applying the Dominant Reaction Pathways approach. We find evidence for the existence of two folding pathways, which differ by the order of formation of the two hairpins. This result is consistent with the analysis of the experimental data on the folding kinetics of WW domains and with the results obtained from large-scale molecular dynamics simulations of this system.