Helix lap-joints as ion-binding sites: DNA-binding motifs and Ca-binding "EF hands" are related by charge and sequence reversal

Hybrid method for representing ions in implicit solvation calculations

Fast and accurate calculations of the electrostatic features of highly charged biomolecules such as DNA, RNA, and highly charged proteins are crucial and challenging tasks. Traditional implicit solvent methods calculate the electrostatic features quickly, but these methods are not able to balance the high net biomolecular charges effectively. Explicit solvent methods add unbalanced ions to neutralize the highly charged biomolecules in molecular dynamic simulations, which require more expensive computing resources. Here we report developing a novel method, Hybridizing Ions Treatment (HIT), which hybridizes the implicit solvent method with an explicit method to realistically calculate the electrostatic potential for highly charged biomolecules. HIT utilizes the ionic distribution from an explicit method to predict the bound ions. The bound ions are then added in the implicit solvent method to perform the electrostatic potential calculations. In this study, two training sets were developed to optimize parameters for HIT. The performance on the testing set demonstrates that HIT significantly improves the electrostatic calculations. Results on molecular motors myosin and kinesin reveal some mechanisms and explain some previous experimental findings. HIT can be widely used to study highly charged biomolecules, including DNA, RNA, molecular motors, and other highly charged biomolecules. The HIT package is available at http://compbio.utep.edu/static/downloads/download_hit.zip.

Structural basis for adPEO-causing mutations in the mitochondrial TWINKLE helicase

TWINKLE is the helicase involved in replication and maintenance of mitochondrial DNA (mtDNA) in mammalian cells. Structurally, TWINKLE is closely related to the bacteriophage T7 gp4 protein and comprises a helicase and primase domain joined by a flexible linker region. Mutations in and around this linker region are responsible for autosomal dominant progressive external ophthalmoplegia (adPEO), a neuromuscular disorder associated with deletions in mtDNA. The underlying molecular basis of adPEO-causing mutations remains unclear, but defects in TWINKLE oligomerization are thought to play a major role. In this study, we have characterized these disease variants by single-particle electron microscopy and can link the diminished activities of the TWINKLE variants to altered oligomeric properties. Our results suggest that the mutations can be divided into those that (i) destroy the flexibility of the linker region, (ii) inhibit ring closure and (iii) change the number of subunits within a helicase ring. Furthermore, we demonstrate that wild-type TWINKLE undergoes large-scale conformational changes upon nucleoside triphosphate binding and that this ability is lost in the disease-causing variants. This represents a substantial advancement in the understanding of the molecular basis of adPEO and related pathologies and may aid in the development of future therapeutic strategies.

The role of key residues in structure, function, and stability of cytochrome-c

Cytochrome-c (cyt-c), a multi-functional protein, plays a significant role in the electron transport chain, and thus is indispensable in the energy-production process. Besides being an important component in apoptosis, it detoxifies reactive oxygen species. Two hundred and eighty-five complete amino acid sequences of cyt-c from different species are known. Sequence analysis suggests that the number of amino acid residues in most mitochondrial cyts-c is in the range 104 ± 10, and amino acid residues at only few positions are highly conserved throughout evolution. These highly conserved residues are Cys14, Cys17, His18, Gly29, Pro30, Gly41, Asn52, Trp59, Tyr67, Leu68, Pro71, Pro76, Thr78, Met80, and Phe82. These are also known as "key residues", which contribute significantly to the structure, function, folding, and stability of cyt-c. The three-dimensional structure of cyt-c from ten eukaryotic species have been determined using X-ray diffraction studies. Structure analysis suggests that the tertiary structure of cyt-c is almost preserved along the evolutionary scale. Furthermore, residues of N/C-terminal helices Gly6, Phe10, Leu94, and Tyr97 interact with each other in a specific manner, forming an evolutionary conserved interface. To understand the role of evolutionary conserved residues on structure, stability, and function, numerous studies have been performed in which these residues were substituted with different amino acids. In these studies, structure deals with the effect of mutation on secondary and tertiary structure measured by spectroscopic techniques; stability deals with the effect of mutation on T m (midpoint of heat denaturation), ∆G D (Gibbs free energy change on denaturation) and folding; and function deals with the effect of mutation on electron transport, apoptosis, cell growth, and protein expression. In this review, we have compiled all these studies at one place. This compilation will be useful to biochemists and biophysicists interested in understanding the importance of conservation of certain residues throughout the evolution in preserving the structure, function, and stability in proteins.

Polyproline and triple helix motifs in host-pathogen recognition

Secondary structure elements often mediate protein-protein interactions. Despite their low abundance in folded proteins, polyproline II (PPII) and its variant, the triple helix, are frequently involved in protein-protein interactions, likely due to their peculiar propensity to be solvent-exposed. We here review the role of PPII and triple helix in mediating host-pathogen interactions, with a particular emphasis to the structural aspects of these processes. After a brief description of the basic structural features of these elements, examples of host-pathogen interactions involving these motifs are illustrated. Literature data suggest that the role played by PPII motif in these processes is twofold. Indeed, PPII regions may directly mediate interactions between proteins of the host and the pathogen. Alternatively, PPII may act as structural spacers needed for the correct positioning of the elements needed for adhesion and infectivity.

Recent investigations have highlighted that collagen triple helix is also a common target for bacterial adhesins. Although structural data on complexes between adhesins and collagen models are rather limited, experimental and theoretical studies have unveiled some interesting clues of the recognition process. Interestingly, very recent data show that not only is the triple helix used by pathogens as a target in the host-pathogen interaction but it may also act as a bait in these processes since bacterial proteins containing triple helix regions have been shown to interact with host proteins.

As both PPII and triple helix expose several main chain non-satisfied hydrogen bond acceptors and donors, both elements are highly solvated. The preservation of the solvation state of both PPII and triple helix upon protein-protein interaction is an emerging aspect that will be here thoroughly discussed.

Doing molecular biophysics: finding, naming, and picturing signal within complexity

A macromolecular structure, as measured data or as a list of coordinates or even on-screen as a full atomic model, is an extremely complex and confusing object. The underlying rules of how it folds, moves, and interacts as a biological entity are even less evident or intuitive to the human mind. To do science on such molecules, or to relate them usefully to higher levels of biology, we need to start with a natural history that names their features in meaningful ways and with multiple representations (visual or algebraic) that show some aspect of their organizing principles. The two of us have jointly enjoyed a highly varied and engrossing career in biophysical research over nearly 50 years. Our frequent changes of emphasis are tied together by two threads: first, by finding the right names, visualizations, and methods to help both ourselves and others to better understand the 3D structures of protein and RNA molecules, and second, by redefining the boundary between signal and noise for complex data, in both directions-sometimes identifying and promoting real signal up out of what seemed just noise, and sometimes demoting apparent signal into noise or systematic error. Here we relate parts of our scientific and personal lives, including ups and downs, influences, anecdotes, and guiding principles such as the title theme.

Alpha-alpha linking motifs and interhelical orientations

While the geometry and sequence preferences of turns that link two beta-strands have been exhaustively explored, the corresponding preferences for sequences that link helical structures have been less well studied. Here we examine the interhelical geometry of two connected helices as a function of their link's length. The interhelical geometry of a helical pair appears to be significantly influenced by the number of linking residues. Furthermore, for relatively short link lengths, a very limited number of predominant conformations are observed, which can be categorized by their phi/psi angles. No more than two predominant linking backbone conformations are observed for a given link length, and some linking backbone conformations correlate strongly with distinctive interhelical geometric parameters. In this study, sequence and hydrogen-bonding patterns were defined for predominant interhelical link motifs. These results should assist in both protein structure prediction and de novo protein design.

Ca2+binding sites in calmodulin and troponin C alter interhelical angle movements

Molecular dynamics analyses were performed to examine conformational changes in the C-domain of calmodulin and the N-domain of troponin C induced by binding of Ca(2+) ions. Analyses of conformational changes in calmodulin and troponin C indicated that the shortening of the distance between Ca(2+) ions and Ca(2+) binding sites of helices caused widening of the distance between Ca(2+) binding sites of helices on opposite sides, while the hydrophobic side chains in the center of helices hardly moved due to their steric hindrance. This conformational change acts as the clothespin mechanism.

Visualizing and quantifying molecular goodness-of-fit: small-probe contact dots with explicit hydrogen atoms

The technique of small-probe contact dot surfaces is described as a method for calculating and displaying the detailed atomic contacts inside or between molecules. It allows one both to measure and to visualize directly the goodness-of-fit of packing interactions. It requires both highly accurate structures and also the explicit inclusion of all hydrogen atoms and their van der Waals interactions. A reference dataset of 100 protein structures was chosen on the basis of resolution (1.7 A or better), crystallographic R-value, non-homology, and the absence of any unusual problems. Hydrogen atoms were added in standard geometry and, where needed, with rotational optimization of OH, SH, and NH+3 positions. Side-chain amide orientations were corrected where required by NH van der Waals clashes, as described in the accompanying paper. It was determined that, in general, methyl groups pack well in the default staggered conformation, except for the terminal methyl groups of methionine residues, which required rotational optimization. The distribution of serious clashes (i.e. non-H-bond overlap of >/=0.4 A) was studied as a function of resolution, alternate conformations, and temperature factor (B), leading to the decision that packing and other structural features would not be analyzed for residues in 'b' alternate conformations or with B-factors of 40 or above. At the level of the fine details analyzed here, structural accuracy improves quite significantly over the range from 1.7 to 1.0 A resolution. These high-resolution structures show impressively well-fitted packing interactions, with some regions thoroughly interdigitated and other regions somewhat sparser. Lower-resolution structures or model structures could undoubtedly be improved in accuracy by the incorporation of this additional information: for example, nucleic acid structures in non-canonical conformations are often very accurate for the bases and much less reliable for the backbone, whose conformation could be specified better by including explicit H atom geometry and contacts. The contact dots are an extremely sensitive method of finding problem areas, and often they can suggest how to make improvements. They can also provide explanations for structural features that have been described only as empirical regularities, which is illustrated by showing that the commonest rotamer of methionine (a left-handed spiral, with all chi values near -60 degrees) is preferred because it provides up to five good H atom van der Waals contacts. This methodology is thus applicable in two different ways: (1) for finding and correcting errors in structure models (either experimental or theoretical); and (2) for analyzing interaction patterns in the molecules themselves.

Structural determinants present in the C-terminal binding protein binding site of adenovirus early region 1A proteins

The C-terminal binding protein (CtBP) has previously been shown to bind to a highly conserved six-amino acid motif very close to the C terminus of adenovirus early region 1A (Ad E1A) proteins. We have developed an enzyme-linked immunosorbent assay that has facilitated the screening of synthetic peptides identical or similar to the binding site on Ad E1A for their ability to bind CtBP and thus inhibit its interaction with Ad12 E1A. It has been shown that amino acids both C-terminal and N-terminal to the original proposed binding site contribute to the interaction of peptides with CtBP. Single amino acid substitutions across the binding site appreciably alter the Kd of the peptide for CtBP, indicative of a marked reduction in the affinity of the peptide for CtBP. The solution structures of synthetic peptides equivalent to the C termini of both Ad5 and Ad12 E1A and two substituted forms of these have been determined by proton NMR spectroscopy. Both the Ad12 and Ad5 peptides dissolved in trifluoroethanol/water mixtures were found to adopt regular secondary structural conformations seen as a series of beta-turns. An Ad12 peptide bearing a substitution that resulted in only very weak binding to CtBP (Ad12 L258G) was found to be random coil in solution. However, a second mutant (Ad12 V256K), which bound to CtBP rather more strongly (although not as well as the wild type), adopted a conformation similar to that of the wild type. We conclude that secondary structure (beta-turns) and an appropriate series of amino acid side chains are necessary for recognition by CtBP.

Helix capping

Helix-capping motifs are specific patterns of hydrogen bonding and hydrophobic interactions found at or near the ends of helices in both proteins and peptides. In an alpha-helix, the first four >N-H groups and last four >C=O groups necessarily lack intrahelical hydrogen bonds. Instead, such groups are often capped by alternative hydrogen bond partners. This review enlarges our earlier hypothesis (Presta LG, Rose GD. 1988. Helix signals in proteins. Science 240:1632-1641) to include hydrophobic capping. A hydrophobic interaction that straddles the helix terminus is always associated with hydrogen-bonded capping. From a global survey among proteins of known structure, seven distinct capping motifs are identified-three at the helix N-terminus and four at the C-terminus. The consensus sequence patterns of these seven motifs, together with results from simple molecular modeling, are used to formulate useful rules of thumb for helix termination. Finally, we examine the role of helix capping as a bridge linking the conformation of secondary structure to supersecondary structure.

Determination of the structure of the N-terminal splice region of the cyclic AMP-specific phosphodiesterase RD1 (RNPDE4A1) by 1H NMR and identification of the membrane association domain using chimeric constructs

A 25-residue peptide representing the membrane targeting N-terminal splice region of the cyclic AMP phosphodiesterase RD1 (RNPDE4A1) was synthesized, and its structure was determined by 1H NMR. Two independently folding helical regions were identified, separated by a highly mobile "hinge" region. The first helical region was formed by an N-terminal amphipathic alpha-helix, and the second consisted of multiple overlapping turns and contained a distinct compact, hydrophobic, tryptophan-rich domain (residues 14-20). Chimeric molecules, formed between the N-terminal region of RD1 and the soluble bacterial protein chloramphenicol acetyltransferase, were used in an in vitro system to determine the features within the splice region that were required for membrane association. The ability of RD1-chloramphenicol acetyltransferase chimera to become membrane-associated was not affected by deletion of any of the following regions: the apolar section (residues 2-7) of the first helical region, the polar part of this region together with the hinge region (residues 8-13), or the polar end of the C-terminal helical region (residues 21-25). In marked contrast, deletion of the compact, hydrophobic tryptophan-rich domain (residues 14-20) found in the second helical region obliterated membrane association. Replacement of this domain with a hydrophobic cassette of seven alanine residues also abolished membrane association, indicating that membrane-association occurred by virtue of specific hydrophobic interactions with residues within the compact, tryptophan-rich domain. The structure of this domain is well defined in the peptide, and although the region is helical, both the backbone and the distribution of side chains are somewhat distorted as compared with an ideal alpha-helix. Hydrophobic interactions, such as the "stacked" rings of residues Pro14 and Trp15, stabilize this domain with the side chain of residue Leu16 adopting a central position, interacting with the side chains of all three tryptophan residues 15, 19, and 20. These bulky side chains thus form a hydrophobic cluster. In contrast, the side chain of residue Val17 is relatively exposed, pointing out from the opposite "face" of the peptide. Although it appears that this compact, tryptophan-rich domain is responsible for membrane association, at present the target site and hence the specific interactions involved in membrane targeting by the RD1 splice region remain unidentified.

Helical fold prediction for the cyclin box

The smooth progression of the eukaryotic cell cycle relies on the periodic activation of members of a family of cell cycle kinases by regulatory proteins called cyclins. Outside of the cell cycle, cyclin homologs play important roles in regulating the assembly of transcription complexes; distant structural relatives of the conserved cyclin core or "box" can also function as general transcription factors (like TFIIB) or survive embedded in the chain of the tumor suppressor, retinoblastoma protein. The present work attempts the prediction of the canonical secondary, supersecondary, and tertiary fold of the minimal cyclin box domain using a combination of techniques that make use of the evolutionary information captured in a multiple alignment of homolog sequences. A tandem set of closely packed, helical modules are predicted to form the cyclin box domain.

Common prevalence of alanine and glycine in mobile reactive centre loops of serpins and viral fusion peptides: do prions possess a fusion peptide?

Serpin reactive centre loops and fusion peptides released by proteolytic cleavage are particularly mobile. Their amino acid compositions reveal a common and unusual abundance of alanine, accompanied by high levels of glycine. These two small residues, which are not simultaneously abundant in stable helices (standard or transmembrane), probably play an important role in mobility. Threonine and valine (also relatively small amino acids) are also abundant in these two kinds of peptides. Moreover, the known 3D structures of an uncleaved serpin reactive centre and a fusion peptide are strikingly similar. Such sequences possess many small residues and are found in several signal peptides and in PrP, a protein associated with spongiform encephalopathies and resembling virus envelope proteins. These properties may be related to the infection mechanisms of these diseases.

Molecular recognition in DNA-binding proteins and enzymes

The molecular basis for the specificity and activity of protein-DNA interactions is currently being established from the combination of results on the structure and biochemistry of DNA-binding proteins and enzymes. Data detailed in the 12 most recent studies on DNA-binding protein and enzyme structures, including the major advances in the elucidation of enzyme-mediated DNA-repair processes, have both increased understanding of DNA recognition and enhanced prospects for the design of novel DNA-binding proteins in the future.

Families and the structural relatedness among globular proteins

Protein structures come in families. Are families "closely knit" or "loosely knit" entities? We describe a measure of relatedness among polymer conformations. Based on weighted distance maps, this measure differs from existing measures mainly in two respects: (1) it is computationally fast, and (2) it can compare any two proteins, regardless of their relative chain lengths or degree of similarity. It does not require finding relative alignments. The measure is used here to determine the dissimilarities between all 12,403 possible pairs of 158 diverse protein structures from the Brookhaven Protein Data Bank (PDB). Combined with minimal spanning trees and hierarchical clustering methods, this measure is used to define structural families. It is also useful for rapidly searching a dataset of protein structures for specific substructural motifs. By using an analogy to distributions of Euclidean distances, we find that protein families are not tightly knit entities.

Atomic structure of the DNA repair [4Fe-4S] enzyme endonuclease III

The crystal structure of the DNA repair enzyme endonuclease III, which recognizes and cleaves DNA at damaged bases, has been solved to 2.0 angstrom resolution with an R factor of 0.185. This iron-sulfur [4Fe-4S] enzyme is elongated and bilobal with a deep cleft separating two similarly sized domains: a novel, sequence-continuous, six-helix domain (residues 22 to 132) and a Greek-key, four-helix domain formed by the amino-terminal and three carboxyl-terminal helices (residues 1 to 21 and 133 to 211) together with the [4Fe-4S] cluster. The cluster is bound entirely within the carboxyl-terminal loop with a ligation pattern (Cys-X6-Cys-X2-Cys-X5-Cys) distinct from all other known [4Fe-4S] proteins. Sequence conservation and the positive electrostatic potential of conserved regions identify a surface suitable for binding duplex B-DNA across the long axis of the enzyme, matching a 46 angstrom length of protected DNA. The primary role of the [4Fe-4S] cluster appears to involve positioning conserved basic residues for interaction with the DNA phosphate backbone. The crystallographically identified inhibitor binding region, which recognizes the damaged base thymine glycol, is a seven-residue beta-hairpin (residues 113 to 119). Location and side chain orientation at the base of the inhibitor binding site implicate Glu112 in the N-glycosylase mechanism and Lys120 in the beta-elimination mechanism. Overall, the structure reveals an unusual fold and a new biological function for [4Fe-4S] clusters and provides a structural basis for studying recognition of damaged DNA and the N-glycosylase and apurinic/apyrimidinic-lyase mechanisms.

Taxonomy and conformational analysis of loops in proteins

We propose a general classification scheme for loops, aperiodic segments of protein structure. In an effort to avoid the geometric complexity created by non-repeating phi psi angles, a morphologic definition that focuses upon the linearity and planarity of loops is utilized. Out of 432 loops (4 to 20 residues in length) extracted from 67 proteins, 205 are classified as linear (straps), 133 as non-linear and planar (omegas), and 86 as non-linear and non-planar (zetas). The remaining 8 are classified as compound loops because they contain a combination of strap, omega, and zeta morphologies. We introduce a structural alphabet as a shorthand notation for describing local conformation. The symbols of this alphabet are based on the virtual dihedral angle joining four consecutive alpha carbons. The notation is used to provide a compact description of loop motifs in phosphate binding and calcium binding proteins. Since similar loop conformations form similar "words", the structural sequence facilitates the search for common structural motifs in a family of loops. Contrary to the view of loops as "random coils", we find loops to have positional preferences for amino acid residues analogous to those previously described for beta-turns.

Constraints on amino acid substitutions in the N-terminal helix of cytochrome c explored by random mutagenesis

The interaction of the N- and C-terminal helices is a hallmark of the cytochrome c family. Oligodeoxyribonucleotide-directed random mutagenesis within the gene encoding the C102T protein variant of Saccharomyces cerevisiae iso-1-cytochrome c was used to generate a library of mutations at the evolutionary invariant residues Gly-6 and Phe-10 in the N-terminal helix. Transformation of this library (contained on a low-copy-number yeast shuttle phagemid) into a yeast strain lacking a functional cytochrome c, followed by selection for cytochrome c function, reveals that 4-10% of the 400 possible amino acid substitutions are compatible with function. DNA sequence analysis of phagemids isolated from transformants exhibiting the functional phenotype elucidates the requirements for a stable helical interface. Basic residues are not tolerated at position 6 or 10. There is a broad volume constraint for amino acids at position 6. The amino acid substitutions observed to be compatible with function at Phe-10 show that the hydrophobic effect alone is sufficient to promote helical association. There are severe constraints that limit the combinations consistent with function, but the number of functionally consistent combinations observed exemplifies the plasticity of proteins.

Structural aspects of protein-DNA recognition

Images

Elucidating the structural chemistry of glycosaminoglycan recognition by protein C inhibitor

Glycosaminoglycans (GAGs) including heparin accelerate the inhibition of serine proteases by serine protease inhibitors (serpins), an essential process in regulating blood coagulation. to analyze the molecular basis for GAG recognition by the plasma serpin protein C inhibitor (PCI; also known as plasminogen activator inhibitor 3), we have constructed a complete, energy-minimized, three-dimensional model of PCI by using the structure of homologous alpha 1-antitrypsin as a template. Sequence analysis, hydrogen-bonding environment, and shape complementarity suggested that the N-terminal residues of PCI, which are not homologous to those of alpha 1-antitrypsin, form an amphipathic alpha-helix, here designated A+ since it precedes the alpha 1-antitrypsin A helix. Electrostatic calculations revealed a single, highly positive surface region arising from both the A+ and H helices, suggesting that this two-helix motif is required for GAG binding by PCI. The dominant role of electrostatic interactions in PCI-heparin binding was confirmed by the strong ionic strength dependence of heparin stimulation. The involvement of the A+ helix in heparin binding was verified by demonstrating that an anti-PCI antibody that specifically binds the A+ peptide blocks heparin binding.

Yeast TATA-box transcription factor gene

The first step in the transcription of most protein-encoding genes in eukaryotes is the binding of a transcription factor to the TATA-box promoter element. This TATA-box transcription factor was purified from extracts of the yeast Saccharomyces cerevisiae by using reconstitution of in vitro transcription reactions as an assay. The activity copurified with a protein whose sodium dodecyl sulfate/polyacrylamide gel mobility is 25 kDa. The sequence of the amino-terminal 21 residues of this protein was determined by sequential Edman degradation. A yeast genomic library was screened with mixed oligonucleotides encoding six residues of the protein sequence. The yeast TATA-box factor gene was cloned, and DNA sequencing revealed a 720-base-pair open reading frame encoding a 27,016-Da protein. The identity of the clone was confirmed by expressing the gene in Escherichia coli and detecting TATA-box factor DNA binding and transcriptional activities in extracts of the recombinant E. coli. The TATA-box factor gene was mapped to chromosome five of S. cerevisiae. RNA blot hybridization and nuclease S1 analysis indicated that the major TATA-box factor mRNA is 1.3 kilobases, including an unusually long 5' untranslated region of 188 +/- 5 nucleotides. Homology searches showed a region of distant similarity to the calcium-binding structures of calpains, a structure that has a conformation similar to the helix-turn-helix motif of DNA binding proteins.

Structural basis of DNA-protein recognition

Recent structure determinations of several repressor-operator complexes have shown how proteins can recognize specific binding sites on DNA. Although each of these repressor proteins belongs to the 'helix-turn-helix' class of DNA-binding proteins, they do not use a simple code for recognition.