Evolutionary conservation of protein dynamics: insights from all-atom molecular dynamics simulations of 'peptidase' domain of Spt16.

Protein function is encoded in its sequence, manifested in its three-dimensional structure, and facilitated by its dynamics. Studies have suggested that protein structures with higher sequence similarity could have more similar patterns of dynamics. However, such studies of protein dynamics within and across protein families typically rely on coarse-grained models, or approximate metrics like crystallographic B-factors. This study uses µs scale molecular dynamics (MD) simulations to explore the conservation of dynamics among homologs of ∼50 kDa N-terminal module of Spt16 (Spt16N). Spt16N from Saccharomyces cerevisiae (Sc-Spt16N) and three of its homologs with 30-40% sequence identities were available in the PDB. To make our data-set more comprehensive, the crystal structure of an additional homolog (62% sequence identity with Sc-Spt16N) was solved at 1.7 Å resolution. Cumulative MD simulations of 6 µs were carried out on these Spt16N structures and on two additional protein structures with varying degrees of similarity to it. The simulations revealed that correlation in patterns of backbone fluctuations vary linearly with sequence identity. This trend could not be inferred using crystallographic B-factors. Further, normal mode analysis suggested a similar pattern of inter-domain (inter-lobe) motions not only among Spt16N homologs, but also in the M24 peptidase structure. On the other hand, MD simulation results highlighted conserved motions that were found unique for Spt16N protein, this along with electrostatics trends shed light on functional aspects of Spt16N.Communicated by Ramaswamy H. Sarma.

Machine learning classifiers aid virtual screening for efficient design of mini-protein therapeutics.

De novo design of mini-proteins (4-12 kDa) has recently been shown to produce new candidates for protein therapeutics. They are temperature stable molecules that bind to the drug target with high affinity for inhibiting its interactions. The development of mini-protein binders requires laboratory screening of tens of thousands of molecules for effective target binding. In this study we trained machine learning classifiers which can distinguish, with 90% accuracy and 80% precision, mini-protein binders from non-binding molecules designed for a particular target; this significantly reduces the number of mini protein candidates for experimental screening. Further, on the basis of our results we propose a multi-stage protocol where a small dataset (few hundred experimentally verified target-specific mini-proteins) can be used to train classifiers for improving the efficiency of mini-protein design for any specific target.

Crystal structure of XCC3289 from Xanthomonas campestris: homology with the N-terminal substrate-binding domain of Lon peptidase.

LonA peptidase is a major component of the protein quality-control mechanism in both prokaryotes and the organelles of eukaryotes. Proteins homologous to the N-terminal domain of LonA peptidase, but lacking its other domains, are conserved in several phyla of prokaryotes, including the Xanthomonadales order. However, the function of these homologous proteins (LonNTD-like proteins) is not known. Here, the crystal structure of the LonNTD-like protein from Xanthomonas campestris (XCC3289; UniProt Q8P5P7) is reported at 2.8 Šresolution. The structure was solved by molecular replacement and contains one polypeptide in the asymmetric unit. The structure was refined to an Rfree of 29%. The structure of XCC3289 consists of two domains joined by a long loop. The N-terminal domain (residues 1-112) consists of an α-helix surrounded by ß-sheets, whereas the C-terminal domain (residues 123-193) is an α-helical bundle. The fold and spatial orientation of the two domains closely resembles those of the N-terminal domains of the LonA peptidases from Escherichia coli and Mycobacterium avium. The structure is also similar to that of cereblon, a substrate-recognizing component of the E3 ubiquitin ligase complex. The N-terminal domains of both LonA and cereblon are known to be involved in specific protein-protein interactions. This structural analysis suggests that XCC3289 and other LonNTD-like proteins might also be capable of such protein-protein interactions.

Structural basis for the unusual substrate specificity of unique two-domain M1 metallopeptidase.

M1 metallopeptidases regulate many important biological processes such as angiogenesis, tumour growth, hormone regulation, and immune cell development. Knowledge of substrate specificity mechanism in this family is valuable. An M1 peptidase from Deinococcus radiodurans (M1dr) with preference for bulky hydrophobic residues at N-terminus of peptide substrates was recently reported. In contrast to Escherichia coli aminopeptidase N, a previously characterized M1 peptidase, M1dr exhibits reduced activity towards peptides with N-terminal Arg or Ala residue. In order to illuminate structural basis of substrate specificity, we report several crystal structures of M1dr with different amino acids bound to the active site. Structural analysis indicated that the enzyme makes subtle adjustments to multiple residues leading to significant volume change of the active site cavity to accommodate residues of varying sizes (Leu to Trp). This study further reveals that the low preference for Arg at N-terminus of peptide substrate arises from a non-productive conformation in which many of the Arg molecules bind where they block the proton donor essential for the peptidase reaction. Hence, this study illuminates the substrate-binding mechanism and also reveals the structural basis for the substrate specificity of M1dr enzyme.

Phylogenetic spread of sequence data affects fitness of consensus enzymes: Insights from triosephosphate isomerase.

The concept of consensus in multiple sequence alignments (MSAs) has been used to design and engineer proteins previously with some success. However, consensus design implicitly assumes that all amino acid positions function independently, whereas in reality, the amino acids in a protein interact with each other and work cooperatively to produce the optimum structure required for its function. Correlation analysis is a tool that can capture the effect of such interactions. In a previously published study, we made consensus variants of the triosephosphate isomerase (TIM) protein using MSAs that included sequences form both prokaryotic and eukaryotic organisms. These variants were not completely native-like and were also surprisingly different from each other in terms of oligomeric state, structural dynamics, and activity. Extensive correlation analysis of the TIM database has revealed some clues about factors leading to the unusual behavior of the previously constructed consensus proteins. Among other things, we have found that the more ill-behaved consensus mutant had more broken correlations than the better-behaved consensus variant. Moreover, we report three correlation and phylogeny-based consensus variants of TIM. These variants were more native-like than the previous consensus mutants and considerably more stable than a wild-type TIM from a mesophilic organism. This study highlights the importance of choosing the appropriate diversity of MSA for consensus analysis and provides information that can be used to engineer stable enzymes.

Two-domain aminopeptidase of M1 family: Structural features for substrate binding and gating in absence of C-terminal domain.

Zinc metallopeptidases of the M1 family (M1 peptidases) with unique metal binding motif HEXXH(X)18E regulate many important biological processes such as tumor growth, angiogenesis, hormone regulation, and immune cell development. Typically, these enzymes exist in three-domain [N-terminal domain (N-domain), catalytic domain, and C-terminal domain (C-domain)] or four-domain (N-domain, catalytic domain, middle domain, and C-domain) format in which N-domain and catalytic domain are more conserved. The C-domain plays important roles in substrate binding and gating. In this study we report the first structure of a two-domain (N-domain and catalytic domain) M1 peptidase at 2.05 Šresolution. Despite the lack of C-domain, the enzyme is active and prefers peptide substrates with large hydrophobic N-terminal residues. Its substrate-bound structure was determined at 1.9 Šresolution. Structural analyses supported by site directed mutagenesis and molecular dynamics simulations reveal structural features that could compensate for the lack of C-domain. A unique loop insertion (loop A) in the N-domain has important roles in gating and desolvation of active site. Three Arg residues of the catalytic domain are involved in substrate-binding roles typically played by positively charged residues of C-domain in other M1 peptidases. Further, its unique exopeptidase sequence motif, LALET, creates a more hydrophobic environment at the S1 subsite (which binds N-terminal residue of the substrate in aminopeptidases) than the more common GXMEN motif in the family. This leads to high affinity for large hydrophobic residues in the S1 subsite, which contributes towards efficient substrate binding in absence of C-domain.

Catalytic triad heterogeneity in S51 peptidase family: Structural basis for functional variability.

Peptidase E (PepE) is a nonclassical serine peptidase with a Ser-His-Glu catalytic triad. It is specific for dipeptides with an N-terminal aspartate residue (Asp-X dipeptidase activity). Its homolog from Listeria monocytogenes (PepElm) has a Ser-His-Asn "catalytic triad." Based on sequence alignment we predicted that the PepE homolog from Deinococcus radiodurans (PepEdr) would have a Ser-His-Asp "catalytic triad." We confirmed this by solving the crystal structure of PepEdr to 2.7 Å resolution. We show that PepElm and PepEdr lack the Asp-X dipeptidase activity. Our analyses suggest that absence of P1 pocket in the active site could be the main reason for this lack of typical activity. Sequence and structural data reveal that the PepE homologs can be divided into long and short PepEs based on presence or absence of a C-terminal tail which adopts a ß-hairpin conformation in the canonical PepE from Salmonella enterica. A long PepE from Bacillus subtilis with Ser-His-Asp catalytic triad exhibits Asp-X dipeptidase activity. Whereas the three long PepEs enzymatically characterized till date have been found to possess the Asp-X dipeptidase activity, the three enzymatically characterized short PepEs lack this activity irrespective of the nature of their catalytic triads. This study illuminates the structural and functional heterogeneity in the S51 family and also provides structural basis for the functional variability among PepE homologs.

Carboxypeptidase in prolyl oligopeptidase family: Unique enzyme activation and substrate-screening mechanisms.

Serine peptidases of the prolyl oligopeptidase (POP) family are of substantial therapeutic importance because of their involvement in diseases such as diabetes, cancer, neurological diseases, and autoimmune disorders. Proper annotation and knowledge of substrate specificity mechanisms in this family are highly valuable. Although endopeptidase, dipeptidyl peptidase, tripeptidyl peptidase, and acylaminoacyl peptidase activities have been reported previously, here we report the first instance of carboxypeptidase activity in a POP family member. We determined the crystal structures of this carboxypeptidase, an S9C subfamily member from Deinococcus radiodurans, in its active and inactive states at 2.3-Å resolution, providing an unprecedented view of assembly and disassembly of the active site mediated by an arginine residue. We observed that this residue is poised to bind substrate in the active structure and disrupts the catalytic triad in the inactive structure. The assembly of the active site is accompanied by the ordering of gating loops, which reduces the effective size of the oligomeric pore. This prevents the entry of larger peptides and constitutes a novel mechanism for substrate screening. Furthermore, we observed structural adaptations that enable its carboxypeptidase activity, with a unique loop and two arginine residues in the active site cavity orienting the peptide substrate for catalysis. Using these structural features, we identified homologs of this enzyme in the POP family and confirmed the presence of carboxypeptidase activity in one of them. In conclusion, we have identified a new type within POP enzymes that exhibits not only unique activity but also a novel substrate-screening mechanism.

Crystal structures and biochemical analyses of intermediate cleavage peptidase: role of dynamics in enzymatic function.

Intermediate cleavage peptidase (Icp55) processes a subset of mitochondrial matrix proteins by removing a bulky residue at their N termini, leaving behind smaller N-terminal residues (icp activity). This contributes towards the stability of the mitochondrial proteome. We report crystal structures of yeast Icp55 including one bound to the apstatin inhibitor. Apart from icp activity, the enzyme was found to exhibit Xaa-Pro aminopeptidase activity in vitro. Structural and biochemical data suggest that the enzyme exists in a rapid equilibrium between monomer and dimer. Furthermore, the dimer, and not the monomer, was found to be the active species with loop dynamics at the dimer interface playing an important role in activity. Based on the new evidence, we propose a model for binding and processing of cellular targets by Icp55. DATABASE: The atomic coordinates and structure factors for the structures of Icp55 (code 6A9T, 6A9U, 6A9V) have been deposited in the Protein Data Bank (PDB) (http://www.pdb.org/).

Structures and activities of widely conserved small prokaryotic aminopeptidases-P clarify classification of M24B peptidases.

M24B peptidases cleaving Xaa-Pro bond in dipeptides are prolidases whereas those cleaving this bond in longer peptides are aminopeptidases-P. Bacteria have small aminopeptidases-P (36-39 kDa), which are diverged from canonical aminopeptidase-P of Escherichia coli (50 kDa). Structure-function studies of small aminopeptidases-P are lacking. We report crystal structures of small aminopeptidases-P from E. coli and Deinococcus radiodurans, and report substrate-specificities of these proteins and their ortholog from Mycobacterium tuberculosis. These are aminopeptidases-P, structurally close to small prolidases except for absence of dipeptide-selectivity loop. We noticed absence of this loop and conserved arginine in canonical archaeal prolidase (Maher et al., Biochemistry. 43, 2004, 2771-2783) and questioned its classification. Our enzymatic assays show that this enzyme is an aminopeptidase-P. Further, our mutagenesis studies illuminate importance of DXRY sequence motif in bacterial small aminopeptidases-P and suggest common evolutionary origin with human XPNPEP1/XPNPEP2. Our analyses reveal sequence/structural features distinguishing small aminopeptidases-P from other M24B peptidases.

Structure of Asp-bound peptidase E from Salmonella enterica: Active site at dimer interface illuminates Asp recognition.

Peptidase-E, a nonclassical serine peptidase, is specific for dipeptides with an N-terminal aspartate. This stringent substrate specificity remains largely unexplained. We report an aspartate-bound structure of peptidase-E at 1.83 Å resolution. In contrast to previous reports, the enzyme forms a dimer, and the active site is located at the dimer interface, well shielded from the solvent. Our findings further suggest that the stringent aspartate specificity of the enzyme is due to electrostatics and molecular complementarity in the active site. The new structural information presented herein may provide insights into the role of functionally important residues in peptidase-E.

Phylogenetic spread of sequence data affects fitness of SOD1 consensus enzymes: Insights from sequence statistics and structural analyses.

Non-natural protein sequences with native-like structures and functions can be constructed successfully using consensus design. This design strategy is relatively well understood in repeat proteins with simple binding function, however detailed studies are lacking in globular enzymes. The SOD1 family is a good model for such studies due to the availability of large amount of sequence and structure data motivated by involvement of human SOD1 in the fatal motor neuron disease amyotrophic lateral sclerosis (ALS). We constructed two consensus SOD1 enzymes from multiple sequence alignments from all organisms and eukaryotic organisms. A significant difference in their catalytic activities shows that the phylogenetic spread of the sequences used affects the fitness of the construct obtained. A mutation in an electrostatic loop and overall design incompatibilities between bacterial and eukaryotic sequences were implicated in this disparity. Based on this analysis, a bioinformatics approach was used to classify mutations thought to cause familial ALS providing a unique high level view of the physical basis of disease-causing aggregation of human SOD1.

Crystal structure of a novel prolidase from Deinococcus radiodurans identifies new subfamily of bacterial prolidases.

Xaa-Pro peptidases (XPP) are dinuclear peptidases of MEROPS M24B family that hydrolyze Xaa-Pro iminopeptide bond with a trans-proline at the second position of the peptide substrate. XPPs specific towards dipeptides are called prolidases while those that prefer longer oligopeptides are called aminopeptidases P. Though XPPs are strictly conserved in bacterial and archaeal species, the structural and sequence features that distinguish between prolidases and aminopeptidases P are not always clear. Here, we report 1.4 Å resolution crystal structure of a novel XPP from Deinococcus radiodurans (XPPdr). XPPdr forms a novel dimeric structure via unique dimer stabilization loops of N-terminal domains such that their C-terminal domains are placed far apart from each other. This novel dimerization is also the consequence of a different orientation of N-terminal domain in XPPdr monomer than those in other known prolidases. The enzymatic assays show that it is a prolidase with broad substrate specificity. Our structural, mutational, and molecular dynamics simulation analyses show that the conserved Arg46 of N-terminal domain is important for the dipeptide selectivity. Our BLAST search found XPPdr orthologs with conserved sequence motifs which correspond to unique structural features of XPPdr, thus identify a new subfamily of bacterial prolidases.

Active site gate of M32 carboxypeptidases illuminated by crystal structure and molecular dynamics simulations.

Enzyme gates are important dynamic features that regulate function. Study of these features is critical for understanding of enzyme mechanism. In this study, the active-site gate of M32 carboxypeptidases (M32CP) is illuminated. Only a handful of members of this family have been structurally and functionally characterized and various aspects of their activity and mechanism are yet not clarified. Here, crystal structure of putative M32CP from Deinococcus radiodurans (M32dr) was solved to 2.4Å resolution. Enzymatic assays confirmed its identity as a carboxypeptidase. Open and relatively closed conformations observed in the structure provided supporting evidence for previously hypothesized hinge motion in this family of enzymes. Molecular dynamics simulations of 1.5µs displayed distinct open and closed conformations revealing amplitude of the motion to be beyond what was observed in the crystal structure. Hinge region and anchoring region of this shell-type gate were identified. A small displacement of 3Å and a helical tilt of 9° propagated by the hinge region translates into a 10Å motion at the top of the gate. The dynamics of the gate was supported by our mutagenesis experiment involving formation of disulphide bond across helices of the gate. The nearly inactive mutant enzyme showed 65-fold increase in the enzymatic activity in presence of reducing agent. Further, while a previously proposed structural basis would have led to its classification in subfamily II, experimentally observed substrate length restriction places M32dr in subfamily I of M32CPs.

Structure of the human aminopeptidase XPNPEP3 and comparison of its in vitro activity with Icp55 orthologs: Insights into diverse cellular processes.

The human aminopeptidase XPNPEP3 is associated with cystic kidney disease and TNF-TNFR2 cellular signaling. Its yeast and plant homolog Icp55 processes several imported mitochondrial matrix proteins leading to their stabilization. However, the molecular basis for the diverse roles of these enzymes in the cell is unknown. Here, we report the crystal structure of human XPNPEP3 with bound apstatin product at 1.65 Å resolution, and we compare its in vitro substrate specificity with those of fungal Icp55 enzymes. In contrast to the suggestions by earlier in vivo studies of mitochondrial processing, we found that these enzymes are genuine Xaa-Pro aminopeptidases, which hydrolyze peptides with proline at the second position (P1'). The mitochondrial processing activity involving cleavage of peptides lacking P1' proline was also detected in the purified enzymes. A wide proline pocket as well as molecular complementarity and capping at the S1 substrate site of XPNPEP3 provide the necessary structural features for processing the mitochondrial substrates. However, this activity was found to be significantly lower as compared with Xaa-Pro aminopeptidase activity. Because of similar activity profiles of Icp55 and XPNPEP3, we propose that XPNPEP3 plays the same mitochondrial role in humans as Icp55 does in yeast. Both Xaa-Pro aminopeptidase and mitochondrial processing activities of XPNPEP3 have implications toward mitochondrial fitness and cystic kidney disease. Furthermore, the presence of both these activities in Icp55 elucidates the unexplained processing of the mitochondrial cysteine desulfurase Nfs1 in yeast. The enzymatic and structural analyses reported here provide a valuable molecular framework for understanding the diverse cellular roles of XPNPEP3.

Crystal structure and biochemical investigations reveal novel mode of substrate selectivity and illuminate substrate inhibition and allostericity in a subfamily of Xaa-Pro dipeptidases.

Xaa-Pro dipeptidase (XPD) catalyzes hydrolysis of iminopeptide bond in dipeptides containing trans-proline as a second residue. XPDs are found in all living organisms and are believed to play an essential role in proline metabolism. Here, we report crystal structures and extensive enzymatic studies of XPD from Xanthomonas campestris (XPDxc), the first such comprehensive study of a bacterial XPD. We also report enzymatic activities of its ortholog from Mycobacterium tuberculosis (XPDmt). These enzymes are strictly dipeptidases with broad substrate specificities. They exhibit substrate inhibition and allostericity, as described earlier for XPD from Lactococcus lactis (XPDll). The structural, mutational and comparative data have revealed a novel mechanism of dipeptide selectivity and substrate binding in these enzymes. Moreover, we have identified conserved sequence motifs that distinguish these enzymes from other prolidases, thus defining a new subfamily. This study provides a suitable structural template for explaining unique properties of this XPDxc subfamily. In addition, we report unique structural features of XPDxc protein like an extended N-terminal tail region and absence of a conserved Tyr residue near the active site.

Gating role of His 72 in TmPurL enzyme uncovered by structural analyses and molecular dynamics simulations.

Histidine is ubiquitous in enzyme active sites but its role is often difficult to ascribe due to ambiguity of protonation state and complex electrostatic and dynamic effects involved. In this study the role of His 72 in TmPurL, a glutamine amidotransferase (GAT) enzyme, is investigated. TmPurL is a large 66kDa enzyme that works as part of an even larger (>100kDa) multi-protein complex. This enzyme complex performs an essential step in the purine biosynthesis pathway by abstracting ammonia from a glutamine molecule and channeling it 30Å away into the active site of TmPurL, incorporating it into a purine biosynthesis intermediate. It is known that His 72 is important for substrate binding and possibly acts as a general base. Comparing apo and holo structural forms of this enzyme has revealed to us a possible gating function of His 72 that could regulate ammonia entry into the active site. Bimodal distribution of the χ1 dihedral angle of this amino acid in molecular dynamics simulations of 2µs supports the hypothesis. Different protonation states of His 72 were found to be conformationally distinct, providing a possible link between catalytic and gating roles of His 72. Ammonia channeling and allostery are discussed for GATs and more specifically for PurL family.

Crystallization and preliminary X-ray crystallographic analysis of an artificial molten-globular-like triosephosphate isomerase protein of mixed phylogenetic origin.

A bioinformatics-based protein-engineering approach called consensus design led to the construction of a chimeric triosephosphate isomerase (TIM) protein called ccTIM (curated consensus TIM) which is as active as Saccharomyces cerevisiae TIM despite sharing only 58% sequence identity with it. The amino-acid sequence of this novel protein is as identical to native sequences from eukaryotes as to those from prokaryotes and shares some biophysical traits with a molten globular protein. Solving its crystal structure would help in understanding the physical implications of its bioinformatics-based sequence. In this report, the ccTIM protein was successfully crystallized using the microbatch-under-oil method and a full X-ray diffraction data set was collected to 2.2 Šresolution using a synchrotron-radiation source. The crystals belonged to space group C2221, with unit-cell parameters a=107.97, b=187.21, c=288.22 Å. Matthews coefficient calculations indicated the presence of six dimers in the asymmetric unit, with an approximate solvent content of 46.2%.

Importance of hydrophobic cavities in allosteric regulation of formylglycinamide synthetase: insight from xenon trapping and statistical coupling analysis.

Formylglycinamide ribonucleotide amidotransferase (FGAR-AT) is a 140 kDa bi-functional enzyme involved in a coupled reaction, where the glutaminase active site produces ammonia that is subsequently utilized to convert FGAR to its corresponding amidine in an ATP assisted fashion. The structure of FGAR-AT has been previously determined in an inactive state and the mechanism of activation remains largely unknown. In the current study, hydrophobic cavities were used as markers to identify regions involved in domain movements that facilitate catalytic coupling and subsequent activation of the enzyme. Three internal hydrophobic cavities were located by xenon trapping experiments on FGAR-AT crystals and further, these cavities were perturbed via site-directed mutagenesis. Biophysical characterization of the mutants demonstrated that two of these three voids are crucial for stability and function of the protein, although being ∼20 Šfrom the active centers. Interestingly, correlation analysis corroborated the experimental findings, and revealed that amino acids lining the functionally important cavities form correlated sets (co-evolving residues) that connect these regions to the amidotransferase active center. It was further proposed that the first cavity is transient and allows for breathing motion to occur and thereby serves as an allosteric hotspot. In contrast, the third cavity which lacks correlated residues was found to be highly plastic and accommodated steric congestion by local adjustment of the structure without affecting either stability or activity.