Calculation of protein ionization equilibria with conformational sampling: pK(a) of a model leucine zipper, GCN4 and barnase.

The use of conformational ensembles provided by nuclear magnetic resonance (NMR) experiments or generated by molecular dynamics (MD) simulations has been regarded as a useful approach to account for protein motions in the context of pK(a) calculations, yet the idea has been tested occasionally. This is the first report of systematic comparison of pK(a) estimates computed from long multiple MD simulations and NMR ensembles. As model systems, a synthetic leucine zipper, the naturally occurring coiled coil GCN4, and barnase were used. A variety of conformational averaging and titration curve-averaging techniques, or combination thereof, was adopted and/or modified to investigate the effect of extensive global conformational sampling on the accuracy of pK(a) calculations. Clustering of coordinates is proposed as an approach to reduce the vast diversity of MD ensembles to a few structures representative of the average electrostatic properties of the system in solution. Remarkable improvement of the accuracy of pK(a) predictions was achieved by the use of multiple MD simulations. By using multiple trajectories the absolute error in pK(a) predictions for the model leucine zipper was reduced to as low as approximately 0.25 pK(a) units. The validity, advantages, and limitations of explicit conformational sampling by MD, compared with the use of an average structure and a high internal protein dielectric value as means to improve the accuracy of pK(a) calculations, are discussed.     

Structure-based design of a leucine zipper protein with new DNA contacting region

We have employed a structure-based design to construct a small folding domain from the F-actin bundling protein villin that contains the amino acids necessary for the DNA binding of the basic leucine zipper protein GCN4 and have compared its DNA binding with GCN4. The monomeric motif folds into a stable domain and binds DNA in a rigid-body mechanism, while its affinity is not higher than that of the basic region peptide. The addition of the leucine zipper region to the folded domain restored its sequence-specific DNA binding comparable to that of GCN4. Unlike the monomeric folded domain, its leucine zipper derivative undergoes a conformational change upon DNA binding. CD spectral and thermodynamic studies indicate that the DNA-contacting region is folded in the presence or absence of DNA and suggest that the junction between the DNA-contacting and the leucine zipper regions transits to a helix in the presence of DNA. These results demonstrate that the structural transition outside the direct-contacting region, which adjusts the precise location of the DNA-contacting region, plays a critical role in the specific complex formation of basic leucine zipper proteins.     

A pH-dependent conformational ensemble mediates proton transport through the influenza A/M2 protein

The influenza A/M2 protein exhibits inwardly rectifying, pH-activated proton transport that saturates at low pH. A comparison of high-resolution structures of the transmembrane domain at high and low pH suggests that pH-dependent conformational changes may facilitate proton conduction by alternately changing the accessibility of the N-terminal and C-terminal regions of the channel as a proton transits through the transmembrane domain. Here, we show that M2 functionally reconstituted in liposomes populates at least three different conformational states over a physiologically relevant pH range, with transition midpoints that are consistent with previously reported His37 pK(a) values. We then develop and test two similar, quantitative mechanistic models of proton transport, where protonation shifts the equilibrium between structural states having different proton affinities and solvent accessibilities. The models account well for a collection of experimental data sets over a wide range of pH values and voltages and require only a small number of adjustable parameters to accurately describe the data. While the kinetic models do not require any specific conformation for the protein, they nevertheless are consistent with a large body of structural information based on high-resolution nuclear magnetic resonance and crystallographic structures, optical spectroscopy, and molecular dynamics calculations.     

An engineered leucine zipper a position mutant with an unusual three-state unfolding pathway

The leucine zipper is a dimeric coiled-coil structural motif consisting of four to six heptad repeats, designated (abcdefg)(n). In the GCN4 leucine zipper, a position 16 in the third heptad is occupied by an Asn residue whereas the other a positions are Val residues. Recently, we have constructed variants of the GCN4 leucine zipper in which the a position Val residues were replaced by Ile. The folding and unfolding of the wild-type GCN4 leucine zipper and the Val to Ile variant both adhere to a simple two-state mechanism. In this study, another variant of the GCN4 leucine zipper was constructed by moving the single Asn residue from a position 16 to a position 9. This switch causes the thermal unfolding of the GCN4 leucine zipper to become three state. The unfolding pathway of this variant was determined by thermal denaturation, limited proteinase K digestion, and sedimentation equilibrium analysis. Our data are consistent with a model in which the variant first unfolds from its N terminus and changes the oligomerization specificity from a native dimer to a partially unfolded intermediate containing a mixture of dimers and trimers and then completely unfolds to unstructured monomers.     

Visualization of energetics and conformations from molecular computer simulations

The quantity of data generated from molecular dynamics simulations and energy minimizations of macromolecules is overwhelming. It is an arduous task to extract the relevant and interesting information from the numerous coordinate sets produced. To help solve this problem, the authors have developed a method to aid the visualization of the relevant information from the simulations. This approach combines animation of the results on a high performance graphics device, such as the PS300, with colour-coded atoms based on changes in energy or conformation. The method will be illustrated using as examples: the molecular mechanics minimization of a nonapeptide, the molecular dynamics simulation of the protein myoglobin, including the analysis of the motion of helices during a 300ps trajectory, and changes in sugar puckering that occur during the molecular dynamics simulation of a DNA oligomer. The method is also applicable for analysing energy components and conformational properties of a fixed conformation.     

Simulating proteins at constant pH: An approach combining molecular dynamics and Monte Carlo simulation

For the structure and function of proteins, the pH of the solution is one of the determining parameters. Current molecular dynamics (MD) simulations account for the solution pH only in a limited way by keeping each titratable site in a chosen protonation state. We present an algorithm that generates trajectories at a Boltzmann distributed ensemble of protonation states by a combination of MD and Monte Carlo (MC) simulation. The algorithm is useful for pH-dependent structural studies and to investigate in detail the titration behavior of proteins. The method is tested on the acidic residues of the protein hen egg white lysozyme. It is shown that small structural changes may have a big effect on the pK(A) values of titratable residues.     

pK(a) calculations along a bacteriorhodopsin molecular dynamics trajectory

Electrostatic calculations of pK(a-values) are reported along a 400 ps molecular dynamics trajectory of bacteriorhodopsin. The sensitivity of calculated pK(a) values to a number of structural factors and factors related to the modelling of the electrostatics are also studied. The results are very sensitive to the choice of internal dielectric constant of the protein (in the interval 2-4). Moreover it is important to include internal water molecules and to average over a long enough portion ( approximately 100 ps) of an equilibrium molecular dynamics trajectory. The internal waters are necessary to get an ion-counter ion complex with the Schiff base and Arg 82 protonated and the aspartic groups (85 and 212) deprotonated. The fluctuations along the MD-trajectory do not change the protonation state of internal residues at neutral pH. However, at other pH values the averaging along a trajectory maybe crucial to get correct protonation states. A relationship is found between the arginine group 82, the aspartic group 85 and the glutamate group 204. Glu 204 is protonated in the ground state but the pK(a) value decreases towards deprotonation when the chromophore isomerizes into the cis state.     

Conformational energy analysis of the leucine repeat regions of C/EBP,GCN4, and the proteins of themyc,jun, andfos oncogenes

It has been recently proposed that certain DNA binding proteins (including C/EBP, GCN4 and themyc, jun, andfos oncogene proteins) share a common structural motif based on helix-promoting regions containing heptad repeat sequences of leucines. It has been suggested that this structure is critical to the biological activity of these proteins, since it facilitates the formation of functional dimers held together by interdigitating leucine side-chains along the hydrophobic interfaces between long α-helical regions of the polypeptide chains in a configuration termed the “leucine zipper.” In this paper, conformational energy analysis is used to determine the preferred three-dimensional structures of the leucine repeat regions of these proteins. The results indicate that, in all cases, the global minimum energy conformation for these regions is an amphipathic α-helix with the leucine side-chains arrayed on one side in such a way to favor “leucine zipper” dimerization. Furthermore, amino acid substitutions in these regions (such as Pro for Leu), that are known to inhibit dimer formation and prevent DNA binding, are found to produce significant conformational changes that disrupt the amphipathic helical structure. Thus, these results provide support for the proposed “leucine zipper” configuration as a critical structural feature of this class of DNA binding proteins.

     

DNA-induced increase in the alpha-helical content of C/EBP and GCN4

Leucine zipper proteins comprise a recently identified class of DNA binding proteins that contain a bipartite structural motif consisting of a "leucine zipper" dimerization domain and a segment rich in basic residues responsible for DNA interaction. Protein fragments encompassing the zipper plus basic region domains (bZip) have previously been used to determine the conformational and dynamic properties of this motif. In the absence of DNA, the coiled-coil portion is alpha-helical and dimeric, whereas the basic region is flexible and partially disordered. Addition of DNA containing a specific recognition sequence induces a fully helical conformation in the basic regions of these fragments. However, the question remained whether the same conformational change would be observed in native bZip proteins where the basic regions might be stabilized in an alpha-helical conformation even in the absence of DNA, through interactions with portions of the protein not included in the bZip motif. We have now examined the DNA-induced conformational transition for an intact bZip protein, GCN4, and for the bZip fragment of C/EBP with two enhancers that are differentially symmetric. Our results are consistent with the induced helical fork model wherein the basic regions are largely flexible in the absence of DNA and become fully helical in the presence of the specific DNA recognition sequence.     

Calculation of symmetric multimer structures from NMR data using a priori knowledge of the monomer structure, co-monomer restraints, and interface mapping: The case of leucine zippers

Summary NMR studies of symmetric multimers are problematic due to the difficulty in distinguishing between intra-, inter-, and co-monomer (mixed) NOE signals. Previously, one of us described a general calculation strategy called dynamic assignment by which this difficulty can be overcome [Nilges, M. (1993) Proteins, 17, 297–309]. Here we describe extensions to the method for handling many co-monomer NOEs and for taking advantege of prior knowledge of monomer structures. The new protocol was developed for the particularly difficult case of leucine zipper (LZ) homodimers, for which the previous protocol proved inefficient. In addition to the problem of dimer symmetry, LZs have a particularly high proportion of co-monomer NOE signals and a high degree of repetition in sequence and structure, leading to significant spectral overlap. Furthermore, the leucine zipper is a rather extended (as opposed to globular) protein domain; accurately determining such a structure based only on the very short distances obtainable by NMR is clearly a challenge to the NMR structure determination method. We have previously shown that, for LZ homodimers, many of the backbone-backbone NOESY cross peaks can be unambiguously assigned as intra-monomer, enabling approximate monomer structures to be calculated. Using model and experimental data sets, we verified that the new protocol converges to the correct dimer structure. The results show that short-range NMR distance data can be sufficient to define accurately the extended LZ. The protocol has been used to derive a novel solution structure of the c-Jun LZ domain. Based on these calculations, we propose the protocol as a prototype for the general case of symmetric multimers where the monomer structure is known.Abbreviations 3D three-dimensional - GCN4-c crystal structure of the GCN4 LZ homodimer - GCN4-s solution structures of GCN4 - GSYM global symmetry - Jun-m model structure of the Jun LZ homodimer - Jun-s solution structure of Jun - LZ leucine zipper - MFP mean force potential - MDSA molecular dynamical simulated annealing - NCS noncrystallographic symmetry - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - vdW van der Waals     

Toward accurate prediction of pKa values for internal protein residues: the importance of conformational relaxation and desolvation energy

Proton uptake or release controls many important biological processes, such as energy transduction, virus replication, and catalysis. Accurate pK(a) prediction informs about proton pathways, thereby revealing detailed acid-base mechanisms. Physics-based methods in the framework of molecular dynamics simulations not only offer pK(a) predictions but also inform about the physical origins of pK(a) shifts and provide details of ionization-induced conformational relaxation and large-scale transitions. One such method is the recently developed continuous constant pH molecular dynamics (CPHMD) method, which has been shown to be an accurate and robust pK(a) prediction tool for naturally occurring titratable residues. To further examine the accuracy and limitations of CPHMD, we blindly predicted the pK(a) values for 87 titratable residues introduced in various hydrophobic regions of staphylococcal nuclease and variants. The predictions gave a root-mean-square deviation of 1.69 pK units from experiment, and there were only two pK(a)'s with errors greater than 3.5 pK units. Analysis of the conformational fluctuation of titrating side-chains in the context of the errors of calculated pK(a) values indicate that explicit treatment of conformational flexibility and the associated dielectric relaxation gives CPHMD a distinct advantage. Analysis of the sources of errors suggests that more accurate pK(a) predictions can be obtained for the most deeply buried residues by improving the accuracy in calculating desolvation energies. Furthermore, it is found that the generalized Born implicit-solvent model underlying the current CPHMD implementation slightly distorts the local conformational environment such that the inclusion of an explicit-solvent representation may offer improvement of accuracy.     

Crystal structure of a super leucine zipper,an extended two‐stranded super long coiled coil

Coiled coil is a ubiquitous structural motif in proteins, with two to seven alpha helices coiled together like the strands of a rope, and coiled coil folding and assembly is not completely understood. A GCN4 leucine zipper mutant with four mutations of K3A, D7A, Y17W, and H18N has been designed, and the crystal structure has been determined at 1.6 Å resolution. The peptide monomer shows a helix trunk with short curved N‐ and C‐termini. In the crystal, two monomers cross in 35° and form an X‐shaped dimer, and each X‐shaped dimer is welded into the next one through sticky hydrophobic ends, thus forming an extended two‐stranded, parallel, super long coiled coil rather than a discrete, two‐helix coiled coil of the wild‐type GCN4 leucine zipper. Leucine residues appear at every seventh position in the super long coiled coil, suggesting that it is an extended super leucine zipper. Compared to the wild‐type leucine zipper, the N‐terminus of the mutant has a dramatic conformational change and the C‐terminus has one more residue Glu 32 determined. The mutant X‐shaped dimer has a large crossing angle of 35° instead of 18° in the wild‐type dimer. The results show a novel assembly mode and oligomeric state of coiled coil, and demonstrate that mutations may affect folding and assembly of the overall coiled coil. Analysis of the formation mechanism of the super long coiled coil may help understand and design self‐assembling protein fibers.     

Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins

Protein stability and function relies on residues being in their appropriate ionization states at physiological pH. In situ residue pK(a)s also provides a sensitive measure of the local protein environment. Multiconformation continuum electrostatics (MCCE) combines continuum electrostatics and molecular mechanics force fields in Monte Carlo sampling to simultaneously calculate side chain ionization and conformation. The response of protein to charges is incorporated both in the protein dielectric constant (epsilon(prot)) of four and by explicit conformational changes. The pK(a) of 166 residues in 12 proteins was determined. The root mean square error is 0.83 pH units, and >90% have errors of <1 pH units whereas only 3% have errors >2 pH units. Similar results are found with crystal and solution structures, showing that the method's explicit conformational sampling reduces sensitivity to the initial structure. The outcome also changes little with protein dielectric constant (epsilon(prot) 4-20). Multiconformation continuum electrostatics titrations show coupling of conformational flexibility and changes in ionization state. Examples are provided where ionizable side chain position (protein G), Asn orientation (lysozyme), His tautomer distribution (RNase A), and phosphate ion binding (RNase A and H) change with pH. Disallowing these motions changes the calculated pK(a).     

Stabilization of coiled-coil peptide domains by introduction of trifluoroleucine

Substitution of leucine residues by 5,5,5-trifluoroleucine at the d-positions of the leucine zipper peptide GCN4-p1d increases the thermal stability of the coiled-coil structure. The midpoint thermal unfolding temperature of the fluorinated peptide is elevated by 13 degrees C at 30 microM peptide concentration. The modified peptide is more resistant to chaotropic denaturants, and the free energy of folding of the fluorinated peptide is 0.5-1.2 kcal/mol larger than that of the hydrogenated form. A similarly fluorinated form of the DNA-binding peptide GCN4-bZip binds to target DNA sequences with affinity and specificity identical to those of the hydrogenated form, while demonstrating enhanced thermal stability. Molecular dynamics simulation on the fluorinated GCN4-p1d peptide using the Surface Generalized Born implicit solvation model revealed that the coiled-coil binding energy is 55% more favorable upon fluorination. These results suggest that fluorination of hydrophobic substructures in peptides and proteins may provide new means of increasing protein stability, enhancing protein assembly, and strengthening receptor-ligand interactions.