Energetic determinants of oligomeric state specificity in coiled coils

Mechanism of negative membrane curvature generation by I-BAR domains

The membrane sculpting ability of BAR domains has been attributed to the intrinsic curvature of their banana-shaped dimeric structure. However, there is often a mismatch between this intrinsic curvature and the diameter of the membrane tubules generated. I-BAR domains are especially mysterious since they are almost flat but generate high negative membrane curvature. Here, we use atomistic implicit-solvent computer modeling to show that the membrane bending of the IRSp53 I-BAR domain is dictated by its higher oligomeric structure, whose curvature is completely unrelated to the intrinsic curvature of the dimer. Two other I-BARs give similar results, whereas a flat F-BAR sheet develops a concave membrane-binding interface, consistent with its observed positive membrane curvature generation. Laterally interacting helical spirals of I-BAR dimers on tube interiors are stable and have an enhanced binding energy that is sufficient for membrane bending to experimentally observed tubule diameters at a reasonable surface density.

Derivatives of the lectin complement pathway in Lophotrochozoa

A plethora of non-overlapping immune molecular mechanisms in metazoans is the most puzzling issue in comparative immunobiology. No valid evolutionary retrospective on these mechanisms has been developed. In this study, we aimed to reveal the origin and evolution of the immune complement-like system in Lophotrochozoa. For this, we analyzed publicly available transcriptomes of prebilaterian and lophotrochozoan species, mapping lineage-specific molecular events on the phylogenetic tree. We found that there were no orthologs of mannose-binding lectin (MBL) and ficolins (FCN) in Lophotrochozoa but C1q-like proteins (C1qL), bearing both a collagen domain and a globular C1q domain, were omnipresent in them. This suggests that among all complement-like activators the C1qL-specific domain architecture was an evolutionarily first. Two novel protostomian MASP-Related Molecules, MReM1 and MReM2, might hypothetically compensate for the loss of a prebilaterian MASP-orthologous gene and act in complex with C1qL and C1qDC as a "proto-activator" of an ancient "proto-complement". We proposed a new model of the complement evolution predicting that numerous lineage-specific complement-like systems should have evolved from a stem "antique" molecular complex. First evolved in the common ancestor of coelomic animals, the "antique" humoral complex consisted of a TEP molecule, the common ancestor of TEP-associated proteases (C2/Bf/Сf/Lf), the common ancestor of MASP-like proteases (MASP/C1r/C1s, MReM1/MReM2) and multimeric recognition proteins (C1q-, MBL- and FCN-homologs). Further evolutionary specialization and expansion of the complex was independent and lineage-specific, examples being the mammalian complement system and the Apogastropoda complement-like complex. The latter includes an impressive array of multimeric recognition proteins, the variable immunoglobulin and lectin domain containing molecules (VIgL), homologous to C1q, MBL, FCN and other lectins. Four novel polymorphic subfamilies of VIgLs were found to be expressed in Apogastropoda: C1q-related proteins (QREP), zona pellucida-related proteins (ZREP), Scavenger Receptor Cys-Rich-related proteins (SREP) and HPA-lectin related proteins (HREP). The transcriptional response of fibrinogen-related proteins of VIgL family (LlFREP), LlQREP and LlSREP to infestation of common periwinkle, Littorina littorea, with digenean parasite Himasthla elongata correlates with that of LlMReM1, supporting the model suggested in this study.

Self-assembled clusters of patchy rod-like molecules

The design of complex self-assembled structures remains a challenging task due to the intricate relationship between the properties of the building blocks and the final morphology of the aggregates. Here, we report such a relationship for rod-like particles with one or two attractive patches based on a combination of computer simulations and analytical theory. We investigated the formation of finite aggregates under various conditions and constructed structure diagrams, which can be used to determine and extrapolate the system composition. The size of the clusters is mainly determined by the size of the attractive patches and their geometrical arrangement. We showed that it is challenging to obtain clusters with more than four particles in high yields, and more complex building blocks or additional molecules would need to be used. Moreover, the particles with patch sizes close to the structure boundaries can switch between the aggregate state by a small change in conditions. These findings can be useful for the development of self-assembling building blocks and for the understanding of protein folds of coiled coils under various conditions.

Computational prediction of the optimal oligomeric state for membrane-inserted β-barrels of protegrin-1 and related mutants

Protegrin-1 is a widely studied 18-residue β-hairpin antimicrobial peptide. Evidence suggests that it acts via a β-barrel pore formation mechanism, but the exact number of peptides comprising the pore state is unknown. In this study, we performed molecular dynamics simulations of β-barrels of protegrin and three related mutants (v14v16l, v14v16a, and r4n) in NCNC parallel topology in implicit membrane pores of varying radius and curvature for oligomeric numbers 6-14. We then identified the optimal pore radius and curvature values for all constructs and determined the total effective energy and the translational and rotational entropic losses. These, along with an estimate of membrane deformation free energy from experimental line tension values, provided an estimate of the overall energetics of formation of each pore state. The results indicated that oligomeric numbers 7-13 are generally stable, allowing the possibility of a heterogeneous pore state. The optimal oligomeric state for protegrin is the nonamer, shifting to higher numbers for the mutants. Protegrin, v14v16l, and r4n are stable as membrane-inserted β-barrels, but v14v16a seems much less so because of its decreased hydrophobicity. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

Automatedde novophasing and model building of coiled-coil proteins

Models generated byde novostructure prediction can be very useful starting points for molecular replacement for systems where suitable structural homologues cannot be readily identified. Protein–protein complexes andde novo-designed proteins are examples of systems that can be challenging to phase. In this study, the potential ofde novomodels of protein complexes for use as starting points for molecular replacement is investigated. The approach is demonstrated using homomeric coiled-coil proteins, which are excellent model systems for oligomeric systems. Despite the stereotypical fold of coiled coils, initial phase estimation can be difficult and many structures have to be solved with experimental phasing. A method was developed for automatic structure determination of homomeric coiled coils from X-ray diffraction data. In a benchmark set of 24 coiled coils, ranging from dimers to pentamers with resolutions down to 2.5 Å, 22 systems were automatically solved, 11 of which had previously been solved by experimental phasing. The generated models contained 71–103% of the residues present in the deposited structures, had the correct sequence and had freeRvalues that deviated on average by 0.01 from those of the respective reference structures. The electron-density maps were of sufficient quality that only minor manual editing was necessary to produce final structures. The method, namedCCsolve, combines methods forde novostructure prediction, initial phase estimation and automated model building into one pipeline.CCsolveis robust against errors in the initial models and can readily be modified to make use of alternative crystallographic software. The results demonstrate the feasibility ofde novophasing of protein–protein complexes, an approach that could also be employed for other small systems beyond coiled coils.

Exploring alternate states and oligomerization preferences of coiled-coils by de novo structure modeling

Homomeric coiled-coils can self-assemble into a wide range of structural states with different helix topologies and oligomeric states. In this study, we have combined de novo structure modeling with stability calculations to simultaneously predict structure and oligomeric states of homomeric coiled-coils. For dimers an asymmetric modeling protocol was developed. Modeling without symmetry constraints showed that backbone asymmetry is important for the formation of parallel dimeric coiled-coils. Collectively, our results demonstrate that high-resolution structure of coiled-coils, as well as parallel and antiparallel orientations of dimers and tetramers, can be accurately predicted from sequence. De novo modeling was also used to generate models of competing oligomeric states, which were used to compare stabilities and thus predict the native stoichiometry from sequence. In a benchmark set of 33 coiled-coil sequences, forming dimers to pentamers, up to 70% of the oligomeric states could be correctly predicted. The calculations demonstrated that the free energy of helix folding could be an important factor for determining stability and oligomeric state of homomeric coiled-coils. The computational methods developed here should be broadly applicable to studies of sequence-structure relationships in coiled-coils and the design of higher order assemblies with improved oligomerization specificity.

H-NST induces LEE expression and the formation of attaching and effacing lesions in enterohemorrhagic Escherichia coli

Background

Enteropathogenic E. coli (EPEC) and enterohemorrhagic E. coli are important causes of morbidity and mortality worldwide. These enteric pathogens contain a type III secretion system (T3SS) responsible for the attaching and effacing (A/E) lesion phenotype. The T3SS is encoded by the locus of enterocyte effacement (LEE) pathogenicity island. The H-NS-mediated repression of LEE expression is counteracted by Ler, the major activator of virulence gene expression in A/E pathogens. A regulator present in EPEC, H-NST, positively affects expression of H-NS regulon members in E. coli K-12, although the effect of H-NST on LEE expression and virulence of A/E pathogens has yet-to-be determined.

Results

We examine the effect of H-NST on LEE expression and A/E lesion formation on intestinal epithelial cells. We find that H-NST positively affects the levels of LEE-encoded proteins independently of ler and induces A/E lesion formation. We demonstrate H-NST binding to regulatory regions of LEE1 and LEE3, the first report of DNA-binding by H-NST. We characterize H-NST mutants substituted at conserved residues including Ala16 and residues Arg60 and Arg63, which are part of a potential DNA-binding domain. The single mutants A16V, A16L, R60Q and the double mutant R60Q/R63Q exhibit a decreased effect on LEE expression and A/E lesion formation. DNA mobility shift assays reveal that these residues are important for H-NST to bind regulatory LEE DNA targets. H-NST positively affects Ler binding to LEE DNA in the presence of H-NS, and thereby potentially helps Ler displace H-NS bound to DNA.

Conclusions

H-NST induces LEE expression and A/E lesion formation likely by counteracting H-NS-mediated repression. We demonstrate that H-NST binds to DNA and identify arginine residues that are functionally important for DNA-binding. Our study suggests that H-NST provides an additional means for A/E pathogens to alleviate repression of virulence gene expression by H-NS to promote virulence capabilities.

A thermodynamic approach to alamethicin pore formation

The structure and energetics of alamethicin Rf30 monomer to nonamer in cylindrical pores of 5 to 11Å radius are investigated using molecular dynamics simulations in an implicit membrane model that includes the free energy cost of acyl chain hydrophobic area exposure. Stable, low energy pores are obtained for certain combinations of radius and oligomeric number. The trimer and the tetramer formed 6Å pores that appear closed while the larger oligomers formed open pores at their optimal radius. The hexamer in an 8Å pore and the octamer in an 11Å pore give the lowest effective energy per monomer. However, all oligomers beyond the pentamer have comparable energies, consistent with the observation of multiple conductance levels. The results are consistent with the widely accepted "barrel-stave" model. The N terminal portion of the molecule exhibits smaller tilt with respect to the membrane normal than the C terminal portion, resulting in a pore shape that is a hybrid between a funnel and an hourglass. Transmembrane voltage has little effect on the structure of the oligomers but enhances or decreases their stability depending on its orientation. Antiparallel bundles are lower in energy than the commonly accepted parallel ones and could be present under certain experimental conditions. Dry aggregates (without an aqueous pore) have lower average effective energy than the corresponding aggregates in a pore, suggesting that alamethicin pores may be excited states that are stabilized in part by voltage and in part by the ion flow itself.

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.

An amino acid packing code for α-helical structure and protein design

This work demonstrates that all packing in α-helices can be simplified to repetitive patterns of a single motif: the knob-socket. Using the precision of Voronoi Polyhedra/Delauney Tessellations to identify contacts, the knob-socket is a four-residue tetrahedral motif: a knob residue on one α-helix packs into the three-residue socket on another α-helix. The principle of the knob-socket model relates the packing between levels of protein structure: the intra-helical packing arrangements within secondary structure that permit inter-helix tertiary packing interactions. Within an α-helix, the three-residue sockets arrange residues into a uniform packing lattice. Inter-helix packing results from a definable pattern of interdigitated knob-socket motifs between two α-helices. Furthermore, the knob-socket model classifies three types of sockets: (1) free, favoring only intra-helical packing; (2) filled, favoring inter-helical interactions; and (3) non, disfavoring α-helical structure. The amino acid propensities in these three socket classes essentially represent an amino acid code for structure in α-helical packing. Using this code, we used a novel yet straightforward approach for the design of α-helical structure to validate the knob-socket model. Unique sequences for three peptides were created to produce a predicted amount of α-helical structure: mostly helical, some helical, and no helix. These three peptides were synthesized, and helical content was assessed using CD spectroscopy. The measured α-helicity of each peptide was consistent with the expected predictions. These results and analysis demonstrate that the knob-socket motif functions as the basic unit of packing and presents an intuitive tool to decipher the rules governing packing in protein structure.

Computational analysis of residue contributions to coiled-coil topology

A variety of features are thought to contribute to the oligomeric and topological specificity of coiled coils. In previous work, we examined the determinants of oligomeric state. Here, we examine the energetic basis for the tendency of six coiled-coil peptides to align their α-helices in antiparallel orientation using molecular dynamics simulations with implicit solvation (EEF1.1). We also examine the effect of mutations known to disrupt the topology of these peptides. In agreement with experiment, ARG or LYS at a or d positions were found to stabilize the antiparallel configuration. The modeling suggests that this is not due to a-a' or d-d' repulsions but due to interactions with e' and g' residues. TRP at core positions also favors the antiparallel configuration. Residues that disfavor parallel dimers, such as ILE at d, are better tolerated in, and thus favor the antiparallel configuration. Salt bridge networks were found to be more stabilizing in the antiparallel configuration for geometric reasons: antiparallel helices point amino acid side chains in opposite directions. However, the structure with the largest number of salt bridges was not always the most stable, due to desolvation and configurational entropy contributions. In tetramers, the extent of stabilization of the antiparallel topology by core residues is influenced by the e' residue on a neighboring helix. Residues at b and c positions in some cases also contribute to stabilization of antiparallel tetramers. This work provides useful rules toward the goal of designing coiled coils with a well-defined and predictable three-dimensional structure.

Multicoil2: predicting coiled coils and their oligomerization states from sequence in the twilight zone

The alpha-helical coiled coil can adopt a variety of topologies, among the most common of which are parallel and antiparallel dimers and trimers. We present Multicoil2, an algorithm that predicts both the location and oligomerization state (two versus three helices) of coiled coils in protein sequences. Multicoil2 combines the pairwise correlations of the previous Multicoil method with the flexibility of Hidden Markov Models (HMMs) in a Markov Random Field (MRF). The resulting algorithm integrates sequence features, including pairwise interactions, through multinomial logistic regression to devise an optimized scoring function for distinguishing dimer, trimer and non-coiled-coil oligomerization states; this scoring function is used to produce Markov Random Field potentials that incorporate pairwise correlations localized in sequence. Multicoil2 significantly improves both coiled-coil detection and dimer versus trimer state prediction over the original Multicoil algorithm retrained on a newly-constructed database of coiled-coil sequences. The new database, comprised of 2,105 sequences containing 124,088 residues, includes reliable structural annotations based on experimental data in the literature. Notably, the enhanced performance of Multicoil2 is evident when tested in stringent leave-family-out cross-validation on the new database, reflecting expected performance on challenging new prediction targets that have minimal sequence similarity to known coiled-coil families. The Multicoil2 program and training database are available for download from http://multicoil2.csail.mit.edu.

Computational design of intermolecular stability and specificity in protein self-assembly

The ability to engineer novel proteins using the principles of molecular structure and energetics is a stringent test of our basic understanding of how proteins fold and maintain structure. The design of protein self-assembly has the potential to impact many fields of biology from molecular recognition to cell signaling to biomaterials. Most progress in computational design of protein self-assembly has focused on α-helical systems, exploring ways to concurrently optimize the stability and specificity of a target state. Applying these methods to collagen self-assembly is very challenging, due to fundamental differences in folding and structure of α- versus triple-helices. Here, we explore various computational methods for designing stable and specific oligomeric systems, with a focus on α-helix and collagen self-assembly.

Computational studies of colicin insertion into membranes: the closed state

Colicins are water-soluble toxins that, upon interaction with membranes, undergo a conformational change, insert, and form pores in them. Pore formation activity is localized in a bundle of 10 α-helices named the pore-forming domain (PFD). There is evidence that colicins attach to the membrane via a hydrophobic hairpin embedded in the core of the PFD. Two main models have been suggested for the membrane-bound state: penknife and umbrella, differing in regard to the orientation of the hydrophobic hairpin with respect to the membrane. The arrangement of the amphipathic helices has been described as either a compact three-dimensional structure or a two-dimensional array of loosely interacting helices on the membrane surface. Using molecular dynamics simulations with an implicit membrane model, we studied the structure and stability of the conformations proposed earlier for four colicins. We find that colicins are initially driven towards the membrane by electrostatic interactions between basic residues and the negatively charged membrane surface. They do not have a unique binding orientation, but in the predominant orientations the central hydrophobic hairpin is parallel to the membrane. In the inserted state, the estimated free energy tends to be lower for the compact arrangements of the amphipathic helix, but the more expanded ones are in better agreement with experimental distance distributions. The difference in energy between penknife and umbrella conformations is small enough for equilibrium to exist between them. Elongation of the hydrophobic hairpin helices and membrane thinning were found unable to produce stabilization of the transmembrane configuration of the hydrophobic hairpin.

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.

Structural specificity in coiled-coil interactions

Coiled coils have a rich history in the field of protein design and engineering. Novel structures, such as the first seven-helix coiled coil, continue to provide surprises and insights. Large-scale datasets quantifying the influence of systematic mutations on coiled-coil stability are a valuable new asset to the area. Scoring methods based on sequence and/or structure can predict interaction preferences in coiled-coil-mediated bZIP transcription factor dimerization. Experimental and computational methods for dealing with the near-degeneracy of many coiled-coil structures appear promising for future design applications.

Variable stability heterodimeric coiled-coils from manipulation of electrostatic interface residue chain length

The design of variable-stability coiled-coil heterodimers is described. The electrostatic interface between helices, formed by contact between side chains in heptad e/g positions, is manipulated to produce complexes ranging in stability from ones that are essentially unstructured to those that cannot be thermally denatured. The tuning is accomplished by incremental extension or contraction of parent glutamic acid and lysine side chains by single methylene units, producing peptides that bear either carboxylic acids or amines separated from the peptide backbone by one to four CH2 groups. Detailed examination of all homodimers and electrostatically compatible heterodimers generates interesting combinations, particularly those in which longer-chain acids are incorporated into peptides paired with lysine-bearing ones. The discovery of very stable dimers allows exchange experiments in which one strand of an original heterodimer is specifically replaced by an added one, even in cases where the original complex features the native-like glutamic acid/lysine pairing. The reported results add significantly to the available design templates for coiled-coil construction and enable the future implementation of various triggered-recognition strategies.