Search

The voltage-gated sodium channel Na

A novel multiple turn conformation has been observed for a segment GPGRAFY in the crystal structure of a complex of HIV-1 gp120 V3 loop peptide with the Fab fragment of a neutralizing antibody [Ghiara ct al. (1994) Science 264, 82-85]. A structural motif has been defined for the peptide segment, employing idealized backbone conformations characterized by ranges of virtual C-alpha torsion angles and bond angles. A search of 122 high-resolution protein crystal structures has permitted identification of 24 examples of similar structural motifs. Two major conformational families have been identified, which differ primarily in the conformation at residue 3. The observed conformation at residue 3 in family 1 is left-handed helical (alpha(L)) and that in family 2 is right-handed helical (alpha(R)). Of the 10 examples in family 1, 9 examples have Gly residues at position 3. Of the 12 examples in family 2, 7 examples have Asn/Asp at position 3. Computer modeling of the V3 loop tip sequence using the two backbone conformational families as starting points leads to minimum-energy conformations in which antigenically important side-chains occupy similar spatial arrangements. This stereochemical analysis of the V3 loop tip sequence suggests a rational basis for the design of synthetic analog peptides for use as viral antagonists or synthetic antigens. (C) Munksgaard 1995.

Most proteins exist in the cell as multi-component assemblies. However, which proteins need to be present simultaneously in order to perform a given function is frequently unknown. The first step toward this goal would be to predict proteins that can function only when in a complexed form. Here, we propose a scheme to distinguish whether the protein components are ordered (stable) or disordered when separated from their complexed partners. We analyze structural characteristics of several types of complexes, such as natively unstructured proteins, ribosomal proteins, two-state and three-state complexes, and crystal-packing dimers. Our analysis makes use of the fact that natively unstructured proteins, which undergo a disorder-to-order transition upon binding their partner, and stable monomeric proteins, which exist as dimers only in their crystal form, provide examples of two vastly different scenarios. We find that ordered monomers can be distinguished from disordered monomers on the basis of the per-residue surface and interface areas, which are significantly smaller for ordered proteins. With this scale, two-state dimers (where the monomers unfold upon dimer separation) and ribosomal proteins are shown to resemble disordered proteins. On the other hand, crystal-packing dimers, whose monomers are stable in solution, fall into the ordered protein category. While there should be a continuum in the distributions, nevertheless, the per-residue scale measures the confidence in the determination of whether a protein can exist as a stable monomer. Further analysis, focusing on the chemical and contact preferences at the interface, interior and exposed surface areas, reveals that disordered proteins lack a strong hydrophobic core and are composed of highly polar surface area. We discuss the implication of our results for de novo design of stable monomeric proteins and peptides.

In this paper we present a detailed atomic model for a protofilament, the most basic organization level, of the amyloid fibre formed by the peptide DFNKF. This pentapeptide is a segment derived from the human calcitonin, a natural amyloidogenic protein. Our model, which represents the outcome of extensive explicit solvent molecular dynamics (MD) simulations of different strand/sheet organizations, is a single beta-sheet filament largely without a hydrophobic core. Nevertheless, this structure is capable of reproducing the main features of the characteristic amyloid fibril organization and provides clues to the molecular basis of its experimental aggregation behaviour. Our results show that the side chains' chemical diversity induces the formation of a complex network of interactions that finally determine the microscopic arrangement of the strands at the protofilament level. This network of interactions, consisting of both side chain-side chain and backbone-side chain interactions, confers on the final single beta-sheet arrangement an unexpected stability, both by enhancing the association of related chemical groups and, at the same time, by shielding the hydrophobic segments from the polar solvent. The chemical physical characterization of this protofilament provides hints to the possible thermodynamical basis of the supra molecular organization that allows the formation of the filaments by lateral association of the preformed protofibrils. Its regular, highly polarized structure shows how other protofilaments can assemble. In terms of structural biology, our results clearly indicate that an amyloid organization implies a degree of complexity far beyond a simple nonspecific association of peptide strands via amide hydrogen bonds.

Proteins in the intracellular lipid-binding protein (iLBP) family show remarkably high structural conservation despite their low-sequence identity. A multiple-sequence alignment using 52 sequences of iLBP family members revealed 15 fully conserved positions, with a disproportionately high number of these (n=7) located in the relatively small helical region. The conserved positions displayed high structural conservation based on comparisons of known iLBP crystal structures. It is striking that the beta-sheet domain had few conserved positions, despite its high structural conservation. This observation prompted us to analyze pair-wise interactions within the beta-sheet region to ask whether structural information was encoded in interacting amino acid pairs. We conducted this analysis on the iLBP family member, cellular retinoic acid-binding protein I (CRABP I), whose folding mechanism is under study in our laboratory. Indeed, an analysis based on a simple classification of hydrophobic and polar amino acids revealed a network of conserved interactions in CRABP I that cluster spatially, suggesting a possible nucleation site for folding. Significantly, a small number of residues participated in multiple conserved interactions, suggesting a key role for these sites in the structure and folding of CRABP I. The results presented here correlate well with available experimental evidence on folding of CRABPs and their family members and suggest future experiments. The analysis also shows the usefulness of considering pair-wise conservation based on a simple classification of amino acids, in analyzing sequences and structures to find common core regions among homologues.

Enzyme function often involves a conformational change. There is a general agreement that loops play a vital role in correctly positioning the catalytically important residues. Nevertheless, predicting the functional loops and most importantly their role in enzyme function remains a difficult task. A major reason for this difficulty is that loops that undergo conformational change are frequently not well conserved in their primary sequence. beta1,4-Galactosyltransferase is one such enzyme. There, the amino acid sequence of a long loop that undergoes a large conformational change upon substrate binding is not well conserved. Our molecular dynamics simulations show that the large conformational change in the long loop is brought about by a second, interacting loop. Interestingly, while the structural change of the second loop is much smaller than that of the long loop, its sequence (particularly glycine residues) is highly conserved. We further examine the generality of the proposition that there are loops that trigger movements but nevertheless show little or no structural changes in crystals. We focus on two other enzymes, enolase and lipase. We chose these enzymes, since they too undergo conformational change upon ligand binding, however, they have different folds and different functions. Through multiple sets of simulations we show that the conformational change of the functional loop(s) is brought about through communication of flexibility by triggering loops that have several glycine residues. We further propose that similar to the conservation of common favorable fold types and structural motifs, evolution has also conserved common "skillful" mechanisms. Mechanisms may be conserved across different folds, sequences and functions, with adaptation to specific enzymatic roles.

Current data suggest that globular domains may form amyloids via different mechanisms. Nevertheless, there are indications that the initiation of the process takes place invariably in the less stable segments of a protein domain. We have studied the sequence and structural conservation of beta(2)-microglobulin that deposits into fibrils in dialysis-related amyloidosis. The dataset includes 51 high-resolution non-redundant structures of the antibody constant domain-like proteins (C1) and 132 related sequences. We describe a set of 30 conserved residues. Among them, 23 are conserved structurally, 16 are conserved sequentially and nine are conserved both sequentially and structurally. Strands A (12-18), G (91-95) and D (45-55) are the less conserved and stable segments of the domain, while strands B (22-28), C (36-41), E (62-70) and F (78-83) are the conserved and stable segments. We find that the conserved residues form a cluster with a network of interactions. The observed pattern of conservation is consistent with experimental data including H/D exchange, urea denaturation and limited proteolysis that suggest that strands A and G do not participate in the amyloid fibril. Additionally, the low conservation of strand D is consistent with the observation that this strand may acquire different conformations as seen in crystal structures of bound and isolated beta(2)-microglobulin. We used a docking technique to suggest a model for a fibril via stacking of beta(2)-microglobulin monomers. Our analysis suggests that the favored monomer building block for fibril elongation is the conformation of the isolated beta(2)-microglobulin, without the beta-bulge on strand D and without strands A and G participating in the fibril beta-sheet structure. This monomer retains all the conserved residues and their network of interactions, increasing the likelihood of its existence in solution. The inter-strand interaction between the two (monomer) building blocks forms a new continuous beta-sheet such that addition of monomers results in a fibril model that has the characteristic cross-beta structure.

It has been estimated that a large fraction of cellular proteins are natively disordered. Current opinion largely holds that natively disordered proteins are more ‘adaptive', leading to advantages in regulation and in binding diverse ligands. Here, we argue for another, simple, physically based reason. Disordered proteins often have large intermolecular interfaces, the size of which is dictated by protein function. For proteins to be stable as monomers with extensive interfaces, protein size would need to be 2–3 times larger. This would either increase cellular crowding or enlarge the size of the cell by 15–30%, owing to the increase in the sequence length. Smaller sizes of cells, proteins, DNA and RNA conserve energy. Thus, disordered proteins provide a simple yet elegant solution to having large intermolecular interfaces, but with smaller protein, genome and cell sizes.

Although 8-anilinonaphthalene-1-sulfonic acid (ANS) is frequently used in protein folding studies, the structural and thermodynamic effects of its binding to proteins are not well understood. Using high-resolution two-dimensional NMR and human interleukin-1 receptor antagonist (IL-1ra) as a model protein, we obtained detailed information on ANS-protein interactions in the absence and presence of urea. The effects of ambient to elevated temperatures on the affinity and specificity of ANS binding were assessed from experiments performed at 25 degrees C and 37 degrees C. Overall, the affinity of ANS was lower at 37 degrees C compared to 25 degrees C, but no significant change in the site specificity of binding was observed from the chemical shift perturbation data. The same site-specific binding was evident in the presence of 5.2 M urea, well within the unfolding transition region, and resulted in selective stabilization of the folded state. Based on the two-state denaturation mechanism, ANS-dependent changes in the protein stability were estimated from relative intensities of two amide resonances specific to the folded and unfolded states of IL-1ra. No evidence was found for any ANS-induced partially denatured or aggregated forms of IL-1ra throughout the experimental conditions, consistent with a cooperative and reversible denaturation process. The NMR results support earlier observations on the tendency of ANS to interact with solvent-exposed positively charged sites on proteins. Under denaturing conditions, ANS binding appears to be selective to structured states rather than unfolded conformations. Interestingly, the binding occurs within a previously identified aggregation-critical region in IL-1ra, thus providing an insight into ligand-dependent protein aggregation.

Dynamic molecular interaction networks underlie biological phenomena. Among the many genes which are involved, p53 plays a central role in networks controlling cellular life and death. It not only operates as a tumor suppressor, but also helps regulate hundreds of genes in response to various types of stress. To accomplish these functions as a guardian of the genome, p53 interacts extensively with both nucleic acids and proteins. This paper examines the physical interfaces of the p53 protein with cellular proteins. Previously, in the analysis of the structures of protein–protein complexes, we have observed that amino acids Trp, Met and Phe are important for protein–protein interactions in general. Here we show that these residues are critical for the many functions of p53. Several clusters of the Trp/Met/Phe residues are involved in the p53 protein–protein interactions. Phe19/Trp23 in the TA1 region extensively binds to the transcriptional factors and the MDM2 protein. Trp53/Phe54 in the TA2 region is crucial for transactivation and DNA replication. Met243 in the core domain interacts with 53BP1, 53BP2 and Rad 51 proteins. Met384/Phe385 in the C-terminal region interacts with the S100B protein and the Bromodomain of the CBP protein. Thus, these residues may assist in elucidating the p53 interactions when structural data are not available.

Understanding and ultimately predicting protein associations is immensely important for functional genomics and drug design. Here, we propose that binding sites have preferred organizations. First, the hot spots cluster within densely packed 'hot regions'. Within these regions, they form networks of interactions. Thus, hot spots located within a hot region contribute cooperatively to the stability of the complex. However, the contributions of separate, independent hot regions are additive. Moreover, hot spots are often already pre-organized in the unbound (free) protein states. Describing a binding site through independent local hot regions has implications for binding site definition, design and parametrization for prediction. The compactness and cooperativity emphasize the similarity between binding and folding. This proposition is grounded in computation and experiment. It explains why summation of the interactions may over-estimate the stability of the complex. Furthermore, statistically, charge-charge coupling of the hot spots is disfavored. However, since within the highly packed regions the solvent is screened, the electrostatic contributions are strengthened. Thus, we propose a new description of protein binding sites: a site consists of (one or a few) self-contained cooperative regions. Since the residue hot spots are those conserved by evolution, proteins binding multiple partners at the same sites are expected to use all or some combination of these regions.

Obtaining atomic resolution structural models of amyloid fibrils is currently impossible, yet crucial for our understanding of the amyloid mechanism. Different pathways in the transformation of a native globular domain to an amyloid fibril invariably involve domain destabilization. Hence, locating the unstable segments of a domain is important for understanding its amyloidogenic transformation and possibly control it. Since relative conservation is suggested to relate to local stability [H. Benyamini, K. Gunasekaran, H. Wolfson, R. Nussinov, Conservation and amyloid formation: a study of the gelsolin-like family, Proteins 51 (2003) 266–282. [24]], we performed an extensive, sequence and structure conservation analysis of the β 2 -microglobulin (β 2 -m) domain. Our dataset include 51 high resolution structures belonging to the “C1 set domain” family and 132 clustered PSI-BLAST search results. Segments of the β 2 -m domain corresponding to strands A (residues 12–18), D (45–55) and G (91–95) were found to be less conserved and stable, while the central strands B (residues 22–28), C (36–41), E (62–70) and F (78–83) were found conserved and stable. Our findings are supported by accumulating observations from various experimental methods, including urea denaturation, limited proteolysis, H/D exchange and structure determination by both NMR and X-ray crystallography. We used our conservation findings together with experimental literature information to suggest a structural model for the polymerized unit of β 2 -m. Pairwise protein docking and subsequent monomer stacking in the same manner suggest a fibril model consistent with the cross-β structure.

Loop flexibility in enzymes plays a vital role in correctly positioning catalytically important residues. This strong relationship between enzyme flexibility and function provides an opportunity to engineer new substrates and inhibitors. It further allows the design of site-directed mutagenesis experiments to explore enzymatic activity through the control of flexibility of a functional loop. Earlier, we described a novel mechanism in which a small loop triggers the motions of a functional loop in three enzymes (beta-1,4-galactosyltransferase, lipase, and enolase) unrelated in sequence, structure, or function. Here, we further address the question of how the interactions between various flexible loops modulate the movements of the functional loop. We examine beta-1,4-galactosyltransferase as a model system in which a Long loop undergoes a large conformational change (moves in space up to 20 A) upon substrate binding in addition to a small loop (Trp loop) that shows a considerably smaller conformational change. Our molecular-dynamics simulations carried out in implicit and explicit solvent show that, in addition to these two loops, two other neighboring loops are also highly flexible. These loops are in contact with either the Long loop or the Trp loop. Analysis of the covariance of the spatial displacement of the residues reveals that coupled motions occur only in one of these two loops. Sequence analysis indicates that loops correlated in their motions also have highly conserved residues involved in the loop-loop interactions. Further, analysis of crystal structures and simulations in explicit water open the possibility that the Trp loop that triggers the movement of the Long loop in the unbound conformation may also play the same role in the substrate-bound conformation through its contact with the conserved and correlated third loop. Our proposition is supported by the observation that four of the five conserved positions in the third loop are at the interface with the Trp loop. Evolution appears to select residues that drive the functional Long loop to a large conformational change. These observations suggest that altering selected loop-loop interactions might modulate the movements of the functional loop.

There is a general agreement that electrostatic interactions play a significant role in the structure and function of biological molecules. However, obtaining quantitative estimation of the electrostatic energy requires computational models that capture the microscopic nature of the heterogeneous environment of macromolecules. This protocol offers elaboration on one of the common methods to calculate the electrostatic energetic contributions using continuum electrostatics. The method involves solving the Poisson-Boltzmann (PB) equation numerically and regarding the solute as having a homogenous dielectric constant. In order to apply this method, a three dimensional structure of the molecule derived from experimental data (crystallography, NMR) or modeling techniques is required. The protocol will focus on the DelPhi program (Accelrys Inc. San Diego), which is one of the most common programs used for the estimation of electrostatic free energy contribution. A simple procedure of assigning criteria and parameters (charge distribution, solvent and solute dielectric constants, iterations, grid resolution, etc) enables one to illustrate an electrostatic potential map and estimate the electrostatic free energy, although with limited accuracy.

Investigators have recently turned to studies of protein families to shed light on the mechanism of protein folding. In small proteins for which detailed analysis has been performed, recent studies show that transition-state structure is generally conserved. The number and structures of populated folding intermediates have been found to vary in homologous families of larger (greater than 100-residue) proteins, reflecting a balance of local and global interactions.