Sequence directionality dramatically affects LCST behavior of elastin

Apr 17, 2018 - In this paper, the thermal responsive properties of the canonical ELP, poly(VPGVG), and its reverse sequence poly(VGPVG) were investiga...
1 downloads 7 Views 2MB Size
Subscriber access provided by UNIV OF DURHAM

Sequence directionality dramatically affects LCST behavior of elastin-like polypeptides Nan K Li, Stefan Roberts, Felipe Garcia Quiroz, Ashutosh Chilkoti, and Yaroslava G Yingling Biomacromolecules, Just Accepted Manuscript • DOI: 10.1021/acs.biomac.8b00099 • Publication Date (Web): 17 Apr 2018 Downloaded from http://pubs.acs.org on April 18, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Sequence directionality dramatically affects LCST behavior of elastin-like polypeptides Nan K. Li†, Stefan Roberts‡, Felipe Garcia Quiroz$, Ashutosh Chilkoti‡, and Yaroslava G. Yingling†* †Department of Materials Science and Engineering, North Carolina State University, 911 Partners Way, Raleigh, NC 27695, USA. ‡ Department of Biomedical Engineering, PO Box 90281, Duke University, Durham, North Carolina 27708, USA $ Howard Hughes Medical Institute, Laboratory of Mammalian Cell Biology & Development, The Rockefeller University, New York, New York 10065, USA

ABSTRACT

Elastin-like polypeptides (ELP) exhibit an inverse temperature transition or lower critical solution temperature (LCST) transition phase behavior in aqueous solutions. In this paper, the thermal responsive properties of the canonical ELP, poly(VPGVG), and its reverse sequence poly(VGPVG) were investigated by turbidity measurements of the cloud point behavior, circular dichroism (CD) measurements, and all-atom molecular dynamics (MD) simulations to gain a molecular understanding of mechanism that controls hysteretic phase behavior. It was shown

ACS Paragon Plus Environment

1

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 46

experimentally that both poly(VPGVG) and poly(VGPVG) undergo a transition from soluble to insoluble in aqueous solution upon heating above the transition temperature (Tt). However, poly(VPGVG) resolubilizes upon cooling below its Tt, whereas the reverse sequence, poly(VGPVG), remains aggregated despite significant undercooling below the Tt. The results from MD simulations indicated that a change in sequence order results in significant differences in the dynamics of the specific residues, especially valines, which lead to extensive changes in the conformations of VPGVG and VGPVG pentamers and, consequently, dissimilar propensities for secondary structure formation and overall structure of polypeptides. These changes affected the relative hydrophilicities of polypeptides above Tt, where poly(VGPVG) is more hydrophilic than poly(VPGVG) with more extended conformation and larger surface area, which led to formation of strong interchain hydrogen bonds responsible for stabilization of the aggregated phase and the observed thermal hysteresis for poly(VGPVG).

KEYWORDS. Molecular dynamics simulations, hysteresis, ELP, VPGVG, VGPVG, LCST, sequence order, syntax

ACS Paragon Plus Environment

2

Page 3 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Introduction Elastin-like polypeptides (ELP) are a class of synthetic polypeptides derived from elastin, which undergo a transition from soluble to insoluble in aqueous solution at a lower critical solution temperature (LCST).1 ELPs, at a given concentration in aqueous solution, are soluble below their LCST cloud point, or transition temperature (Tt), and undergo hydrophobic collapse and form insoluble aggregates at temperatures greater than the Tt.2,

3

ELPs are characterized by their

temperature responsiveness that can be tuned at the sequence level and their biocompatibility, which make them interesting as a class of synthetic biopolymers for a variety of applications including drug delivery engineering scaffolds

4-6

11, 12

, protein purification7-9, responsive nanoparticles10, and tissue

. The most widely studied ELP is poly(VPGVG) whose LCST

behavior as a function of concentration13-15 and ionic strength16 in aqueous solutions has been extensively characterized. This pentapeptide motif was also generalized into the VPGXG repeat unit by Urry, where the fourth residue (X) in the pentapeptide repeat is often termed the “guest’ residue that can be substituted with other amino acids except proline.14, 17 ELPs that exhibit LCST behavior may also display a hysteresis loop in their temperaturedependent turbidimetry profiles. For example, it has been long known that poly(VPAVG) displays significant thermal hysteresis in its heating-cooling cycle,18,

19

so that considerable

undercooling below the Tt is needed to resolubilize poly(VPAVG). We also recently investigated the LCST phase behavior of a large library of polypeptides composed of LCST motifs, which uncovered additional forms of hysteretic phase behaviors and large differences in thermal hysteresis for polypeptides with very subtle changes in amino acid sequence20. While this thermal hysteresis behavior opens exciting possibilities for the design of materials with new

ACS Paragon Plus Environment

3

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 46

properties, the molecular determinants of how amino acid sequence governs hysteretic phase behavior remain poorly understood. The amino acid composition of an ELP is not the only sequence dependent factor that can affect its LCST behavior. Changes in amino acid order are likely to disrupt native secondary structure propensities of well-folded polypeptides, which can have a major effect on their function. However, for intrinsically disordered ELPs that exist in dynamic and flexible conformations, the rationalization of how changes in amino acid order influence their functions and phase behaviors are far from obvious. To provide insights into how the amino acid order controls LCST behavior of ELPs, we examined a unique pair of polypeptides20, poly(VPGVG) and poly(VGPVG), whose sequences are essentially mirror images or “reverse” motifs and hence identical if read in the opposite directions, N- to C-terminus or vice versa. A comparison of the LCST phase behavior of poly(VPGVG) and poly(VGPVG) allows us to examine the dependence between the arrangement of amino acids and LCST phase behavior with no change in the overall composition. The thermal behavior of aqueous solutions of these two polymers were examined using turbidity measurements and all-atom molecular dynamics (MD) simulations. Although the turbidity measurements can effectively describe the macroscopic “cloud point” behavior of an ELP, changes in the molecular properties associated with sequence reversal are less easy to characterize experimentally. In our previous study, the effect of temperature on the structure, dynamics and association of (VPGVG)18 in aqueous solution was investigated using atomistic MD simulations. And from that study we concluded that the LCST phase behavior of poly(VPGVG) is a collective phenomenon that originates from the correlated gradual changes in single polypeptide structure and the abrupt change in properties of hydration water around the peptide and is a result of a competition between peptide−peptide and peptide−water

ACS Paragon Plus Environment

4

Page 5 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

interactions.21 We hence carried out MD simulations to reveal how sequence reversal impacts the thermal behaviors of ELPs (VPGVG)18 and (VGPVG)18. Moreover, we attempt to also elucidate the mechanism underlying the hysteretic behavior observed in the LCST phase transition of (VGPVG)18. Materials and Methods Molecular Dynamics Simulations In this study, we employed the Amber 1122 and 1223 programs to run MD simulations. The ff99SB force field was employed for peptides and water was modeled using the TIP3P model24. The force field chosen here was demonstrated to be able to reproduce the inverse temperature transition properties of poly(VPGVG) in our previous studies.21, 25 The non-bonded interactions were truncated at 9 Å cutoff with a 0.00001 tolerance for Ewald convergence and the long-range electrostatic interactions were taken into account by Particle Mesh Ewald (PME) summation26. Simulations of a single poly(VPGVG) and poly(VGPVG) peptides in water were performed at ten different temperatures between 290 K and 350 K, which span the temperature range relevant to the experimental study of LCST phase transitions. The initial structures for the models of ELPs were built based on dihedral angles for each residue in Urry’s β-spiral model27, which was proposed based on pioneering work in the development of ELP by chemically synthesizing poly(VPGVG) by Urry’s group and widely used as the initial structure by previous MD studies on ELPs.21, 25, 28, 29 The initial structures were then solvated in explicit water with the closest distance between any polypeptide atom and the edge of the periodic box to be at least 8 Å to ensure that the polypeptide did not cross periodic boundaries and interact with its own image. However, the initial structure was completely distorted at the very early stage of MD trajectory. The simulation box for (VPGVG)18 contained 6856 water molecules with a size of 5.9 nm 7.7

ACS Paragon Plus Environment

5

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 46

nm 5.0 nm in the equilibrium state and the simulation box for (VGPVG)18 contained around 8651 water molecules with a size of 6.9 nm 8.7 nm

5.7 nm in the equilibrium state. The

equilibrium box size varied slightly at different temperatures. The system was then equilibrated in eight stages, including several energy minimization, heating and MD run cycles. The equilibration protocol was the same as described in our previous study21. The system temperature was maintained using the Berendsen thermostat.30 The SHAKE algorithm was used to constrain the position of the hydrogen atoms.31 The production simulations were performed for at least 70 ns with a 2 fs time step. With a sampling rate of one frame per 1 ps (every 500 time steps) of the MD trajectories, we would expect that not only enough samples were collected, but also the simulation efficiency was enhanced. To ensure the convergence of all simulations, we calculated the time autocorrelation function for in-plane and out-of-plane backbone rotations, from which the relaxation time of the peptide backbone was estimated to be less than 25 ns (Figure S5, S6). The statistical and clustering analyses were carried out on MD trajectories for the last 40 ns using in-house scripts along with the tool suite accompanying Amber12.0. The interaction energies were calculated using the molecular mechanics energy function in NAMD 2.732. The hydrogen bond analyses were performed using an angle cut-off of 30 and a distance cutoff of 3.5 Å. To identify the most probable peptide structures from the MD simulations trajectory, the hierarchical RMSD-based clustering algorithm33 was used. Specifically, the last 40 ns of a trajectory from the single-peptide MD simulations were clustered to produce three structural clusters using the pairwise RMSD between frames as a metric and comparing the carbon alpha atoms with a critical distance of 12 Å. Representative snapshots of the ELP structures from the most populated clusters at 290 K and 350 K are shown in Figures S3 and S4. The lowest-energy

ACS Paragon Plus Environment

6

Page 7 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

representative single polypeptide structures from these clusters were then taken to be the initial structures in the simulations of the interaction between two polypeptides. The two selected single polypeptide structures were then placed side-by-side with a 32.3~32.4 Å distance between the centers of mass of the polypeptides and a 7 Å distance between the two closest surfaces. The structures were then solvated in explicit solvent with the closest distance between any solute atom and the edge of the periodic box specified as 12 Å and simulated at 350 K. The simulation box contained 10493 water molecules with the size of 9.2 nm 7.9 nm

5.8 nm for double-

poly(VPGVG) system and 9738 water molecules with size of 6.6 nm 9.4 nm 6.4 nm for double-poly(VGPVG) system. Equilibration protocols and MD simulations were the same as described for single-peptide simulations21. The production steps of the simulations were carried out at 350 K for more than 15 ns with a 2 fs time step and performed three times. Experimental protocol Recombinant synthesis and characterization of polypeptides: We synthesized genes encoding polymers (VPGVG)40, (VPGVG)80, (VGPVG)40 and (VGPVG)80 using recursive direction ligation by plasmid reconstruction and expressed them in E. coli BL21(DL3, Edgebio) as previously described34,35. Polymers of VPGVG were purified by inverse transition cycling (ITC)36. Polymers of VGPVG were purified from the insoluble fraction after sonication and centrifugation using 6 M guanidinium chloride followed by dialysis into water and centrifugation to remove insoluble proteins. Polymers of VGPVG, which remain soluble upon dialysis, were further purified by cycles of ITC using NaCl instead of temperature to trigger the phase transition37. The LCST phase behaviors of the purified polymers were studied by temperaturedependent turbidimetry using a Cary UV-vis spectrophotometer (Agilent) at a wavelength of 350 nm and heating/cooling rates of 1 C/min, except where otherwise specified. CD experiments

ACS Paragon Plus Environment

7

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 46

were carried out using an Aviv Model 202 instrument and 1 mm quartz cells (Hellma USA). Scans were carried out in H2O with a polymer concentration of 10 µM. Polymers were scanned in triplicate from 190 -260 nm in 1 nm steps with a 1 s averaging time. Data points with a dynode voltage above 500 V were removed from analysis.

ACS Paragon Plus Environment

8

Page 9 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Figure 1: Experimental comparison of LCST phase behavior of poly(VPGVG) and poly(VGPVG). (a,b) Temperature dependent turbidimetry for (VPGVG)40/80 and (VGPVG)40/80 (50 µM, PBS). Polymers with identical composition but sequencereversed orientation (reverse sequence) have different LCST phase behaviors. (VPGVG)40 and (VPGVG)80 —canonical ELP sequences—shows reversible phase behavior, but the reverse sequence, poly(VGPVG), shows irreversible aggregation as a function of solution temperature. (c) The reverse sequence VGPVG has a slightly higher transition temperature across all measured concentrations despite the identical amino acid composition. (d) CD measurements of the soluble polymers (10 µM at 15 C in H2O) reveal a significant difference in their secondary structure. Poly(VPGVG) has a negative “shelf” at 218 nm, indicating the presence of beta turns 35, 36, and this peak is absent in the CD spectrum of the reverse sequence, poly(VGPVG). Both polymers otherwise appear highly disordered. Results Experimental results Turbidity measurements for poly(VPGVG) and poly(VGPVG) reveal sharp differences in their LCST phase behavior (Figures 1a,b and Figures S1a,b). For both the 40-mer and 80-mer lengths,

ACS Paragon Plus Environment

9

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 46

the two polymers have a similar Tt upon heating; however, poly(VPGVG) is fully reversible whereas the reverse sequence poly(VGPVG) remains aggregated despite significant undercooling below its Tt. Above its Tt, poly(VPGVG)40 and poly(VPGVG)80 aggregates also show signs of sedimentation (the optical density drops gradually as the coacervates settle below the UV light path). This sedimentation is not observed for poly(VGPVG) polymers on the time scale of our experiments, indicating that the aggregates are either smaller or less dense than those

ACS Paragon Plus Environment

10

Page 11 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

of poly(VPGVG). The irreversible aggregation observed for poly(VGPVG) is also independent of both concentration (Fig. S2), and heating and cooling rates, with a nearly identical Tt and complete irreversibility observed for heating and cooling rates ranging from 0.5– 5

C/min (Fig.

S1c-f). Interestingly, the Tt for poly(VPGVG) at both 40 and 80-mers tends to be slightly higher than that for poly(VGPVG), despite identical amino acid composition (Fig. 1c). This difference in transition temperatures, consistent for all concentrations and for both polymer chain lengths, is

ACS Paragon Plus Environment

11

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 46

likely due to differences in polypeptide conformation. CD measurements for both polymers reveal slightly different molar ellipticity peaks (Fig. 1d and Fig. S1g). Poly(VPGVG) appears largely disordered with a slight shelf at 218 nm which has historically been attributed to unstable Pro-Gly beta turns in the polymer

38, 39

. Poly(VGPVG), which has a reversed Gly-Pro sequence,

is also largely disordered but does not exhibit a peak at 218 nm in its CD spectrum, indicating a

ACS Paragon Plus Environment

12

Page 13 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

slight shift in the secondary structure likely produced from disruption of the Pro-Gly beta-turn structure.

Figure 2. Temperature dependence of the properties of a single (VPGVG)18 and (VGPVG)18 ELP chain: (a) radius of gyration (Rg), (b) npw, the number of peptide-water hydrogen bonds, (c) solvent accessible surface area (SASA, (d) peptide-water interaction energy. Error bars represent the standard deviation. Structural properties and water-peptide interactions To explain the differences observed in the experiments, we examined dynamics of forward sequence (VPGVG)18 and reverse sequence (VGPVG)18 in water as a function of temperature using MD simulations. We chose 18 repeats because the forward sequence (VPGVG)18 has been

ACS Paragon Plus Environment

13

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 46

investigated systematically and has been shown to exhibit inverse transition thermal behaviors.21 The evolution of single chain Rg for (VPGVG)18 and (VGPVG)18 show that both peptides adopt a collapsed state over a temperature range between 290 K and 350 K (Fig. 2, S3, S4).21 However, while the Rg of (VPGVG)18 gradually decreases, the Rg for the reverse sequence (VGPVG)18 stays approximately constant as the temperature increases. The comparison of the number of peptide-peptide (npp) and water mediated peptide-water-peptide (npwp) intrachain hydrogen bonds between the two polymers (Figure S7a) reveals that for (VPGVG)18, npwp monotonically decreases while npp increases as a function of temperature. However, for (VGPVG)18 system all observables vary only slightly upon heating. The temperature induced changes in the solvent accessible surface area (SASA) and number of water-peptide hydrogen bonds (npw) from the MD simulations (Fig 1b.c) indicated that for (VPGVG)18, SASA and npw decrease gradually upon heating until a sudden drop at around 330 K, which we previously connected to the Tt21. In contrast, for (VGPVG)18, the dependence of SASA and npw on the temperature is gradual with no abrupt change (Fig. 1 b,c). At high temperatures (VGPVG)18 has a greater solvent accessible surface area (SASA) and forms more hydrogen bonds with water molecules than (VPGVG)18, which indicate that the reversal of the amino acid order makes (VPGVG)18 more hydrophilic at elevated temperatures than (VGPVG)18. In our previous studies, we identified that the number of water molecules in the first hydration shell is an important parameter for characterization of the phase behavior, which is defined as the number of water molecules within a 2.23 Å distance of any atom on the polypeptide21, 25. This distance is taken as the first minima from the radial distribution function (RDF) between the oxygen of water and an atom on the polypeptide.21 Below 330 K, the number of first layer hydration water molecules around (VPGVG)18 and (VGPVG)18 is very similar (Fig. S7b). Above

ACS Paragon Plus Environment

14

Page 15 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

330 K, however, the differences between the two polymers become more significant, as (VPGVG)18 exhibits a sharp decrease which is correlated with a similar decrease in the SASA, unlike (VGPVG)18. The interaction energy between polypeptides and water provides information for the hydrophobicity change of a single peptide with temperature. As shown on Figure 2d, polypeptide-water interactions become energetically more unfavorable as the temperature increases. Moreover, at high temperatures the difference in the interaction energies between polypeptides indicates that the surface hydrophobicity of a single molecule of (VGPVG)18 is less than that of (VPGVG)18. Overall, reversing the amino acid order appears to cause minor structural and energetic changes in ELPs at low temperatures. However, the poly(VGPVG) is larger (higher Rg) and more hydrophilic (larger SASA, npw, and lower peptide-water interaction energy) than poly(VPGVG) at high temperatures, although their overall amino acid compositions are identical. Dihedral angles, secondary structures and role of individual residues. To explore the structural differences associated with different temperature-responsive behaviors of (VPGVG)18 and (VGPVG)18 in aqueous solution, we calculated the torsion angles of each amino acid (Figure 3 and S8). The high intensity region in the Ramachandran plot indicates the most favorable, low-energy torsion angle pairs. In our previous study, we discussed the features of torsion angles in VPGVG pentamer at low and high temperatures and compared our simulation results with multiple experimental studies in details.21 Overall, our simulations show that for (VPGVG)18, the change in dihedral angle as a function of temperature agrees well with the experimental observations.

40, 41, 42

Here, we found that structural dynamics of backward

sequence, (VGPVG)18 showed similarities but also significant differences to the forward

ACS Paragon Plus Environment

15

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 46

sequence, (VPGVG)18, especially for valine residues (Figure 3). Val1 in VPGVG is preceded by proline and shows torsion angle distribution that is distinctive and different from the other valines in both VPGVG and VGPVG pentamer. The ψ angle of Val1 in VPGVG only populates around 150° and φ is around −130° or −75°, whereas Val4’ torsion angles are distributed between four regions where ψ is around 150° or −2° and φ is around −130° or −90°. However, Val1 and Val4 in VGPVG occupy similar conformational spaces in the left side of the Ramachandran plots. The ψ angle display the bimodal distribution with (ψ ~150° or −5°, respectively) while ϕ ranges from ~−45° to ~−170° with one peak around −126° for a4_1, a4_2, and d4_1, and −107.7° for d4_2. The torsion angles of Val4 residues in VGPVG also clustered at

ACS Paragon Plus Environment

16

Page 17 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

the left-handed alpha helical conformational space marked as d4_3 (52°, 30.3°). Glycine and proline are usually observed to have very different ϕ, ψ distribution outlines than the other amino acids, since glycine is conformationally less constrained and proline is conformationally more constrained. 43 In VPGVG, the Pro2 and Gly3 primarily adopt the torsion angle pairs (φ = −60°, ψ = 120°), (φ = −75°, ψ = 150°) as the i+1 and i+2 residues of a type-II βturn structure. Gly2 in VGPVG is also a pre-proline residue; however, Gly2 does not have a side chain which allows high flexibility as well as the torsion angles which are normally not allowed for other amino acid residues. Moreover, the most populated region for Gly3 in VPGVG is missing in the plots for Gly2 in VGPVG, indicating missing turn structures in the segment

Figure 3. (a-d) Ramachandran plot (ϕ, ψ distributions) for Val1 and Val4 in (VPGVG)18 and (VGPVG)18, respectively, at 350 K, colored by intensities. The most populated regions on the Ramachandran plots are labelled. (e-f) Representative structures for (e) G(VPGVG)V segment and (f) G(VGPVG)V segment at 350 K built with the dihedral angles obtained from the highest intensity points in Ramachandran plots (Figure 3 and Figure S8).

ACS Paragon Plus Environment

17

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 46

VGPV. The most populated angles adopted by Gly2 in VGPVG correspond to the high intensity regions of the Ramachandran plot labeled by c3_1 (−76.5°, 150.6°), c3_2 (70°, −141.1°) and c4_1 (98.27°, −140.0°) as shown in Figure S8. The Gly5 in VGPVG are scattered over various conformational spaces of the Ramachandran plot, which may adopt various secondary structural motifs such as Helix, Turn and β-strand. Interestingly, the Ramachandran plot pattern for Gly5 in poly(VGPVG) at 350 K is similar to Gly5 in poly(VPGVG) at 290 K; the Ramachandran plot pattern for Gly5 in poly(VGPVG) at 290K is similar to Gly5 in poly(VPGVG) at 350 K. Gly5 residue in poly(VGPVG) at 290 K and 350 K occupy two areas with relatively high number of occurrences, such as e3_1 (−80.76°,135°) and e4_2 (−107.7°, −2.0°), respectively. In contrast to glycine amino acids, ψ angle for Pro3 in VGPVG varies in a wider range from −62° to 173° which indicates a larger flexibility of proline residues in VGPVG than that in VPGVG, where φ is fixed around −65° due to its backbone stereochemistry. To indicate the difference in the most probable conformation we build the representative structures for pentamers of VPGVG and VGPVG based on the most populated torsion angles obtained from Ramachandran plots, as shown in Figure 3e and f. The VGPVG segment has a more stretched out conformation than VPGVG at 350K (Fig. 3 e,f), which can explain the observed differences in radius of gyration. The secondary structure motifs were determined by a common structure recognition algorithm, the DSSP method, which is based on H-bonding patterns44. In the DSSP algorithm, the residues bracketed by the hydrogen bond between residue i and residue i+n (n=3, 4, 5) are considered to adopt a turn structure. The helical structure includes residues from α- , 3-10 and π helix; the βstrand structure includes residues from β-bridges and extended strands involved in parallel and anti- parallel β-sheets. The secondary structural propensities of each single polypeptide at various temperatures are listed in Table S2. Distinct changes were observed due to the reversing

ACS Paragon Plus Environment

18

Page 19 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

of the amino acid order in the secondary structural propensities, which is in agreement with the experimental observations (Figure 1). In general, there is a larger propensity for helix structures in the reverse sequence, (VGPVG)18 chain than in the forward sequence (VPGVG)18 at all temperatures, which is consistent with the observations from the Ramachandran plot. In addition, (VPGVG)18 adopts more turn structures than (VGPVG)18 for most groups at different temperatures. The propensity for β-strands are very close for the two polypeptides at low temperature; however, there is a smaller propensity for β-strands for (VGPVG)18 chain at high temperature than (VPGVG)18, which is probably due to a higher flexibility of prolines and valines in poly(VGPVG). The occurrence of each residue in each type of secondary structural motifs (Turn, β-strands, and Helix) at 290 K and 350 K is presented in Table S3 and plotted in Figure S9.

In the forward sequence, VPGVG, Pro2-Gly3-Val4 segment mostly formed Turns and Gly5Val1 mostly formed β-strands. In reverse VGPVG sequence, residues adopt diverse secondary structures. For example, more than 8% of each residue in Gly2-Pro3-Val4 segment adopt Helix structures, while less than 1.1% of each residue in Pro2-Gly3-Val4 segment adopt Helix structures. This could be due to proline position within the sequence which has asymmetrical structure and restrains the conformational freedom of residues before it, which can break the symmetry of peptide chains. The secondary structures differ significantly due to the proximity of residues to proline, which demonstrate the importance of the sequence direction and the local amino acid order on the formation of the secondary structure of polypeptides. The reversing of amino acid order from poly(VPGVG) to poly(VGPVG) can be described at the molecular level as the mutation of Pro2-Gly3 to Gly2-Pro3. Statistical analysis was conducted on

ACS Paragon Plus Environment

19

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 46

Table 1. Distance between Cα of Val1 and Cα of Val4 in pentamer (VPGVG) and (VGPVG) at 290 K and 350 K, respectively. The numbers listed in the third column are the probability of the distances below 7 Å. Average distance, Å

Probability (distances < 7 Å)

(VPGVG)18, 290 K

8.29±1.24

11.11%

(VPGVG)18, 350 K

8.41±1.49

33.33%

(VGPVG)18, 290 K

7.53±1.22

33.33%

(VGPVG)18, 350 K

8.60±1.33

11.11%

the distances (D1, D2) between the α-carbon atoms in Val1 and Val4 in forward and reverse sequences based on the last 40 ns trajectories. The segment VPGV or VGPV adopt β-turn structures when the distance between the α-carbon atoms in Val1 and Val4 is below 7 Å based on the distance criterion of β-turn structures.45 As listed in Table 1, 11.1% and 33.3% of VPGV segments form β-turns (D1< 7 Å) in forward sequence at 290 K and 350 K, respectively. Interestingly, these probabilities are reversed as 33.3% and 11.1% of VGPV segments form βturns (D2< 7 Å) in the backward sequence at 290 K and 350 K, respectively. The difference is due to the fact that the segment VPGV can intrinsically adopt type-II β-turn structures which increased as the temperature increased from 290 K and 350 K. Therefore, the “mutation” of Pro2-Gly3 to Gly2-Pro3 cause significant differences in the local residue interactions and consequently secondary structure formation, which is consistent with CD measurement (Figure 1).

ACS Paragon Plus Environment

20

Page 21 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Interactions between each amino acid residue with their aqueous environment, determine the conformation of polypeptide and mediate further intermolecular interactions (Figure 4a). In general, the interactions between each amino acid and water molecules decreases as the temperature increased, except for Val4 in the reverse sequence VGPVG. The energy plots clearly confirm that the interactions between each residue type and water in the two sequences are not identical, although Pro2/3 and Val4-Gly5 in these two sequences can be argued to show very similar hydrophobicity. Val1 in the VPGVG interacts with water weaker than Val1 in VGPVG, possibly due to the steric restriction of Val1 in VPGVG by Pro2. The number of hydrogen bonds formed between each residue and water molecules reveals less structural heterogeneity in VGPVG than VPGVG (Figure 4b). The plot of SASA (Figure 4c) for each residue shows very little difference between Pro2, Val4 and Gly5 in VPGVG and Pro3, Val4 and Gly5 in VGPVG. However, the Val1 in VGPVG is more accessible by water molecules than Val1 in VPGVG and Gly2 in VGPVG is less accessible than Gly3 in VPGVG. The compensation between the SASA of these two residues explains the overall SASA trends of both ELP chains.

ACS Paragon Plus Environment

21

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 46

Figure 4. Temperature dependence of properties of a single ELP chain: (a) Interaction energy of each residue in an ELP pentamer with water. (b) Number of hydrogen bonds formed between each residue in an ELP pentamer and water. (c) SASA of each residue in an ELP pentamer. The error bars represent the standard deviation.

Hydration water network We further characterized the structure of hydration water around the polypeptide chains by measuring the size of the largest H-bonded water network around each polypeptide (Figure 5a),

ACS Paragon Plus Environment

22

Page 23 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

which represents how polypeptide influence the water structure and dynamics46, the effect of orientational ordering of water molecules and water-water interactions near the polypeptide surface.21 In the hydration water network, all water molecules are considered to form hydrogen bonds when ROO ≤ 3.5 Å and φO1···O2-H2 ≤ 30° (Figure 5b). The probability distribution of the largest water network around the (VPGVG)18 and (VGPVG)18 at T = 290 K, 310 K, 330 K and 350 K are displayed in Figure 5a. For both polypeptides, the water-water H-bonds in the hydration water shell surrounding each polypeptide generally decreases with increasing temperature, indicating a more disordered water network, especially for (VGPVG)18. We observed a higher thermal stability of the H-bonded water network around (VPGVG)18 than that around (VGPVG)18, whereas the hydration water of (VGPVG)18 has a higher thermal-sensitivity than that of (VPGVG)18. This contrasts with a greater number of hydration water molecules in the water shell of (VGPVG)18 which also has a higher SASA; however, the water molecules form more hydrogen bonds with (VGPVG)18, thereby decreasing the water network in the hydration shell.

ACS Paragon Plus Environment

23

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 46

Figure 5. (a) Probability distribution of the size Nmax of the largest water network around polypeptide. (b) A snapshot of (VPGVG)18 and the largest water network around it at 290 K Polypeptide interactions

ACS Paragon Plus Environment

24

Page 25 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

To examine the effect of temperature on the early stages of polypeptide aggregation, we performed simulations of two temperature-equilibrated polypeptides at 350 K. The simulation

ACS Paragon Plus Environment

25

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

set-up and MD simulation protocols were the same as described in our previous paper.

Page 26 of 46

21

The

distance between the center-of-mass of two chains (Dpp) (Figure 6a) shows that the both

ACS Paragon Plus Environment

26

Page 27 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

polypeptides aggregate during the simulations. The interaction energy (Epp) and number of interchain hydrogen bonds (Npp-double) profiles (Figure 6b and c) between the polypeptides for both

ACS Paragon Plus Environment

27

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 46

cases are correlated with each other and show that the reverse sequence poly(VGPVG) form more compact aggregates with more inter-chain hydrogen bonds and a lower interaction energy

ACS Paragon Plus Environment

28

Page 29 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

between two chains than the forward sequence poly(VPGVG).

Figure 6. Interaction between two ELPs at 350 K. (a) Distance between the center-of-mass of two polypeptides, (b) peptide-peptide interaction energy and (c) the number of inter-chain hydrogen bonds.

ACS Paragon Plus Environment

29

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 46

Discussions The “Pro-Gly” sequence motif has been demonstrated to be important for the elastomeric mechanical behaviors observed within a number of native protein-based polymers.54 It is reasonable to assume that it is also a recurrent if not dominant structural feature associated with the reversibility in the LCST behaviors of ELPs. This is because the type II β-turn per pentamer VPGVG involves proline-glycine pair as the corner residues. The presence of a hydrogen bond between Val1 and Val4 in VPGV is important for stabilization of the type-II β-turn structure formation. The peptide with Pro-Gly to Gly-Pro mutation still exhibits LCST behavior, but loses the propensity for the β-turn formation, which results in a more extended conformation (larger Rg and smaller npp) and a larger solvent exposure (larger SASA) of a peptide. The CD measurements of both soluble ELPs reveal a significant difference in the secondary structure between VPGVG and VGPVG peptides. Poly(VPGVG) has a negative “shelf” at 218 nm, which has historically been attributed to Pro-Gly β-turns in the polymer

38, 39

, and this peak

is absent in the CD spectrum of the reverse sequence, poly(VGPVG). Simulation results showed that, at all temperatures, (VGPVG)18 tends to form more Helix-structures than (VPGVG)18 and at higher temperatures a smaller number of β-strands and β-turns than (VPGVG)18. Such a dependence between amino acid order and secondary structures is consistent with a previous study that examined 1288 PDB structural fragments, and reported that about 21% of the protein sequence segments exhibits different secondary structure propensities in the forward and reverse directions, as assessed by the PREDATOR algorithm. 52, 53 The distinct coacervation process, which is exhibited by ELPs upon heating, is considered to be a complex and multistep transition in thermodynamics point of view.19, 55, 56 It was proposed that the coacervation process can be simplified into intra-molecular interaction and inter-molecular

ACS Paragon Plus Environment

30

Page 31 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

aggregation stages, and intra-molecular interactions during the initial chain folding prevail.

57, 58

Thus, it is reasonable to assume that the dissociation in the cooling process would occur in the reverse order of the two interaction stages; i.e. the order of intrachain interaction and interchain interaction is reversed. Upon cooling, the dissociation of the inter-chain aggregation occurs first. Collectively, the consecutive inter-chain hydrogen bonding arising in the poly(VGPVG) appeared to perturb the dissociation by the additional interchain interactions. In other words, each poly(VPGVG) aggregate directly and then quickly dissociates into individual chain when the solution is cooled to below the LCST, which is different from the dissociation mechanism of poly(VGPVG) aggregates where some aggregates are not able to completely dissociate upon cooling due to formation of additional strong inter-chain hydrogen bonds. Previously, the consistent but very small degree of hysteresis in the heating and cooling cycle was also observed in the PNIPAM water solution. 47-49 Based on FTIR data, it has been proposed that some of inter-chain hydrogen bonds could persist upon cooling, especially when the temperature is not far away from the phase transition temperature.50 Furthermore, a comparison between PNIPAM and PDEAM indicates that the hysteresis originates from additional hydrogen bonds formed in their collapsed state at temperatures higher than the LCST.51 We propose that the hysteretic thermal behavior observed for poly(VGPVG) (Figure 1) is because at high temperatures, a single poly(VGPVG) chain is in a loose conformation (Figure 2) that permits strong inter-molecular interactions within the aggregate (Figure 6) with formation of more interchain hydrogen bonds as compared to poly(VPGVG). These inter-chain hydrogen bonds can stabilize the aggregated phase by providing stable inter-chain contacts that are not easily disrupted by cooling. Moreover, more extended conformation and larger surface area in (VGPVG)18 increase the chances for the formation of intermolecular hydrogen bonds.

ACS Paragon Plus Environment

31

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 46

Conclusions This paper presents results from a combined study using atomistic MD simulations and experimental observations (turbidity and CD measurements) on the effect of amino acid order on the structural and dynamics changes of ELP chains as a function of temperature. Our experimental observations show different thermal hysteresis behavior of the forward sequence, poly(VPGVG), and its reverse sequence, poly(VGPVG); whereas poly(VPGVG)s showed no thermal hysteresis in its LCST phase behavior, while the reverse sequence poly(VGPVG)s exhibits significant thermal hysteresis as the polypeptides remain aggregated despite significant undercooling below their Tt. The irreversible aggregation observed for poly(VGPVG)s is independent of concentration. Their behavior is also independent of both heating and cooling rates with nearly identical Tt and complete irreversibility observed for heating and cooling rates ranging from 0.5 – 5

C/min.

Our MD simulations show that both (VPGVG)18 and (VGPVG)18 show an increase in hydrophobicity with increased temperature, which is the common characteristic feature of neutral thermo-responsive polymers that exhibit LCST phase behavior. However, (VGPVG)18 is more hydrophilic than (VGPVG)18 at high temperatures despite the fact that both sequences have the same amino acid composition and molecular weight. (VPGVG)18 also assumes a more compact structure as indicated by a smaller Rg at high temperature then (VGPVG)18. The examination of the association kinetics revealed additional details about interchain association. The aggregation of two poly(VGPVG) molecules above the LCST can be regarded as the consequence of the competition between hydrophilic peptide–water and hydrophobic peptide-peptide interactions. We found that a pair of (VGPVG)18 peptides show stronger aggregation with a larger number of inter-chain hydrogen bonds as compared to (VPGVG)18,

ACS Paragon Plus Environment

32

Page 33 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

which is probably due to a more extended conformation and larger solvent exposure of a (VGPVG)18 peptide. Thus, thermal hysteresis exhibited by poly(VGPVG) may originate from formation of additional inter-chain hydrogen bonds that stabilize the aggregated phase by providing stable inter-chain contacts that are not easily disrupted by cooling. In summary, the conformational preferences of amino acids, which are impacted by sequence reversal, control the thermal hysteresis of ELPs. The change in rotational degrees of freedom within the sequence in (VPGVG)18 and (VGPVG)18 leads to the differences in the overall compaction as a function of temperature, which in turn control intermolecular interactions and hydrogen bonds formation within the aggregate, which consequently results in a thermal hysteresis. Supporting Information Available: Additional information on the simulation systems, convergence of simulations and analysis. This material is available free of charge via the Internet at http://pubs.acs.org. AUTHOR INFORMATION Corresponding Author *[email protected] Funding Sources This work was supported by the NSF's Research Triangle MRSEC (DMR-1121107) and by a grant from the NIH TO A.C. (GM-61232) ACKNOWLEDGMENT

ACS Paragon Plus Environment

33

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 46

This work was supported by the NSF's Research Triangle MRSEC (DMR-1121107) and by the NIH (GM-61232). Computer support was provided by the High Performance Computing Center at North Carolina State University.

REFERENCES 1. van Eldijk, M. B.; McGann, C. L.; Kiick, K. L.; van Hest, J. C. M., Elastomeric Polypeptides. Top Curr Chem 2012, 310, 71-116. 2. Amiram, M.; Quiroz, F. G.; Callahan, D. J.; Chilkoti, A., A highly parallel method for synthesizing DNA repeats enables the discovery of 'smart' protein polymers. Nat Mater 2011, 10, (2), 141-148. 3. Tamburro, A. M., A never-ending love story with elastin: a scientific autobiography. Nanomedicine-Uk 2009, 4, (4), 469-487. 4. Chilkoti, A.; Dreher, M. R.; Meyer, D. E., Design of thermally responsive, recombinant polypeptide carriers for targeted drug delivery. Adv Drug Deliver Rev 2002, 54, (8), 1093-1111. 5. McDaniel, J. R.; Bhattacharyya, J.; Vargo, K. B.; Hassouneh, W.; Hammer, D. A.; Chilkoti, A., Self-Assembly of Thermally Responsive Nanoparticles of a Genetically Encoded Peptide Polymer by Drug Conjugation. Angewandte Chemie-International Edition 2013, 52, (6), 1683-1687. 6. Dreher, M. R.; Raucher, D.; Balu, N.; Colvin, O. M.; Ludeman, S. M.; Chilkoti, A., Evaluation of an elastin-like polypeptide-doxorubicin conjugate for cancer therapy. J Control Release 2003, 91, (1-2), 31-43. 7. Meyer, D. E.; Chilkoti, A., Purification of recombinant proteins by fusion with thermallyresponsive polypeptides. Nat Biotechnol 1999, 17, (11), 1112-1115. 8. Ge, X.; Yang, D. S. C.; Trabbic-Carlson, K.; Kim, B.; Chilkoti, A.; Filipe, C. D. M., Selfcleavable stimulus responsive tags for protein purification without chromatography. Journal of the American Chemical Society 2005, 127, (32), 11228-11229. 9. Trabbic-Carlson, K.; Meyer, D. E.; Liu, L.; Piervincenzi, R.; Nath, N.; LaBean, T.; Chilkoti, A., Effect of protein fusion on the transition temperature of an environmentally responsive elastin-like polypeptide: a role for surface hydrophobicity? Protein Eng Des Sel 2004, 17, (1), 57-66. 10. Ziebarth, J.; Wang, Y., Coarse-Grained Molecular Dynamics Simulations of DNA Condensation by Block Copolymer and Formation of Core−Corona Structures. J. Phys. Chem. B 2010, 114, (19), 6225-6232. 11. Lim, D. W.; Nettles, D. L.; Setton, L. A.; Chilkoti, A., Rapid cross-linking of elastin-like polypeptides with (hydroxymethyl)phosphines in aqueous solution. Biomacromolecules 2007, 8, (5), 1463-1470. 12. McHale, M. K.; Setton, L. A.; Chilkoti, A., Synthesis and in vitro evaluation of enzymatically cross-linked elastin-like polypeptide gels for cartilaginous tissue repair. Tissue Eng 2005, 11, (11-12), 1768-1779. 13. Vrhovski, B.; Jensen, S.; Weiss, A. S., Coacervation characteristics of recombinant human tropoelastin. Eur J Biochem 1997, 250, (1), 92-98.

ACS Paragon Plus Environment

34

Page 35 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

14. Urry, D. W.; Luan, C. H.; Parker, T. M.; Gowda, D. C.; Prasad, K. U.; Reid, M. C.; Safavy, A., Temperature of Polypeptide Inverse Temperature Transition Depends on Mean Residue Hydrophobicity. Journal of the American Chemical Society 1991, 113, (11), 4346-4348. 15. Luan, C. H.; Parker, T. M.; Prasad, K. U.; Urry, D. W., Differential Scanning Calorimetry Studies of Nacl Effect on the Inverse Temperature Transition of Some Elastin-Based Polytetrapeptides, Polypentapeptides, and Polynonapeptides. Biopolymers 1991, 31, (5), 465475. 16. Nuhn, H.; Klok, H. A., Secondary Structure Formation and LCST Behavior of Short Elastin-Like Peptides. Biomacromolecules 2008, 9, (10), 2755-2763. 17. Urry, D. W., Physical chemistry of biological free energy transduction as demonstrated by elastic protein-based polymers. J Phys Chem B 1997, 101, (51), 11007-11028. 18. Rincón, A. C.; Molina-Martinez, I. T.; de las Heras, B.; Alonso, M.; Baílez, C.; Rodríguez-Cabello, J. C.; Herrero-Vanrell, R., Biocompatibility of elastin-like polymer poly(VPAVG) microparticles: in vitro and in vivo studies. Journal of Biomedical Materials Research Part A 2006, 78A, (2), 343-351. 19. Reguera, J.; Lagaron, J. M.; Alonso, M.; Reboto, V.; Calvo, B.; Rodriguez-Cabello, J. C., Thermal behavior and kinetic analysis of the chain unfolding and refolding and of the concomitant nonpolar solvation and desolvation of two elastin-like polymers. Macromolecules 2003, 36, (22), 8470-8476. 20. Quiroz, F. G.; Chilkoti, A., Sequence heuristics to encode phase behaviour in intrinsically disordered protein polymers. Nat Mater 2015, 14, (11), 1164-1171. 21. Li, N. K.; Quiroz, F. G.; Hall, C. K.; Chilkoti, A.; Yingling, Y. G., Molecular Description of the LCST Behavior of an Elastin-Like Polypeptide. Biomacromolecules 2014, 15, (10), 35223530. 22. Case, D. A.; Darden, T. A.; T.E. Cheatham, I.; Simmerling, C. L.; Wang, J.; Duke, R. E.; Luo, R.; Walker, R. C.; Zhang, W.; Merz, K. M.; Roberts, B.; Hayik, S.; Roitberg, A.; Seabra, G.; Swails, J.; Goetz, A. W.; Kolossváry, I.; Wong, K. F.; Paesani, F.; Vanicek, J.; Wolf, R. M.; Liu, J.; Wu, X.; Brozell, S. R.; Steinbrecher, T.; Gohlke, H.; Cai, Q.; Ye, X.; Wang, J.; Hsieh, M.-J.; Cui, G.; Roe, D. R.; Mathews, D. H.; Seetin, M. G.; Salomon-Ferrer, R.; Sagui, C.; Babin, V.; Luchko, T.; Gusarov, S.; Kovalenko, A.; Kollman, P. A. AMBER 11, University of California: San Francisco, 2011. 23. Case, D. A.; Darden, T. A.; T.E. Cheatham, I.; Simmerling, C. L.; Wang, J.; Duke, R. E.; Luo, R.; Walker, R. C.; Zhang, W.; Merz, K. M.; Roberts, B.; Hayik, S.; Roitberg, A.; Seabra, G.; Swails, J.; Goetz, A. W.; Kolossváry, I.; Wong, K. F.; Paesani, F.; Vanicek, J.; Wolf, R. M.; Liu, J.; Wu, X.; Brozell, S. R.; Steinbrecher, T.; Gohlke, H.; Cai, Q.; Ye, X.; Wang, J.; Hsieh, M.-J.; Cui, G.; Roe, D. R.; Mathews, D. H.; Seetin, M. G.; Salomon-Ferrer, R.; Sagui, C.; Babin, V.; Luchko, T.; Gusarov, S.; Kovalenko, A.; Kollman, P. A. AMBER 12, University of California: San Francisco, 2012. 24. Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L., Comparison of Simple Potential Functions for Simulating Liquid Water. J Chem Phys 1983, 79, (2), 926-935. 25. Zhao, B.; Li, N. K.; Yingling, Y. G.; Hall, C. K., LCST Behavior is Manifested in a Single Molecule: Elastin-Like polypeptide (VPGVG)n. Biomacromolecules 2016, 17, (1), 111118. 26. Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G., A Smooth Particle Mesh Ewald Method. J Chem Phys 1995, 103, (19), 8577-8593.

ACS Paragon Plus Environment

35

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 46

27. Urry, D. W., Entropic Elastic Processes in Protein Mechanisms .1. Elastic Structure Due to an Inverse Temperature Transition and Elasticity Due to Internal Chain Dynamics. J Protein Chem 1988, 7, (1), 1-34. 28. Krukau, A.; Brovchenko, I.; Geiger, A., Temperature-induced conformational transition of a model elastin-like peptide GVG(VPGVG)(3) in water. Biomacromolecules 2007, 8, (7), 2196-2202. 29. Li, B.; Alonso, D. O. V.; Daggett, V., The molecular basis for the inverse temperature transition of elastin. J Mol Biol 2001, 305, (3), 581-592. 30. Berendsen, H. J. C.; Postma, J. P. M.; Vangunsteren, W. F.; Dinola, A.; Haak, J. R., Molecular-Dynamics with Coupling to an External Bath. J Chem Phys 1984, 81, (8), 3684-3690. 31. Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H. J. C., Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. Journal of Computational Physics 1977, 23, (3), 327-341. 32. Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kale, L.; Schulten, K., Scalable molecular dynamics with NAMD. J Comput Chem 2005, 26, (16), 1781-1802. 33. Shao, J. Y.; Tanner, S. W.; Thompson, N.; Cheatham, T. E., Clustering molecular dynamics trajectories: 1. Characterizing the performance of different clustering algorithms. J Chem Theory Comput 2007, 3, (6), 2312-2334. 34. McDaniel, J. R.; Mackay, J. A.; Quiroz, F. G.; Chilkoti, A., Recursive directional ligation by plasmid reconstruction allows rapid and seamless cloning of oligomeric genes. Biomacromolecules 2010, 11, (4), 944-52. 35. Quiroz, F. G.; Chilkoti, A., Sequence heuristics to encode phase behaviour in intrinsically disordered protein polymers. Nat Mater 2015, 14, (11), 1164-71. 36. McDaniel, J. R.; Callahan, D. J.; Chilkoti, A., Drug delivery to solid tumors by elastinlike polypeptides. Adv Drug Deliv Rev 2010, 62, (15), 1456-67. 37. Cho, Y.; Zhang, Y.; Christensen, T.; Sagle, L. B.; Chilkoti, A.; Cremer, P. S., Effects of Hofmeister anions on the phase transition temperature of elastin-like polypeptides. J Phys Chem B 2008, 112, (44), 13765-71. 38. Reiersen, H.; Clarke, A. R.; Rees, A. R., Short elastin-like peptides exhibit the same temperature-induced structural transitions as elastin polymers: implications for protein engineering. J Mol Biol 1998, 283, (1), 255-64. 39. Bochicchio, B.; Pepe, A.; Tamburro, A. M., Investigating by CD the molecular mechanism of elasticity of elastomeric proteins. Chirality 2008, 20, (9), 985-94. 40. Gross, P. C.; Possart, W.; Zeppezauer, M., An alternative structure model for the polypentapeptide in elastin. Z Naturforsch C 2003, 58, (11-12), 873-878. 41. Yao, X. L.; Hong, M., Structure distribution in an elastin-mimetic peptide (VPGVG)(3) investigated by solid-state NMR. Journal of the American Chemical Society 2004, 126, (13), 4199-4210. 42. Ohgo, K.; Ashida, J.; Kumashiro, K. K.; Asakura, T., Structural determination of an elastin-mimetic model peptide, (Val-Pro-Gly-Val-Gly)(6), studied by C-13 CP/MAS NMR chemical shifts, two-dimensional off magic angle spinning spin-diffusion NMR, rotational echo double resonance, and statistical distribution of torsion angles from protein data bank. Macromolecules 2005, 38, (14), 6038-6047.

ACS Paragon Plus Environment

36

Page 37 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

43. Lovell, S. C.; Davis, I. W.; Adrendall, W. B.; de Bakker, P. I. W.; Word, J. M.; Prisant, M. G.; Richardson, J. S.; Richardson, D. C., Structure validation by C alpha geometry: phi,psi and C beta deviation. Proteins 2003, 50, (3), 437-450. 44. Kabsch, W.; Sander, C., Dictionary of Protein Secondary Structure - Pattern-Recognition of Hydrogen-Bonded and Geometrical Features. Biopolymers 1983, 22, (12), 2577-2637. 45. Bornot, A.; de Brevern, A. G., Protein beta-turn assignments. Bioinformation 2006, 1, (5), 153-155. 46. Abseher, R.; Schreiber, H.; Steinhauser, O., The influence of a protein on water dynamics in its vicinity investigated by molecular dynamics simulation. Proteins 1996, 25, (3), 366-378. 47. Ding, Y. W.; Ye, X. D.; Zhang, G. Z., Microcalorimetric investigation on aggregation and dissolution of poly(N-isopropylacrylamide) chains in water. Macromolecules 2005, 38, (3), 904-908. 48. Maeda, Y.; Yamamoto, H.; I, I., Effects of ionization on the phase behavior of poly(Nisopropylacrylamide-co-acrylic acid) and poly(N,N-diethylacrylamide-co-acrylic acid) in water. Colloid Polym Sci 2004, 282, (11), 1268-1273. 49. Maeda, Y.; Nakamura, T.; Ikeda, I., Change in solvation of poly(N,N-diethylacrylamide) during phase transition in aqueous solutions as observed by IR spectroscopy. Macromolecules 2002, 35, (27), 10172-10177. 50. Cheng, H.; Shen, L.; Wu, C., LLS and FTIR studies on the hysteresis in association and dissociation of poly(N-isopropylacrylamide) chains in water. Macromolecules 2006, 39, (6), 2325-2329. 51. Lu, Y. J.; Zhou, K. J.; Ding, Y. W.; Zhang, G. Z.; Wu, C., Origin of hysteresis observed in association and dissociation of polymer chains in water. Physical Chemistry Chemical Physics 2010, 12, (13), 3188-3194. 52. Frishman, D.; Argos, P., Seventy-five percent accuracy in protein secondary structure prediction. Proteins 1997, 27, (3), 329-35. 53. Park, J.; Dietmann, S.; Heger, A.; Holm, L., Estimating the significance of sequence order in protein secondary structure and prediction. Bioinformatics 2000, 16, (11), 978-987. 54. Kim, W.; Conticello, V. P., Protein Engineering Methods for Investigation of StructureFunction Relationships in Protein-Based Elastomeric Materials. Polymer Reviews 2007, 47, (1), 93-119. 55. Kaibara, K.; Akinari, Y.; Okamoto, K.; Uemura, Y.; Yamamoto, S.; Kodama, H.; Kondo, M., Characteristic interaction of Ca2+ ions with elastin coacervate: Ion transport study across coacervate layers of alpha-elastin and elastin model polypeptide, (Val-Pro-Gly-Val-Gly)(n). Biopolymers 1996, 39, (2), 189-198. 56. Reguera, J.; Urry, D. W.; Parker, T. M.; McPherson, D. T.; Rodriguez-Cabello, J. C., Effect of NaCl on the exothermic and endothermic components of the inverse temperature transition of a model elastin-like polymer. Biomacromolecules 2007, 8, (2), 354-8. 57. Ribeiro, A.; Arias, F. J.; Reguera, J.; Alonso, M.; Rodriguez-Cabello, J. C., Influence of the Amino-Acid Sequence on the Inverse Temperature Transition of Elastin-Like Polymers. Biophys. J. 2009, 97, (1), 312-320. 58. Teeuwen, R. L. M.; de Wolf, F. A.; Zuilhof, H.; van Hest, J. C. M., Elastin-like polypeptides of different molecular weights show independent transition temperatures when mixed. Soft Matter 2009, 5, (21), 4305-4310.

ACS Paragon Plus Environment

37

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 46

ACS Paragon Plus Environment

38

Page 39 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Table of Contents Graphic

ACS Paragon Plus Environment

39

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1: Experimental comparison of LCST phase behavior of poly(VPGVG) and poly(VGPVG). (a,b) Temperature dependent turbidimetry for (VPGVG)40/80 and (VGPVG)40/80 (50 µM, PBS). Polymers with identical composition but sequence-reversed orientation (reverse sequence) have different LCST phase behaviors. (VPGVG)40 and (VPGVG)80 —canonical ELP sequences—shows reversible phase behavior, but the reverse sequence, poly(VGPVG), shows irreversible aggregation as a function of solution temperature. (c) The reverse sequence VGPVG has a slightly higher transition temperature across all measured concentrations despite the identical amino acid composition. (d) CD measurements of the soluble polymers (10 µM at 15⁰C in H2O) reveal a significant difference in their secondary structure. Poly(VPGVG) has a negative “shelf” at 218 nm, indicating the presence of beta turns 35, 36, and this peak is absent in the CD spectrum of the reverse sequence, poly(VGPVG). Both polymers otherwise appear highly disordered. 180x144mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 40 of 46

Page 41 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Figure 2. Temperature dependence of the properties of a single (VPGVG)18 and (VGPVG)18 ELP chain: (a) radius of gyration (Rg), (b) npw, the number of peptide-water hydrogen bonds, (c) solvent accessible surface area (SASA, (d) peptide-water interaction energy. Error bars represent the standard deviation. 159x128mm (300 x 300 DPI)

ACS Paragon Plus Environment

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. (a-d) Ramachandran plot (ϕ, ψ distributions) for Val1 and Val4 in (VPGVG)18 and (VGPVG)18, respectively, at 350 K, colored by intensities. The most populated regions on the Ramachandran plots are labelled. (e-f) Representative structures for (e) G(VPGVG)V segment and (f) G(VGPVG)V segment at 350 K built with the dihedral angles obtained from the highest intensity points in Ramachandran plots (Figure 3 and Figure S8). 142x95mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 42 of 46

Page 43 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Figure 4. Temperature dependence of properties of a single ELP chain: (a) Interaction energy of each residue in an ELP pentamer with water. (b) Number of hydrogen bonds formed between each residue in an ELP pentamer and water. (c) SASA of each residue in an ELP pentamer. The error bars represent the standard deviation. 94x241mm (300 x 300 DPI)

ACS Paragon Plus Environment

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. (a) Probability distribution of the size Nmax of the largest water network around polypeptide. (b) A snapshot of (VPGVG)18 and the largest water network around it at 290 K 91x133mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 44 of 46

Page 45 of 46 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biomacromolecules

Figure 6. Interaction between two ELPs at 350 K. (a) Distance between the center-of-mass of two polypeptides, (b) peptide-peptide interaction energy and (c) the number of inter-chain hydrogen bonds. 98x232mm (300 x 300 DPI)

ACS Paragon Plus Environment

Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

88x35mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 46 of 46