What is the relationship between the global structures of apo and holo proteins?

It is well known that ligand binding and release may induce a wide range of structural changes in a receptor protein, varying from small movements of loops or side chains in the binding pocket to large‐scale domain hinge‐bending and shear motions or even partial unfolding that facilitates the capture and release of a ligand. An interesting question is what in general are the conformational changes triggered by ligand binding? The aim of this work is analyze the magnitude of structural changes in a protein resulting from ligand binding to assess if the state of ligand binding needs to be included in template‐based protein structure prediction algorithms. To address this issue, a nonredundant dataset of 521 paired protein structures in the ligand‐free and ligand‐bound form was created and used to estimate the degree of both local and global structure similarity between the apo and holo forms. In most cases, the proteins undergo relatively small conformational rearrangements of their tertiary structure upon ligand binding/release (most root‐mean‐square‐deviations from native, RMSD, are <1 Å). However, a clear difference was observed between single‐ and multiple‐domain proteins. For the latter, RMSD changes greater than 1 Å and sometimes larger were found for almost 1/3 of the cases; these are mainly associated with large‐scale hinge‐bending movements of entire domains. The changes in the mutual orientation of individual domains in multiple‐domain proteins upon ligand binding were investigated using a mechanistic model based on mass‐weighted principal axes as well as interface buried surface calculations. Some preferences toward the anticipated mechanism of protein domain movements are predictable based on the examination of just the ligand‐free structural form. These results have applications to protein structure prediction, particularly in the context of protein domain assembly, if additional information concerning ligand binding is exploited. Proteins 2008. © 2007 Wiley‐Liss, Inc.     

Model-free methods of analyzing domain motions in proteins from simulation: A comparison of normal mode analysis and molecular dynamics simulation of lysozyme

Model-free methods are introduced to determine quantities pertaining to protein domain motions from normal mode analyses and molecular dynamics simulations. For the normal mode analysis, the methods are based on the assumption that in low frequency modes, domain motions can be well approximated by modes of motion external to the domains. To analyze the molecular dynamics trajectory, a principal component analysis tailored specifically to analyze interdomain motions is applied. A method based on the curl of the atomic displacements is described, which yields a sharp discrimination of domains, and which defines a unique interdomain screw-axis. Hinge axes are defined and classified as twist or closure axes depending on their direction. The methods have been tested on lysozyme. A remarkable correspondence was found between the first normal mode axis and the first principal mode axis, with both axes passing within 3 Å of the alpha-carbon atoms of residues 2, 39, and 56 of human lysozyme, and near the interdomain helix. The axes of the first modes are overwhelmingly closure axes. A lesser degree of correspondence is found for the second modes, but in both cases they are more twist axes than closure axes. Both analyses reveal that the interdomain connections allow only these two degrees of freedom, one more than provided by a pure mechanical hinge. Proteins 27:425–437, 1997. © 1997 Wiley-Liss, Inc.     

The S1 ribosomal protein family contains a unique conservative domain

The number of amino acid residues contained in the S1 ribosomal protein of various bacteria varies in a wide range: from 111 to 863 residues in Spiroplasma kunkelii and Treponema pallidum, respectively. The architecture of this protein is traditionally (in particular, because of unknown spatial structure) represented as repeated S1 domains, the copy number of which depends on the protein length. The data on the copy number and boundaries of these domains is available in specialized databases, such as SMART, Pfam, and PROSITE; however, these data can be rather different for the same object. In this work, we used the approach utilizing analysis of predicted secondary structure (PsiPred program). This allowed us to detect the structural domains in S1 protein sequences; their copy number varied from one to six. Alignment of the S1 proteins containing different numbers of domains with the S1 RNA-binding domain of Escherichia coli polynucleotide phosphorylase provided for discovering a domain within this family displaying the maximal homology to the E. coli domain. This conservative domain migrates along the chain, and its location in the proteins with different numbers of domains follows a certain pattern. Similar to the S1 domain of polynucleotide phosphorylase, residues Phe19, Phe22, His34, Asp64, and Arg68 in this conservative domain are clustered on the surface to form an RNA-binding site.     

Concerted structural changes in the peptidase and the propeller domains of prolyl oligopeptidase are required for substrate binding

Prolyl oligopeptidase contains a peptidase domain and its catalytic triad is covered by the central tunnel of a seven-bladed beta-propeller. This domain makes the enzyme an oligopeptidase by excluding large structured peptides from the active site. The apparently rigid crystal structure does not explain how the substrate can approach the catalytic groups. Two possibilities of substrate access were investigated: either blades 1 and 7 of the propeller domain move apart, or the peptidase and/or propeller domains move to create an entry site at the domain interface. Engineering disulfide bridges to the expected oscillating structures prevented such movements, which destroyed the catalytic activity and precluded substrate binding. This indicated that concerted movements of the propeller and the peptidase domains are essential for the enzyme action.     

The dynamic cycle of bacterial translation initiation factor IF3

Initiation factor IF3 is an essential protein that enhances the fidelity and speed of bacterial mRNA translation initiation. Here, we describe the dynamic interplay between IF3 domains and their alternative binding sites using pre-steady state kinetics combined with molecular modelling of available structures of initiation complexes. Our results show that IF3 accommodates its domains at velocities ranging over two orders of magnitude, responding to the binding of each 30S ligand. IF1 and IF2 promote IF3 compaction and the movement of the C-terminal domain (IF3C) towards the P site. Concomitantly, the N-terminal domain (IF3N) creates a pocket ready to accept the initiator tRNA. Selection of the initiator tRNA is accompanied by a transient accommodation of IF3N towards the 30S platform. Decoding of the mRNA start codon displaces IF3C away from the P site and rate limits translation initiation. 70S initiation complex formation brings IF3 domains in close proximity to each other prior to dissociation and recycling of the factor for a new round of translation initiation. Altogether, our results describe the kinetic spectrum of IF3 movements and highlight functional transitions of the factor that ensure accurate mRNA translation initiation.     

Ligand-free and -bound structures of the binding protein (LivJ) of the Escherichia coli ABC leucine/isoleucine/valine transport system: trajectory and dynamics of the interdomain rotation and ligand specificity

The leucine/isoleucine/valine-binding protein (LIVBP or LivJ) serves as the primary high-affinity receptor of the Escherichia coli ABC-type transporter for the three aliphatic amino acids. The first structure of LIVBP determined previously without bound ligand showed a molecule comprised of two domains which are far apart and bisected by a wide open, solvent-accessible cleft. Here we report the crystal structures of another ligand-free state and three complexes with the aliphatic amino acids. In the present ligand-free structure, the two domains are farther apart. In the three very similar complex structures, the two domains are in close proximity, and each desolvated ligand is completely engulfed in the cleft and bound by both domains. The two different ligand-free structures, combined with those of the very similar ligand-bound structures, indicate the trajectory and backbone torsion angle changes of the hinges that accompany domain closure and play crucial functional roles. The amino acids are bound by polar and nonpolar interactions, occurring predominantly in one domain. Consistent with the protein specificity, the aliphatic side chains of the ligands lie in a hydrophobic pocket fully formed following domain or cleft closure. Comparison of the structures of LIVBP with several different binding proteins indicates no correlations between the magnitudes of the hinge-bending angles and the protein masses, the ligand sizes, or the number of segments connecting the two domains. Results of normal-mode analysis and molecular dynamics simulations are consistent with the trajectory and intrinsic flexibility of the interdomain hinges and the dominance of one domain in ligand binding in the open state.     

Essential dynamics of helices provide a functional classification of EF-hand proteins

Low energy modes have been calculated for the largest possible number of available representatives (>150) of EF-hand domains belonging to different members of the calcium-binding EF-hand protein superfamily. These proteins are the major actors in signal transduction. The latter, in turn, relies on the dynamical properties of the systems, in particular on the relative movements of the four helices characterizing each EF-hand domain upon calcium binding. The peculiar structural and dynamical features of this protein superfamily are systematically investigated by a novel approach, where the lowest energy (essential) modes are described in the space of the six interhelical angles among the four helices constituting the EF-hand domain. The modes, obtained through a general and transferable coarse-graining scheme, identify the easy directions of helical motions. It is found that, for most proteins, the two lowest energy modes are sufficient to capture most of the helices' fluctuation dynamics. Strikingly, the comparison of such modes for all possible pairs of EF-hand domain representatives reveals that only few easy directions are preferred within this large protein superfamily. This enables us to introduce a novel dynamics-based classification of EF-hand domains that complements existing structure-based characterizations from an unexplored biological perspective.     

Proteome Evolution and the Metabolic Origins of Translation and Cellular Life

The origin of life has puzzled molecular scientists for over half a century. Yet fundamental questions remain unanswered, including which came first, the metabolic machinery or the encoding nucleic acids. In this study we take a protein-centric view and explore the ancestral origins of proteins. Protein domain structures in proteomes are highly conserved and embody molecular functions and interactions that are needed for cellular and organismal processes. Here we use domain structure to study the evolution of molecular function in the protein world. Timelines describing the age and function of protein domains at fold, fold superfamily, and fold family levels of structural complexity were derived from a structural phylogenomic census in hundreds of fully sequenced genomes. These timelines unfold congruent hourglass patterns in rates of appearance of domain structures and functions, functional diversity, and hierarchical complexity, and revealed a gradual build up of protein repertoires associated with metabolism, translation and DNA, in that order. The most ancient domain architectures were hydrolase enzymes and the first translation domains had catalytic functions for the aminoacylation and the molecular switch-driven transport of RNA. Remarkably, the most ancient domains had metabolic roles, did not interact with RNA, and preceded the gradual build-up of translation. In fact, the first translation domains had also a metabolic origin and were only later followed by specialized translation machinery. Our results explain how the generation of structure in the protein world and the concurrent crystallization of translation and diversified cellular life created further opportunities for proteomic diversification.     

A method for the analysis of domain movements in large biomolecular complexes

A new method for the analysis of domain movements in large, multichain, biomolecular complexes is presented. The method is applicable to any molecule for which two atomic structures are available that represent a conformational change indicating a possible domain movement. The method is blind to atomic bonding and atom type and can, therefore, be applied to biomolecular complexes containing different constituent molecules such as protein, RNA, or DNA. At the heart of the method is the use of blocks located at grid points spanning the whole molecule. The rotation vector for the rotation of atoms from each block between the two conformations is calculated. Treating components of these vectors as coordinates means that each block is associated with a point in a “rotation space” and that blocks with atoms that rotate together, perhaps as part of the same rigid domain, will have colocated points. Thus a domain can be identified from the clustering of points from blocks that span it. Domain pairs are accepted for analysis of their relative movements in terms of screw axes based upon a set of reasonable criteria. Here, we report on the application of the method to biomolecules covering a considerable size range: hemoglobin, liver alcohol dehydrogenase, S‐Adenosylhomocysteine hydrolase, aspartate transcarbamylase, and the 70S ribosome. The results provide a depiction of the conformational change within each molecule that is easily understood, giving a perspective that is expected to lead to new insights. Of particular interest is the allosteric mechanism in some of these molecules. Results indicate that common boundaries between subunits and domains are good regions to focus on as movement in one subunit can be transmitted to another subunit through such interfaces. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Fold Evolution before LUCA: Common Ancestry of SH3 Domains and OB Domains

SH3 and OB are the simplest, oldest, and most common protein domains within the translation system. SH3 and OB domains are β-barrels that are structurally similar but are topologically distinct. To transform an OB domain to a SH3 domain, β-strands must be permuted in a multistep and evolutionarily implausible mechanism. Here, we explored relationships between SH3 and OB domains of ribosomal proteins, initiation, and elongation factors using a combined sequence- and structure-based approach. We detect a common core of SH3 and OB domains, as a region of significant structure and sequence similarity. The common core contains four β-strands and a loop, but omits the fifth β-strand, which is variable and is absent from some OB and SH3 domain proteins. The structure of the common core immediately suggests a simple permutation mechanism for interconversion between SH3 and OB domains, which appear to share an ancestor. The OB domain was formed by duplication and adaptation of the SH3 domain core, or vice versa, in a simple and probable transformation. By employing the folding algorithm AlphaFold2, we demonstrated that an ancestral reconstruction of a permuted SH3 sequence folds into an OB structure, and an ancestral reconstruction of a permuted OB sequence folds into a SH3 structure. The tandem SH3 and OB domains in the universal ribosomal protein uL2 share a common ancestor, suggesting that the divergence of these two domains occurred before the last universal common ancestor.     

The carboxy-terminal domain of the DExDH protein YxiN is sufficient to confer specificity for 23S rRNA

DE x DH proteins are believed to modulate the structures of RNAs and ribonucleoprotein complexes by disrupting RNA helices and RNA-protein interactions. All DE x DH proteins contain a two-domain catalytic core that enables their RNA-dependent ATPase and RNA helicase activities. The catalytic core may be flanked by ancillary domains that are proposed to confer substrate specificity and facilitate the unique functions of individual proteins. The Escherichia coli DE x DH protein DbpA and its Bacillus subtilis ortholog YxiN have similar 75aa carboxy-terminal domains, and both proteins are specifically targeted to 23S rRNA. Here we demonstrate that the carboxy-terminal domain of YxiN is sufficient to confer RNA specificity by characterizing a chimera in which this domain is appended to the core domains of E.coli SrmB, a DE x DH protein with no apparent substrate specificity. Both the RNA-dependent ATPase and RNA helicase activities of the chimera are specifically activated by 23S rRNA and abolished by sequence changes within hairpin 92, a critical recognition element for Y x iN. These data support a model in which the carboxy-terminal domain binds hairpin 92 to target the protein to 23S rRNA.     

Intron/exon structure of the human gene for the muscle isozyme of glycogen phosphorylase

The intron/exon organization of the human gene for glycogen phosphorylase has been determined. The segments of the polypeptide chain that corresponds to the 19 exons of the gene are examined for relationships between the three-dimensional structure to the protein and gene structure. Only weak correlations are observed between domains of phosphorylase and exons. The nucleotide binding domains that are found in phosphorylase and other glycolytic enzymes are examined for relationships between exons of the genes and structures of the domains. When mapped to the three-dimensional structures, the intron/exon boundaries are shown to be widely distributed in this family of protein domains.     

Computer-aided NMR assay for detecting natively folded structural domains

Structural genomics projects require strategies for rapidly recognizing protein sequences appropriate for routine structure determination. For large proteins, this strategy includes the dissection of proteins into structural domains that form stable native structures. However, protein dissection essentially remains an empirical and often a tedious process. Here, we describe a simple strategy for rapidly identifying structural domains and assessing their structures. This approach combines the computational prediction of sequence regions corresponding to putative domains with an experimental assessment of their structures and stabilities by NMR and biochemical methods. We tested this approach with nine putative domains predicted from a set of 108 Thermus thermophilus HB8 sequences using PASS, a domain prediction program we previously reported. To facilitate the experimental assessment of the domain structures, we developed a generic 6-hour His-tag-based purification protocol, which enables the sample quality evaluation of a putative structural domain in a single day. As a result, we observed that half of the predicted structural domains were indeed natively folded, as judged by their HSQC spectra. Furthermore, two of the natively folded domains were novel, without related sequences classified in the Pfam and SMART databases, which is a significant result with regard to the ability of structural genomics projects to uniformly cover the protein fold space.     

Prediction of protein domain boundaries from sequence alone

We present here a simple approach to identify domain boundaries in proteins of an unknown three-dimensional structure. Our method is based on the hypothesis that a high-side chain entropy of a region in a protein chain must be compensated by a high-residue interaction energy within the region, which could correlate with a well-structured part of the globule, that is, with a domain unit. For protein domains, this means that the domain boundary is conditioned by amino acid residues with a small value of side chain entropy, which correlates with the side chain size. On the one hand, relatively high Ala and Gly content on the domain boundary results in high conformational entropy of the backbone chain between the domains. On the other hand, the presence of Pro residues leads to the formation of hinges for a relative orientation of domains. The method was applied to 646 proteins with two contiguous domains extracted from the SCOP database with a success rate of 63%. We also report the prediction of domain boundaries for CASP5 targets obtained with the same method.     

Multiple phosphate positions in the catalytic site of glycogen phosphorylase: structure of the pyridoxal-5''-pyrophosphate coenzyme-substrate analog.

The three-dimensional structure of an R-state conformer of glycogen phosphorylase containing the coenzyme-substrate analog pyridoxal-5'-diphosphate at the catalytic site (PLPP-GPb) has been refined by X-ray crystallography to a resolution of 2.87 A. The molecule comprises four subunits of phosphorylase related by approximate 222 symmetry. Whereas the quaternary structure of R-state PLPP-GPb is similar to that of phosphorylase crystallized in the presence of ammonium sulfate (Barford, D. & Johnson, L.N., 1989, Nature 340, 609-616), the tertiary structures differ in that the two domains of the PLPP-GPb subunits are rotated apart by 5 degrees relative to the T-state conformation. Global differences among the four subunits suggest that the major domains of the phosphorylase subunit are connected by a flexible hinge. The two different positions observed for the terminal phosphate of the PLPP are interpreted as distinct phosphate subsites that may be occupied at different points along the reaction pathway. The structural basis for the unique ability of R-state dimers to form tetramers results from the orientation of subunits with respect to the dyad axis of the dimer. Residues in opposing dimers are in proper registration to form tetramers only in the R-state.     

The role of internal duplication in the evolution of multi-domain proteins

Many proteins consist of several structural domains. These multi-domain proteins have likely been generated by selective genome growth dynamics during evolution to perform new functions as well as to create structures that fold on a biologically feasible time scale. Domain units frequently evolved through a variety of genetic shuffling mechanisms. Here we examine the protein domain statistics of more than 1000 organisms including eukaryotic, archaeal and bacterial species. The analysis extends earlier findings on asymmetric statistical laws for proteome to a wider variety of species. While proteins are composed of a wide range of domains, displaying a power-law decay, the computation of domain families for each protein reveals an exponential distribution, characterizing a protein universe composed of a thin number of unique families. Structural studies in proteomics have shown that domain repeats, or internal duplicated domains, represent a small but significant fraction of genome. In spite of its importance, this observation has been largely overlooked until recently. We model the evolutionary dynamics of proteome and demonstrate that these distinct distributions are in fact rooted in an internal duplication mechanism. This process generates the contemporary protein structural domain universe, determines its reduced thickness, and tames its growth. These findings have important implications, ranging from protein interaction network modeling to evolutionary studies based on fundamental mechanisms governing genome expansion.     

Structure of PA1221, a nonribosomal peptide synthetase containing adenylation and peptidyl carrier protein domains

Many bacteria use large modular enzymes for the synthesis of polyketide and peptide natural products. These multidomain enzymes contain integrated carrier domains that deliver bound substrates to multiple catalytic domains, requiring coordination of these chemical steps. Nonribosomal peptide synthetases (NRPSs) load amino acids onto carrier domains through the activity of an upstream adenylation domain. Our lab recently determined the structure of an engineered two-domain NRPS containing fused adenylation and carrier domains. This structure adopted a domain-swapped dimer that illustrated the interface between these two domains. To continue our investigation, we now examine PA1221, a natural two-domain protein from Pseudomonas aeruginosa. We have determined the amino acid specificity of this new enzyme and used domain specific mutations to demonstrate that loading the downstream carrier domain within a single protein molecule occurs more quickly than loading of a nonfused carrier domain intermolecularly. Finally, we have determined crystal structures of both apo- and holo-PA1221 proteins, the latter using a valine-adenosine vinylsulfonamide inhibitor that traps the adenylation domain-carrier domain interaction. The protein adopts an interface similar to that seen with the prior adenylation domain-carrier protein construct. A comparison of these structures with previous structures of multidomain NRPSs suggests that a large conformational change within the NRPS adenylation domains guides the carrier domain into the active site for thioester formation.