On the liquid demixing of water + elastin-like polypeptide mixtures: bimodal re-entrant phase behaviour.

Water + elastin-like polypeptides (ELPs) exhibit a transition temperature below which the chains transform from collapsed to expanded states, reminiscent of the cold denaturation of proteins. This conformational change coincides with liquid-liquid phase separation. A statistical-thermodynamics theory is used to model the fluid-phase behavior of ELPs in aqueous solution and to extrapolate the behavior at ambient conditions over a range of pressures. At low pressures, closed-loop liquid-liquid equilibrium phase behavior is found, which is consistent with that of other hydrogen-bonding solvent + polymer mixtures. At pressures evocative of deep-sea conditions, liquid-liquid immiscibility bounded by two lower critical solution temperatures (LCSTs) is predicted. As pressure is increased further, the system exhibits two separate regions of closed-loop of liquid-liquid equilibrium (LLE). The observation of bimodal LCSTs and two re-entrant LLE regions herald a new type of binary global phase diagram: Type XII. At high-ELP concentrations the predicted phase diagram resembles a protein pressure denaturation diagram; possible "molten-globule"-like states are observed at low concentration.

Navigating in foldonia: Using accelerated molecular dynamics to explore stability, unfolding and self-healing of the ß-solenoid structure formed by a silk-like polypeptide.

The ß roll molecules with sequence (GAGAGAGQ)10 stack via hydrogen bonding to form fibrils which have been themselves been used to make viral capsids of DNA strands, supramolecular nanotapes and pH-responsive gels. Accelerated molecular dynamics (aMD) simulations are used to investigate the unfolding of a stack of two ß roll molecules, (GAGAGAGQ)10, to shed light on the folding mechanism by which silk-inspired polypeptides form fibrils and to identify the dominant forces that keep the silk-inspired polypeptide in a ß roll configuration. Our study shows that a molecule in a stack of two ß roll molecules unfolds in a step-wise fashion mainly from the C terminal. The bottom template is found to play an important role in stabilizing the ß roll structure of the molecule on top by strengthening the hydrogen bonds in the layer that it contacts. Vertical hydrogen bonds within the ß roll structure are considerably weaker than lateral hydrogen bonds, signifying the importance of lateral hydrogen bonds in stabilizing the ß roll structure. Finally, an intermediate structure was found containing a ß hairpin and an anti-parallel ß sheet consisting of strands from the top and bottom molecules, revealing the self-healing ability of the ß roll stack.

Predicting the Fluid-Phase Behavior of Aqueous Solutions of ELP (VPGVG) Sequences Using SAFT-VR.

The statistical associating fluid theory for potentials of variable range (SAFT-VR) is used to predict the fluid phase behavior of elastin-like polypeptide (ELP) sequences in aqueous solution with special focus on the loci of lower critical solution temperatures (LCSTs). A SAFT-VR model for these solutions is developed following a coarse-graining approach combining information from atomistic simulations and from previous SAFT models for previously reported relevant systems. Constant-pressure temperature-composition phase diagrams are determined for solutions of (VPGVG)n sequences + water with n = 1 to 300. The SAFT-VR equation of state lends itself to the straightforward calculation of phase boundaries so that complete fluid-phase equilibria can be calculated efficiently. A broad range of thermodynamic conditions of temperature and pressure are considered, and regions of vapor-liquid and liquid-liquid coexistence, including LCSTs, are found. The calculated phase boundaries at low concentrations match those measured experimentally. The temperature-composition phase diagrams of the aqueous ELP solutions at low pressure (0.1 MPa) are similar to those of types V and VI phase behavior in the classification of Scott and van Konynenburg. An analysis of the high-pressure phase behavior confirms, however, that a closed-loop liquid-liquid immiscibility region, separate from the gas-liquid envelope, is present for aqueous solutions of (VPGVG)30; such a phase diagram is typical of type VI phase behavior. ELPs with shorter lengths exhibit both liquid-liquid and gas-liquid regions, both of which become less extensive as the chain length of the ELP is decreased. The strength of the hydrogen-bonding interaction is also found to affect the phase diagram of the (VPGVG)30 system in that the liquid-liquid and gas-liquid regions expand as the hydrogen-bonding strength is decreased and shrink as it is increased. The LCSTs of the mixtures are seen to decrease as the ELP chain length is increased.

The Composition of the Mobile Phase Affects the Dynamic Chiral Recognition of Drug Molecules by the Chiral Stationary Phase.

More than half of all pharmaceuticals are chiral compounds. Although the enantiomers of chiral compounds have the same chemical structure, they can exhibit marked differences in physiological activity; therefore, it is important to remove the undesirable enantiomer. Chromatographic separation of chiral enantiomers is one of the best available methods to get enantio-pure substances, but the optimization of the experimental conditions can be very time-consuming. One of the most widely used chiral stationary phases, amylose tris(3,5-dimethylphenyl carbamate) (ADMPC), has been extensively investigated using both experimental and computational methods; however, the dynamic nature of the interaction between enantiomers and ADMPC, as well as the solvent effects on the ADMPC-enantiomer interaction, are currently absent from models of the chiral recognition mechanism. Here we use QM/MM and molecular dynamics (MD) simulations to model the enantiomers of flavanone on ADMPC in either methanol or heptane/2-propanol (IPA) (90/10) to elucidate the chiral recognition mechanism from a new dynamic perspective. In atomistic MD simulations, the 12-mer model of ADMPC is found to hold the 4/3 left-handed helical structure in both methanol and heptane/IPA (90/10); however, the ADMPC polymer is found to have a more extended average structure in heptane/IPA (90/10) than in methanol. This results from the differences in the distribution of solvent molecules close to the backbone of ADMPC leads to changes in the distribution of the (φ, ψ) dihedral angles of the glycoside bond (between adjacent monomers) that define the structure of the polymer. Our simulations have shown that the lifetime of hydrogen bonds formed between ADMPC and flavanone enantiomers in the MD simulations are able to reproduce the elution order observed in experiments for both the methanol and the heptane/IPA solvent systems. Furthermore, the ratios of hydrogen-bonding-lifetime-related properties also capture the solvent effects, in that heptane/IPA (90/10) is found to make the separation between the two enantiomers of flavanone less effective than methanol, which agrees with the experimental separation factors of 0.9 versus 0.4 for R/S, respectively.

LCST Behavior is Manifested in a Single Molecule: Elastin-Like polypeptide (VPGVG)n.

The physical origin of the lower critical solution temperature (LCST) behavior of a variety of fluids, including elastin-like polypeptides (ELPs), has been studied for the past few decades. As is the case for polymer solutions, LCST behavior of ELPs is invariably reported for large systems of molecules and is considered evidence for collective behavior. In contrast, we find evidence for properties changes associated with LCST behavior in a single molecule by performing long atomic-level molecular dynamics simulation on the ELP sequences (Val-Pro-Gly-Val-Gly)n for four different length peptides over a wide range of temperatures. We observe a sharp transition in the number of hydrogen bonds between peptide and water and in the number of water molecules within the first hydration shell as temperature rises; this is used to locate the transition temperature. The dependence of the transition temperatures of ELPs on their lengths agrees well with experiments in that both have the same power law exponents. Our simulations reveal that the tendency for pentamers (VPGVG) in ELPs of all lengths to lose H-bonds with water or to gain H-bonds with themselves as temperature rises is independent of the length of the chain in which they are embedded. Thus, the transition temperature of ELPs in pure water is determined by two factors: the hydrogen bonding tendency of the pentamers and the number of pentamers per ELP. Moreover, the hydrogen bonding tendency of pentamers depends only on their sequences, not on the ELP chain length.

Dock 'n roll: folding of a silk-inspired polypeptide into an amyloid-like beta solenoid.

Polypeptides containing the motif ((GA)mGX)n occur in silk and have a strong tendency to self-assemble. For example, polypeptides containing (GAGAGAGX)n, where X = G or H have been observed to form filaments; similar sequences but with X = Q have been used in the design of coat proteins (capsids) for artificial viruses. The structure of the (GAGAGAGX)m filaments has been proposed to be a stack of peptides in a ß roll structure with the hydrophobic side chains pointing outwards (hydrophobic shell). Another possible configuration, a ß roll or ß solenoid structure which has its hydrophobic side chains buried inside (hydrophobic core) was, however, overlooked. We perform ground state analysis as well as atomic-level molecular dynamics simulations, both on single molecules and on two-molecule stacks of the silk-inspired sequence (GAGAGAGQ)10, to decide whether the hydrophobic core or the hydrophobic shell configuration is the most stable one. We find that a stack of two hydrophobic core molecules is energetically more favorable than a stack of two hydrophobic shell molecules. A shell molecule initially placed in a perfect ß roll structure tends to rotate its strands, breaking in-plane hydrogen bonds and forming out-of-plane hydrogen bonds, while a core molecule stays in the ß roll structure. The hydrophobic shell structure has type II' ß turns whereas the core configuration has type II ß turns; only the latter secondary structure agrees well with solid-state NMR experiments on a similar sequence (GA)15. We also observe that the core stack has a higher number of intra-molecular hydrogen bonds and a higher number of hydrogen bonds between stack and water than the shell stack. Hence, we conclude that the hydrophobic core configuration is the most likely structure. In the stacked state, each peptide has more intra-molecular hydrogen bonds than a single folded molecule, which suggests that stacking provides the extra stability needed for molecules to reach the folded state.

Simulation study of the ability of a computationally-designed peptide to recognize target tRNALys3 and other decoy tRNAs.

In this paper, we investigate the ability of our computationally-designed peptide, Pept10 (PNWNGNRWLNNCLRG), to recognize the anticodon stem and loop (ASL) domain of the hypermodified tRNALys3 (mcm5 s2 U34 ,ms2 t6 A37 ), a reverse transcription primer of HIV replication. Five other ASLs, the singly modified ASLLys3 (ms2 t6 A37 ), ASLLys3 (s2 U34 ), ASLLys3 (Ψ39 ), ASLLys1,2 (t6 A37 ), and ASLGlu (s2 U34 ), were used as decoys. Explicit-solvent atomistic molecular dynamics simulations were performed to examine the process of binding of Pept10 with the target ASLLys3 (mcm5 s2 U34 ,ms2 t6 A37 ) and the decoy ASLs. Simulation results demonstrated that Pept10 is capable of recognizing the target ASLLys3 (mcm5 s2 U34 ,ms2 t6 A37 ) as well as one of the decoys, ASLLys3 (Ψ39 ), but screens out the other four decoy ASLs. The interchain van der Waals (VDW) and charge-charge (ELE + EGB) energies for the two best complexes were evaluated to shed light on the molecular recognition mechanism between Pept10 and ASLs. The results indicated that Pept10 recognizes and binds to the target ASLLys3 (mcm5 s2 U34 ,ms2 t6 A37 ) through residues W3 and R7 which interact with the nucleotides mcm5 s2 U34 , U35 , and ms2 t6 A37 via the interchain VDW energy. Pept10 also recognizes the decoy ASLLys3 (Ψ39 ) through residue R14 which contacts the nucleotide U36 via the interchain VDW energy. Regardless of the type of ASL, the positively charged arginines on Pept10 are attracted to the negatively charged phosphate linkages on the ASL via the interchain ELE + EGB energy, thereby enhancing the binding affinity.