Unfolding of a leucine zipper is not a simple two-state transition

MEILB2-BRME1 forms a V-shaped DNA clamp upon BRCA2-binding in meiotic recombination

DNA double-strand break repair by homologous recombination has a specialised role in meiosis by generating crossovers that enable the formation of haploid germ cells. This requires meiosis-specific MEILB2-BRME1, which interacts with BRCA2 to facilitate loading of recombinases onto resected DNA ends. Here, we report the crystal structure of the MEILB2-BRME1 2:2 core complex, revealing a parallel four-helical assembly that recruits BRME1 to meiotic double-strand breaks in vivo. It forms an N-terminal β-cap that binds to DNA, and a MEILB2 coiled-coil that bridges to C-terminal ARM domains. Upon BRCA2-binding, MEILB2-BRME1 2:2 complexes dimerize into a V-shaped 2:4:4 complex, with rod-like MEILB2-BRME1 components arranged at right-angles. The β-caps located at the tips of the MEILB2-BRME1 limbs are separated by 25 nm, allowing them to bridge between DNA molecules. Thus, we propose that BRCA2 induces MEILB2-BRME1 to function as a DNA clamp, connecting resected DNA ends or homologous chromosomes to facilitate meiotic recombination.

Thermodynamics and S-Palmitoylation Dependence of Interactions between Human Aquaporin-4 M1 Tetramers in Model Membranes

Aquaporin-4 (AQP4) is a water channel protein found primarily in the central nervous system (CNS) that helps to regulate water-ion homeostasis. AQP4 exists in two major isoforms: M1 and M23. While both isoforms have a homotetrameric quaternary structure and are functionally identical when transporting water, the M23 isoform forms large protein aggregates known as orthogonal arrays of particles (OAPs). In contrast, the M1 isoform creates a peripheral layer around the outside of these OAPs, suggesting a thermodynamically stable interaction between the two. Structurally, the M1 isoform has an N-terminal tail that is 22 amino acids longer than the M23 isoform and contains two solvent-accessible cysteines available for S-palmitoylation at cysteine-13 (Cys-13) and cysteine-17 (Cys-17) in the amino acid sequence. Earlier work suggests that the palmitoylation of these cysteines might aid in regulating AQP4 assemblies. This work discusses the thermodynamic driving forces for M1 protein-protein interactions and how the palmitoylation state of M1 affects them. Using temperature-dependent single-particle tracking, the standard state free energies, enthalpies, and entropies were measured for these interactions. Furthermore, we present a binding model based on measured thermodynamics and a structural modeling study. The results of this study demonstrate that the M1 isoform will associate with itself according to the following expressions: 2[AQP4-M1]4 ↔ [[AQP4-M1]4]2 when palmitoylated and 3[AQP4-M1]4 ↔ [AQP4-M1]4 + [[AQP4-M1]4]2 ↔ [[AQP4-M1]4]3 when depalmitoylated. This is primarily due to a conformational change induced by adding the palmitic acid groups at Cys-13 and Cys-17 in the N-terminal tails of the homotetramers. In addition, a statistical mechanical model was developed to estimate the Gibbs free energy, enthalpy, and entropy for forming dimers and trimers. These results were in good agreement with experimental values.

Heat sensitive E-helix cut ferritin nanocages for facile and high-efficiency loading of doxorubicin

Ferritin possesses a stable and uniform cage structure, along with tumor-targeting properties and excellent biocompatibility, making it a promising drug delivery vehicle. However, the current ferritin drug loading strategy involves complex steps and harsh reaction conditions, resulting in low yield and recovery of drug loading, which limits the clinical application prospects of ferritin nanomedicine. In this study, we utilized the high-efficiency heat-sensitivity of the multiple channel switch structures of the E-helix-cut ferritin mutant (Ecut-HFn) and Cu2+ assistance to achieve high-efficiency loading of chemotherapeutic drugs in a one-step process at low temperatures. This method features mild reaction conditions (45 °C), high loading efficiency (about 110 doxorubicin (Dox) per Ecut-HFn), and improved protein and Dox recovery rates (with protein recovery rate around 94 % and Dox recovery rate reaching up to 45 %). The prepared ferritin-Dox particles (Ecut-HFn-Cu-Dox) exhibit a uniform size distribution, good stability, and retain the natural tumor targeting ability of ferritin. Overall, this temperature-controlled drug loading strategy utilizing heat-sensitivity ferritin mutants is energy-saving, environmentally friendly, efficient, and easy to operate, offering a new perspective for scaling up the industrial production of ferritin drug carriers.

Decoding the role of coiled-coil motifs in human prion-like proteins

Prions are self-propagating proteins that cause fatal neurodegenerative diseases in humans. However, increasing evidence suggests that eukaryotic cells exploit prion conformational conversion for functional purposes. A recent study delineated a group of twenty prion-like proteins in humans, characterized by the presence of low-complexity glutamine-rich sequences with overlapping coiled-coil (CCs) motifs. This is the case of Mediator complex subunit 15 (MED15), which is overexpressed in a wide range of human cancers. Biophysical studies demonstrated that the prion-like domain (PrLD) of MED15 forms homodimers in solution, sustained by CCs interactions. Furthermore, the same coiled-coil (CC) region plays a crucial role in the PrLD structural transition to a transmissible β-sheet amyloid state. In this review, we discuss the role of CCs motifs and their contribution to amyloid transitions in human prion-like domains (PrLDs), while providing a comprehensive overview of six predicted human prion-like proteins involved in transcription, gene expression, or DNA damage response and associated with human disease, whose PrLDs contain or overlap with CCs sequences. Finally, we try to rationalize how these molecular signatures might relate to both their function and involvement in disease.

Combining Constrained Heptapeptide Cassettes with Computational Design To Create Coiled-Coil Targeting Helical Peptides

A total of 32 heptapeptides have been synthesized and characterized to establish the effect of K → D (i → i + 4) lactamization upon their ability to adopt a helical conformation. Because most parallel and dimeric coiled-coil sequences can be deconvoluted into gabcdef repeats, we have introduced fixed solvent exposed b → f (K → D) constraints into this design scaffold. Interfacial " a" hydrophobic (L/I/V/N) and " e/g" electrostatic (E/K) options (4 × 2 × 2 = 16 cassettes) were introduced as core drivers of coiled-coil stability and specificity. All present as random coils when linear but adopt a helical conformation upon lactamization. Helicity varied in magnitude from 34 to 68%, indicating different levels of constraint tolerance within the context of a sequence required to be helical for function. Using the oncogenic transcription factor cJun as an exemplar, we next utilized our bCIPA coiled-coil screening engine to select four cassettes of highest predicted affinity when paired with four gabcdef cassettes within the full-length cJun target counterpart (164 = 65 536 combinations). This information was coupled with observed helicity for each constrained cassette to select for the best balance of predicted affinity when linear and experimentally validated helicity when constrained. As a control, the same approach was taken using cassettes of high predicted target affinity but with lower experimentally validated helicity. The approach provides a novel platform of modular heptapeptide cassettes experimentally validated and separated by helical content. Appropriate cassettes can be selected and conjugated to produce longer peptides, in which constraints impart appropriate helicity such that a wide range of targets can be engaged with high affinity and selectivity.

Atomistic simulations indicate the functional loop-to-coiled-coil transition in influenza hemagglutinin is not downhill

Influenza hemagglutinin (HA) mediates viral entry into host cells through a large-scale conformational rearrangement at low pH that leads to fusion of the viral and endosomal membranes. Crystallographic and biochemical data suggest that a loop-to-coiled-coil transition of the B-loop region of HA is important for driving this structural rearrangement. However, the microscopic picture for this proposed "spring-loaded" movement is missing. In this study, we focus on understanding the transition of the B loop and perform a set of all-atom molecular dynamics simulations of the full B-loop trimeric structure with the CHARMM36 force field. The free-energy profile constructed from our simulations describes a B loop that stably folds half of the postfusion coiled coil in tens of microseconds, but the full coiled coil is unfavorable. A buried hydrophilic residue, Thr59, is implicated in destabilizing the coiled coil. Interestingly, this conserved threonine is the only residue in the B loop that strictly differentiates between the group 1 and 2 HA molecules. Microsecond-scale constant temperature simulations revealed that kinetic traps in the structural switch of the B loop can be caused by nonnative, intramonomer, or intermonomer β-sheets. The addition of the A helix stabilized the postfusion state of the B loop, but introduced the possibility for further β-sheet structures. Overall, our results do not support a description of the B loop in group 2 HAs as a stiff spring, but, rather, it allows for more structural heterogeneity in the placement of the fusion peptides during the fusion process.

Downsizing Proto-oncogene cFos to Short Helix-Constrained Peptides That Bind Jun

The oncogenic transcription factor activator protein-1 (AP-1) is a DNA-binding protein that assembles through dimerization of Fos and Jun protein subunits, their leucine-rich helical sequences entwining into a coiled-coil structure. This study reports on downsizing the proto-oncogene cFos protein (380 residues) to shorter peptides (37-25 residues) modified with helix-inducing constraints to enhance binding to Jun. A crystal structure is reported for a 37-residue Fos-derived peptide (FosW) bound to Jun. This guided iterative downsizing of FosW to shorter peptide sequences that were constrained into stable water-soluble α-helices by connecting amino acid side chains to form cyclic pentapeptide components. Structural integrity in the presence and absence of Jun was assessed by circular dichroism spectroscopy, while the thermodynamics of binding to cFos was measured by isothermal titration calorimetry. A 25-residue constrained peptide, one-third shorter yet 25% more helical than the structurally characterized 37-residue Fos-derived peptide, retained 80% of the binding free energy as a result of preorganization in a Jun-binding helix conformation, with the entropy gain (TΔS = +3.2 kcal/mol) compensating for the enthalpy loss. Attaching a cell-penetrating peptide (TAT_48-57) and a nuclear localization signal (SV40) promoted cell uptake, localization to the nucleus, and inhibition of the proliferation of two breast cancer cell lines.

A Coiled-Coil Peptide Shaping Lipid Bilayers upon Fusion

A system based on two designed peptides, namely the cationic peptide K, (KIAALKE)3, and its complementary anionic counterpart called peptide E, (EIAALEK)3, has been used as a minimal model for membrane fusion, inspired by SNARE proteins. Although the fact that docking of separate vesicle populations via the formation of a dimeric E/K coiled-coil complex can be rationalized, the reasons for the peptides promoting fusion of vesicles cannot be fully explained. Therefore it is of significant interest to determine how the peptides aid in overcoming energetic barriers during lipid rearrangements leading to fusion. In this study, investigations of the peptides' interactions with neutral PC/PE/cholesterol membranes by fluorescence spectroscopy show that tryptophan-labeled K∗ binds to the membrane (KK∗ ∼6.2 103 M-1), whereas E∗ remains fully water-solvated. 15N-NMR spectroscopy, depth-dependent fluorescence quenching, CD-spectroscopy experiments, and MD simulations indicate a helix orientation of K∗ parallel to the membrane surface. Solid-state 31P-NMR of oriented lipid membranes was used to study the impact of peptide incorporation on lipid headgroup alignment. The membrane-immersed K∗ is found to locally alter the bilayer curvature, accompanied by a change of headgroup orientation relative to the membrane normal and of the lipid composition in the vicinity of the bound peptide. The NMR results were supported by molecular dynamics simulations, which showed that K reorganizes the membrane composition in its vicinity, induces positive membrane curvature, and enhances the lipid tail protrusion probability. These effects are known to be fusion relevant. The combined results support the hypothesis for a twofold role of K in the mechanism of membrane fusion: 1) to bring opposing membranes into close proximity via coiled-coil formation and 2) to destabilize both membranes thereby promoting fusion.

Coiled coils as possible models of protein structure evolution

Coiled coils are formed by two or more α-helices wrapped around one another. This structural motif often guides di-, tri- or multimerization of proteins involved in diverse biological processes such as membrane fusion, signal transduction and the organization of the cytoskeleton. Although coiled coil motifs seem conceptually simple and their existence was proposed in the early 1950s, the high variability of the motif makes coiled coil prediction from sequence a difficult task. They might be confused with intrinsically disordered sequences and even more with a recently described structural motif, the charged single α-helix. By contrast, the versatility of coiled coil structures renders them an ideal candidate for protein (re)design and many novel variants have been successfully created to date. In this paper, we review coiled coils in the light of protein evolution by putting our present understanding of the motif and its variants in the context of structural interconversions. We argue that coiled coils are ideal subjects for studies of subtle and large-scale structural changes because of their well-characterized and versatile nature.

Interplay between Lipid Interaction and Homo-coiling of Membrane-Tethered Coiled-Coil Peptides

The designed coiled-coil-forming peptides E [(EIAALEK)3] and K [(KIAALKE)3] are known to trigger efficient membrane fusion when they are tethered to lipid vesicles in the form of lipopeptides. Knowledge of their secondary structure is a key element in understanding their role in membrane fusion. Special conditions can be found at the interface of the membrane, where the peptides are confined in close proximity to other peptide molecules as well as to the lipid interface. Consequently, different structural states were proposed for the peptides when tethered to this interface. Due to the multitude of possible states, determining the structure solely on the basis of circular dichroism (CD) spectra at a single temperature can be misleading. In addition, it has not yet been possible to unambiguously distinguish between the membrane-bound and the coiled-coil states of these peptides by means of infrared (IR) spectroscopy due to their very similar amide I' bands. Here, the molecular basis of this similarity is investigated by means of site-specific (13)C-labeled FTIR spectroscopy. Structural similarities between the membrane-interacting helix of K and the homo-coiled-coil-forming helix of E are shown to cause the similar spectroscopic properties. Furthermore, the peptide structure is investigated using temperature-dependent CD and IR spectroscopy, and it is shown that the different states can be distinguished on the basis of their thermal behavior. It is shown that the two peptides behave fundamentaly differently when tethered to the lipid membrane, which implies that their role during membrane fusion is different and the mechanism of this process is asymmetric.

A stable trimeric influenza hemagglutinin stem as a broadly protective immunogen

The identification of human broadly neutralizing antibodies (bnAbs) targeting the hemagglutinin (HA) stem revitalized hopes of developing a universal influenza vaccine. Using a rational design and library approach, we engineered stable HA stem antigens ("mini-HAs") based on an H1 subtype sequence. Our most advanced candidate exhibits structural and bnAb binding properties comparable to those of full-length HA, completely protects mice in lethal heterologous and heterosubtypic challenge models, and reduces fever after sublethal challenge in cynomolgus monkeys. Antibodies elicited by this mini-HA in mice and nonhuman primates bound a wide range of HAs, competed with human bnAbs for HA stem binding, neutralized H5N1 viruses, and mediated antibody-dependent effector activity. These results represent a proof of concept for the design of HA stem mimics that elicit bnAbs against influenza A group 1 viruses.

Dihedral angle entropy measures for intrinsically disordered proteins

Protein stability is based on a delicate balance between energetic and entropic factors. Intrinsically disordered proteins (IDPs) interacting with a folded partner protein in the act of binding can order the IDP to form the correct functional interface by decrease in the overall free energy. In this work, we evaluate the part of the entropic cost of ordering an IDP arising from their dihedral states. The IDP studied is a leucine zipper dimer that we simulate with molecular dynamics and find that it does show disorder in six phi and psi dihedral angles of the N terminal sequence of one monomer. Essential to ascertain is the degree of disorder in the IDP, and we do so by considering the entire, discretized probability distribution function of N dihedrals with M conformers per dihedral. A compositional clustering method is introduced, whereby the NS = N(M) states are formed from the Cartesian product of each dihedral's conformational space. Clustering is carried out with a version of a k-means algorithm that accounts for the circular nature of dihedral angles. For the 12 dihedrals each found to have three conformers, among the resulting 531441 states, their populations show that the first 100 (500) most populated states account for ∼65% (∼90%) of the entire population, indicating that there are strong dependencies among the dihedrals' conformations. These state populations are used to evaluate a Kullback-Leibler divergence entropy measure and obtain the dihedral configurational entropy S. At 300 K, TS ∼ 3 kcal/mol, showing that IDP entropy, while roughly half that would be expected from independently distributed dihedrals, can be a decisive contributor to the free energy of this IDP binding and ordering.

pH-jump induced leucine zipper folding beyond the diffusion limit

The folding of a pH-sensitive leucine zipper, that is, a GCN4 mutant containing eight glutamic acid residues, has been investigated. A pH-jump induced by a caged proton (o-nitrobenzaldehyde, oNBA) is employed to initiate the process, and time-resolved IR spectroscopy of the amide I band is used to probe it. The experiment has been carefully designed to minimize the buffer capacity of the sample solution so that a large pH jump can be achieved, leading to a transition from a completely unfolded to a completely folded state with a single laser shot. In order to eliminate the otherwise rate-limiting diffusion-controlled step of the association of two peptides, they have been covalently linked. The results for the folding kinetics of the cross-linked peptide are compared with those of an unlinked peptide, which reveals a detailed picture of the folding mechanism. That is, folding occurs in two steps, one on an ∼1-2 μs time scale leading to a partially folded α-helix even in the monomeric case and a second one leading to the final coiled-coil structure on distinctively different time scales of ∼30 μs for the cross-linked peptide and ∼200 μs for the unlinked peptide. By varying the initial pH, it is found that the folding mechanism is consistent with a thermodynamic two-state model, despite the fact that a transient intermediate is observed in the kinetic experiment.

Unique features of the anti-parallel, heterodimeric coiled-coil interaction between methyl-cytosine binding domain 2 (MBD2) homologues and GATA zinc finger domain containing 2A (GATAD2A/p66α)

The methyl-cytosine binding domain 2 (MBD2)-nucleosome remodeling and deacetylase (NuRD) complex recognizes methylated DNA and silences expression of associated genes through histone deacetylase and nucleosome remodeling functions. Our previous structural work demonstrated that a coiled-coil interaction between MBD2 and GATA zinc finger domain containing 2A (GATAD2A/p66α) proteins recruits the chromodomain helicase DNA-binding protein (CHD4/Mi2β) to the NuRD complex and is necessary for MBD2-mediated DNA methylation-dependent gene silencing in vivo (Gnanapragasam, M. N., Scarsdale, J. N., Amaya, M. L., Webb, H. D., Desai, M. A., Walavalkar, N. M., Wang, S. Z., Zu Zhu, S., Ginder, G. D., and Williams, D. C., Jr. (2011) p66α-MBD2 coiled-coil interaction and recruitment of Mi-2 are critical for globin gene silencing by the MBD2-NuRD complex. Proc. Natl. Acad. Sci. U.S.A. 108, 7487-7492). The p66α-MBD2 interaction differs from most coiled-coils studied to date by forming an anti-parallel heterodimeric complex between two peptides that are largely monomeric in isolation. To further characterize unique features of this complex that drive heterodimeric specificity and high affinity binding, we carried out biophysical analyses of MBD2 and the related homologues MBD3, MBD3-like protein 1 (MBD3L1), and MBD3-like protein 2 (MBD3L2) as well as specific mutations that modify charge-charge interactions and helical propensity of the coiled-coil domains. Analytical ultracentrifugation analyses show that the individual peptides remain monomeric in isolation even at 300 μM in concentration for MBD2. Circular dichroism analyses demonstrate a direct correlation between helical content of the coiled-coil domains in isolation and binding affinity for p66α. Furthermore, complementary electrostatic surface potentials and inherent helical content of each peptide are necessary to maintain high-affinity association. These factors lead to a binding affinity hierarchy of p66α for the different MBD2 homologues (MBD2 ≈ MBD3 > MBD3L1 ≈ MBD3L2) and suggest a hierarchical regulatory model in tissue and life cycle stage-specific silencing by NuRD complexes.

Energetics of oligomeric protein folding and association

In nature, proteins most often exist as complexes, with many of these consisting of identical subunits. Understanding of the energetics governing the folding and misfolding of such homooligomeric proteins is central to understanding their function and misfunction, in disease or biotechnology. Much progress has been made in defining the mechanisms and thermodynamics of homooligomeric protein folding. In this review, we outline models as well as calorimetric and spectroscopic methods for characterizing oligomer folding, and describe extensive results obtained for diverse proteins, ranging from dimers to octamers and higher order aggregates. To our knowledge, this area has not been reviewed comprehensively in years, and the collective progress is impressive. The results provide evolutionary insights into the development of subunit interfaces, mechanisms of oligomer folding, and contributions of oligomerization to protein stability, function and regulation. Thermodynamic analyses have also proven valuable for understanding protein misfolding and aggregation mechanisms, suggesting new therapeutic avenues. Successful recent designs of novel, functional proteins demonstrate increased understanding of oligomer folding. Further rigorous analyses using multiple experimental and computational approaches are still required, however, to achieve consistent and accurate prediction of oligomer folding energetics. Modeling the energetics remains challenging but is a promising avenue for future advances.

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.

Structural analysis of the human SYCE2-TEX12 complex provides molecular insights into synaptonemal complex assembly

The successful completion of meiosis is essential for all sexually reproducing organisms. The synaptonemal complex (SC) is a large proteinaceous structure that holds together homologous chromosomes during meiosis, providing the structural framework for meiotic recombination and crossover formation. Errors in SC formation are associated with infertility, recurrent miscarriage and aneuploidy. The current lack of molecular information about the dynamic process of SC assembly severely restricts our understanding of its function in meiosis. Here, we provide the first biochemical and structural analysis of an SC protein component and propose a structural basis for its function in SC assembly. We show that human SC proteins SYCE2 and TEX12 form a highly stable, constitutive complex, and define the regions responsible for their homotypic and heterotypic interactions. Biophysical analysis reveals that the SYCE2-TEX12 complex is an equimolar hetero-octamer, formed from the association of an SYCE2 tetramer and two TEX12 dimers. Electron microscopy shows that biochemically reconstituted SYCE2-TEX12 complexes assemble spontaneously into filamentous structures that resemble the known physical features of the SC central element (CE). Our findings can be combined with existing biological data in a model of chromosome synapsis driven by growth of SYCE2-TEX12 higher-order structures within the CE of the SC.

A basis set of de novo coiled-coil peptide oligomers for rational protein design and synthetic biology

Protein engineering, chemical biology, and synthetic biology would benefit from toolkits of peptide and protein components that could be exchanged reliably between systems while maintaining their structural and functional integrity. Ideally, such components should be highly defined and predictable in all respects of sequence, structure, stability, interactions, and function. To establish one such toolkit, here we present a basis set of de novo designed α-helical coiled-coil peptides that adopt defined and well-characterized parallel dimeric, trimeric, and tetrameric states. The designs are based on sequence-to-structure relationships both from the literature and analysis of a database of known coiled-coil X-ray crystal structures. These give foreground sequences to specify the targeted oligomer state. A key feature of the design process is that sequence positions outside of these sites are considered non-essential for structural specificity; as such, they are referred to as the background, are kept non-descript, and are available for mutation as required later. Synthetic peptides were characterized in solution by circular-dichroism spectroscopy and analytical ultracentrifugation, and their structures were determined by X-ray crystallography. Intriguingly, a hitherto widely used empirical rule-of-thumb for coiled-coil dimer specification does not hold in the designed system. However, the desired oligomeric state is achieved by database-informed redesign of that particular foreground and confirmed experimentally. We envisage that the basis set will be of use in directing and controlling protein assembly, with potential applications in chemical and synthetic biology. To help with such endeavors, we introduce Pcomp, an on-line registry of peptide components for protein-design and synthetic-biology applications.

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.

Structural and functional insights into DNA-end processing by the archaeal HerA helicase-NurA nuclease complex

Helicase–nuclease systems dedicated to DNA end resection in preparation for homologous recombination (HR) are present in all kingdoms of life. In thermophilic archaea, the HerA helicase and NurA nuclease cooperate with the highly conserved Mre11 and Rad50 proteins during HR-dependent DNA repair. Here we show that HerA and NurA must interact in a complex with specific subunit stoichiometry to process DNA ends efficiently. We determine crystallographically that NurA folds in a toroidal dimer of intertwined RNaseH-like domains. The central channel of the NurA dimer is too narrow for double-stranded DNA but appears well suited to accommodate one or two strands of an unwound duplex. We map a critical interface of the complex to an exposed hydrophobic epitope of NurA abutting the active site. Based upon the presented evidence, we propose alternative mechanisms of DNA end processing by the HerA-NurA complex.

A hamiltonian replica exchange method for building protein-protein interfaces applied to a leucine zipper

Leucine zippers consist of alpha helical monomers dimerized (or oligomerized) into alpha superhelical structures known as coiled coils. Forming the correct interface of a dimer from its monomers requires an exploration of configuration space focused on the side chains of one monomer that must interdigitate with sites on the other monomer. The aim of this work is to generate good interfaces in short simulations starting from separated monomers. Methods are developed to accomplish this goal based on an extension of a previously introduced [Su and Cukier, J. Phys. Chem. B 113, 9595, (2009)] hamiltonian temperature replica exchange method (HTREM), which scales the hamiltonian in both potential and kinetic energies that was used for the simulation of dimer melting curves. The new method, HTREM_MS (MS designates mean square), focused on interface formation, adds restraints to the hamiltonians for all but the physical system, which is characterized by the normal molecular dynamics force field at the desired temperature. The restraints in the nonphysical systems serve to prevent the monomers from separating too far, and have the dual aims of enhancing the sampling of close in configurations and breaking unwanted correlations in the restrained systems. The method is applied to a 31-residue truncation of the 33-residue leucine zipper (GCN4-p1) of the yeast transcriptional activator GCN4. The monomers are initially separated by a distance that is beyond their capture length. HTREM simulations show that the monomers oscillate between dimerlike and monomerlike configurations, but do not form a stable interface. HTREM_MS simulations result in the dimer interface being faithfully reconstructed on a 2 ns time scale. A small number of systems (one physical and two restrained with modified potentials and higher effective temperatures) are sufficient. An in silico mutant that should not dimerize because it lacks charged residues that provide electrostatic stabilization of the dimer does not with HTREM_MS, giving confidence in the method. The interface formation time scale is sufficiently short that using HTREM_MS as a screening tool to validate leucine zipper design methods may be feasible.

Enhancing β3 -peptide bundle stability by design

We reported recently that certain β(3) -peptides self-assemble in aqueous solution into discrete bundles of unique structure and defined stoichiometry. The first β-peptide bundle reported was the octameric Zwit-1F, whose fold is characterized by a well-packed, leucine-rich core and a salt-bridge-rich surface. Close inspection of the Zwit-1F structure revealed four nonideal interhelical salt-bridge interactions whose heavy atom-heavy atom distances were longer than found in natural proteins of known structure. Here we demonstrate that the thermodynamic stability of a β-peptide bundle can be enhanced by optimizing the length of these four interhelical salt bridges. Combined with previous work on the role of internal packing residues, these results provide another critical step in the "bottom-up" formation of β-peptide assemblies with defined sizes, reproducible structures, and sophisticated function.

Conserved noncanonical residue Gly-126 confers instability to the middle part of the tropomyosin molecule

Tropomyosin (Tm) is a two-stranded α-helical coiled-coil protein with a well established role in regulation of actin cytoskeleton and muscle contraction. It is believed that many Tm functions are enabled by its flexibility whose nature has not been completely understood. We hypothesized that the well conserved non-canonical residue Gly-126 causes local destabilization of Tm. To test this, we substituted Gly-126 in skeletal muscle α-Tm either with an Ala residue, which should stabilize the Tm α-helix, or with an Arg residue, which is expected to stabilize both α-helix and coiled-coil structure of Tm. We have shown that both mutations dramatically reduce the rate of Tm proteolysis by trypsin at Asp-133. Differential scanning calorimetry was used for detailed investigation of thermal unfolding of the Tm mutants, both free in solution and bound to F-actin. It was shown that a significant part of wild type Tm unfolds in a non-cooperative manner at low temperature, and both mutations confer cooperativity to this part of the Tm molecule. The size of the flexible middle part of Tm is estimated to be 60-70 amino acid residues, about a quarter of the Tm molecule. Thus, our results show that flexibility is unevenly distributed in the Tm molecule and achieves the highest extent in its middle part. We conclude that the highly conserved Gly-126, acting in concert with the previously identified non-canonical Asp-137, destabilizes the middle part of Tm, resulting in a more flexible region that is important for Tm function.

Analyzing protein folding cooperativity by differential scanning calorimetry and NMR spectroscopy

Some marginally stable proteins undergo microsecond time scale folding reactions that involve significant populations of partly ordered forms, making it difficult to discern individual steps in their folding pathways. It has been suggested that many of these proteins fold non-cooperatively, with no significant barriers to separate the energy landscape into distinct thermodynamic states. Here we present an approach for studying the cooperativity of rapid protein folding with a combination of differential scanning calorimetry (DSC), nuclear magnetic resonance (NMR) relaxation dispersion experiments, and an analysis of the temperature dependence of amide (1)H and (15)N chemical shifts. We applied this method to the PBX homeodomain (PBX-HD), which folds on the microsecond time scale and produces a broad DSC thermogram with an elevated and steeply sloping native-state heat capacity baseline, making it a candidate for barrierless folding. However, by globally fitting the NMR thermal melt and DSC data, and by comparing these results to those obtained from the NMR relaxation dispersion experiments, we show that the native form of the protein undergoes two-state exchange with a small population of the thermally denatured form, well below the melting temperature. This result directly demonstrates the coexistence of distinct folded and unfolded forms and firmly establishes that folding of PBX-HD is cooperative. Further, we see evidence of large-scale structural and dynamical changes within the native state by NMR, which helps to explain the broad and shallow DSC profile. This study illustrates the potential of combining calorimetry with NMR dynamics experiments to dissect mechanisms of protein folding.

Critical interactions in the stability control region of tropomyosin

Our laboratory has recently described a stability control region in the two-stranded alpha-helical coiled-coil alpha-tropomyosin that accounts for overall protein stability but is not required for folding (Hodges et al., 2009). We have used a synthetic peptide approach to investigate three stability control sites within the stability control region (residues 97-118). Two of the sites, electrostatic cluster 1 (97-104, EELDRAQE) and electrostatic cluster 2 (112-118, KLEEAEK), feature sequences with unusually high charge density and the potential to form multiple intrachain and interchain salt bridges (ionic attractions). A third site (105-111, RLATALQ) features an e position Leu residue, an arrangement known previously to enhance coiled-coil stability modestly. A native peptide and seven peptide analogs of the tropomyosin sequence 85-119 were prepared by Fmoc solid-phase peptide synthesis. Thermal stability measurements by circular dichroism (CD) spectroscopy revealed the following T(m) values for the native peptide and three key analogs: 52.9 degrees C (Native), 46.0 degrees C (R101A), 45.3 degrees C (K112A/K118A), and 27.9 degrees C (L110A). The corresponding DeltaT(m) values for the analogs, relative to the native peptide, are -6.9 degrees C, -7.6 degrees C, and -25.0 degrees C, respectively. The dramatic contribution to stability made by L110e is three times greater than the contribution of either electrostatic cluster 1 or 2, likely resulting from a novel hydrophobic interaction not previously observed. These thermal stability results were corroborated by temperature profiling analyses using reversed-phase high-performance liquid chromatography (RP-HPLC). We believe that the combined contributions of the interactions within the three stability control sites are responsible for the effect of the stability control region in tropomyosin, with the Leu110e contribution being most critical.

The importance of being flexible: the case of basic region leucine zipper transcriptional regulators

Large volumes of protein sequence and structure data acquired by proteomic studies led to the development of computational bioinformatic techniques that made possible the functional annotation and structural characterization of proteins based on their primary structure. It has become evident from genome-wide analyses that many proteins in eukaryotic cells are either completely disordered or contain long unstructured regions that are crucial for their biological functions. The content of disorder increases with evolution indicating a possibly important role of disorder in the regulation of cellular systems. Transcription factors are no exception and several proteins of this class have recently been characterized as premolten/molten globules. Yet, mammalian cells rely on these proteins to control expression of their 30,000 or so genes. Basic region:leucine zipper (bZIP) DNA-binding proteins constitute a major class of eukaryotic transcriptional regulators. This review discusses how conformational flexibility "built" into the amino acid sequence allows bZIP proteins to interact with a large number of diverse molecular partners and to accomplish their manifold cellular tasks in a strictly regulated and coordinated manner.

Microcalorimetry of proteins and their complexes

Ultrasensitive microcalorimetric techniques for measuring the heat capacities of proteins in dilute solutions over a broad temperature range (DSC) and the heats of protein reactions at fixed temperatures (ITC) are described and the methods of working with these instruments are considered. Particular attention is paid to analyzing the thermal properties of individual proteins, their stability, the energetics of their folding, and their association with specific macromolecular partners. Use of these calorimetric methods is illustrated with examples of small compact globular proteins, small proteins having loose noncompact structure, multidomain proteins, and protein complexes, particularly with DNA.

Microsecond melting of a folding intermediate in a coiled-coil peptide, monitored by T-jump/UV Raman spectroscopy

A truncated version of the GCN4 coiled-coil peptide has been studied by ultraviolet resonance Raman (UVRR) spectroscopy with 197 nm excitation, where amide modes are optimally enhanced. Although the CD melting curve could be satisfactorily described with a two-state transition having Tm = 30 degrees C, singular value decomposition of the UVRR data yielded three principal components, whose temperature dependence implicates an intermediate form between the folded and unfolded forms, with formation and melting temperatures of 10 and 40 degrees C. Two alpha-helical amide III bands, at 1340 and 1300 cm(-1), melted out selectively at 10 and 40 degrees C, respectively, and are assigned to hydrated and unhydrated helical regions. The hydrated regions are proposed to be melted in the intermediate form, while the unhydrated regions are intact. Time-resolved UVRR spectra following laser-induced temperature jumps revealed two relaxations, with time constants of 0.2 and 15 mus. These are suggested to reflect the melting times of hydrated and unhydrated helices. The unhydrated helical region may be associated with a 14-residue "trigger" sequence that has been identified in the C-terminal half of GCN4. Dehydration of helices may be a key step in the folding of coiled-coils.