Molecular-dynamics simulations of C- and N-terminal peptide derivatives of GCN4-p1 in aqueous solution

The native GCN4 leucine-zipper domain does not uniquely specify a dimeric oligomerization state

The dimerization domain of the yeast transcription factor GCN4, one of the first coiled-coil proteins to be structurally characterized at high resolution, has served as the basis for numerous fundamental studies on α-helical folding. Mutations in the GCN4 leucine zipper are known to change its preferred oligomerization state from dimeric to trimeric or tetrameric; however, the wild-type sequence has been assumed to encode a two-chain assembly exclusively. Here we demonstrate that the GCN4 coiled-coil domain can populate either a dimer or trimer fold, depending on environment. We report high-resolution crystal structures of the wild-type sequence in dimeric and trimeric assemblies. Biophysical measurements suggest populations of both oligomerization states under certain experimental conditions in solution. We use parallel tempering molecular dynamics simulations on the microsecond time scale to compare the stability of the dimer and trimer folded states in isolation. In total, our results suggest that the folding behavior of the well-studied GCN4 leucine-zipper domain is more complex than was previously appreciated. Our results have implications in ongoing efforts to establish predictive algorithms for coiled-coil folds and the selection of coiled-coil model systems for design and mutational studies where oligomerization state specificity is an important consideration.

Simulations of temperature and salt concentration effects on bZIP, a basic region leucine zipper

Basic region leucine zipper (bZIP) transcription factors are dimeric proteins that recognize DNA. The monomers consist of a leucine zipper subdomain responsible for dimerization and a highly basic DNA recognition subdomain. Twelve explicit solvent molecular dynamics (MD) trajectories were run on the GCN4 bZIP transcriptional factor in the absence of DNA at three temperatures and two ion concentrations (0 mM with Cl(-) ions to neutralize the bZIP and 200 mM with additional Na(+) and Cl(-) ions) to probe the conformational ensemble that the basic region samples. In most trajectories, the basic region exhibits an alligator-jaw-like opening and closing (only one monomer moves), versus scissor-like motion, by a mainly rigid body, hinge motion centered on three "fork" residues that span the basic region to the coiled coil. In this motion, the α-helical character of the basic region monomers is mostly maintained. A broad range of distances is accessed, consistent with the absence of particular interactions for the basic region monomers. In two of the trajectories, the basic region monomers "collapse" to form a stable state. The coiled coil, leucine zipper subdomain is very stable for all of the trajectories. Ion solvation of the charged residue side chains is transient, on the scale of a few picoseconds. There is no evidence for persistent specific ion salt bridges to charged residues. For 0 mM, only certain basic region positively charged residues are substantially Cl(-) ion salt bridged. For 200 mM, in addition, some basic region negatively (positively) charged residues are salt bridged to Na(+) (Cl(-)) ions. The different ion solvation patterns at the two ion concentrations are not greatly temperature sensitive, and the conformational sampling found in the MD is remarkably unperturbed by ion concentration and/or temperature.

Methods of NMR structure refinement: molecular dynamics simulations improve the agreement with measured NMR data of a C-terminal peptide of GCN4-p1

The C-terminal trigger sequence is essential in the coiled-coil formation of GCN4-p1; its conformational properties are thus of importance for understanding this process at the atomic level. A solution NMR model structure of a peptide, GCN4p16-31, encompassing the GCN4-p1 trigger sequence was proposed a few years ago. Derived using a standard single-structure refinement protocol based on 172 nuclear Overhauser effect (NOE) distance restraints, 14 hydrogen-bond and 11 phi torsional-angle restraints, the resulting set of 20 NMR model structures exhibits regular alpha-helical structure. However, the set slightly violates some measured NOE bounds and does not reproduce all 15 measured (3)J(H(N)-H(Calpha))-coupling constants, indicating that different conformers of GCN4p16-31 might be present in solution. With the aim to resolve structures compatible with all NOE upper distance bounds and (3)J-coupling constants, we executed several structure refinement protocols employing unrestrained and restrained molecular dynamics (MD) simulations with two force fields. We find that only configurational ensembles obtained by applying simultaneously time-averaged NOE distance and (3)J-coupling constant restraining with either force field reproduce all the experimental data. Additionally, analyses of the simulated ensembles show that the conformational variability of GCN4p16-31 in solution admitted by the available set of 187 measured NMR data is larger than represented by the set of the NMR model structures. The conformations of GCN4p16-31 in solution differ in the orientation not only of the side-chains but also of the backbone. The inconsistencies between the NMR model structures and the measured NMR data are due to the neglect of averaging effects and the inclusion of hydrogen-bond and torsional-angle restraints that have little basis in the primary, i.e. measured NMR data.

Computational Studies of a Protein-based Nanoactuator for Nanogripping Applications

The design hypothesis, architectures, and computational modeling of a novel peptide-based nanoactuator are presented in this paper. We engineered the α-helical coiled-coil portion of the yeast transcriptional activator peptide called GCN4 to obtain an environmentally responsive nanoactuator. The dimeric coiled-coil peptide consists of two identical approximately 4.5 nm long and approximately 3 nm wide polypeptide chains. The actuation mechanism depends on the modification of electrostatic charges along the peptide by varying the pH of the solution resulting in the reversible movement of helices and, therefore, creating the motion of an actuator. Using molecular dynamics simulations we showed that pH changes led to a reversible opening of up to 1.5 nm which is approximately 150% of the initial separation of the nanoactuator. We also investigated the forces generated by the nanoactuator upon pH actuation, using a new method based on a modified steered molecular dynamics technique. Owing to its open and close motion resembling that of tweezers, the new nanoactuator can potentially be used as a nanogripper in various nanomanipulation tasks such as detection and removal of heavy metal ions during nanofabrication processes or as a molecular switch.

Molecular dynamics guided study of salt bridge length dependence in both fluorinated and non-fluorinated parallel dimeric coiled-coils

The alpha-helical coiled-coil is one of the most common oligomerization motifs found in both native and engineered proteins. To better understand the stability and dynamics of the coiled-coil motifs, including those modified by fluorination, several fluorinated and nonfluorinated parallel dimeric coiled-coil protein structures were designed and modeled. We also attempt to investigate how changing the length and geometry of the important stabilizing salt bridges influences the coiled-coil protein structure. Molecular dynamics (MD) and free energy simulations with AMBER used a particle mesh Ewald treatment of the electrostatics in explicit TIP3P solvent with balanced force field treatments. Preliminary studies with legacy force fields (ff94, ff96, and ff99) show a profound instability of the coiled-coil structures in short MD simulation. Significantly, better behavior is evident with the more balanced ff99SB and ff03 protein force fields. Overall, the results suggest that the coiled-coil structures can readily accommodate the larger acidic arginine or S-2,7-diaminoheptanedoic acid mutants in the salt bridge, whereas substitution of the smaller L-ornithine residue leads to rapid disruption of the coiled-coil structure on the MD simulation time scale. This structural distortion of the secondary structure allows both the formation of large hydration pockets proximal to the charged groups and within the hydrophobic core. Moreover, the increased structural fluctuations and movement lead to a decrease in the water occupancy lifetimes in the hydration pockets. In contrast, analysis of the hydration in the stable dimeric coiled-coils shows high occupancy water sites along the backbone residues with no water occupancy in the hydrophobic core, although transitory water interactions with the salt bridge residues are evident. The simulations of the fluorinated coiled-coils suggest that in some cases fluorination electrostatically stabilizes the intermolecular coiled-coil salt bridges. Structural analyses also reveal different side chain rotamer preferences for leucine when compared with 5,5,5,5',5',5'-hexafluoroleucine mutants. These observed differences in the side chain rotamer populations suggest differential changes in the side chain conformational entropy upon coiled-coil formation when the protein is fluorinated. The free energy of hydration of the isolated 5,5,5,5',5',5'-hexafluoroleucine amino acid is calculated to be 1.1 kcal/mol less stable than leucine; this hydrophobic penalty in the monomer may provide a driving force for coiled-coil dimer formation. Estimation of the ellipticity at 222 nm from a series of snapshots from the MD simulations with DicroCalc shows distinct increases in the ellipticity when the coiled-coil is fluorinated, which suggests that the helicity in the folded coiled-coils is greater when fluorinated.

Reversible pH-controlled DNA-binding peptide nanotweezers: an in-silico study

We describe the molecular dynamics (MD)-aided engineering design of mutant peptides based on the α-helical coiled-coil GCN4 leucine zipper peptide (GCN4-p1) in order to obtain environmentally-responsive nanotweezers. The actuation mechanism of the nanotweezers depends on the modification of electrostatic charges on the residues along the length of the coiled coil. Modulating the solution pH between neutral and acidic values results in the reversible movement of helices toward and away from each other and creates a complete closed-open-closed transition cycle between the helices. Our results indicate that the mutants show a reversible opening of up to 15 Å (1.5 nm; approximately 150% of the initial separation) upon pH actuation. Investigation on the physicochemical phenomena that influence conformational properties, structural stability, and reversibility of the coiled-coil peptide-based nanotweezers revealed that a rationale- and design-based approach is needed to engineer stable peptide or macromolecules into stimuli-responsive devices. The efficacy of the mutant that demonstrated the most significant reversible actuation for environmentally responsive modulation of DNA-binding activity was also demonstrated. Our results have significant implications in bioseparations and in the engineering of novel transcription factors.

Configurational entropy elucidates the role of salt-bridge networks in protein thermostability

Detailed knowledge of how networks of surface salt bridges contribute to protein thermal stability is essential not only to understand protein structure and function but also to design thermostable proteins for industrial applications. Experimental studies investigating thermodynamic stability through measurements of free energy associated with mutational alterations in proteins provide only macroscopic evidence regarding the structure of salt-bridge networks and assessment of their contribution to protein stability. Using explicit-solvent molecular dynamics simulations to provide insight on the atomic scale, we investigate here the structural stability, defined in terms of root-mean-square fluctuations, of a short polypeptide designed to fold into a stable trimeric coiled coil with a well-packed hydrophobic core and an optimal number of intra- and interhelical surface salt bridges. We find that the increase of configurational entropy of the backbone and side-chain atoms and decreased pair correlations of these with increased temperature are consistent with nearly constant atom-positional root-mean-square fluctuations, increased salt-bridge occupancies, and stronger electrostatic interactions in the coiled coil. Thus, our study of the coiled coil suggests a mechanism in which well-designed salt-bridge networks could accommodate stochastically the disorder of increased thermal motion to produce thermostability.

Molecular basis of coiled-coil formation

Coiled coils have attracted considerable interest as design templates in a wide range of applications. Successful coiled-coil design strategies therefore require a detailed understanding of coiled-coil folding. One common feature shared by coiled coils is the presence of a short autonomous helical folding unit, termed "trigger sequence," that is indispensable for folding. Detailed knowledge of trigger sequences at the molecular level is thus key to a general understanding of coiled-coil formation. Using a multidisciplinary approach, we identify and characterize here the molecular determinants that specify the helical conformation of the monomeric early folding intermediate of the GCN4 coiled coil. We demonstrate that a network of hydrogen-bonding and electrostatic interactions stabilize the trigger-sequence helix. This network is rearranged in the final dimeric coiled-coil structure, and its destabilization significantly slows down GCN4 leucine zipper folding. Our findings provide a general explanation for the molecular mechanism of coiled-coil formation.