Structures Of Related Proteins Research Articles

Alignment of distantly related protein structures: algorithm, bound and implications to homology modeling

Building an accurate alignment of a large set of distantly related protein structures is still very challenging. This article presents a novel method 3DCOMB that can generate a multiple structure alignment (MSA) with not only as many conserved cores as possible, but also high-quality pairwise alignments. 3DCOMB is unique in that it makes use of both local and global structure environments, combined by a statistical learning method, to accurately identify highly similar fragment blocks (HSFBs) among all proteins to be aligned. By extending the alignments of these HSFBs, 3DCOMB can quickly generate an accurate MSA without using progressive alignment. 3DCOMB significantly excels others in aligning distantly related proteins. 3DCOMB can also generate correct alignments for functionally similar regions among proteins of very different structures while many other MSA tools fail. 3DCOMB is useful for many real-world applications. In particular, it enables us to find out that there is still large improvement room for multiple template homology modeling while several other MSA tools fail to do so. 3DCOMB is available at http://ttic.uchicago.edu/~jinbo/software.htm. jinboxu@gmail.com Supplementary data are available at Bioinformatics online.

Structural Analysis of Papain-Like NlpC/P60 Superfamily Enzymes with a Circularly Permuted Topology Reveals Potential Lipid Binding Sites

NlpC/P60 superfamily papain-like enzymes play important roles in all kingdoms of life. Two members of this superfamily, LRAT-like and YaeF/YiiX-like families, were predicted to contain a catalytic domain that is circularly permuted such that the catalytic cysteine is located near the C-terminus, instead of at the N-terminus. These permuted enzymes are widespread in virus, pathogenic bacteria, and eukaryotes. We determined the crystal structure of a member of the YaeF/YiiX-like family from Bacillus cereus in complex with lysine. The structure, which adopts a ligand-induced, “closed” conformation, confirms the circular permutation of catalytic residues. A comparative analysis of other related protein structures within the NlpC/P60 superfamily is presented. Permutated NlpC/P60 enzymes contain a similar conserved core and arrangement of catalytic residues, including a Cys/His-containing triad and an additional conserved tyrosine. More surprisingly, permuted enzymes have a hydrophobic S1 binding pocket that is distinct from previously characterized enzymes in the family, indicative of novel substrate specificity. Further analysis of a structural homolog, YiiX (PDB 2if6) identified a fatty acid in the conserved hydrophobic pocket, thus providing additional insights into possible function of these novel enzymes.

A tetrameric structure is not essential for activity in dihydrodipicolinate synthase (DHDPS) from Mycobacterium tuberculosis

Dihydrodipicolinate synthase (DHDPS) is a validated antibiotic target for which a new approach to inhibitor design has been proposed: disrupting native tetramer formation by targeting the dimer–dimer interface. In this study, rational design afforded a variant of Mycobacterium tuberculosis, Mtb-DHDPS-A204R, with disrupted quaternary structure. X-ray crystallography (at a resolution of 2.1Å) revealed a dimeric protein with an identical fold and active-site structure to the tetrameric wild-type enzyme. Analytical ultracentrifugation confirmed the dimeric structure in solution, yet the dimeric mutant has similar activity to the wild-type enzyme. Although the affinity for both substrates was somewhat decreased, the high catalytic competency of the enzyme was surprising in the light of previous results showing that dimeric variants of the Escherichia coli and Bacillus anthracis DHDPS enzymes have dramatically reduced activity compared to their wild-type tetrameric counterparts. These results suggest that Mtb-DHDPS-A204R is similar to the natively dimeric enzyme from Staphylococcus aureus, and highlight our incomplete understanding of the role played by oligomerisation in relating protein structure and function.

The Structural Biology Knowledgebase: a portal to protein structures, sequences, functions, and methods

The Protein Structure Initiative’s Structural Biology Knowledgebase (SBKB, URL: http://sbkb.org) is an open web resource designed to turn the products of the structural genomics and structural biology efforts into knowledge that can be used by the biological community to understand living systems and disease. Here we will present examples on how to use the SBKB to enable biological research. For example, a protein sequence or Protein Data Bank (PDB) structure ID search will provide a list of related protein structures in the PDB, associated biological descriptions (annotations), homology models, structural genomics protein target status, experimental protocols, and the ability to order available DNA clones from the PSI:Biology-Materials Repository. A text search will find publication and technology reports resulting from the PSI’s high-throughput research efforts. Web tools that aid in research, including a system that accepts protein structure requests from the community, will also be described. Created in collaboration with the Nature Publishing Group, the Structural Biology Knowledgebase monthly update also provides a research library, editorials about new research advances, news, and an events calendar to present a broader view of structural genomics and structural biology.

Identification of Local Conformational Similarity in Structurally Variable Regions of Homologous Proteins Using Protein Blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to α-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the “structurally variable” regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of ‘variable’ regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.

Packing of the Extracellular Domain Hydrophobic Core Has Evolved to Facilitate Pentameric Ligand-gated Ion Channel Function

Protein function depends on conformational flexibility and folding stability. Loose packing of hydrophobic cores is not infrequent in proteins, as the enhanced flexibility likely contributes to their biological function. Here, using experimental and computational approaches, we show that eukaryotic pentameric ligand-gated ion channels are characterized by loose packing of their extracellular domain β-sandwich cores, and that loose packing contributes to their ability to rapidly switch from closed to open channel states in the presence of ligand. Functional analyses of GABA(A) receptors show that increasing the β-core packing disrupted GABA-mediated currents, with impaired GABA efficacy and slowed GABA current activation and desensitization. We propose that loose packing of the hydrophobic β-core developed as an evolutionary strategy aimed to facilitate the allosteric mechanisms of eukaryotic pentameric ligand-gated ion channels.

The structural basis of function in Cys-loop receptors

Cys-loop receptors are membrane-spanning neurotransmitter-gated ion channels that are responsible for fast excitatory and inhibitory transmission in the peripheral and central nervous systems. The best studied members of the Cys-loop family are nACh, 5-HT3, GABAA and glycine receptors. All these receptors share a common structure of five subunits, pseudo-symmetrically arranged to form a rosette with a central ion-conducting pore. Some are cation selective (e.g. nACh and 5-HT3) and some are anion selective (e.g. GABAA and glycine). Each receptor has an extracellular domain (ECD) that contains the ligand-binding sites, a transmembrane domain (TMD) that allows ions to pass across the membrane, and an intracellular domain (ICD) that plays a role in channel conductance and receptor modulation. Cys-loop receptors are the targets for many currently used clinically relevant drugs (e.g. benzodiazepines and anaesthetics). Understanding the molecular mechanisms of these receptors could therefore provide the catalyst for further development in this field, as well as promoting the development of experimental techniques for other areas of neuroscience.In this review, we present our current understanding of Cys-loop receptor structure and function. The ECD has been extensively studied. Research in this area has been stimulated in recent years by the publication of high-resolution structures of nACh receptors and related proteins, which have permitted the creation of many Cys loop receptor homology models of this region. Here, using the 5-HT3 receptor as a typical member of the family, we describe how homology modelling and ligand docking can provide useful but not definitive information about ligand interactions. We briefly consider some of the many Cys-loop receptors modulators. We discuss the current understanding of the structure of the TMD, and how this links to the ECD to allow channel gating, and consider the roles of the ICD, whose structure is poorly understood. We also describe some of the current methods that are beginning to reveal the differences between different receptor states, and may ultimately show structural details of transitions between them.

Structure of Pleiotrophin- and Hepatocyte Growth Factor-binding Sulfated Hexasaccharide Determined by Biochemical and Computational Approaches

Endogenous pleiotrophin and hepatocyte growth factor (HGF) mediate the neurite outgrowth-promoting activity of chondroitin sulfate (CS)/dermatan sulfate (DS) hybrid chains isolated from embryonic pig brain. CS/DS hybrid chains isolated from shark skin have a different disaccharide composition, but also display these activities. In this study, pleiotrophin- and HGF-binding domains in shark skin CS/DS were investigated. A high affinity CS/DS fraction was isolated using a pleiotrophin-immobilized column. It showed marked neurite outgrowth-promoting activity and strong inhibitory activity against the binding of pleiotrophin to immobilized CS/DS chains from embryonic pig brain. The inhibitory activity was abolished by chondroitinase ABC or B, and partially reduced by chondroitinase AC-I. A pentasulfated hexasaccharide with a novel structure was isolated from the chondroitinase AC-I digest using pleiotrophin affinity and anion exchange chromatographies. It displayed a potent inhibitory effect on the binding of HGF to immobilized shark skin CS/DS chains, suggesting that the pleiotrophin- and HGF-binding domains at least partially overlap in the CS/DS chains involved in the neuritogenic activity. Computational chemistry using molecular modeling and calculations of the electrostatic potential of the hexasaccharide and two pleiotrophin-binding octasaccharides previously isolated from CS/DS hybrid chains of embryonic pig brain identified an electronegative zone potentially involved in the molecular recognition of the oligosaccharides by pleiotrophin. Homology modeling of pleiotrophin based on a related midkine protein structure predicted the binding pocket of pleiotrophin for the oligosaccharides and provided new insights into the molecular mechanism of the interactions between the oligosaccharides and pleiotrophin.

Docking and scoring protein interactions: CAPRI 2009

Protein docking algorithms are assessed by evaluating blind predictions performed during 2007-2009 in Rounds 13-19 of the community-wide experiment on critical assessment of predicted interactions (CAPRI). We evaluated the ability of these algorithms to sample docking poses and to single out specific association modes in 14 targets, representing 11 distinct protein complexes. These complexes play important biological roles in RNA maturation, G-protein signal processing, and enzyme inhibition and function. One target involved protein-RNA interactions not previously considered in CAPRI, several others were hetero-oligomers, or featured multiple interfaces between the same protein pair. For most targets, predictions started from the experimentally determined structures of the free (unbound) components, or from models built from known structures of related or similar proteins. To succeed they therefore needed to account for conformational changes and model inaccuracies. In total, 64 groups and 12 web-servers submitted docking predictions of which 4420 were evaluated. Overall our assessment reveals that 67% of the groups, more than ever before, produced acceptable models or better for at least one target, with many groups submitting multiple high- and medium-accuracy models for two to six targets. Forty-one groups including four web-servers participated in the scoring experiment with 1296 evaluated models. Scoring predictions also show signs of progress evidenced from the large proportion of correct models submitted. But singling out the best models remains a challenge, which also adversely affects the ability to correctly rank docking models. With the increased interest in translating abstract protein interaction networks into realistic models of protein assemblies, the growing CAPRI community is actively developing more efficient and reliable docking and scoring methods for everyone to use.

Crystal structure of a fructokinase homolog from Halothermothrix orenii

Fructokinase (FRK; EC 2.7.1.4) catalyzes the phosphorylation of d-fructose to d-fructose 6-phosphate (F6P). This irreversible and near rate-limiting step is a central and regulatory process in plants and bacteria, which channels fructose into a metabolically active state for glycolysis. Towards understanding the mechanism of FRK, here we report the crystal structure of a FRK homolog from a thermohalophilic bacterium Halothermothrixorenii (Hore_18220 in sequence databases). The structure of the Hore_18220 protein reveals a catalytic domain with a Rossmann-like fold and a beta-sheet "lid" for dimerization. Based on comparison of Hore_18220 to structures of related proteins, we propose its mechanism of action, in which the lid serves to regulate access to the substrate binding sites. Close relationship of Hore_18220 and plant FRK enzymes allows us to propose a model for the structure and function of FRKs.

Structural classification by the Lipase Engineering Database: a case study of Candida antarctica lipase A

BackgroundThe Lipase Engineering Database (LED) integrates information on sequence, structure and function of lipases, esterases and related proteins with the α/β hydrolase fold. A new superfamily for Candida antarctica lipase A (CALA) was introduced including the recently published crystal structure of CALA. Since CALA has a highly divergent sequence in comparison to other α/β hydrolases, the Lipase Engineering Database was used to classify CALA in the frame of the already established classification system. This involved the comparison of CALA to similar structures as well as sequence-based comparisons against the content of the LED.ResultsThe new release 3.0 (December 2009) of the Lipase Engineering Database contains 24783 sequence entries for 18585 proteins as well as 656 experimentally determined protein structures, including the structure of CALA. In comparison to the previous release [1] with 4322 protein and 167 structure entries this update represents a significant increase in data volume. By comparing CALA to representative structures from all superfamilies, a structure from the deacetylase superfamily was found to be most similar to the structure of CALA. While the α/β hydrolase fold is conserved in both proteins, the major difference is found in the cap region. Sequence alignments between both proteins show a sequence similarity of only 15%. A multisequence alignment of both protein families was used to create hidden Markov models for the cap region of CALA and showed that the cap region of CALA is unique among all other proteins of the α/β hydrolase fold. By specifically comparing the substrate binding pocket of CALA to other binding pockets of α/β hydrolases, the binding pocket of Candida rugosa lipase was identified as being highly similar. This similarity also applied to the lid of Candida rugosa lipase in comparison to the potential lid of CALA.ConclusionThe LED serves as a valuable tool for the systematic analysis of single proteins or protein families. The updated release 3.0 was used for the evaluation of α/β hydrolases. The HTML version of the database with new features is available at http://www.led.uni-stuttgart.de and provides sequences, structures and a set of analysis tools including phylogenetic trees and HMM profiles

IMGT/3Dstructure-DB and IMGT/DomainGapAlign: a database and a tool for immunoglobulins or antibodies, T cell receptors, MHC, IgSF and MhcSF

IMGT/3Dstructure-DB is the three-dimensional (3D) structure database of IMGT®, the international ImMunoGenetics information system® that is acknowledged as the global reference in immunogenetics and immunoinformatics. IMGT/3Dstructure-DB contains 3D structures of immunoglobulins (IG) or antibodies, T cell receptors (TR), major histocompatibility complex (MHC) proteins, antigen receptor/antigen complexes (IG/Ag, TR/peptide/MHC) of vertebrates; 3D structures of related proteins of the immune system (RPI) of vertebrates and invertebrates, belonging to the immunoglobulin and MHC superfamilies (IgSF and MhcSF, respectively) and found in complexes with IG, TR or MHC. IMGT/3Dstructure-DB data are annotated according to the IMGT criteria, using IMGT/DomainGapAlign, and based on the IMGT-ONTOLOGY concepts and axioms. IMGT/3Dstructure-DB provides IMGT gene and allele identification (CLASSIFICATION), region and domain delimitations (DESCRIPTION), amino acid positions according to the IMGT unique numbering (NUMEROTATION) that are used in IMGT/3Dstructure-DB cards, results of contact analysis and renumbered flat files. In its Web version, the IMGT/DomainGapAlign tool analyses amino acid sequences, per domain. Coupled to the IMGT/Collier-de-Perles tool, it provides an invaluable help for antibody engineering and humanization design based on complementarity determining region (CDR) grafting as it precisely defines the standardized framework regions (FR-IMGT) and CDR-IMGT. IMGT/3Dstructure-DB and IMGT/DomainGapAlign are freely available at http://www.imgt.org.

Prediction of the spatial structure of proteins: Emphasis on membrane targets

Knowledge of the spatial structure of proteins is a prerequisite for both awareness of their functional mechanisms and the framework for rational drug discovery and design. Meanwhile, direct structural determination is often hampered or impractical due to the complexity, expensiveness, and limited capabilities of experimental techniques. These issues are especially pronounced for integral membrane proteins. On numerous occasions, the theoretical prediction of protein structures may facilitate the process by exploiting physical or empirical principles. This paper surveys modern techniques for the prediction of the spatial structure of proteins using computer algorithms, and the main emphasis is placed on the most "complex" targets - membrane proteins (MPs). The first part of the review describes de novo methods based on empirical physical principles; in the second part, a comparative modeling philosophy, which accounts for the structure of related proteins, is described. Special focus is made regarding pharmacologically relevant classes of G-coupled receptors, receptor tyrosine ki-nases, and other MPs. Algorithms for the assessment of the models quality and potential fields of application of computer models are discussed.

Structural Basis of Activation of Cys-Loop Receptors: the Extracellular–Transmembrane Interface as a Coupling Region

Cys-loop receptors mediate rapid transmission throughout the nervous system by converting a chemical signal into an electric one. They are pentameric proteins with an extracellular domain that carries the transmitter binding sites and a transmembrane region that forms the ion pore. Their essential function is to couple the binding of the agonist at the extracellular domain to the opening of the ion pore. How the structural changes elicited by agonist binding are propagated through a distance of 50 A to the gate is therefore central for the understanding of the receptor function. A step forward toward the identification of the structures involved in gating has been given by the recently elucidated high-resolution structures of Cys-loop receptors and related proteins. The extracellular-transmembrane interface has attracted attention because it is a structural transition zone where beta-sheets from the extracellular domain merge with alpha-helices from the transmembrane domain. Within this zone, several regions form a network that relays structural changes from the binding site toward the pore, and therefore, this interface controls the beginning and duration of a synaptic response. In this review, the most recent findings on residues and pairwise interactions underlying channel gating are discussed, the main focus being on the extracellular-transmembrane interface.

Building and assessing atomic models of proteins from structural templates: learning and benchmarks.

One approach to predict a protein fold from a sequence (a target) is based on structures of related proteins that are used as templates. We present an algorithm that examines a set of candidates for templates, builds from each of the templates an atomically detailed model, and ranks the models. The algorithm performs a hierarchical selection of the best model using a diverse set of signals. After a quick and suboptimal screening of template candidates from the protein data bank, the current method fine-tunes the selection to a few models. More detailed signals test the compatibility of the sequence and the proposed structures, and are merged to give a global fitness measure using linear programming. This algorithm is a component of the prediction server LOOPP (http://www.loopp.org). Large-scale training and tests sets were designed and are presented. Recent results of the LOOPP server in CASP8 are discussed.

Heat-induced protein structure and subfractions in relation to protein degradation kinetics and intestinal availability in dairy cattle

The objectives of this study were to reveal protein structures of feed tissues affected by heat processing at a cellular level, using the synchrotron-based Fourier transform infrared microspectroscopy as a novel approach, and quantify protein structure in relation to protein digestive kinetics and nutritive value in the rumen and intestine in dairy cattle. The parameters assessed included 1) protein structure α-helix to β-sheet ratio; 2) protein subfractions profiles; 3) protein degradation kinetics and effective degradability; 4) predicted nutrient supply using the intestinally absorbed protein supply (DVE)/degraded protein balance (OEB) system for dairy cattle. In this study, Vimy flaxseed protein was used as a model feed protein and was autoclave-heated at 120°C for 20, 40, and 60min in treatments T1, T2, and T3, respectively. The results showed that using the synchrotron-based Fourier transform infrared microspectroscopy revealed and identified the heat-induced protein structure changes. Heating at 120°C for 40 and 60min increased the protein structure α-helix to β-sheet ratio. There were linear effects of heating time on the ratio. The heating also changed chemical profiles, which showed soluble CP decreased upon heating with concomitant increases in nonprotein nitrogen, neutral, and acid detergent insoluble nitrogen. The protein subfractions with the greatest changes were PB1, which showed a dramatic reduction, and PB2, which showed a dramatic increase, demonstrating a decrease in overall protein degradability. In situ results showed a reduction in rumen-degradable protein and in rumen-degradable dry matter without differences between the treatments. Intestinal digestibility, determined using a 3-step in vitro procedure, showed no changes to rumen undegradable protein. Modeling results showed that heating increased total intestinally absorbable protein (feed DVE value) and decreased degraded protein balance (feed OEB value), but there were no differences between the treatments. There was a linear effect of heating time on the DVE and a cubic effect on the OEB value. Our results showed that heating changed chemical profiles, protein structure α-helix to β-sheet ratio, and protein subfractions; decreased rumen-degradable protein and rumen-degradable dry matter; and increased potential nutrient supply to dairy cattle. The protein structure α-helix to β-sheet ratio had a significant positive correlation with total intestinally absorbed protein supply and negative correlation with degraded protein balance.

Improving protein structure similarity searches using domain boundaries based on conserved sequence information

BackgroundThe identification of protein domains plays an important role in protein structure comparison. Domain query size and composition are critical to structure similarity search algorithms such as the Vector Alignment Search Tool (VAST), the method employed for computing related protein structures in NCBI Entrez system. Currently, domains identified on the basis of structural compactness are used for VAST computations. In this study, we have investigated how alternative definitions of domains derived from conserved sequence alignments in the Conserved Domain Database (CDD) would affect the domain comparisons and structure similarity search performance of VAST.ResultsAlternative domains, which have significantly different secondary structure composition from those based on structurally compact units, were identified based on the alignment footprints of curated protein sequence domain families. Our analysis indicates that domain boundaries disagree on roughly 8% of protein chains in the medium redundancy subset of the Molecular Modeling Database (MMDB). These conflicting sequence based domain boundaries perform slightly better than structure domains in structure similarity searches, and there are interesting cases when structure similarity search performance is markedly improved.ConclusionStructure similarity searches using domain boundaries based on conserved sequence information can provide an additional method for investigators to identify interesting similarities between proteins with known structures. Because of the improvement in performance of structure similarity searches using sequence domain boundaries, we are in the process of implementing their inclusion into the VAST search and MMDB resources in the NCBI Entrez system.

Predicting the Accuracy of Ligand Overlay Methods with Random Forest Models

The accuracy of binding mode prediction using standard molecular overlay methods (ROCS, FlexS, Phase, and FieldCompare) is studied. Previous work has shown that simple decision tree modeling can be used to improve accuracy by selection of the best overlay template. This concept is extended to the use of Random Forest (RF) modeling for template and algorithm selection. An extensive data set of 815 ligand-bound X-ray structures representing 5 gene families was used for generating ca. 70,000 overlays using four programs. RF models, trained using standard measures of ligand and protein similarity and Lipinski-related descriptors, are used for automatically selecting the reference ligand and overlay method maximizing the probability of reproducing the overlay deduced from X-ray structures (i.e., using rmsd < or = 2 A as the criteria for success). RF model scores are highly predictive of overlay accuracy, and their use in template and method selection produces correct overlays in 57% of cases for 349 overlay ligands not used for training RF models. The inclusion in the models of protein sequence similarity enables the use of templates bound to related protein structures, yielding useful results even for proteins having no available X-ray structures.

CONSERVED WATER MOLECULES IN X-RAY STRUCTURES HIGHLIGHT THE ROLE OF WATER IN INTRAMOLECULAR AND INTERMOLECULAR INTERACTIONS

Water molecules immobilized on a protein or DNA surface are known to play an important role in intramolecular and intermolecular interactions. Comparative analysis of related three-dimensional (3D) structures allows to predict the locations of such water molecules on the protein surface. We have developed and implemented the algorithm WLAKE detecting "conserved" water molecules, i.e. those located in almost the same positions in a set of superimposed structures of related proteins or macromolecular complexes. The problem is reduced to finding maximal cliques in a certain graph. Despite exponential algorithm complexity, the program works appropriately fast for dozens of superimposed structures. WLAKE was used to predict functionally significant water molecules in enzyme active sites (transketolases) as well as in intermolecular (ETS-DNA complexes) and intramolecular (thiol-disulfide interchange protein) interactions. The program is available online at http://monkey.belozersky.msu.ru/~evgeniy/wLake/wLake.html.

Hydrophobic Core Mutations in CI2 Globally Perturb Fast Side-Chain Dynamics Similarly without Regard to Position

Protein dynamics is currently an area of intense research because of its importance as complementary information to the huge quantity of available data relating protein structure and function. Because it is usually the influence of dynamics on function that is studied, the physical determinants of the distribution of flexibility in proteins have not been explored as thoroughly. In the present NMR study, an expanded suite of five (2)H relaxation experiments was used to characterize the picosecond-to-nanosecond side-chain dynamics of chymotrypsin inhibitor 2 (CI2) and five hydrophobic core mutants, some of which are members of the folding nucleus. Because CI2 is a homologue of the serine protease inhibitor eglin c, which has already been extensively characterized in terms of its dynamics, it was possible to compare not only side-chain dynamics but also the responses of these dynamics to analogous mutations. Remarkably, each of the five core mutations in CI2 led to similar and reproducible increases in side-chain flexibility throughout the entire structure. Although the expanded suite of (2)H relaxation experiments does not affect model selection for the vast majority of residues, it did enable the detection of increasing levels of nanosecond-scale motions in CI2's reactive site binding loop as the L68 side chain was progressively shortened by mutation. Collectively, we observed that the CI2 mutants are more dynamically similar to each other than to the more rigid wild-type CI2, from which we propose that wild-type CI2 has been optimized to a specific level of rigidity which may aid in its function as a serine protease inhibitor. We also observed that the pattern of side-chain dynamics of CI2 is quantitatively similar to eglin c, but that this similarity is lost upon mutating both proteins at an equivalent position. Finally, (15)N relaxation was used to characterize the backbone dynamics of wild-type and mutant CI2. Interestingly, mutation at folding nucleus positions led to widespread increases in backbone flexibility, whereas non-folding-nucleus positions led to increases in flexibility in the C-terminal half of the protein only.