Slipknot or Crystallographic Error: A Computational Analysis of the Plasmodium falciparum DHFR Structural Folds

The presence of protein structures with atypical folds in the Protein Data Bank (PDB) is rare and may result from naturally occurring knots or crystallographic errors. Proper characterisation of such folds is imperative to understanding the basis of naturally existing knots and correcting crystallographic errors. If left uncorrected, such errors can frustrate downstream experiments that depend on the structures containing them. An atypical fold has been identified in P. falciparum dihydrofolate reductase (PfDHFR) between residues 20–51 (loop 1) and residues 191–205 (loop 2). This enzyme is key to drug discovery efforts in the parasite, necessitating a thorough characterisation of these folds. Using multiple sequence alignments (MSA), a unique insert was identified in loop 1 that exacerbates the appearance of the atypical fold-giving it a slipknot-like topology. However, PfDHFR has not been deposited in the knotted proteins database, and processing its structure failed to identify any knots within its folds. The application of protein homology modelling and molecular dynamics simulations on the DHFR domain of P. falciparum and those of two other organisms (E. coli and M. tuberculosis) that were used as molecular replacement templates in solving the PfDHFR structure revealed plausible unentangled or open conformations of these loops. These results will serve as guides for crystallographic experiments to provide further insights into the atypical folds identified.     

Importance of an extreme C-terminal motif of a family I.3 lipase for stability

A five-residue sequence motif (VTLVG) located at positions 15-19 from the C-terminus of family I.3 lipase from Pseudomonas sp. MIS38 (PML) and an extreme C-terminal motif (DGIVIA) located at the C-terminus of PML are relatively well conserved in the passenger proteins of type 1 secretion system (T1SS). To analyze the role of these motifs, four mutant proteins of PML (PMLΔ5, PMLΔ10, 3A-PML and 2A-PML) were constructed. PMLΔ5 and PMLΔ10 lack the C-terminal 5 and 10 residues of PML, respectively. 3A-PML has triple mutations within an extreme C-terminal motif and 2A-PML has double mutations within a five-residue sequence motif. Secretion of these proteins was analyzed using Escherichia coli DH5 cells carrying Lip system (T1SS for family I.3 lipase). The secretion level of 2A-PML was dramatically reduced when compared with that of PML, whereas the secretion level of 3A-PML was comparable to that of PML, indicating that a five-residue sequence motif, instead of an extreme C-terminal motif, is required for secretion of PML. None of the mutations and truncations seriously affects the enzymatic activity of PML. However, 3A-PML, PMLΔ5 and PMLΔ10 were less stable than PML by 2.1, 7.6 and 7.6°C in T(1/2), respectively, and by 5.0, 21.3 and 17.9 kJ/mol in ΔG(H(2)O), respectively. These results indicate that an extreme C-terminal motif of PML is important for stability.     

Rational design and selection of bivalent peptide ligands of thrombin incorporating P4-P1 tetrapeptide sequences: from good substrates to potent inhibitors

The tetrapeptide Phe-Asn-Pro-Arg is a structurally optimized sequence for binding to the active site of thrombin. By conjugating this tetrapeptide or some variants to a C-terminal fragment of hirudin, we were able to generate a series of new bivalent inhibitors of thrombin containing only genetically encodable natural amino acids. We found that synergistic binding to both the active site and an exosite of thrombin can be enhanced through substitutions of amino acid residues at the P3 and P3' sites of the active-site directed sequence, Phe(P4)-Xaa(P3)-Pro(P2)-Arg(P1)-Pro(P1')-Gln(P2')-Yaa(P3'). Complementary to rational design, a phage library was constructed to explore further the residue requirements at the P4, P3 and P3' sites for bivalent and optimized two-site binding. Very significantly, panning of the phage library has led to thrombin-inhibitory peptides possessing strong anti-clotting activities in the low nanomolar range and yet interfering only partially the catalytic active site of thrombin. Modes of action of the newly discovered bivalent inhibitors are rationalized in light of the allosteric properties of thrombin, especially the interplay between the proteolytic action and regulatory binding occurring at thrombin surfaces remote from the catalytic active site.     

Aggregation Condition–Structure Relationship of Mouse Prion Protein Fibrils

Prion diseases are associated with conformational conversion of cellular prion protein into a misfolded pathogenic form, which resembles many properties of amyloid fibrils. The same prion protein sequence can misfold into different conformations, which are responsible for variations in prion disease phenotypes (prion strains). In this work, we use atomic force microscopy, FTIR spectroscopy and magic-angle spinning NMR to devise structural models of mouse prion protein fibrils prepared in three different denaturing conditions. We find that the fibril core region as well as the structure of its N- and C-terminal parts is almost identical between the three fibrils. In contrast, the central part differs in length of β-strands and the arrangement of charged residues. We propose that the denaturant ionic strength plays a major role in determining the structure of fibrils obtained in a particular condition by stabilizing fibril core interior-facing glutamic acid residues.     

Role of paired basic residues of protein C-termini in phospholipid binding

It is a well known phenomenon that the occurrence of severaldistinct amino acids at the C-terminus of proteins is non-random.We have analysed all Saccharomyces cerevisiae proteins predictedby computer databases and found lysine to be the most frequentresidue both at the last (–1) and at the penultimate aminoacid (–2) positions. To test the hypothesis that C-terminalbasic residues efficiently bind to phospholipids we randomlyexpressed GST-fusion proteins from a yeast genomic library.Fifty-four different peptide fragments were found to bind phospholipidsand 40% of them contained lysine/arginine residues at the (–1)or (–2) positions. One peptide showed high sequence similaritywith the yeast protein Sip18p. Mutational analysis revealedthat both C-terminal lysine residues of Sip18p are essentialfor phospholipid-binding in vitro. We assume that basic aminoacid residues at the (–1) and (–2) positions inC-termini are suitable to attach the C-terminus of a given proteinto membrane components such as phospholipids, thereby stabilizingthe spatial structure of the protein or contributing to itssubcellular localization. This mechanism could be an additionalexplanation for the C-terminal amino acid bias observed in proteinsof several species.     

Intersubunit disulfide-bonded {lambda}-Cro protein

Site-directed mutagenesis has been employed to substitute cysteinefor valine at position 55, which is located on the dimer interfaceof the Cro protein of bacteriophage . It has been found thatthe Cys55 Cro protein (Cro VC55) spontaneously forms a stabledisulfide-bonded dimer in the absence of a reducing agent. UV—CDand NMR data showed that the mutant protein retains the conformationof the wild Cro protein and has acquired significant heat-stability.However, its specific DNA-binding activity is reduced severaltimes compared with that of the wild Cro. Photochemically induceddynamic nuclear polarization (CIDNP) spectra demonstrated thata conformational change of Cro VC55 did not take place uponthe formation of a complex with O_R3, in contrast to the caseof the wild Cro. These data suggest that the induced fitting,like loosening, of the two subunits of the wild Cro dimer contributesto the enhancement of its affinity to its operator DNA, whichresults in a specific interaction between Cro and O_R3.     

Sequence and structural analysis of two designed proteins with 88% identity adopting different folds

Protein folding is a natural phenomenon by which a sequence of amino acids folds into a unique functional three-dimensional structure. Although the sequence code that governs folding remains a mystery, one can identify key inter-residue contacts responsible for a given topology. In nature, there are many pairs of proteins of a given length that share little or no sequence identity. Similarly, there are many proteins that share a common topology but lack significant evidence of homology. In order to tackle this problem, protein engineering studies have been used to determine the minimal number of amino acid residues that codes for a particular fold. In recent years, the coupling of theoretical models and experiments in the study of protein folding has resulted in providing some fruitful clues. He et al. have designed two proteins with 88% sequence identity, which adopt different folds and functions. In this work, we have systematically analysed these two proteins by performing pentapeptide search, secondary structure predictions, variation in inter-residue interactions and residue-residue pair preferences, surrounding hydrophobicity computations, conformational switching and energy computations. We conclude that the local secondary structural preference of the two designed proteins at the Nand C-terminal ends to adopt either coil or strand conformation may be a crucial factor in adopting the different folds. Early on during the process of folding, both proteins may choose different energetically favourable pathways to attain the different folds.     

Unfolding events of Chymotrypsin Inhibitor 2 (CI2) revealed by Monte Carlo (MC) simulations and their consistency from structure-based analysis of conformations

We investigated the behavior of the conformations of Chymotrypsin Inhibitor (CI2) from the native to the denatured states, obtained in Monte Carlo (MC)/Metropolis simulations, where a low-resolution model is used together with knowledge-based potentials. New conformations starting from the X-ray native structure are generated by random perturbations along with a constraint to increase the radius of gyration. Unfolding is also simulated by unrestrained simulations at a higher temperature. All simulations yield a similar sequence of unfolding events. The preferred pathway starts with loss of native contacts between (N-terminal)-β₃ and continues with β₂-β₃. The persistence of the contacts between β₁ and β₂ at intermediate values of the fraction of native contacts (Q); whereas, highly unfolded conformations with only some helical contacts persisting at low values of Q, are observed. Structure-based analysis of the fluctuations of the unfolded conformations by Gaussian Network Model (GNM) reveals that the termini of the chain—C terminus being more mobile—depict relatively higher flexibility with a native-like hinge near β₂ that divides the structure into two domains. The fluctuations of the two domains are negatively correlated, with partly folded α-helix and a small hydrophobic cluster in the middle of the chain displaying positively correlated fluctuations. The most persistent short-range rotational bond correlations are observed between the residues of α-helix, C terminus of the β₁-part of the reactive site loop, and around the C terminus of the β₂. The latter regions also appear as hot spots; i.e. high frequency fluctuating regions, of the structure surviving in unfolded conformations. The results imply that the unfolded CI2 has an intrinsic ability to undergo correlated fluctuations along with some residual native structure specifically induced by its sequence, consisting at the lowest level of a single hinge.     

The structural plasticity of the C terminus of p21Cip1 is a determinant for target protein recognition

The cyclin-dependent kinase inhibitory protein p21(Cip1) might play multiple roles in cell-cycle regulation through interaction of its C-terminal domain with a defined set of cellular proteins such as proliferating cell nuclear antigen (PCNA), calmodulin (CaM), and the oncoprotein SET. p21(Cip1) could be described as an intrinsically unstructured protein in solution although the C-terminal domain adopts a well-defined extended conformation when bound to PCNA. However, the molecular mechanism of the interaction with CaM and the oncoprotein SET is not well understood, partly because of the lack of structural information. In this work, a peptide derived from the C-terminal domain of p21(Cip1) that covers the binding domain of the three above-mentioned proteins was used to demonstrate that the C-terminal domain of p21 recognizes multiple ligands through its ability to adopt multiple conformations. The conformation is dictated by tertiary contacts rather than by the primary sequence of the protein. Our results suggest that the C-terminal domain of p21(Cip1) adopts an extended structure when bound to PCNA and probably when bound to the oncoprotein SET, but an alpha helix when bound to CaM.     

Computational analysis of alpha-helical membrane protein structure: implications for the prediction of 3D structural models

Relatively little has been known about the structure of alpha-helical membrane proteins, since until recently few structures had been crystallized. These limited data have restricted structural analyses to the prediction of secondary structure, rather than tertiary folds. In order to address this, this paper describes an analysis of the 23 available membrane protein structures. A number of findings are made that are of particular relevance to transmembrane helix packing: (1) on average lipid-tail-accessible transmembrane residues are significantly more hydrophobic, less conserved and contain different residue types to buried residues; (2) charged residues are not always buried and, when accessible to membrane lipid tails, few are paired with another charge and instead they often interact with phospholipid head-groups or with other residue types; (3) a significant proportion of lipid-tail-accessible charged and polar residues form hydrogen bonds only with residues one turn away in the same helix (intra-helix); (4) pore-lining residues are usually hydrophobic and it is difficult to distinguish them from buried residues in terms of either residue type or conservation; and (5) information was gained about the proportion of helices that tend to contribute to lining a pore and the resulting pore diameter. These findings are discussed with relevance to the prediction of membrane protein 3D structure.     

Sequence-structure specificity--how does an inverse folding approach work?

The inverse folding approach is a powerful tool in protein structure prediction when the native state of a sequence adopts one of the known protein folds. This is because some proteins show strong sequence- structure specificity in inverse folding experiments that allow gaps and insertions in the sequence-structure alignment. In those cases when structures similar to their native folds are included in the structure database, the z-scores (which measure the sequence-structure specificity) of these folds are well separated from those of other alternative structures. In this paper, we seek to understand the origin of this sequence-structure specificity and to identify how the specificity arises on passing from a short peptide chain to the entire protein sequence. To accomplish this objective, a simplified version of inverse folding, gapless inverse folding, is performed using sequence fragments of different sizes from 53 proteins. The results indicate that usually a significant portion of the entire protein sequence is necessary to show sequence-structure specificity, but there are regions in the sequence that begin to show this specificity at relatively short fragment size (15-20 residues). An island picture, in which the regions in the sequence that recognize their own native structure grow from some seed fragments, is observed as the fragment size increases. Usually, more similar structures to the native states are found in the top-scoring structural fragments in these high-specificity regions.      

The crystal structure of the ubiquitin-like (UbL) domain of human homologue A of Rad23 (hHR23A) protein

The human homologue of the yeast Rad23 protein, hHR23A, plays dual roles in DNA repair as well as in translocating polyubiquitinated proteins to the proteasome. We determined the three-dimensional structure of its ubiquitin-like (UbL) domain by X-ray crystallography. It has the same overall structure and fold characteristics as ubiquitin and other members of the UbL domain family, with overall root mean square deviations in Cα positions in the range of 1.0-1.3 ?. There are local differences in the α1-β3 loop where hHR23A UbL domain has three more residues constituting a bigger loop. Analysis of the crystal packing revealed a possible dimeric arrangement mediated by the three residues (Leu10, Ile49 and Met75) that are known to be critical for molecular interactions. In contrast to the overall well-defined structure, these three residues are either disordered or have multiple conformations, suggesting that conformation variability is an important property of the binding surface. The electrostatic potentials at the binding surface are conserved among the family, with the hHR23B domain being the most similar to this structure. The intra-molecular complexes formed by the UbL domain of hHR23A with its UbA1 or UbA2 domains was studied by comparative homology modelling, which suggests these two interactions are structurally similar and are mutually exclusive.     

Expression, purification and characterization of the homeodomain of rat ISL-1 protein

Isl-1 is a member of a family of Homeodomains containing proteins that possess an N-terminal pair of zinc binding LIM domains. The Isl-1 gene in rat codes for a protein that binds to the insulin gene enhancer and is also involved in regulation of amylin and proglucagon genes. A DNA sequence coding for 66 amino acid residues containing the C-terminal homeodomain fragment of Isl-1 was expressed as a soluble protein in Escherichia coli. Here, we describe a procedure which allows the rapid native purification of recombinant homeodomain protein fused to an N- terminal tag of six histidines. The purified homeodomain showed DNA- binding activity to its cognate DNA sequence. An enhanced binding activity is observed in the presence of a reducing agent in electrophoretic mobility shift assays. The DNA binding was further characterized by circular dichroism spectroscopy. Addition of DNA to the homeodomain did not change the overall secondary structure content, but the thermal and chemical denaturing profiles were altered. A stabilization of the secondary structure was observed upon DNA binding. The free energy of unfolding at 23 degrees C was 7 kJ mol(-1) in absence of DNA and 29 kJ mol(-1) in the presence of DNA.      

Anchoring Intrinsically Disordered Proteins to Multiple Targets: Lessons from N-Terminus of the p53 Protein

Anchor residues, which are deeply buried upon binding, play an important role in protein-protein interactions by providing recognition specificity and facilitating the binding kinetics. Up to now, studies on anchor residues have been focused mainly on ordered proteins. In this study, we investigated anchor residues in intrinsically disordered proteins (IDPs) which are flexible in the free state. We identified the anchor residues of the N-terminus of the p53 protein (Glu17-Asn29, abbreviated as p53N) which are involved in binding with two different targets (MDM2 and Taz2), and analyzed their side chain conformations in the unbound states. The anchor residues in the unbound p53N were found to frequently sample conformations similar to those observed in the bound complexes (i.e., Phe19, Trp23, and Leu26 in the p53N-MDM2 complex, and Leu22 in the p53N-Taz2 complex). We argue that the bound-like conformations of the anchor residues in the unbound state are important for controlling the specific interactions between IDPs and their targets. Further, we propose a mechanism to account for the binding promiscuity of IDPs in terms of anchor residues and molecular recognition features (MoRFs).     

Recognition of analogous and homologous protein folds--assessment of prediction success and associated alignment accuracy using empirical substitution matrices

Fold recognition methods aim to use the information in the known protein structures (the targets) to identify that the sequence of a protein of unknown structure (the probe) will adopt a known fold. This paper highlights that the structural similarities sought by these methods can be divided into two types: remote homologues and analogues. Homologues are the result of divergent evolution and often share a common function. We define remote homologues as those that are not easily detectable by sequence comparison methods alone. Analogues do not have a common ancestor and generally do not have a common function. Several sets of empirical matrices for residue substitution, secondary structure conservation and residue accessibility conservation have previously been derived from aligned pairs of remote homologues and analogues (Russell et al., J. Mol. Biol., 1997, 269, 423-439). Here a method for fold recognition, FOLDFIT, is introduced that uses these matrices to match the sequences, secondary structures and residue accessibilities of the probe and target. The approach is evaluated on distinct datasets of analogous and remotely homologous folds. The accuracy of FOLDFIT with the different matrices on the two datasets is contrasted to results from another fold recognition method (THREADER) and to searches using mutation matrices in the absence of any structural information. FOLDFIT identifies at top rank 12 out of 18 remotely homologous folds and five out of nine analogous folds. The average alignment accuracies for residue and secondary structure equivalencing are much higher for homologous folds (residue approximately 42%, secondary structure approximately 78%) than for analogues folds (approximately 12%, approximately 47%). Sequence searches alone can be successful for several homologues in the testing sets but nearly always fail for the analogues. These results suggest that the recognition of analogous and remotely homologous folds should be assessed separately. This study has implications for the development and comparative evaluation of fold recognition algorithms.      

Prediction of protein residue contacts with a PDB-derived likelihood matrix

Proteins with similar folds often display common patterns ofresidue variability. A widely discussed question is how thesepatterns can be identified and deconvoluted to predict proteinstructure. In this respect, correlated mutation analysis (CMA)has shown considerable promise. CMA compares multiple membersof a protein family and detects residues that remain constantor mutate in tandem. Often this behavior points to structuralor functional interdependence between residues. CMA has beenused to predict pairs of amino acids that are distant in theprimary sequence but likely to form close contacts in the nativethree-dimensional structure. Until now these methods have usedevolutionary or biophysical models to score the fit betweenresidues. We wished to test whether empirical methods, derivedfrom known protein structures, would provide useful predictivepower for CMA. We analyzed 672 known protein structures, derivedcontact likelihood scores for all possible amino acid pairs,and used these scores to predict contacts. We then tested themethod on 118 different protein families for which structureshave been solved to atomic resolution. The mean performancewas almost seven times better than random prediction. Used inconcert with secondary structure prediction, the new CMA methodcould supply restraints for predicting still undetermined structures.     

Utilization of AlphaFold2 to Predict MFS Protein Conformations after Selective Mutation

The major facilitator superfamily (MFS) is the largest secondary transporter family and is responsible for transporting a broad range of substrates across the biomembrane. These proteins are involved in a series of conformational changes during substrate transport. To decipher the transport mechanism, it is necessary to obtain structures of these different conformations. At present, great progress has been made in predicting protein structure based on coevolutionary information. In this study, AlphaFold2 was used to predict different conformational structures for 69 MFS transporters of E. coli after the selective mutation of residues at the interface between the N- and C-terminal domains. The predicted structures for these mutants had small RMSD values when compared to structures obtained using X-ray crystallography, which indicates that AlphaFold2 predicts the structure of MSF transporters with high accuracy. In addition, different conformations of other transporter family proteins have been successfully predicted based on mutation methods. This study provides a structural basis to study the transporting mechanism of the MFS transporters and a method to probe dynamic conformation changes of transporter family proteins when performing their function.