Structural analysis of prolyl oligopeptidases using molecular docking and dynamics: insights into conformational changes and ligand binding

Prolyl oligopeptidase (POP) is considered as an important pharmaceutical target for the treatment of numerous diseases. Despite enormous studies on various aspects of POPs structure and function still some of the questions are intriguing like conformational dynamics of the protein and interplay between ligand entry/egress. Here, we have used molecular modeling and docking based approaches to unravel questions like differences in ligand binding affinities in three POP species (porcine, human and A. thaliana). Despite high sequence and structural similarity, they possess different affinities for the ligands. Interestingly, human POP was found to be more specific, selective and incapable of binding to a few planar ligands which showed extrapolation of porcine POP in human context is more complicated. Possible routes for substrate entry and product egress were also investigated by detailed analyses of molecular dynamics (MD) simulations for the three proteins. Trajectory analysis of bound and unbound forms of three species showed differences in conformational dynamics, especially variations in β-propeller pore size, which was found to be hidden by five lysine residues present on blades one and seven. During simulation, β-propeller pore size was increased by ∼2 Å in porcine ligand-bound form which might act as a passage for smaller product movement as free energy barrier was reduced, while there were no significant changes in human and A. thaliana POPs. We also suggest that these differences in pore size could lead to fundamental differences in mode of product egress among three species. This analysis also showed some functionally important residues which can be used further for in vitro mutagenesis and inhibitor design. This study can help us in better understanding of the etiology of POPs in several neurodegenerative diseases.     

Evolutionary relationships of the prolyl oligopeptidase family enzymes.

The prolyl oligopeptidase (POP) family of serine proteases includes prolyl oligopeptidase, dipeptidyl peptidase IV, acylaminoacyl peptidase and oligopeptidase B. The enzymes of this family specifically hydrolyze oligopeptides with less than 30 amino acids. Many of the POP family enzymes have evoked pharmaceutical interest as they have roles in the regulation of peptide hormones and are involved in a variety of diseases such as dementia, trypanosomiasis and type 2 diabetes. In this study we have clarified the evolutionary relationships of these four POP family enzymes and analyzed POP sequences from different sources. The phylogenetic trees indicate that the four enzymes were present in the last common ancestor of all life forms and that the beta-propeller domain has been part of the family for billions of years. There are striking differences in the mutation rates between the enzymes and POP was found to be the most conserved enzyme of this family. However, the localization of this enzyme has changed throughout evolution, as three archaeal POPs seem to be membrane bound and one third of the bacterial as well as two eukaryotic POPs were found to be secreted out of the cell. There are also considerable distinctions between the mutation rates of the different substrate binding subsites of POP. This information may help in the development of species-specific POP inhibitors.     

Distribution,classification, domain architectures and evolution of prolyl oligopeptidases in prokaryotic lineages

Background

Prolyl oligopeptidases (POPs) are proteolytic enzymes, widely distributed in all the kingdoms of life. Bacterial POPs are pharmaceutically important enzymes, yet their functional and evolutionary details are not fully explored. Therefore, current analysis is aimed at understanding the distribution, domain architecture, probable biological functions and gene family expansion of POPs in bacterial and archaeal lineages.

Results

Exhaustive sequence analysis of 1,202 bacterial and 91 archaeal genomes revealed ~3,000 POP homologs, with only 638 annotated POPs. We observed wide distribution of POPs in all the analysed bacterial lineages. Phylogenetic analysis and co-clustering of POPs of different phyla suggested their common functions in all the prokaryotic species. Further, on the basis of unique sequence motifs we could classify bacterial POPs into eight subtypes. Analysis of coexisting domains in POPs highlighted their involvement in protein-protein interactions and cellular signaling. We proposed significant extension of this gene family by characterizing 39 new POPs and 158 new α/β hydrolase members.

Conclusions

Our study reflects diversity and functional importance of POPs in bacterial species. Many genomes with multiple POPs were identified with high sequence variations and different cellular localizations. Such anomalous distribution of POP genes in different bacterial genomes shows differential expansion of POP gene family primarily by multiple horizontal gene transfer events.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-985) contains supplementary material, which is available to authorized users.     

E. coli 5'-nucleotidase undergoes a hinge-bending domain rotation resembling a ball-and-socket motion

Structures of nine independent conformers of E. coli 5'-nucleotidase (5'-NT) have been analyzed using four different crystal forms. These data show that the two-domain protein undergoes an unusual 96 degrees hinge-bending domain rotation. Structures of the open and closed forms with substrates and inhibitors reveal that the substrate moves by approximately 25 A with the large domain rotation into the catalytic site. The domain motions derived from a comparison of the nine conformations agree well with motions obtained from a normal mode analysis in that all independent domain rotations are around axes that are roughly located in the plane which includes the domain centers and the hinge. Two residues, Lys355 and Gly356, form the core of the hinge region and undergo a large change of the main-chain torsion angles. The hinge-bending movement observed for 5'-nucleotidase differs markedly from a classical hinge-bending closure motion which involves an opening of the substrate or ligand-binding cleft between two domains. In contrast, the movement observed in 5'-nucleotidase resembles that of a ball-and-socket joint. The smaller C-terminal domain rotates approximately around its center such that the residues at the domain interface move in a sliding motion along the interface. Few direct interdomain contacts and a layer of water molecules between the two domains facilitate the sliding motion.     

The loops facing the active site of prolyl oligopeptidase are crucial components in substrate gating and specificity

Prolyl oligopeptidase (POP) has emerged as a drug target for neurological diseases. A flexible loop structure comprising loop A (res. 189–209) and loop B (res. 577–608) at the domain interface is implicated in substrate entry to the active site. Here we determined kinetic and structural properties of POP with mutations in loop A, loop B, and in two additional flexible loops (the catalytic His loop, propeller Asp/Glu loop). POP lacking loop A proved to be an inefficient enzyme, as did POP with a mutation in loop B (T590C). Both variants displayed an altered substrate preference profile, with reduced ligand binding capacity. Conversely, the T202C mutation increased the flexibility of loop A, enhancing the catalytic efficiency beyond that of the native enzyme. The T590C mutation in loop B increased the preference for shorter peptides, indicating a role in substrate gating. Loop A and the His loop are disordered in the H680A mutant crystal structure, as seen in previous bacterial POP structures, implying coordinated structural dynamics of these loops. Unlike native POP, variants with a malfunctioning loop A were not inhibited by a 17-mer peptide that may bind non-productively to an exosite involving loop A. Biophysical studies suggest a predominantly closed resting state for POP with higher flexibility at the physiological temperature. The flexible loop A, loop B and His loop system at the active site is the main regulator of substrate gating and specificity and represents a new inhibitor target.     

Crystal structure of hen apo-ovotransferrin. Both lobes adopt an open conformation upon loss of iron.

The three-dimensional crystal structure of hen apo-ovotransferrin has been solved by molecular replacement and refined by simulated annealing and restrained least squares to a 3.0-A resolution. The final model, which comprises 5312 protein atoms (residues 1 to 686) and 28 carbohydrate atoms (from two monosaccharides attached to Asn(473)), gives an R-factor of 0.231 for the 11,989 observed reflections between 20.0- and 3.0-A resolution. In the structure, both empty iron binding clefts are in the open conformation, lending weight to the theory that Fe(3+) binding or release in transferrin proceeds via a mechanism that involves domain opening and closure. Upon opening, the domains rotate essentially as rigid bodies. The two domains of the N-lobe rotate away from one another by 53 degrees, whereas the C-lobe domains rotate away each another by 35 degrees. These rotations take place about an axis that passes through the two beta-strands, linking the domains. The domains of each lobe make different contacts with one another in the open and closed forms. These contacts form two interdomain interfaces on either side of the rotation axis, and domain opening or closing produces a see-saw motion between these two alternative close-packed interfaces. The interdomain disulfide bridge (Cys(478)-Cys(671)), found only in the C-lobe, may restrict domain opening but does not completely prevent it.     

Domain motion and interdomain hot spots in a multidomain enzyme

The aim of this article is to analyze conformational changes by comparing 10 different structures of Pseudomonas aeruginosa phosphomannomutase/phosphoglucomutase (PMM/PGM), a four‐domain enzyme in which both substrate binding and catalysis require substantial movement of the C‐terminal domain. We focus on changes in interdomain and active site crevices using a method called computational solvent mapping rather than superimposing the structures. The method places molecular probes (i.e., small organic molecules containing various functional groups) around the protein to find hot spots. One of the most important hot spots is in the active site, consistent with the ability of the enzyme to bind both glucose and mannose phosphosugar substrates. The protein has eight additional hot spots at domain‐domain interfaces and hinge regions. The locations and nature of six of these hot spots vary between the open, half‐open, and closed conformers of the enzyme, in good agreement with the ligand‐induced conformational changes. In the closed structures the number of probe clusters at the hinge region significantly depends on the position of the phosphorylated oxygen in the substrate (e.g., glucose 1‐phosphate versus glucose 6‐phosphate), but the protein remains almost unchanged in terms of the overall RMSD, indicating that computational solvent mapping is a more sensitive approach to detect changes in binding sites and interdomain crevices. Focusing on multidomain proteins we show that the subresolution conformational differences revealed by the mapping are in fact significant, and present a general statistical method of analysis to determine the significance of rigid body domain movements in X‐ray structures.     

Interdomain mobility in di-ubiquitin revealed by NMR

Domain orientation and dynamics can play an essential role in the function of multidomain proteins. Lys48-linked polyubiquitin chains, the principal signal for proteasomal protein degradation, adopt a closed conformation at physiological conditions, in which the functionally important residues Leu8, Ile44, and Val70 are sequestered at the interdomain interface. This interface must open in order for these groups to become available for interactions with various chain-recognition factors. Knowledge of the mechanism of domain motion leading to the opening of the interdomain interface in polyubiqutin is, therefore, essential for the understanding of the processes controlling molecular recognition events in polyubiquitin signaling. Here we use NMR to characterize the interdomain dynamics that open the interface in a di-ubiquitin chain. This process occurs via domain reorientations on a 10-ns time scale and with the amplitudes that are sufficient for making functionally important hydrophobic residues in polyubiquitin available for direct interactions with various ubiquitin-binding factors. The analysis revealed the structures of the interconverting conformational states of di-ubiquitin and the rates and amplitudes of this process at near-physiological and acidic pH. The proposed mechanism of domain reorientation is quite general and could serve as a paradigm of interdomain mobility in other multidomain systems.     

Kinetic and mechanistic studies of prolyl oligopeptidase from the hyperthermophile Pyrococcus furiosus

Prolyl oligopeptidase (POP) is widely distributed in mammals, where it is implicated in neuropeptide processing. It is also present in some bacteria and archaea. Because POP is found in mesophilic and hyperthermophilic organisms, and is distributed among all three phylogenetic domains, studies of its function and structure could lead to new insights about the evolution of enzyme mechanisms and thermostability. Kinetic studies were conducted on the POP of the hyperthermophilic archaeon Pyrococcus furiosus (Pfu) 85 degrees C in both H(2)O and D(2)O. Pfu POP displayed many similarities to mammalian POPs, however the solvent isotope effect (k(0)/k(1)) was 2.2 at both high and low pH, indicating that general base/acid catalysis is the rate-limiting step. The pH-rate profiles indicated a three-deprotonation process with pK(a) values of 4.3, 7.2, and 9.1. The temperature dependence of these values revealed a heat of ionization of 4.7 kJ/mol for pK(es1) and 22 kJ/mol for pK(es2), suggesting the catalytic involvement of a carboxyl group and an imidazole group, respectively. Temperature dependence of the catalytic rate was assessed at pH 6.0 and 7.6. Entropy values of -119 and -143 Jmol(-1)K(-1) were calculated at the respective pH values, with a corresponding difference in enthalpy of 8.5 kJ/mol. These values suggest that two or three hydrogen bonds are broken during the transition state of the acidic enzyme form, whereas only one or two are broken during the transition state of the basic enzyme form. A model has been constructed for Pfu POP based on the crystal structure of porcine POP and the sequence alignment. The similarities demonstrated for POPs from these two organisms reflect the most highly conserved characteristics of this class of serine protease, whereas the differences between these enzymes highlights the large evolutionary distance between them. Such fundamental information is crucial to our understanding of the function of proteins at high temperature.     

Distinguishing unfolding and functional conformational transitions of calmodulin using ultraviolet resonance Raman spectroscopy

Calmodulin (CaM) is a ubiquitous moderator protein for calcium signaling in all eukaryotic cells. This small calcium‐binding protein exhibits a broad range of structural transitions, including domain opening and folding–unfolding, that allow it to recognize a wide variety of binding partners in vivo. While the static structures of CaM associated with its various binding activities are fairly well‐known, it has been challenging to examine the dynamics of transition between these structures in real‐time, due to a lack of suitable spectroscopic probes of CaM structure. In this article, we examine the potential of ultraviolet resonance Raman (UVRR) spectroscopy for clarifying the nature of structural transitions in CaM. We find that the UVRR spectral change (with 229 nm excitation) due to thermal unfolding of CaM is qualitatively different from that associated with opening of the C‐terminal domain in response to Ca²⁺ binding. This spectral difference is entirely due to differences in tertiary contacts at the interdomain tyrosine residue Tyr138, toward which other spectroscopic methods are not sensitive. We conclude that UVRR is ideally suited to identifying the different types of structural transitions in CaM and other proteins with conformation‐sensitive tyrosine residues, opening a path to time‐resolved studies of CaM dynamics using Raman spectroscopy.     

Interdomain interactions in hinge-bending transitions.

BACKGROUND: The mechanisms that allow or constrain protein movement have not been understood. Here we study interdomain interactions in proteins to investigate hinge-bending motions. RESULTS: We find a limited number of salt bridges and hydrogen bonds at the interdomain interface, in both the "closed" and the "open" conformations. Consistently, analysis of 222 salt bridges in an independently selected database indicates that most salt bridges form within rather than between independently folding hydrophobic units. Calculations show that these interdomain salt bridges either destabilize or only marginally stabilize the closed conformation in most proteins. In contrast, the nonpolar buried surface area between the moving parts can be extensive in the closed conformations. However, when the nonpolar buried surface area is large, we find that at the interdomain interface in the open conformation it may be as large or larger than in the closed conformation. Hence, the energetic penalty of opening the closed conformation is overcome. Consistently, a large nonpolar surface area buried in the closed interdomain interface accompanies limited opening of the domains, yielding a larger interface. CONCLUSIONS: Short-range electrostatic interactions are largely absent between moving domains. Interdomain nonpolar buried surface area may be large in the closed conformation, but it is largely offset by the area buried in the open conformation. In such cases the opening of the domains appears to be relatively small. This may allow prediction of the extent of domain opening. Such predictions may have implications for the shape and size of the binding pockets in drug/protein design.     

A possible role of prolyl oligopeptidase during Linum usitatissimum (flax) seed development

Involvement of prolyl oligopeptidases (POPs) in the control of several mammalian peptide hormone signalling pathways has been studied extensively in recent years. POPs are ubiquitous enzymes, but little attention has been paid to understanding their function in plants. Using a cDNA-AFLP approach, two flax (Linum usitatissimum) POP ESTs were identified as being specifically expressed in the early stages of flax seed development. This specific expression was confirmed using real time RT-PCR and in situ hybridisation approaches. Seed expression of Arabidopsis POP genes was measured and showed no specificity. Comparison between results obtained with flax and Arabidopsis is discussed in order to address a hypothetic function for POPs during seed formation. These results provide the first insights into POP gene expression and hypothetical function in plants.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Mechanisms of intramolecular communication in a hyperthermophilic acylaminoacyl peptidase: a molecular dynamics investigation

Protein dynamics and the underlying networks of intramolecular interactions and communicating residues within the three-dimensional (3D) structure are known to influence protein function and stability, as well as to modulate conformational changes and allostery. Acylaminoacyl peptidase (AAP) subfamily of enzymes belongs to a unique class of serine proteases, the prolyl oligopeptidase (POP) family, which has not been thoroughly investigated yet. POPs have a characteristic multidomain three-dimensional architecture with the active site at the interface of the C-terminal catalytic domain and a β-propeller domain, whose N-terminal region acts as a bridge to the hydrolase domain. In the present contribution, protein dynamics signatures of a hyperthermophilic acylaminoacyl peptidase (AAP) of the prolyl oligopeptidase (POP) family, as well as of a deletion variant and alanine mutants (I12A, V13A, V16A, L19A, I20A) are reported. In particular, we aimed at identifying crucial residues for long range communications to the catalytic site or promoting the conformational changes to switch from closed to open ApAAP conformations. Our investigation shows that the N-terminal α1-helix mediates structural intramolecular communication to the catalytic site, concurring to the maintenance of a proper functional architecture of the catalytic triad. Main determinants of the effects induced by α1-helix are a subset of hydrophobic residues (V16, L19 and I20). Moreover, a subset of residues characterized by relevant interaction networks or coupled motions have been identified, which are likely to modulate the conformational properties at the interdomain interface.     

Unique GMP-binding site in Mycobacterium tuberculosis guanosine monophosphate kinase

Bacterial nucleoside monophosphate (NMP) kinases, which convert NMPs to nucleoside diphosphates (NDP), are investigated as potential antibacterial targets against pathogenic bacteria. Herein, we report the biochemical and structural characterization of GMP kinase from Mycobacterium tuberculosis (GMPKMt). GMPKMt is a monomer with an unusual specificity for ATP as a phosphate donor, a lower catalytic efficiency compared with eukaryotic GMPKs, and it carries two redox-sensitive cysteines in the central CORE domain. These properties were analyzed in the light of the high-resolution crystal structures of unbound, GMP-bound, and GDP-bound GMPKMt. The latter structure was obtained in both an oxidized form, in which the cysteines form a disulfide bridge, and a reduced form which is expected to correspond to the physiological enzyme. GMPKMt has a modular domain structure as most NMP kinases. However, it departs from eukaryotic GMPKs by the unusual conformation of its CORE domain, and by its partially open LID and GMP-binding domains which are the same in the apo-, GMP-bound, and GDP-bound forms. GMPKMt also features a unique GMP binding site which is less close-packed than that of mammalian GMPKs, and in which the replacement of a critical tyrosine by a serine removes a catalytic interaction. In contrast, the specificity of GMPKMt for ATP may be a general feature of GMPKs because of an invariant structural motif that recognizes the adenine base. Altogether, differences in domain dynamics and GMP binding between GMPKMt and mammalian GMPKs should reveal clues for the design of GMPKMt-specific inhibitors.     

Crystal structure of novel NADP-dependent 3-hydroxyisobutyrate dehydrogenase from Thermus thermophilus HB8

3-Hydroxyisobutyrate, a central metabolite in the valine catabolic pathway, is reversibly oxidized to methylmalonate semialdehyde by a specific dehydrogenase belonging to the 3-hydroxyacid dehydrogenase family. To gain insight into the function of this enzyme at the atomic level, we have determined the first crystal structures of the 3-hydroxyisobutyrate dehydrogenase from Thermus thermophilus HB8: holo enzyme and sulfate ion complex. The crystal structures reveal a unique tetrameric oligomerization and a bound cofactor NADP+. This bacterial enzyme may adopt a novel cofactor-dependence on NADP, whereas NAD is preferred in eukaryotic enzymes. The protomer folds into two distinct domains with open/closed interdomain conformations. The cofactor NADP+ with syn nicotinamide and the sulfate ion are bound to distinct sites located at the interdomain cleft of the protomer through an induced-fit domain closure upon cofactor binding. From the structural comparison with the crystal structure of 6-phosphogluconate dehydrogenase, another member of the 3-hydroxyacid dehydrogenase family, it is suggested that the observed sulfate ion and the substrate 3-hydroxyisobutyrate share the same binding pocket. The observed oligomeric state might be important for the catalytic function through forming the active site involving two adjacent subunits, which seems to be conserved in the 3-hydroxyacid dehydrogenases. A kinetic study confirms that this enzyme has strict substrate specificity for 3-hydroxyisobutyrate and serine, but it cannot distinguish the chirality of the substrates. Lys165 is likely the catalytic residue of the enzyme.     

Plant odour plumes as mediators of plant–insect interactions

Insect olfactory orientation along odour plumes has been studied intensively with respect to pheromonal communication, whereas little knowledge is available on how plant odour plumes (POPs) affect olfactory searching by an insect for its host plants. The primary objective of this review is to examine the role of POPs in the attraction of insects. First, we consider parameters of an odour source and the environment which determine the size, shape and structure of an odour plume, and we apply that knowledge to POPs. Second, we compare characteristics of insect pheromonal plumes and POPs. We propose a ‘POP concept’ for the olfactory orientation of insects to plants. We suggest that: (i) an insect recognises a POP by means of plant volatile components that are encountered in concentrations higher than a threshold detection limit and that occur in a qualitative and quantitative blend indicating a resource; (ii) perception of the fine structure of a POP enables an insect to distinguish a POP from an unspecific odorous background and other interfering plumes; and (iii) an insect can follow several POPs to their sources, and may leave the track of one POP and switch to another one if this conveys a signal with higher reliability or indicates a more suitable resource. The POP concept proposed here may be a useful tool for research in olfactory‐mediated plant–insect interactions.     

Interdomain salt‐bridges in the Ebola virus protein VP40 and their role in domain association and plasma membrane localization

The Ebola virus protein VP40 is a transformer protein that possesses an extraordinary ability to accomplish multiple functions by transforming into various oligomeric conformations. The disengagement of the C‐terminal domain (CTD) from the N‐terminal domain (NTD) is a crucial step in the conformational transformations of VP40 from the dimeric form to the hexameric form or octameric ring structure. Here, we use various molecular dynamics (MD) simulations to investigate the dynamics of the VP40 protein and the roles of interdomain interactions that are important for the domain–domain association and dissociation, and report on experimental results of the behavior of mutant variants of VP40. The MD studies find that various salt‐bridge interactions modulate the VP40 domain dynamics by providing conformational specificity through interdomain interactions. The MD simulations reveal a novel salt‐bridge between D45‐K326 when the CTD participates in a latch‐like interaction with the NTD. The D45‐K326 salt‐bridge interaction is proposed to help domain–domain association, whereas the E76‐K291 interaction is important for stabilizing the closed‐form structure. The effects of the removal of important VP40 salt‐bridges on plasma membrane (PM) localization, VP40 oligomerization, and virus like particle (VLP) budding assays were investigated experimentally by live cell imaging using an EGFP‐tagged VP40 system. It is found that the mutations K291E and D45K show enhanced PM localization but D45K significantly reduced VLP formation.     

Methionine residues lining the substrate pathway in prolyl oligopeptidase from Pleurotus eryngii play an important role in substrate recognition

Family S9 prolyl oligopeptidases (POPs) are of interest as pharmacological targets. We recently found that an S9 POP from Pleurotus eryngii showed altered substrate specificity following H₂O₂ treatment. Oxidation of Met203 on the non-catalytic β-propeller domain resulted in decreased activity toward non-aromatic aminoacyl-para-nitroanilides (pNAs) while maintaining its activity toward aromatic aminoacyl-pNAs. Given that the other Met residues should also be oxidized by H₂O₂ treatment, we constructed mutants in which all the Met residues were substituted with other amino acids. Analysis of the mutants showed that Met570 in the catalytic domain is another potent residue for the altered substrate specificity following oxidation. Met203 and Met570 lie on the surfaces of two different domains and form part of a funnel from the surface to the active center. Our findings indicate that the funnel forms the substrate pathway and plays a role in substrate recognition.     

Segment swapping aided the evolution of enzyme function: The case of uroporphyrinogen III synthase

In an earlier study, we showed that two‐domain segment‐swapped proteins can evolve by domain swapping and fusion, resulting in a protein with two linkers connecting its domains. We proposed that a potential evolutionary advantage of this topology may be the restriction of interdomain motions, which may facilitate domain closure by a hinge‐like movement, crucial for the function of many enzymes. Here, we test this hypothesis computationally on uroporphyrinogen III synthase, a two‐domain segment‐swapped enzyme essential in porphyrin metabolism. To compare the interdomain flexibility between the wild‐type, segment‐swapped enzyme (having two interdomain linkers) and circular permutants of the same enzyme having only one interdomain linker, we performed geometric and molecular dynamics simulations for these species in their ligand‐free and ligand‐bound forms. We find that in the ligand‐free form, interdomain motions in the wild‐type enzyme are significantly more restricted than they would be with only one interdomain linker, while the flexibility difference is negligible in the ligand‐bound form. We also estimated the entropy costs of ligand binding associated with the interdomain motions, and find that the change in domain connectivity due to segment swapping results in a reduction of this entropy cost, corresponding to ～20% of the total ligand binding free energy. In addition, the restriction of interdomain motions may also help the functional domain‐closure motion required for catalysis. This suggests that the evolution of the segment‐swapped topology facilitated the evolution of enzyme function for this protein by influencing its dynamic properties. Proteins 2016; 85:46–53. © 2016 Wiley Periodicals, Inc.