A structural study for the optimisation of functional motifs encoded in protein sequences

Background  

A large number of PROSITE patterns select false positives and/or miss known true positives. It is possible that – at least in some cases – the weak specificity and/or sensitivity of a pattern is due to the fact that one, or maybe more, functional and/or structural key residues are not represented in the pattern. Multiple sequence alignments are commonly used to build functional sequence patterns. If residues structurally conserved in proteins sharing a function cannot be aligned in a multiple sequence alignment, they are likely to be missed in a standard pattern construction procedure.     

A structural database for k-turn motifs in RNA

The kink-turn (k-turn) is a common structural motif in RNA that introduces a tight kink into the helical axis. k-turns play an important architectural role in RNA structures and serve as binding sites for a number of proteins. We have created a database of known and postulated k-turn sequences and three-dimensional (3D) structures, available via the internet. This site provides (1) a database of sequence and structure, as a resource for the RNA community, and (2) a tool to enable the manipulation and comparison of 3D structures where known.     

A computational tool for designing FRET protein biosensors by rigid-body sampling of their conformational space

Biosensors relying on the fluorescence resonance energy transfer (FRET) between fluorescent proteins have been used for live-cell imaging of cellular events including Ca(2+) signaling. The efficiency of energy transfer between the donor and acceptor fluorescent proteins depends on the relative distance and orientation between them, which become altered by conformational changes of a fused sensory protein caused by a cellular event. In this way, changes in FRET efficiency of Ca(2+) biosensors can be correlated with Ca(2+) concentrations. The design of these FRET biosensors can be improved by modeling conformational changes before and after a cellular event. Hence, a computational tool called FPMOD was developed to predict FRET efficiency changes by constructing FRET biosensors and sampling their conformational space through rigid-body rotation. We showed with FPMOD that our computational modeling approach can qualitatively predict the FRET efficiencies of a range of biosensors, which had strong agreement with experimental results.     

Gibbs sampling and helix-cap motifs

Protein backbones have characteristic secondary structures, including α-helices and β-sheets. Which structure is adopted locally is strongly biased by the local amino acid sequence of the protein. Accurate (probabilistic) mappings from sequence to structure are valuable for both secondary-structure prediction and protein design. For the case of α-helix caps, we test whether the information content of the sequence–structure mapping can be self-consistently improved by using a relaxed definition of the structure. We derive helix-cap sequence motifs using database helix assignments for proteins of known structure. These motifs are refined using Gibbs sampling in competition with a null motif. Then Gibbs sampling is repeated, allowing for frameshifts of ±1 amino acid residue, in order to find sequence motifs of higher total information content. All helix-cap motifs were found to have good generalization capability, as judged by training on a small set of non-redundant proteins and testing on a larger set. For overall prediction purposes, frameshift motifs using all training examples yielded the best results. Frameshift motifs using a fraction of all training examples performed best in terms of true positives among top predictions. However, motifs without frameshifts also performed well, despite a roughly one-third lower total information content.     

An approach to systematic detection of protein structural motifs

A procedure to detect similar local structures of proteins fromC coordinates is presented. First, the conformations of seven-residuepeptide segments are approximated by a limited number of representatives,each of which is assigned a symbol. Thus, the overall conformationof a protein is represented by a symbol string. The comparisonof these symbol strings using a sequence alignment techniquethen gives pairs of similar local structures. These pairs areconsidered candidates of structural motifs. The applicationof the procedure to the analysis of 93 proteins gave 858 pairsof similar local structures, which included several well-knownstructural motifs such as the nucleotide-binding ßß-unitand the calcium-binding EF hand. The characterization of aminoacid patterns of similar local structures given by the procedureshould be useful for the development of protein structure predictionbased on the acquisition of empirical rules from a large-scaledatabase.     

A new structural framework for integrating replication protein A into DNA processing machinery

By coupling the protection and organization of single-stranded DNA (ssDNA) with recruitment and alignment of DNA processing factors, replication protein A (RPA) lies at the heart of dynamic multi-protein DNA processing machinery. Nevertheless, how RPA coordinates biochemical functions of its eight domains remains unknown. We examined the structural biochemistry of RPA’s DNA-binding activity, combining small-angle X-ray and neutron scattering with all-atom molecular dynamics simulations to investigate the architecture of RPA’s DNA-binding core. The scattering data reveal compaction promoted by DNA binding; DNA-free RPA exists in an ensemble of states with inter-domain mobility and becomes progressively more condensed and less dynamic on binding ssDNA. Our results contrast with previous models proposing RPA initially binds ssDNA in a condensed state and becomes more extended as it fully engages the substrate. Moreover, the consensus view that RPA engages ssDNA in initial, intermediate and final stages conflicts with our data revealing that RPA undergoes two (not three) transitions as it binds ssDNA with no evidence for a discrete intermediate state. These results form a framework for understanding how RPA integrates the ssDNA substrate into DNA processing machinery, provides substrate access to its binding partners and promotes the progression and selection of DNA processing pathways.     

SMoS: a database of structural motifs of protein superfamilies

The Structural Motifs of Superfamilies (SMoS) database provides information about the structural motifs of aligned protein domain superfamilies. Such motifs among structurally aligned multiple members of protein superfamilies are recognized by the conservation of amino acid preference and solvent inaccessibility and are examined for the conservation of other features like secondary structural content, hydrogen bonding, non-polar interaction and residue packing. These motifs, along with their sequence and spatial orientation, represent the conserved core structure of each superfamily and also provide the minimal requirement of sequence and structural information to retain each superfamily fold.     

Human EGF-derived direct and reverse short linear motifs: conformational dynamics insight into the receptor-binding residues

Short linear motifs (SLiMs) have been recognized to perform diverse functions in a variety of regulatory proteins through the involvement in protein–protein interactions, signal transduction, cell cycle regulation, protein secretion, etc. However, detailed molecular mechanisms underlying their functions including roles of definite amino acid residues remain obscure. In our previous studies, we demonstrated that conformational dynamics of amino acid residues in oligopeptides derived from regulatory proteins such as alpha-fetoprotein (AFP), carcino-embryonic antigen (CEA), and pregnancy specific β1-glycoproteins (PSGs) contributes greatly to their biological activities. In the present work, we revealed the 22-member linear modules composed of direct and reverse AFP_14–20-like heptapeptide motifs linked by CxxGY/FxGx consensus motif within epidermal growth factor (EGF), growth factors of EGF family and numerous regulatory proteins containing EGF-like modules. We showed, first, the existence of similarity in amino acid signatures of both direct and reverse motifs in terms of their physicochemical properties. Second, molecular dynamics (MD) simulation study demonstrated that key receptor-binding residues in human EGF in the aligned positions of the direct and reverse motifs may have similar distribution of conformational probability densities and dynamic behavior despite their distinct physicochemical properties. Third, we found that the length of a polypeptide chain (from 7 to 53 residues) has no effect, while disulfide bridging and backbone direction significantly influence the conformational distribution and dynamics of the residues. Our data may contribute to the atomic level structure–function analysis and protein structure decoding; additionally, they may provide a basis for novel protein/peptide engineering and peptide-mimetic drug design.     

Force field based conformational analysis of RNA structural motifs: GNRA tetraloops and their pyrimidine relatives

The protocol of conformational analysis applied here to ribonucleotide oligomers combines conformational search in the space of torsion angles and energy minimization using the AMBER4.1 force field with a continuum treatment of electrostatic solute-solvent interactions. RNA fragments with 5′-GGGCGNNAGCCU-3′ sequences commonly fold into hairpins with four-membered loops. The combinatorial search for acceptable conformations using the MC-SYM program was restricted to loop nucleotides and yielded roughly 1500 structures being compatible with a double-stranded stem. After energy minimization by the JUMNA program (without applying any experimental constraints), these structures converged into an ensemble of 74 different conformers including 26 structures which contained the sheared G-A base pair observed in experimental studies of GNRA tetraloops. Energetic analysis shows that inclusion of solvent electrostatic effects is critically important for the selection of conformers that agree with experimentally determined structures. The continuum model accounts for solvent polarization by means of the electrostatic reaction field. In the case of GNRA loop sequences, the contributions of the reaction field shift relative stabilities towards conformations showing most of the structural features derived from NMR studies. The agreement of computed conformations with the experimental structures of GAAA, GCAA, and GAGA tetraloops suggests that the continuum treatment of the solvent represents a definitive improvement over methods using simple damping models in electrostatic energy calculations. Application of the procedure described here to the evaluation of the relative stabilities of conformers resulting from searching the conformational space of RNA structural motifs provides some progress in (non-homology based) RNA 3D-structure prediction. Received: 20 January 1999 / Revised version: 4 June 1999 / Accepted: 10 June 1999     

LMProt: an efficient algorithm for Monte Carlo sampling of protein conformational space

A new and efficient Monte Carlo algorithm for sampling protein configurations in the continuous space is presented; the efficiency of this algorithm, named Local Moves for Proteins (LMProt), was compared to other alternative algorithms. For this purpose, we used an intrachain interaction energy function that is proportional to the root mean square deviation (rmsd) with respect to alpha-carbons from native structures of real proteins. For phantom chains, the LMProt method is approximately 10(4) and 20 times faster than the algorithms Thrashing (no local moves) and Sevenfold Way (local moves), respectively. Additionally, the LMProt was tested for real chains (excluded-volume all-atoms model); proteins 5NLL (138 residues) and 1BFF (129 residues) were used to determine the folding success xi as a function of the number eta of residues involved in the chain movements, and as a function of the maximum amplitude of atomic displacement delta r(max). Our results indicate that multiple local moves associated with relative chain flexibility, controlled by appropriate adjustments for eta and delta r(max), are essential for configurational search efficiency.     

The fragment transformation method to detect the protein structural motifs

To identify functional structural motifs from protein structures of unknown function becomes increasingly important in recent years due to the progress of the structural genomics initiatives. Although certain structural patterns such as the Asp-His-Ser catalytic triad are easy to detect because of their conserved residues and stringently constrained geometry, it is usually more challenging to detect a general structural motifs like, for example, the betabetaalpha-metal binding motif, which has a much more variable conformation and sequence. At present, the identification of these motifs usually relies on manual procedures based on different structure and sequence analysis tools. In this study, we develop a structural alignment algorithm combining both structural and sequence information to identify the local structure motifs. We applied our method to the following examples: the betabetaalpha-metal binding motif and the treble clef motif. The betabetaalpha-metal binding motif plays an important role in nonspecific DNA interactions and cleavage in host defense and apoptosis. The treble clef motif is a zinc-binding motif adaptable to diverse functions such as the binding of nucleic acid and hydrolysis of phosphodiester bonds. Our results are encouraging, indicating that we can effectively identify these structural motifs in an automatic fashion. Our method may provide a useful means for automatic functional annotation through detecting structural motifs associated with particular functions.     

Strategies for dealing with conformational sampling in structural calculations of flexible or kinked transmembrane peptides.

Peptides corresponding to transmembrane (TM) segments from membrane proteins provide a potential route for the determination of membrane protein structure. We have determined that 2 functionally critical TM segments from the mammalian Na+/H+ exchanger display well converged structure in regions separated by break points. The flexibility of these break points results in conformational sampling in solution. A brief review of available NMR structures of helical membrane proteins demonstrates that there are a number of published structures showing similar properties. Such flexibility is likely indicative of kinks in the full-length protein. This minireview focuses on methods and protocols for NMR structure calculation and analysis of peptide structures under conditions of conformational sampling. The methods outlined allow the identification and analysis of structured peptides containing break points owing to conformational sampling and the differentiation between oligomerization and ensemble-averaged observation of multiple peptide conformations.     

A structural insight into the C-terminal RNA recognition motifs of T-cell intracellular antigen-1 protein

T-cell intracellular antigen-1 (TIA-1) plays a pleiotropic role in cell homeostasis through the regulation of alternative pre-mRNA splicing and mRNA translation by recognising uridine-rich sequences of RNAs. TIA-1 contains three RNA recognition motifs (RRMs) and a glutamine-rich domain. Here, we characterise its C-terminal RRM2 and RRM3 domains. Notably, RRM3 contains an extra novel N-terminal α-helix (α(1)) which protects its single tryptophan from the solvent exposure, even in the two-domain RRM23 context. The α(1) hardly affects the thermal stability of RRM3. On the contrary, RRM2 destabilises RRM3, indicating that both modules are tumbling together, which may influence the RNA binding activity of TIA-1.