Regular Secondary Structure Elements Research Articles

Domino Effect in Allosteric Signaling of Peptide Binding

While being a thoroughly studied model of dynamic allostery in a small protein, the pathway of signal transduction in the PDZ3 domain has not been fully determined. Here, we investigate peptide binding to the PDZ3 domain by conventional and fully data-driven analyses of molecular dynamics simulations. First, we identify isoleucine 37 as a key residue by widely used computational procedures such as cross-correlation and community network analysis. Simulations of the Ile37Ala mutant show disruption of the coordinated movements of spatially close regular elements of secondary structure. Then, we employ a recently developed unsupervised, data-driven procedure to determine an optimized reaction coordinate (slowest-relaxation eigenvector) of peptide binding. We use this reaction coordinate to improve sampling by restarting additional simulations from the transition state region. Significant differences in the distributions of some of the pairwise residue distances in the bound and unbound states emerge from the projection onto the optimized reaction coordinate. The unsupervised analysis shows that allosteric signaling is transduced from the β2 strand, which forms part of the peptide binding site, to the spatially adjacent β3 and β4 strands, and from there to the α3 helix. The domino-like transmission of a (peptide binding) signal along β strands and α helices that are close in three-dimensional space is likely to be a general mechanism of allostery in single-domain proteins.

MrParse: finding homologues in the PDB and the EBI AlphaFold database for molecular replacement and more.

Crystallographers have an array of search-model options for structure solution by molecular replacement (MR). The well established options of homologous experimental structures and regular secondary-structure elements or motifs are increasingly supplemented by computational modelling. Such modelling may becarried out locally or may use pre-calculated predictions retrieved from databases such as the EBI AlphaFold database. MrParse is a new pipeline to help to streamline the decision process in MR by consolidating bioinformatic predictions in one place. When reflection data are provided, MrParse can rank any experimental homologues found using eLLG, which indicates the likelihood that a given search model will work in MR. Inbuilt displays of predicted secondary structure, coiled-coil and transmembrane regions further inform the choice of MR protocol. MrParse can also identify and rank homologues in the EBI AlphaFold database, a function that will also interest other structural biologists and bioinformaticians.

NMR assignments and secondary structure distribution of emfourin, a novel proteinaceous protease inhibitor.

Emfourin (M4in) from Serratia proteamaculans is a new proteinaceous inhibitor of protealysin-like proteases (PLPs), a subgroup of the well-known and widely represented metallopeptidase M4 family. Although the biological role of PLPs is debatable, data published indicate their involvement in pathogenesis, including bacterial invasion into eukaryotic cells, suppression of immune defense of some animals, and destruction of plant cell walls. Gene colocalization into a bicistronic operon observed for some PLPs and their inhibitors (as in the case of M4in) implies a mutually consistent functioning of both entities. The originality of the amino acid sequence of M4in suggests it belongs to a previously unknown protein family and this encourages structural studies. In this work, we report a near-complete assignment of 1H, 13C, and 15N resonances of recombinant M4in and its structural-dynamic properties derived from the chemical shifts. According the NMR data analysis, the M4in molecule comprises 3-5 helical elements and 4-6 β-strands, at least two of which are apparently antiparallel, ascribing this obviously globular protein to the α + β structural class. Besides, two disordered regions also exist in the central loops between the regular secondary structural elements. The obtained data provide the basis for determining the high-resolution structure as well as functioning mechanism of M4in that can be used for development of new antibacterial therapeutic strategies.

Protein secondary structure motifs: A kinematic construction.

The kinematic geometry of protein backbone structures, constrained by either single or multiple hydrogen bonds (H-bonds), possibly in a periodic array, is discussed. These structures include regular secondary structure elements α-helices and β-sheets but also include other short H-bond stabilized irregular structural elements like β-turns. The work here shows that the variations observed in such structures have simple geometrical correlations consistent with constrained motion kinematics. A new classification of the ideal helices is given, in terms of the parameter α, the angle at a Cα atom to its two neighboring Cα 's along the helix, and shown how it can be generalized to include nonideal helices. Specifically, we derive an analytical expression of the backbone dihedrals, (ϕ, ψ), in terms of the parameter α subject to the constraint that the peptide planes are parallel to the helical axis. Helices constructed in this way exhibit near-vertical alignment of the C = O and N - H units and are the canonical objects of this study. These expressions are easily modifiable to include perturbations of parameters relevant to nonplanar peptide units and noncanonical angles. The addition of a second parameter, ε0 , inclination of successive peptide planes along a helix with respect to the helical axis leads to a generalization of the previous expression and provides an efficient parametrization of such structures in terms of coordinates consistent with H-bond parameters. An analogs parametrization of β-turns, using inverse kinematic methods, is also given. Besides offering a unifying viewpoint, our results may find useful applications to protein and peptide design.

Interactions of p53 with poly(ADP-ribose) and DNA induce distinct changes in protein structure as revealed by ATR-FTIR spectroscopy.

Due to multiple domains and in part intrinsically disordered regions, structural analyses of p53 remain a challenging task, particularly in complex with DNA and other macromolecules. Here, we applied a novel attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopic approach to investigate changes in secondary structure of full-length p53 induced by non-covalent interactions with DNA and poly(ADP-ribose) (PAR). To validate our approach, we confirmed a positive regulatory function of p53’s C-terminal domain (CTD) with regard to sequence-specific DNA binding and verified that the CTD mediates p53–PAR interaction. Further, we demonstrate that DNA and PAR interactions result in distinct structural changes of p53, indicating specific binding mechanisms via different domains. A time-dependent analysis of the interplay of DNA and PAR binding to p53 revealed that PAR represents p53’s preferred binding partner, which efficiently controls p53–DNA interaction. Moreover, we provide infrared spectroscopic data on PAR pointing to the absence of regular secondary structural elements. Finally, temperature-induced melting experiments via CD spectroscopy show that DNA binding stabilizes the structure of p53, while PAR binding can shift the irreversible formation of insoluble p53 aggregates to higher temperatures. In conclusion, this study provides detailed insights into the dynamic interplay of p53 binding to DNA and PAR at a formerly inaccessible molecular level.

\u03b2-Ginkgotides: Hyperdisulfide-constrained peptides from Ginkgo biloba

Hyperdisulfide-constrained peptides are distinguished by their high stability and diverse functions. Thus far, these peptides have been reported from animals only but their occurrence in plants are rare. Here, we report the discovery, synthesis and characterization of a hyperdisulfide-constrained peptides family of approximately 2 kDa, β-ginkgotides (β-gB1 and β-gB2) from Ginkgo biloba. Proteomic analysis showed β-ginkgotides contain 18‒20 amino acids, of which 16 residues form a conserved six-cysteine core with a highly clustered cysteine spacing of C‒CC‒C‒CC, an arrangement that has not been reported in cysteine-rich peptides. Disulfide mapping revealed a novel disulfide connectivity of CysI‒IV, CysII‒VI and CysIII‒V. Oxidative folding of synthetic β-gB1 to the native form was obtained in 70% yield. The synthetic β-gB1 displays a compact structure with no regular secondary structural elements, as determined by NMR spectroscopy. Transcriptomic analysis showed precursor βgb1 has a four-domain architecture and revealed an additional 76 β-ginkgotide-like peptides in 59 different gymnosperms, but none in angiosperms. Phylogenetic clustering analysis demonstrated β-ginkgotides belong to a new cysteine-rich peptide family. β-Ginkgotide is resistant to thermal, chemical and proteolytic degradation. Together, β-ginkgotides represent the first-in-class hyperdisulfide-constrained peptide family from plants with a novel scaffold that could be useful for engineering metabolically stable peptidyl therapeutics.

Protein Domain Mimics as Modulators of Protein-Protein Interactions.

Protein-protein interactions (PPIs) are ubiquitous in biological systems and often misregulated in disease. As such, specific PPI modulators are desirable to unravel complex PPI pathways and expand the number of druggable targets available for therapeutic intervention. However, the large size and relative flatness of PPI interfaces make them challenging molecular targets. This Account describes our systematic approach using secondary and tertiary protein domain mimics (PDMs) to specifically modulate PPIs. Our strategy focuses on mimicry of regular secondary and tertiary structure elements from one of the PPI partners to inspire rational PDM design. We have compiled three databases (HIPPDB, SIPPDB, and DIPPDB) of secondary and tertiary structures at PPI interfaces to guide our designs and better understand the energetics of PPI secondary and tertiary structures. Our efforts have focused on three of the most common secondary and tertiary structures: α-helices, β-strands, and helix dimers (e.g., coiled coils). To mimic α-helices, we designed the hydrogen bond surrogate (HBS) as an isosteric PDM and the oligooxopiperazine helix mimetic (OHM) as a topographical PDM. The nucleus of the HBS approach is a peptide macrocycle in which the N-terminal i, i + 4 main-chain hydrogen bond is replaced with a covalent carbon-carbon bond. In mimicking a main-chain hydrogen bond, the HBS approach stabilizes the α-helical conformation while leaving all helical faces available for functionalization to tune binding affinity and specificity. The OHM approach, in contrast, envisions a tetrapeptide to mimic one face of a two-turn helix. We anticipated that placement of ethylene bridges between adjacent amides constrains the tetrapeptide backbone to mimic the i, i + 4, and i + 7 side chains on one face of an α-helix. For β-strands, we developed triazolamers, a topographical PDM where the peptide bonds are replaced by triazoles. The triazoles simultaneously stabilize the extended, zigzag conformation of β-strands and transform an otherwise ideal protease substrate into a stable molecule by replacement of the peptide bonds. We turned to a salt bridge surrogate (SBS) approach as a means for stabilizing very short helix dimers. As with the HBS approach, the SBS strategy replaces a noncovalent interaction with a covalent bond. Specifically, we used a bis-triazole linkage that mimics a salt bridge interaction to drive helix association and folding. Using this approach, we were able to stabilize helix dimers that are less than half of the length required to form a coiled coil from two independent strands. In addition to demonstrating the stabilization of desired structures, we have also shown that our designed PDMs specifically modulate target PPIs in vitro and in vivo. Examples of PPIs successfully targeted include HIF1α/p300, p53/MDM2, Bcl-xL/Bak, Ras/Sos, and HIV gp41. The PPI databases and designed PDMs created in these studies will aid development of a versatile set of molecules to probe complex PPI functions and, potentially, PPI-based therapeutics.

Predictive approach for protein aggregation: Correlation of protein surface characteristics and conformational flexibility to protein aggregation propensity.

The aggregation of proteins became one of the major challenges in the development of biopharmaceu-ticals since the formation of aggregates can affect drug quality and immunogenicity. However, aggregation mechanisms are highly complex and the investigation requires cost, time, and material intensive experi-mental effort. In the present work, the predictive power of protein characteristics for the phase behavior of three different proteins which are very similar in size and structure was studied. In particular, the surface hydrophobicity, zeta potential, and conformational flexibility of human lysozyme, lysozyme from chicken egg white, and α-lactalbumin at pH 3, 5, 7, and 9 were assessed and examined for correlation with experimental stability studies focusing on protein phase behavior induced by sodium chloride and ammonium sulfate. The molecular dynamics (MD) simulation based study of the conformational flexibility without precipitants was able to identify highly flexible protein regions which could be associated to the less regular secondary structure elements and random coiled and terminal regions in particular. Conformational flex-ibility of the entire protein structure and protein surface hydrophobicity could be correlated to differing aggregation propensities among the studied proteins and could be identified to be applicable for predic-tion of protein phase behavior in aqueous solution without precipitants. For prediction of protein phase behavior and aggregation propensity in aqueous solution with precipitants, protein flexibility was further studied in dependency of salt concentration and species by means of human lysozyme. Even though the results of the salt dependent MD simulations could not be shown to be sufficient for prediction of salt depending phase behavior, this study revealed a more pronounced destabilizing effect of ammonium sulfate in comparison to sodium chloride and thus, was found to be in good agreement with theoretical considerations along the Hofmeister series as well as experimentally evaluated phase behavior. Biotechnol. Bioeng. 2017;114: 1170-1183. © 2016 Wiley Periodicals, Inc.

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.

Deducing the functional characteristics of the human selenoprotein SELK from the structural properties of its intrinsically disordered C-terminal domain.

The intrinsically disordered proteins (IDPs) cannot be described by a single structural representation but, due to their high structural fluctuation, through conformational ensembles. Certainly, molecular dynamics (MD) simulations represent a useful tool to study their different conformations capturing the conformational distribution. Our group is focusing on the structural characterization of proteins belonging to the seleno-proteome due to their involvement in cancer. They present disordered domains central for their biological function, and, in particular, SELK is a single-pass transmembrane protein that resides in the endoplasmic reticulum membrane (ER) with a C-terminal domain exposed to the cytoplasm that is known to interact with different components of the endoplasmic reticulum associated to the protein degradation (ERAD) pathway. This protein is found to be up-expressed in hepatocellular carcinoma and in other cancers. In this work we performed a detailed analysis of the C-terminal domain sequence of SELK and discovered that it is characterized by many prolines, and four negatively and eleven positively charged residues, which are crucial for its biological activity. This region can be considered as a weak polyelectrolyte and, specifically, a polycation, with high disordered propensity and different phosphorylation sites dislocated along the sequence. Then, we modeled its three-dimensional structure by performing MD simulations in water at neutral pH to analyze the structural stability as well as to identify the presence of HUB residues that play a key structural role as evidenced by the residue-residue interaction network analysis. Through this approach, we demonstrate that the C-terminal domain of SELK (i) presents a poor content of regular secondary structure elements, (ii) is dynamically stabilized by a network of intra-molecular H-bonds and H-bonds with water molecules, (iii) is highly fluctuating and, therefore, can be described only through a conformational ensemble, where we evidenced a distribution of equilibrium conformers which continuously inter-change their conformations. Finally to verify the specific role played by the negative charges, we also performed MD simulations at acidic pH. Overall, all the obtained results evidenced that SELK has the dynamic structural features to be defined as a HUB protein able to interact with multiple members. Therefore, considering the possible role that this protein can have in cancer development and progression, it can represent a target for drug design studies.

Integration of QUARK and I-TASSER for Ab Initio Protein Structure Prediction in CASP11.

We tested two pipelines developed for template-free protein structure prediction in the CASP11 experiment. First, the QUARK pipeline constructs structure models by reassembling fragments of continuously distributed lengths excised from unrelated proteins. Five free-modeling (FM) targets have the model successfully constructed by QUARK with a TM-score above 0.4, including the first model of T0837-D1, which has a TM-score = 0.736 and RMSD = 2.9 Å to the native. Detailed analysis showed that the success is partly attributed to the high-resolution contact map prediction derived from fragment-based distance-profiles, which are mainly located between regular secondary structure elements and loops/turns and help guide the orientation of secondary structure assembly. In the Zhang-Server pipeline, weakly scoring threading templates are re-ordered by the structural similarity to the ab initio folding models, which are then reassembled by I-TASSER based structure assembly simulations; 60% more domains with length up to 204 residues, compared to the QUARK pipeline, were successfully modeled by the I-TASSER pipeline with a TM-score above 0.4. The robustness of the I-TASSER pipeline can stem from the composite fragment-assembly simulations that combine structures from both ab initio folding and threading template refinements. Despite the promising cases, challenges still exist in long-range beta-strand folding, domain parsing, and the uncertainty of secondary structure prediction; the latter of which was found to affect nearly all aspects of FM structure predictions, from fragment identification, target classification, structure assembly, to final model selection. Significant efforts are needed to solve these problems before real progress on FM could be made. Proteins 2016; 84(Suppl 1):76-86. © 2015 Wiley Periodicals, Inc.

Prediction of temperature factors from protein sequence

Protein flexibility is useful in structural and functional aspect of proteins. We have analyzed the local primary protein sequence features that in combination can predict the B-value of amino acid residues directly from the protein sequence. We have also analyzed the distribution of B-value in different regions of protein three dimensional structures. On an average, the normalized Bvalue decreases by 0.1055 with every 0.5Å increase in the distance of the residue from protein surface. The residues in the loop regions have higher B-values as compared to the residues present in other regular secondary structural elements. Buried residues which are present in the protein core are more rigid (lower B-values) than the residues present on the protein surface. Similarly, the hydrophobic residues which tend to be present in the protein core have lower average B-value than the polar residues. Finally, we have proposed the method based on Support Vector Regression (SVR) to predict the B-value from protein primary sequence. Our result shows that, the SVR model achieved the correlation coefficient of 0.47 which is comparable to existing methods.

Use of (113)Cd NMR to probe the native metal binding sites in metalloproteins: an overview.

Our laboratories have actively published in this area for several years and the objective of this chapter is to present as comprehensive an overview as possible. Following a brief review of the basic principles associated with (113)Cd NMR methods, we will present the results from a thorough literature search for (113)Cd chemical shifts from metalloproteins. The updated (113)Cd chemical shift figure in this chapter will further illustrate the excellent correlation of the (113)Cd chemical shift with the nature of the coordinating ligands (N, O, S) and coordination number/geometry, reaffirming how this method can be used not only to identify the nature of the protein ligands in uncharacterized cases but also the dynamics at the metal binding site. Specific examples will be drawn from studies on alkaline phosphatase, Ca(2+) binding proteins, and metallothioneins.In the case of Escherichia coli alkaline phosphatase, a dimeric zinc metalloenzyme where a total of six metal ions (three per monomer) are involved directly or indirectly in providing the enzyme with maximal catalytic activity and structural stability, (113)Cd NMR, in conjunction with (13)C and (31)P NMR methods, were instrumental in separating out the function of each class of metal binding sites. Perhaps most importantly, these studies revealed the chemical basis for negative cooperativity that had been reported for this enzyme under metal deficient conditions. Also noteworthy was the fact that these NMR studies preceded the availability of the X-ray crystal structure.In the case of the calcium binding proteins, we will focus on two proteins: calbindin D(9k) and calmodulin. For calbindin D(9k) and its mutants, (113)Cd NMR has been useful both to follow actual changes in the metal binding sites and the cooperativity in the metal binding. Ligand binding to calmodulin has been studied extensively with (113)Cd NMR showing that the metal binding sites are not directly involved in the ligand binding. The (113)Cd chemical shifts are, however, exquisitely sensitive to minute changes in the metal ion environment.In the case of metallothionein, we will reflect upon how (113)Cd substitution and the establishment of specific Cd to Cys residue connectivity by proton-detected heteronuclear (1)H-(113)Cd multiple-quantum coherence methods (HMQC) was essential for the initial establishment of the 3D structure of metallothioneins, a protein family deficient in the regular secondary structural elements of α-helix and β-sheet and the first native protein identified with bound Cd. The (113)Cd NMR studies also enabled the characterization of the affinity of the individual sites for (113)Cd and, in competition experiments, for other divalent metal ions: Zn, Cu, and Hg.

The Interconversion between a Flexible β-Sheet and a Fibril β-Arrangement Constitutes the Main Conformational Event during Misfolding of PSD95-PDZ3 Domain

The temperature-induced misfolding pathway of PDZ3, the third PDZ domain of the PSD95 neuronal protein, is populated by a trimeric β-sheet-rich intermediate state that leads to a stepwise and reversible formation of supramacromolecular structures. Using FTIR, we have found that misfolding of this pathway is not due to different ensembles of a variety of precursors, but comes mainly from the interconversion of a flexible β-sheet of the domain to wormlike fibrils. The appearance of the wormlike fibril FTIR component is also accompanied by a slight decrease of the band that corresponds to loops in the native state, whereas the rest of the regular elements of secondary structure are fairly well maintained upon misfolding. Transmission electron microscope micrographs have confirmed the presence of wormlike fibrils upon heating at 60°C, where the trimeric intermediate is maximally populated. Toxicity assays in the human neuroblastoma cell line SH-SY5Y show that cytotoxicity increases as the aggregation pathway proceeds. NMR analysis of chemical shifts as a function of temperature has revealed, as one of the main conformational aspects of such an interconversion at the residue level, that the β-sheet arrangement around strand β3 promotes the change that drives misfolding of the PDZ3 domain.

Quantifying protein modularity and evolvability: A comparison of different techniques

Modularity increases evolvability by reducing constraints on adaptation and by allowing preexisting parts to function in new contexts for novel uses. Protein evolution provides an excellent context to study the causes and consequences of biological modularity. In order to address such questions, however, an index for protein modularity is necessary. This paper proposes a simple index for protein modularity—“module density”—which is the number of evolutionarily independent modules that compose a protein divided by the number of amino acids in the protein. The decomposition of proteins into constituent modules can be accomplished by either of two classes of methods. The first class of methods relies on “suppositional” criteria to assign amino acids to modules, whereas the second class of methods relies on “coevolutionary” criteria for this task. One simple and practical method from the first class consists of approximating the number of modules in a protein as the number of regular secondary structure elements (i.e., helices and sheets). Methods based on coevolutionary criteria require more elaborate data, but they have the advantage of being able to specify modules without prior assumptions about why they exist. Given the increasing availability of datasets sampling protein mutational spectra (e.g., from comparative genomics, experimental evolution, and computational prediction), methods based on coevolutionary criteria will likely become more promising in the near future. The ability to meaningfully quantify protein modularity via simple indices has the potential to aid future efforts to understand protein evolutionary rate determinants, improve molecular evolution models and engineer novel proteins.

Terminal Amino Acids Disturb Xylanase Thermostability and Activity

Protein structure is composed of regular secondary structural elements (α-helix and β-strand) and non-regular region. Unlike the helix and strand, the non-regular region consists of an amino acid defined as a disordered residue (DR). When compared with the effect of the helix and strand, the effect of the DR on enzyme structure and function is elusive. An Aspergillus niger GH10 xylanase (Xyn) was selected as a model molecule of (β/α)(8) because the general structure consists of ~10% enzymes. The Xyn has five N-terminal DRs and one C-terminal DR, respectively, which were deleted to construct three mutants, XynΔN, XynΔC, and XynΔNC. Each mutant was ~2-, 3-, or 4-fold more thermostable and 7-, 4-, or 4-fold more active than the Xyn. The N-terminal deletion decreased the xylanase temperature optimum for activity (T(opt)) 6 °C, but the C-terminal deletion increased its T(opt) 6 °C. The N- and C-terminal deletions had opposing effects on the enzyme T(opt) but had additive effects on its thermostability. The five N-terminal DR deletions had more effect on the enzyme kinetics but less effect on its thermo property than the one C-terminal DR deletion. CD data showed that the terminal DR deletions increased regular secondary structural contents, and hence, led to slow decreased Gibbs free energy changes (ΔG(0)) in the thermal denaturation process, which ultimately enhanced enzyme thermostabilities.

Measurement of amide hydrogen exchange rates with the use of radiation damping

A simple method for measuring amide hydrogen exchange rates is presented, which is based on the selective inversion of water magnetization with the use of radiation damping. Simulations show that accurate exchange rates can be measured despite the complications of radiation damping and cross relaxation to the exchange process between amide and water protons. This method cannot eliminate the contributions of the exchange-relayed NOE and direct NOE to the measured exchange rates, but minimize the direct NOE contribution. In addition, the amides with a significant amount of such indirect contributions are possible to be identified from the shape of the exchange peak intensity profiles or/and from the apparent relaxation rates of amide protons which are extracted from fitting the intensity profiles to an equation established here for our experiment. The method was tested on ubiquitin and also applied to an acyl carrier protein. The amide exchange rates for the acyl carrier protein at two pHs indicate that the entire protein is highly dynamic on the second timescale. Low protection factors for the residues in the regular secondary structural elements also suggest the presence of invisible unfolded species. The highly dynamic nature of the acyl carrier protein may be crucial for its interactions with its substrate and enzymes.

Recognition and Specificity Determinants of the Human Cbx Chromodomains

The eight mammalian Cbx proteins are chromodomain-containing proteins involved in regulation of heterochromatin, gene expression, and developmental programs. They are evolutionarily related to the Drosophila HP1 (dHP1) and Pc (dPc) proteins that are key components of chromatin-associated complexes capable of recognizing repressive marks such as trimethylated Lys-9 and Lys-27, respectively, on histone H3. However, the binding specificity and function of the human homologs, Cbx1–8, remain unclear. To this end we employed structural, biophysical, and mutagenic approaches to characterize the molecular determinants of sequence contextual methyllysine binding to human Cbx1–8 proteins. Although all three human HP1 homologs (Cbx1, -3, -5) replicate the structural and binding features of their dHP counterparts, the five Pc homologs (Cbx2, -4, -6, -7, -8) bind with lower affinity to H3K9me3 or H3K27me3 peptides and are unable to distinguish between these two marks. Additionally, peptide permutation arrays revealed a greater sequence tolerance within the Pc family and suggest alternative nonhistone sequences as potential binding targets for this class of chromodomains. Our structures explain the divergence of peptide binding selectivity in the Pc subfamily and highlight previously unrecognized features of the chromodomain that influence binding and specificity.

Protein Structural Modularity and Robustness Are Associated with Evolvability

Theory suggests that biological modularity and robustness allow for maintenance of fitness under mutational change, and when this change is adaptive, for evolvability. Empirical demonstrations that these traits promote evolvability in nature remain scant however. This is in part because modularity, robustness, and evolvability are difficult to define and measure in real biological systems. Here, we address whether structural modularity and/or robustness confer evolvability at the level of proteins by looking for associations between indices of protein structural modularity, structural robustness, and evolvability. We propose a novel index for protein structural modularity: the number of regular secondary structure elements (helices and strands) divided by the number of residues in the structure. We index protein evolvability as the proportion of sites with evidence of being under positive selection multiplied by the average rate of adaptive evolution at these sites, and we measure this as an average over a phylogeny of 25 mammalian species. We use contact density as an index of protein designability, and thus, structural robustness. We find that protein evolvability is positively associated with structural modularity as well as structural robustness and that the effect of structural modularity on evolvability is independent of the structural robustness index. We interpret these associations to be the result of reduced constraints on amino acid substitutions in highly modular and robust protein structures, which results in faster adaptation through natural selection.

High‐Resolution Crystal Structure of a Lasso Peptide

Lasso peptides are a group of a growing class of bioactive peptides of microbial origin. The first crystal structure of a member of this family, the glucagon receptor antagonist BI-32169, shows that the fold is built predominantly by regular secondary structural elements and a tight network of hydrogen bonds that are partially shielded from solvent by hydrophobic amino acid side chains. This results in an extraordinary stable structure that is resistant to thermal unfolding or proteolytic digestion, which facilitates its biological function.