An amino acid code to define a protein's tertiary packing surface

One difficult aspect of the protein‐folding problem is characterizing the nonspecific interactions that define packing in protein tertiary structure. To better understand tertiary structure, this work extends the knob‐socket model by classifying the interactions of a single knob residue packed into a set of contiguous sockets, or a pocket made up of 4 or more residues. The knob‐socket construct allows for a symbolic two‐dimensional mapping of pockets. The two‐dimensional mapping of pockets provides a simple method to investigate the variety of pocket shapes to understand the geometry of protein tertiary surfaces. The diversity of pocket geometries can be organized into groups of pockets that share a common core, which suggests that some interactions in pockets are ancillary to packing. Further analysis of pocket geometries displays a preferred configuration that is right‐handed in α‐helices and left‐handed in β‐sheets. The amino acid composition of pockets illustrates the importance of nonpolar amino acids in packing as well as position specificity. As expected, all pocket shapes prefer to pack with hydrophobic knobs; however, knobs are not selective for the pockets they pack. Investigating side‐chain rotamer preferences for certain pocket shapes uncovers no strong correlations. These findings allow a simple vocabulary based on knobs and sockets to describe protein tertiary packing that supports improved analysis, design, and prediction of protein structure. Proteins 2016; 84:201–216. © 2015 Wiley Periodicals, Inc.     

Molecular dynamics simulations of a calmodulin-peptide complex in solution

We have performed an 4-ns MD simulation of calmodulin complexed with a target peptide in explicit water, under realistic conditions of constant temperature and pressure, in the presence of a physiological concentration of counterions and using Ewald summation to avoid truncation of long-range electrostatic forces. During the simulation the system tended to perform small fluctuations around a structure similar to, but somewhat looser than the starting crystal structure. The calmodulin-peptide complex was quite rigid and did not exhibit any large amplitude domain motions such as previously seen in apo- and calcium-bound calmodulin. We analyzed the calmodulin-peptide interactions by calculating buried surface areas, CHARMM interaction energies and continuum model interaction free energies. In the trajectory, the protein surface area buried by contact with the peptide is 1373 A(2) approximately evenly divided between the calmodulin N-terminal, C-terminal and central linker regions. A majority of this buried surface, 803 A(2), comes from nonpolar residues, in contrast to the protein as a whole, for which the surface is made up of mostly polar and charged groups. Our continuum calculations indicate that the largest favorable contribution to peptide binding comes from burial of molecular surface upon complex formation. Electrostatic contributions are favorable but smaller in the trajectory structures, and actually unfavorable for binding in the crystal structure. Since nonpolar groups make up most of buried surface of the protein, our calculations suggest that the hydrophobic effect is the main driving force for binding the helical peptide to calmodulin, consistent with thermodynamic analysis of experimental data. Besides the burial of nonpolar surface area, secondary contributions to peptide binding come from burial of polar surface and electrostatic interactions. In the nonpolar interactions a crucial role is played by the nine methionines of calmodulin. In the electrostatic interactions the negatively charged protein residues and positively charged peptide residues play a dominant role.     

A fluorescent reporter for mapping cellular protein‐protein interactions in time and space

We introduce a fluorescent reporter for monitoring protein–protein interactions in living cells. The method is based on the Split‐Ubiquitin method and uses the ratio of two auto‐fluorescent reporter proteins as signal for interaction (SPLIFF). The mating of two haploid yeast cells initiates the analysis and the interactions are followed online by two‐channel time‐lapse microscopy of the diploid cells during their first cell cycle. Using this approach we could with high spatio‐temporal resolution visualize the differences between the interactions of the microtubule binding protein Stu2p with two of its binding partners, monitor the transient association of a Ran‐GTPase with its receptors at the nuclear pore, and distinguish between protein interactions at the polar cortical domain at different phases of polar growth. These examples further demonstrate that protein–protein interactions identified from large‐scale screens can be effectively followed up by high‐resolution single‐cell analysis.     

Dissection of the protein G B1 domain binding site for human IgG Fc fragment.

The contribution to the free energy of binding of each of the residues forming the binding site for a human IgG Fc fragment on the surface of the B1 domain of protein G was determined by alanine-scanning mutagenesis. The interface between these two proteins is atypical in that it is smaller than usual, polar in character, and involves two well-defined "knobs-into-holes" interactions. The bulk of the free energy of binding is contributed by three central residues, which make hydrogen bonds across the interface. Of these, the most critical interaction is formed by Glu27, which acts as a charged knob on the surface of the B1 domain, inserting into a polar hole on the Fc fragment. A single alanine mutation of this residue virtually abolishes stable complex formation. Formation of a stable interface between these two proteins is therefore dominated by a small, polar "hot spot."     

Molecular Dynamics Simulations of a Calmodulin-peptide Complex in Solution

Abstract

We have performed an 4-ns MD simulation of calmodulin complexed with a target peptide in explicit water, under realistic conditions of constant temperature and pressure, in the presence of a physiological concentration of counterions and using Ewald summation to avoid truncation of long-range electrostatic forces. During the simulation the system tended to perform small fluctuations around a structure similar to, but somewhat looser than the starting crystal structure. The calmodulin-peptide complex was quite rigid and did not exhibit any large amplitude domain motions such as previously seen in apo- and calcium-bound calmodulin. We analyzed the calmodulin-peptide interactions by calculating buried surface areas, CHARMM interaction energies and continuum model interaction free energies. In the trajectory, the protein surface area buried by contact with the peptide is 1373 Ų, approximately evenly divided between the calmodulin N-terminal, C-terminal and central linker regions. A majority of this buried surface, 803 ·A², comes from nonpolar residues, in contrast to the protein as a whole, for which the surface is made up of mostly polar and charged groups. Our continuum calculations indicate that the largest favorable contribution to pep- tide binding comes from burial of molecular surface upon complex formation. Electrostatic contributions are favorable but smaller in the trajectory structures, and actually unfavorable for binding in the crystal structure. Since nonpolar groups make up most of buried surface of the protein, our calculations suggest that the hydrophobic effect is the main driving force for binding the helical peptide to calmodulin, consistent with thermodynamic analysis of experimental data. Besides the burial of nonpolar surface area, secondary contributions to peptide binding come from burial of polar surface and electrostatic interactions. In the nonpolar interactions a crucial role is played by the nine methionines of calmodulin. In the electrostatic interactions the negatively charged protein residues and positively charged peptide residues play a dominant role.     

Carbohydrate binding module recognition of xyloglucan defined by polar contacts with branching xyloses and CH‐Π interactions

Engineering of novel carbohydrate‐binding proteins that can be utilized in various biochemical and biotechnical applications would benefit from a deeper understanding of the biochemical interactions that determine protein‐carbohydrate specificity. In an effort to understand further the basis for specificity we present the crystal structure of the multi‐specific carbohydrate‐binding module (CBM) X‐2 L110F bound to a branched oligomer of xyloglucan (XXXG). X‐2 L110F is an engineered CBM that can recognize xyloglucan, xylans and β‐glucans. The structural observations of the present study compared with previously reported structures of X‐2 L110F in complex with linear oligomers, show that the π‐surface of a phenylalanine, F110, allows for interactions with hydrogen atoms on both linear (xylopentaose and cellopentaose) and branched ligands (XXXG). Furthermore, X‐2 L110F is shown to have a relatively flexible binding cleft, as illustrated in binding to XXXG. This branched ligand requires a set of reorientations of protein side chains Q72, N31, and R142, although these residues have previously been determined as important for binding to xylose oligomers by mediating polar contacts. The loss of these polar contacts is compensated for in binding to XXXG by polar interactions mediated by other protein residues, T74, R115, and Y149, which interact mainly with the branching xyloses of the xyloglucan oligomer. Taken together, the present study illustrates in structural detail how CH‐π interactions can influence binding specificity and that flexibility is a key feature for the multi‐specificity displayed by X‐2 L110F, allowing for the accommodation of branched ligands. Proteins 2014; 82:3466–3475. © 2014 Wiley Periodicals, Inc.     

Statistical analysis of protein structures suggests that buried ionizable residues in proteins are hydrogen bonded or form salt bridges

It is well known that nonpolar residues are largely buried in the interior of proteins, whereas polar and ionizable residues tend to be more localized on the protein surface where they are solvent exposed. Such a distribution of residues between surface and interior is well understood from a thermodynamic point: nonpolar side chains are excluded from the contact with the solvent water, whereas polar and ionizable groups have favorable interactions with the water and thus are preferred at the protein surface. However, there is an increasing amount of information suggesting that polar and ionizable residues do occur in the protein core, including at positions that have no known functional importance. This is inconsistent with the observations that dehydration of polar and in particular ionizable groups is very energetically unfavorable. To resolve this, we performed a detailed analysis of the distribution of fractional burial of polar and ionizable residues using a large set of ?2600 nonhomologous protein structures. We show that when ionizable residues are fully buried, the vast majority of them form hydrogen bonds and/or salt bridges with other polar/ionizable groups. This observation resolves an apparent contradiction: the energetic penalty of dehydration of polar/ionizable groups is paid off by favorable energy of hydrogen bonding and/or salt bridge formation in the protein interior. Our conclusion agrees well with the previous findings based on the continuum models for electrostatic interactions in proteins. Proteins 2011; © 2011 Wiley‐Liss, Inc.     

Exploring the mechanism how Marburg virus VP35 recognizes and binds dsRNA by molecular dynamics simulations and free energy calculations

Filoviruses often cause terrible infectious disease which has not been successfully dealt with pharmacologically. All filoviruses encode a unique protein termed VP35 which can mask doubled‐stranded RNA to deactivate interferon. The interface of VP35–dsRNA would be a feasible target for structure‐based antiviral agent design. To explore the essence of VP35–dsRNA interaction, molecular dynamics simulation combined with MM‐GBSA calculations were performed on Marburg virus VP35–dsRNA complex and several mutational complexes. The energetic analysis indicates that nonpolar interactions provide the main driving force for the binding process. Although the intermolecular electrostatic interactions play important roles in VP35–dsRNA interaction, the whole polar interactions are unfavorable for binding which result in a low binding affinity. Compared with wild type VP35, the studied mutants F228A, R271A, and K298A have obviously reduced binding free energies with dsRNA reflecting in the reduction of polar or nonpolar interactions. The results also indicate that the loss of binding affinity for one dsRNA strand would abolish the total binding affinity. Three important residues Arg271, Arg294, and Lys298 which makes the largest contribution for binding in VP35 lose their binding affinity significantly in mutants. The uncovering of VP35–dsRNA recognition mechanism will provide some insights for development of antiviral drug. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 849–860, 2014.     

Structural evidence for a dominant role of nonpolar interactions in the binding of a transport/chemosensory receptor to its highly polar ligands.

The receptor, a maltose/maltooligosaccharide-binding protein, has been found to be an excellent system for the study of molecular recognition because its polar and nonpolar binding functions are segregated into two globular domains. The X-ray structures of the "closed" and "open" forms of the protein complexed with maltose and maltotetraitol have been determined. These sugars have approximately 3 times more accessible polar surface (from OH groups) than nonpolar surface (from small clusters of sugar ring CH bonds). In the closed structures, the oligosaccharides are buried in the groove between the two domains of the protein and bound by extensive hydrogen bonding interactions of the OH groups with the polar residues confined mostly in one domain and by nonpolar interactions of the CH clusters with four aromatic residues lodged in the other domain. Substantial contacts between the sugar hydroxyls and aromatic residues are also formed. In the open structures, the oligosaccharides are bound almost exclusively in the domain rich in aromatic residues. This finding, along with the analysis of buried surface area due to complex formations in the open and closed structures, supports a major role for nonpolar interactions in initial ligand binding even when the ligands have significantly greater potential for highly specific polar interactions.     

The interplay between effector binding and allostery in an engineered protein switch

The protein design rules for engineering allosteric regulation are not well understood. A fundamental understanding of the determinants of ligand binding in an allosteric context could facilitate the design and construction of versatile protein switches and biosensors. Here, we conducted extensive in vitro and in vivo characterization of the effects of 285 unique point mutations at 15 residues in the maltose‐binding pocket of the maltose‐activated β‐lactamase MBP317‐347. MBP317‐347 is an allosteric enzyme formed by the insertion of TEM‐1 β‐lactamase into the E. coli maltose binding protein (MBP). We find that the maltose‐dependent resistance to ampicillin conferred to the cells by the MBP317‐347 switch gene (the switch phenotype) is very robust to mutations, with most mutations slightly improving the switch phenotype. We identified 15 mutations that improved switch performance from twofold to 22‐fold, primarily by decreasing the catalytic activity in the absence of maltose, perhaps by disrupting interactions that cause a small fraction of MBP in solution to exist in a partially closed state in the absence of maltose. Other notable mutations include K15D and K15H that increased maltose affinity 30‐fold and Y155K and Y155R that compromised switching by diminishing the ability of maltose to increase catalytic activity. The data also provided insights into normal MBP physiology, as select mutations at D14, W62, and F156 retained high maltose affinity but abolished the switch's ability to substitute for MBP in the transport of maltose into the cell. The results reveal the complex relationship between ligand binding and allostery in this engineered switch.     

Mapping disease‐related missense mutations in the immunoglobulin‐like fold domain of lamin A/C reveals novel genotype–phenotype associations for laminopathies

Mutations in A‐type nuclear lamins cause laminopathies. However, genotype–phenotype correlations using the 340 missense mutations within the LMNA gene are unclear: partially due to the limited availability of three‐dimensional structure. The immunoglobulin (Ig)‐like fold domain has been solved, and using bioinformatics tools (including Polyphen‐2, Fold X, Parameter OPtimized Surfaces, and PocketPicker) we characterized 56 missense mutations for position, surface exposure, change in charge and effect on Ig‐like fold stability. We find that 21 of the 27 mutations associated with a skeletal muscle phenotype are distributed throughout the Ig‐like fold, are nonsurface exposed and predicted to disrupt overall stability of the Ig‐like fold domain. Intriguingly, the remaining 6 mutations clustered, had higher surface exposure, and did not affect stability. The majority of 9 lipodystrophy or 10 premature aging syndrome mutations also did not disrupt Ig‐like fold domain stability and were surface exposed and clustered in distinct regions that overlap predicted binding pockets. Although buried, the 10 cardiac mutations had no other consistent properties. Finally, most lipodystrophy and premature aging mutations resulted in a ‐1 net charge change, whereas skeletal muscle mutations caused no consistent net charge changes. Since premature aging, lipodystrophy and the subset of 6 skeletal muscle mutations cluster tightly in distinct, charged regions, they likely affect lamin A/C –protein/DNA/RNA interactions: providing a consistent genotype–phenotype relationship for mutations in this domain. Thus, this subgroup of skeletal muscle laminopathies that we term the ‘Skeletal muscle cluster’, may have a distinct pathological mechanism. These novel associations refine the ability to predict clinical features caused by certain LMNA missense mutations. Proteins 2014; 82:904–915. © 2013 Wiley Periodicals, Inc.     

Proteome‐wide analysis of phospho‐regulated PDZ domain interactions

A key function of reversible protein phosphorylation is to regulate protein–protein interactions, many of which involve short linear motifs (3–12 amino acids). Motif‐based interactions are difficult to capture because of their often low‐to‐moderate affinities. Here, we describe phosphomimetic proteomic peptide‐phage display, a powerful method for simultaneously finding motif‐based interaction and pinpointing phosphorylation switches. We computationally designed an oligonucleotide library encoding human C‐terminal peptides containing known or predicted Ser/Thr phosphosites and phosphomimetic variants thereof. We incorporated these oligonucleotides into a phage library and screened the PDZ (PSD‐95/Dlg/ZO‐1) domains of Scribble and DLG1 for interactions potentially enabled or disabled by ligand phosphorylation. We identified known and novel binders and characterized selected interactions through microscale thermophoresis, isothermal titration calorimetry, and NMR. We uncover site‐specific phospho‐regulation of PDZ domain interactions, provide a structural framework for how PDZ domains accomplish phosphopeptide binding, and discuss ligand phosphorylation as a switching mechanism of PDZ domain interactions. The approach is readily scalable and can be used to explore the potential phospho‐regulation of motif‐based interactions on a large scale.     

Mcl‐1–Bim complexes accommodate surprising point mutations via minor structural changes

Mcl‐1 is an antiapoptotic Bcl‐2‐family protein that protects cells against death. Structures of Mcl‐1, and of other anti‐apoptotic Bcl‐2 proteins, reveal a surface groove into which the α‐helical BH3 regions of certain proapoptotic proteins can bind. Despite high overall structural conservation, differences in this groove afford binding specificity that is important for the mechanism of Bcl‐2 family function. We report the crystal structure of human Mcl‐1 bound to a BH3 peptide derived from human Bim and the structures for three complexes that accommodate large physicochemical changes at conserved Bim sites. The mutations had surprisingly modest effects on complex stability, and the structures show that Mcl‐1 can undergo small changes to accommodate the mutant ligands. For example, a shift in a leucine side chain fills a hole left by an isoleucine‐to‐alanine mutation at the first hydrophobic buried position of Bim BH3. Larger changes are also observed, with shifting of helix α3 accommodating an isoleucine‐to‐tyrosine mutation at this same position. We surveyed the variation in available Mcl‐1 and Bcl‐x_L structures and observed moderate flexibility that is likely critical for facilitating interactions of diverse BH3‐only proteins with Mcl‐1. With the antiapoptotic Bcl‐2 family members attracting significant attention as therapeutic targets, these structures contribute to our growing understanding of how specificity is achieved and can help to guide the design of novel inhibitors that target Mcl‐1.     

A polymetamorphic protein

Arc repressor is a homodimeric protein with a ribbon‐helix–helix fold. A single polar‐to‐hydrophobic substitution (N11L) at a solvent‐exposed position leads to population of an alternate dimeric fold in which 3₁₀ helices replace a β‐sheet. Here we find that the variant Q9V/N11L/R13V (S‐VLV), with two additional polar‐to‐hydrophobic surface mutations in the same β‐sheet, forms a highly stable, reversibly folded octamer with approximately half the?α‐helical content of wild‐type Arc. At low protein concentration and low ionic strength, S‐VLV also populates both dimeric topologies previously observed for N11L, as judged by NMR chemical shift comparisons. Thus, accumulation of simple hydrophobic mutations in Arc progressively reduces fold specificity, leading first to a sequence with two folds and then to a manifold bridge sequence with at least three different topologies. Residues 9–14 of S‐VLV form a highly hydrophobic stretch that is predicted to be amyloidogenic, but we do not observe aggregates of higher order than octamer. Increases in sequence hydrophobicity can promote amyloid aggregation but also exert broader and more complex effects on fold specificity. Altered native folds, changes in fold coupled to oligomerization, toxic pre‐amyloid oligomers, and amyloid fibrils may represent a near continuum of accessible alternatives in protein structure space.     

Organic dyes as small molecule protein–protein interaction inhibitors for the CD40–CD154 costimulatory interaction

It is becoming increasingly clear that small molecules can often act as effective protein–protein interaction (PPI) inhibitors, an area of increasing interest for its many possible therapeutic applications. We have identified several organic dyes and related small molecules that (i) concentration‐dependently inhibit the important CD40–CD154 costimulatory interaction with activities in the low micromolar (µM) range, (ii) show selectivity toward this particular PPI, (iii) seem to bind on the surface of CD154, and (iv) concentration‐dependently inhibit the CD154‐induced B cell proliferation. They were identified through an iterative activity screening/structural similarity search procedure starting with suramin as lead, and the best smaller compounds, the main focus of the present work, achieved an almost 3‐fold increase in ligand efficiency (ΔG⁰/nonhydrogen atom = 0.8 kJ/N_nHa) approaching the average of known promising small‐molecule PPI inhibitors (～1.0 kJ/N_nHa). Since CD154 is a member of the tumor necrosis factor (TNF) superfamily of cell surface interaction molecules, inhibitory activities on the TNF‐R1–TNF‐α interactions were also determined to test for specificity, and the compounds selected here all showed more than 30‐fold selectivity toward the CD40–CD154 interaction. Because of their easy availability in various structural scaffolds and because of their good protein‐binding ability, often explored for tissue‐specific staining and other purposes, such organic dyes can provide a valuable addition to the chemical space searched to identify small molecule PPI inhibitors in general. Copyright © 2009 John Wiley & Sons, Ltd.     

Different combinations of atomic interactions predict protein‐small molecule and protein‐DNA/RNA affinities with similar accuracy

Interactions between proteins and other molecules play essential roles in all biological processes. Although it is widely held that a protein's ligand specificity is determined primarily by its three‐dimensional structure, the general principles by which structure determines ligand binding remain poorly understood. Here we use statistical analyses of a large number of protein?ligand complexes with associated binding‐affinity measurements to quantitatively characterize how combinations of atomic interactions contribute to ligand affinity. We find that there are significant differences in how atomic interactions determine ligand affinity for proteins that bind small chemical ligands, those that bind DNA/RNA and those that interact with other proteins. Although protein‐small molecule and protein‐DNA/RNA binding affinities can be accurately predicted from structural data, models predicting one type of interaction perform poorly on the others. Additionally, the particular combinations of atomic interactions required to predict binding affinity differed between small‐molecule and DNA/RNA data sets, consistent with the conclusion that the structural bases determining ligand affinity differ among interaction types. In contrast to what we observed for small‐molecule and DNA/RNA interactions, no statistical models were capable of predicting protein?protein affinity with >60% correlation. We demonstrate the potential usefulness of protein‐DNA/RNA binding prediction as a possible tool for high‐throughput virtual screening to guide laboratory investigations, suggesting that quantitative characterization of diverse molecular interactions may have practical applications as well as fundamentally advancing our understanding of how molecular structure translates into function. Proteins 2015; 83:2100–2114. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.