Ensemble refinement shows conformational flexibility in crystal structures of human complement factor D

Structural dynamics and functional cooperativity of human NQO1 by ambient temperature serial crystallography and simulations

The human NQO1 (hNQO1) is a flavin adenine nucleotide (FAD)-dependent oxidoreductase that catalyzes the two-electron reduction of quinones to hydroquinones, being essential for the antioxidant defense system, stabilization of tumor suppressors, and activation of quinone-based chemotherapeutics. Moreover, it is overexpressed in several tumors, which makes it an attractive cancer drug target. To decipher new structural insights into the flavin reductive half-reaction of the catalytic mechanism of hNQO1, we have carried serial crystallography experiments at new ID29 beamline of the ESRF to determine, to the best of our knowledge, the first structure of the hNQO1 in complex with NADH. We have also performed molecular dynamics simulations of free hNQO1 and in complex with NADH. This is the first structural evidence that the hNQO1 functional cooperativity is driven by structural communication between the active sites through long-range propagation of cooperative effects across the hNQO1 structure. Both structural results and MD simulations have supported that the binding of NADH significantly decreases protein dynamics and stabilizes hNQO1 especially at the dimer core and interface. Altogether, these results pave the way for future time-resolved studies, both at x-ray free-electron lasers and synchrotrons, of the dynamics of hNQO1 upon binding to NADH as well as during the FAD cofactor reductive half-reaction. This knowledge will allow us to reveal unprecedented structural information of the relevance of the dynamics during the catalytic function of hNQO1.

Conformational dynamics of adenylate kinase in crystals

Adenylate kinase is a ubiquitous enzyme in living systems and undergoes dramatic conformational changes during its catalytic cycle. For these reasons, it is widely studied by genetic, biochemical, and biophysical methods, both experimental and theoretical. We have determined the basic crystal structures of three differently liganded states of adenylate kinase from Methanotorrus igneus, a hyperthermophilic organism whose adenylate kinase is a homotrimeric oligomer. The multiple copies of each protomer in the asymmetric unit of the crystal provide a unique opportunity to study the variation in the structure and were further analyzed using advanced crystallographic refinement methods and analysis tools to reveal conformational heterogeneity and, thus, implied dynamic behaviors in the catalytic cycle.

MD simulations of macromolecular crystals: Implications for the analysis of Bragg and diffuse scattering

Some of our most detailed information about structure and dynamics of macromolecules comes from X-ray-diffraction studies in crystalline environments. More than 170,000 atomic models have been deposited in the Protein Data Bank, and the number of observations (typically of intensities of Bragg diffraction peaks) is generally quite large, when compared to other experimental methods. Nevertheless, the general agreement between calculated and observed intensities is far outside the experimental precision, and the majority of scattered photons fall between the sharp Bragg peaks, and are rarely taken into account. This chapter considers how molecular dynamics simulations can be used to explore the connections between microscopic behavior in a crystalline lattice and observed scattering intensities, and point the way to new atomic models that could more faithfully recapitulate Bragg intensities and extract useful information from the diffuse scattering that lies between those peaks.

Refinement of multiconformer ensemble models from multi-temperature X-ray diffraction data

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363 K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

Refinement of Multiconformer Ensemble Models from Multi-temperature X-ray Diffraction Data

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

New TSPO Crystal Structures of Mutant and Heme-Bound Forms with Altered Flexibility, Ligand Binding, and Porphyrin Degradation Activity

The ancient protein TSPO (translocator protein 18kD) is found in all kingdoms and was originally identified as a binding site of benzodiazepine drugs. Its physiological function remains unclear, although porphyrins are conserved ligands. Several crystal structures of bacterial TSPO and nuclear magnetic resonance structures of a mouse form have revealed monomer and dimer configurations, but there have been no reports of structures with a physiological ligand. Here, we present the first X-ray structures of Rhodobacter sphaeroides TSPO with a physiological ligand bound. Two different variants (substituting threonine for alanine at position 139 (A139T) and phenylalanine for alanine at position 138 (A138F)) yielded well-diffracting crystals giving structures of both apo- and heme-containing forms. Both variants have wild-type micromolar affinity for heme and protoporphyrin IX, but A139T has very low ability to accelerate the breakdown of porphyrin in the presence of light and oxygen. The binding of heme to one protomer of the dimer of either mutant induces a more rigid structure, both in the heme-binding protomer and the protomer without heme bound, demonstrating an allosteric response. Ensemble refinement of the X-ray data reveals distinct regions of altered flexibility in response to single heme binding to the dimer. The A139T variant shows a more rigid structure overall, which may relate to extra hydrogen bonding of waters captured in the heme crevice. As TSPO has been suggested to have a role in heme delivery from mitochondria to the cytoplasm, the new structures provide potential clues regarding the structural basis of such activity.

Functional Identification of Complement Factor D and Analysis of Its Expression during GCRV Infection in Grass Carp (Ctenopharyngodon idella)

Complement factor D (Df) is a serine protease well known for activating the alternative pathway (AP) in mammals by promoting the cleavage of complement component 3 (C3), thus becoming involved in innate defense. In teleost fish, however, the functional mechanisms of Df in the AP and against pathogen infection are far from clear. In the present study, we cloned and characterized the Df gene, CiDf, from grass carp (Ctenopharyngodon idella) and analyzed its function in promoting C3 cleavage and expression changes after grass carp reovirus (GCRV) infection. The open reading frame of CiDf was found to be 753 bp, encoding 250 amino acids with a molecular mass of 27.06 kDa. CiDf harbors a conserved Tryp_SPc domain, with three conserved residues representing the catalytic triad and three conserved binding sites in the substrate specificity pocket. Pairwise alignment showed that CiDf shares the highest identity (96%) and similarity (98%) with Df from Anabarilius grahami. Phylogenetic analysis indicated that CiDf and other fish Dfs formed a distinct evolutionary branch. Similar to most Dfs from other vertebrates, the CiDf gene structure is characterized by four introns and five exons. The incubation of recombinant CiDf protein with grass carp serum significantly increased the C3b content, demonstrating the conserved function of CiDf in the AP in promoting C3 cleavage, similar to Dfs in mammals. CiDf mRNA expression was widely detected in various tissues and levels were relatively higher in the liver, spleen, and intestine of grass carp. During GCRV infection over a 168-hour period, a high level of CiDf mRNA expression in the liver, spleen, and intestine was maintained at 144 and 168 h, suggesting AP activity at the late stage of GCRV infection. Collectively, the above results reveal the conserved structure and function of CiDf and its distinct expression patterns after GCRV infection, which provide a key basis for studying the roles of Df and AP during GCRV infection in the grass carp C. idella.

Complement Factor D as a Strategic Target for Regulating the Alternative Complement Pathway

The complement system is central to first-line defense against invading pathogens. However, excessive complement activation and/or the loss of complement regulation contributes to the development of autoimmune diseases, systemic inflammation, and thrombosis. One of the three pathways of the complement system, the alternative complement pathway, plays a vital role in amplifying complement activation and pathway signaling. Complement factor D, a serine protease of this pathway that is required for the formation of C3 convertase, is the rate-limiting enzyme. In this review, we discuss the function of factor D within the alternative pathway and its implication in both healthy physiology and disease. Because the alternative pathway has a role in many diseases that are characterized by excessive or poorly mediated complement activation, this pathway is an enticing target for effective therapeutic intervention. Nonetheless, although the underlying disease mechanisms of many of these complement-driven diseases are quite well understood, some of the diseases have limited treatment options or no approved treatments at all. Therefore, in this review we explore factor D as a strategic target for advancing therapeutic control of pathological complement activation.

Exploring ligand dynamics in protein crystal structures with ensemble refinement

Understanding the dynamics of ligands bound to proteins is an important task in medicinal chemistry and drug design. However, the dominant technique for determining protein-ligand structures, X-ray crystallography, does not fully account for dynamics and cannot accurately describe the movements of ligands in protein binding sites. In this article, an alternative method, ensemble refinement, is used on six protein-ligand complexes with the aim of understanding the conformational diversity of ligands in protein crystal structures. The results show that ensemble refinement sometimes indicates that the flexibility of parts of the ligand and some protein side chains is larger than that which can be described by a single conformation and atomic displacement parameters. However, since the electron-density maps are comparable and Rfree values are slightly increased, the original crystal structure is still a better model from a statistical point of view. On the other hand, it is shown that molecular-dynamics simulations and automatic generation of alternative conformations in crystallographic refinement confirm that the flexibility of these groups is larger than is observed in standard refinement. Moreover, the flexible groups in ensemble refinement coincide with groups that give high atomic displacement parameters or non-unity occupancy if optimized in standard refinement. Therefore, the conformational diversity indicated by ensemble refinement seems to be qualitatively correct, indicating that ensemble refinement can be an important complement to standard crystallographic refinement as a tool to discover which parts of crystal structures may show extensive flexibility and therefore are poorly described by a single conformation. However, the diversity of the ensembles is often exaggerated (probably partly owing to the rather poor force field employed) and the ensembles should not be trusted in detail.

Uncovering Membrane-Bound Models of Coagulation Factors by Combined Experimental and Computational Approaches

In the life sciences, including hemostasis and thrombosis, methods of structural biology have become indispensable tools for shedding light on underlying mechanisms that govern complex biological processes. Advancements of the relatively young field of computational biology have matured to a point where it is increasingly recognized as trustworthy and useful, in part due to their high space–time resolution that is unparalleled by most experimental techniques to date. In concert with biochemical and biophysical approaches, computational studies have therefore proven time and again in recent years to be key assets in building or suggesting structural models for membrane-bound forms of coagulation factors and their supramolecular complexes on membrane surfaces where they are activated. Such endeavors and the proposed models arising from them are of fundamental importance in describing and understanding the molecular basis of hemostasis under both health and disease conditions. We summarize the body of work done in this important area of research to drive forward both experimental and computational studies toward new discoveries and potential future therapeutic strategies.

Unique active-site and subsite features in the arabinogalactan-degrading GH43 exo-β-1,3-galactanase from Phanerochaete chrysosporium

Arabinogalactan proteins (AGPs) are plant proteoglycans with functions in growth and development. However, these functions are largely unexplored, mainly because of the complexity of the sugar moieties. These carbohydrate sequences are generally analyzed with the aid of glycoside hydrolases. The exo-β-1,3-galactanase is a glycoside hydrolase from the basidiomycete Phanerochaete chrysosporium (Pc1,3Gal43A), which specifically cleaves AGPs. However, its structure is not known in relation to its mechanism bypassing side chains. In this study, we solved the apo and liganded structures of Pc1,3Gal43A, which reveal a glycoside hydrolase family 43 subfamily 24 (GH43_sub24) catalytic domain together with a carbohydrate-binding module family 35 (CBM35) binding domain. GH43_sub24 is known to lack the catalytic base Asp conserved among other GH43 subfamilies. Our structure in combination with kinetic analyses reveals that the tautomerized imidic acid group of Gln263 serves as the catalytic base residue instead. Pc1,3Gal43A has three subsites that continue from the bottom of the catalytic pocket to the solvent. Subsite -1 contains a space that can accommodate the C-6 methylol of Gal, enabling the enzyme to bypass the β-1,6-linked galactan side chains of AGPs. Furthermore, the galactan-binding domain in CBM35 has a different ligand interaction mechanism from other sugar-binding CBM35s, including those that bind galactomannan. Specifically, we noted a Gly → Trp substitution, which affects pyranose stacking, and an Asp → Asn substitution in the binding pocket, which recognizes β-linked rather than α-linked Gal residues. These findings should facilitate further structural analysis of AGPs and may also be helpful in engineering designer enzymes for efficient biomass utilization.

Discovery of allosteric binding sites by crystallographic fragment screening

Understanding allosteric regulation of proteins is fundamental to our study of protein structure and function. Moreover, allosteric binding pockets have become a major target of drug discovery efforts in recent years. However, even though the function of almost every protein can be influenced by allostery, it remains a challenge to discover, rationalise and validate putative allosteric binding pockets. This review examines how the discovery and analysis of putative allosteric binding sites have been influenced by the availability of centralised facilities for crystallographic fragment screening, along with newly developed computational methods for modelling low occupancy features. We discuss the experimental parameters required for success, and how new methods could influence the field in the future. Finally, we reflect on the general problem of how to translate these findings into actual ligand development programs.

Substrate-recognition mechanism of tomato β-galactosidase 4 using X-ray crystallography and docking simulation

TBG4 recognize multiple linkage types substrates due to having a spatially wide subsite + 1. This feature allows the degradation of AGI, AGII, and AGP leading to the fruit ripening. β-galactosidase (EC 3. 2. 1. 23) catalyzes the hydrolysis of β-galactan and release of D-galactose. Tomato has at least 17 β-galactosidases (TBGs), of which, TBG 4 is responsible for fruit ripening. TBG4 hydrolyzes not only β-1,4-bound galactans, but also β-1,3- and β-1,6-galactans. In this study, we compared each enzyme-substrate complex using X-ray crystallography, ensemble refinement, and docking simulation to understand the broad substrate-specificity of TBG4. In subsite - 1, most interactions were conserved across each linkage type of galactobioses; however, some differences were seen in subsite + 1, owing to the huge volume of catalytic pocket. In addition to this, docking simulation indicated TBG4 to possibly have more positive subsites to recognize and hydrolyze longer galactans. Taken together, our results indicated that during tomato fruit ripening, TBG4 plays an important role by degrading arabinogalactan I (AGI), arabinogalactan II (AGII), and the carbohydrate moiety of arabinogalactan protein (AGP).

Mechanistic and Crystallographic Studies of Azoreductase AzoA from Bacillus wakoensis A01

The azoreductase AzoA from the alkali-tolerant Bacillus wakoensis A01 has been studied to reveal its structural and mechanistic details. For this, a recombinant expression system was developed which yields impressive amounts of fully active enzyme. The purified holo enzyme is remarkably solvent-tolerant and thermostable with an apparent melting temperature of 71 °C. The dimeric enzyme contains FMN as a prosthetic group and is strictly NADH dependent. While AzoA shows a negligible ability to use molecular oxygen as an electron acceptor, it is efficient in reducing various azo dyes and quinones. The kinetic and catalytic mechanism has been studied in detail using steady state kinetic analyses and stopped-flow studies. The data show that AzoA performs quinone and azo dye reductions via a two-electron transfer. Moreover, quinones were shown to be much better substrates (kcat values of 100-400 s-1 for several naphtoquinones) when compared with azo dyes. This suggests that the physiological role of AzoA and sequence-related microbial reductases is linked to quinone reductions and that they can better be annotated as quinone reductases. The structure of AzoA has been determined in complex with FMN at 1.8 Å resolution. AzoA displays unique features in the active site providing clues for explaining its catalytic and thermostability features. An uncommon loop, when compared with sequence-related reductases, forms an active site lid with Trp60 acting as an anchor. Several Trp60 mutants have been analyzed disclosing an important role of this residue in the stability of AzoA, while they retained activity. Structural details are discussed in relation to other azo and quinone reductases. This study provides new insights into the molecular functioning of AzoA and sequence-related reductases.

Allostery in Coagulation Factor VIIa Revealed by Ensemble Refinement of Crystallographic Structures

A critical step in injury-induced initiation of blood coagulation is the formation of the complex between the trypsin-like protease coagulation factor VIIa (FVIIa) and its cofactor tissue factor (TF), which converts FVIIa from an intrinsically poor enzyme to an active protease capable of activating zymogens of downstream coagulation proteases. Unlike its constitutively active ancestor trypsin, FVIIa is allosterically activated (by TF). Here, ensemble refinement of crystallographic structures, which uses multiple copies of the entire structure as a means of representing structural flexibility, is applied to explore the impacts of inhibitor binding to trypsin and FVIIa, as well as cofactor binding to FVIIa. To assess the conformational flexibility and its role in allosteric pathways in these proteases, main-chain hydrogen bond networks are analyzed by calculating the hydrogen-bond propensity. Mapping pairwise propensity differences between relevant structures shows that binding of the inhibitor benzamidine to trypsin has a minor influence on the protease flexibility. For FVIIa, in contrast, the protease domain is "locked" into the catalytically competent trypsin-like configuration upon benzamidine binding as indicated by the stabilization of key structural features: the nonprime binding cleft and the oxyanion hole are stabilized, and the effect propagates from the active site region to the calcium-binding site and to the vicinity of the disulphide bridge connecting with the light chain. TF binding to FVIIa furthermore results in stabilization of the 170 loop, which in turn propagates an allosteric signal from the TF-binding region to the active site. Analyses of disulphide bridge energy and flexibility reflect the striking stability difference between the unregulated enzyme and the allosterically activated form after inhibitor or cofactor binding. The ensemble refinement analyses show directly, for the first time to our knowledge, whole-domain structural footprints of TF-induced allosteric networks present in x-ray crystallographic structures of FVIIa, which previously only have been hypothesized or indirectly inferred.

Integrating Experiment and Theory to Understand TCR-pMHC Dynamics

The conformational dynamism of proteins is well established. Rather than having a single structure, proteins are more accurately described as a conformational ensemble that exists across a rugged energy landscape, where different conformational sub-states interconvert. The interaction between αβ T cell receptors (TCR) and cognate peptide-MHC (pMHC) is no exception, and is a dynamic process that involves substantial conformational change. This review focuses on technological advances that have begun to establish the role of conformational dynamics and dynamic allostery in TCR recognition of the pMHC and the early stages of signaling. We discuss how the marriage of molecular dynamics (MD) simulations with experimental techniques provides us with new ways to dissect and interpret the process of TCR ligation. Notably, application of simulation techniques lags behind other fields, but is predicted to make substantial contributions. Finally, we highlight integrated approaches that are being used to shed light on some of the key outstanding questions in the early events leading to TCR signaling.

Discovery and Design of First Benzylamine-Based Ligands Binding to an Unlocked Conformation of the Complement Factor D

Complement Factor D, a serine protease of the S1 family and key component of the alternative pathway amplification loop, represents a promising target for the treatment of several prevalent and rare diseases linked to the innate immune system. Previously reported FD inhibitors have been shown to bind to the FD active site in its self-inhibited conformation characterized by the presence of a salt bridge at the bottom of the S1 pocket between Asp189 and Arg218. We report herein a new set of small-molecule FD ligands that harbor a basic S1 binding moiety directly binding to the carboxylate of Asp189, thereby displacing the Asp189-Arg218 ionic interaction and significantly changing the conformation of the self-inhibitory loop.

The maximum penalty criterion for ridge regression: application to the calibration of the force constant in elastic network models

Tikhonov regularization, or ridge regression, is a popular technique to deal with collinearity in multivariate regression. We unveil a formal analogy between ridge regression and statistical mechanics, where the objective function is comparable to a free energy, and the ridge parameter plays the role of temperature. This analogy suggests two novel criteria for selecting a suitable ridge parameter: specific-heat (C_v) and maximum penalty (MP). We apply these fits to evaluate the relative contributions of rigid-body and internal fluctuations, which are typically highly collinear, to crystallographic B-factors. This issue is particularly important for computational models of protein dynamics, such as the elastic network model (ENM), since the amplitude of the predicted internal motion is commonly calibrated using B-factor data. After validation on simulated datasets, our results indicate that rigid-body motions account on average for more than 80% of the amplitude of B-factors. Furthermore, we evaluate the ability of different fits to reproduce the amplitudes of internal fluctuations in X-ray ensembles from the B-factors in the corresponding single X-ray structures. The new ridge criteria are shown to be markedly superior to the commonly used two-parameter fit that neglects rigid-body rotations and to the full fits regularized under generalized cross-validation. In conclusion, the proposed fits ensure a more robust calibration of the ENM force constant and should prove valuable in other applications.

Previously Hidden Dynamics at the TCR-Peptide-MHC Interface Revealed

A structural characterization of the interaction between αβ TCRs and cognate peptide-MHC (pMHC) is central to understanding adaptive T cell-mediated immunity. X-ray crystallography, although the source of much structural data, traditionally provides only a static snapshot of the protein. Given the emerging evidence for the important role of conformational dynamics in protein function, we interrogated 309 crystallographic structures of pMHC complexes using ensemble refinement, a technique that can extract dynamic information from the x-ray data. Focusing on a subset of human pMHC class I systems, we found that in many cases, ensemble methods were able to uncover previously hidden evidence of significant conformational plasticity, thereby revealing additional information that can build upon and significantly enhance functional interpretations that are based on a single static structure. Notable examples include the interpretation of differences in the disease association of HLA subtypes, the relationship between peptide prominence and TCR recognition, the role of conformational flexibility in vaccine design, and the discrimination between induced fit and conformational selection models of TCR binding. We show that the currently widespread practice of analyzing pMHC interactions via the study of a single crystallographic structure does not make use of pertinent and easily accessible information from x-ray data concerning alternative protein conformations. This new analysis therefore not only highlights the capacity for ensemble methods to significantly enrich the interpretation of decades of structural data but also provides previously missing information concerning the dynamics of existing characterized TCR-pMHC interactions.

Protein crystal lattices are dynamic assemblies: the role of conformational entropy in the protein condensed phase

Until recently, the occurrence of conformational entropy in protein crystal contacts was considered to be a very unlikely event. A study based on the most accurately refined protein structures demonstrated that side-chain conformational entropy and static disorder might be common in protein crystal lattices. The present investigation uses structures refined using ensemble refinement to show that although paradoxical, conformational entropy is likely to be the major factor in the emergence and integrity of the protein condensed phase. This study reveals that the role of shape entropy and local entropic forces expands beyond the onset of crystallization. For the first time, the complete pattern of intermolecular interactions by protein atoms in crystal lattices is presented, which shows that van der Waals interactions dominate in crystal formation.

Molecular Mechanisms of Drug Action: X-ray Crystallography at the Basis of Structure-based and Ligand-based Drug Design

Biological systems are recognized for their complexity and diversity and yet we sometimes manage to cure disease via the administration of small chemical drug molecules. At first, active ingredients were found accidentally and at that time there did not seem a need to understand the molecular mechanism of drug functioning. However, the urge to develop new drugs, the discovery of multipurpose characteristics of some drugs, and the necessity to remove unwanted secondary drug effects, incited the pharmaceutical sector to rationalize drug design. This did not deliver success in the years directly following its conception, but it drove the evolution of biochemical and biophysical techniques to enable the characterization of molecular mechanisms of drug action. Functional and structural data generated by biochemists and structural biologists became a valuable input for computational biologists, chemists and bioinformaticians who could extrapolate in silico, based on variations in the structural aspects of the drug molecules and their target. This opened up new avenues with much improved predictive power because of a clearer perception of the role and impact of structural elements in the intrinsic affinity and specificity of the drug for its target. In this chapter, we review how crystal structures can initiate structure-based drug design in general.

Determinants of Macromolecular Specificity from Proteomics-Derived Peptide Substrate Data

BACKGROUND

Recent advances in proteomics methodologies allow for high throughput profiling of proteolytic cleavage events. The resulting substrate peptide distributions provide deep insights in the underlying macromolecular recognition events, as determinants of biomolecular specificity identified by proteomics approaches may be compared to structure-based analysis of corresponding protein-protein interfaces.

METHOD

Here, we present an overview of experimental and computational methodologies and tools applied in the area and provide an outlook beyond the protein class of proteases.

RESULTS AND CONCLUSION

We discuss here future potential, synergies and needs of the emerging overlap disciplines of proteomics and structure-based modelling.

Identifying and Visualizing Macromolecular Flexibility in Structural Biology

Structural biology comprises a variety of tools to obtain atomic resolution data for the investigation of macromolecules. Conventional structural methodologies including crystallography, NMR and electron microscopy often do not provide sufficient details concerning flexibility and dynamics, even though these aspects are critical for the physiological functions of the systems under investigation. However, the increasing complexity of the molecules studied by structural biology (including large macromolecular assemblies, integral membrane proteins, intrinsically disordered systems, and folding intermediates) continuously demands in-depth analyses of the roles of flexibility and conformational specificity involved in interactions with ligands and inhibitors. The intrinsic difficulties in capturing often subtle but critical molecular motions in biological systems have restrained the investigation of flexible molecules into a small niche of structural biology. Introduction of massive technological developments over the recent years, which include time-resolved studies, solution X-ray scattering, and new detectors for cryo-electron microscopy, have pushed the limits of structural investigation of flexible systems far beyond traditional approaches of NMR analysis. By integrating these modern methods with powerful biophysical and computational approaches such as generation of ensembles of molecular models and selective particle picking in electron microscopy, more feasible investigations of dynamic systems are now possible. Using some prominent examples from recent literature, we review how current structural biology methods can contribute useful data to accurately visualize flexibility in macromolecular structures and understand its important roles in regulation of biological processes.

Releasing the brakes in coagulation Factor IXa by co-operative maturation of the substrate-binding site

Coagulation Factor IX is positioned at the merging point of the intrinsic and extrinsic blood coagulation cascades. Factor IXa (activated Factor IX) serves as the trigger for amplification of coagulation through formation of the so-called Xase complex, which is a ternary complex of Factor IXa, its substrate Factor X and the cofactor Factor VIIIa on the surface of activated platelets. Within the Xase complex the substrate turnover by Factor IXa is enhanced 200000-fold; however, the mechanistic and structural basis for this dramatic enhancement remains only partly understood. A multifaceted approach using enzymatic, biophysical and crystallographic methods to evaluate a key set of activity-enhanced Factor IXa variants has demonstrated a delicately balanced bidirectional network. Essential molecular interactions across multiple regions of the Factor IXa molecule co-operate in the maturation of the active site. This maturation is specifically facilitated by long-range communication through the Ile212–Ile213 motif unique to Factor IXa and a flexibility of the 170-loop that is further dependent on the conformation in the Cys168–Cys182 disulfide bond. Ultimately, the network consists of compensatory brakes (Val16 and Ile213) and accelerators (Tyr99 and Phe174) that together allow for a subtle fine-tuning of enzymatic activity.

E pluribus unum, no more: from one crystal, many conformations

Several distinct computational approaches have recently been implemented to represent conformational heterogeneity from X-ray crystallography datasets that are averaged in time and space. As these modeling methods mature, newly discovered alternative conformations are being used to derive functional protein mechanisms. Room temperature X-ray data collection is emerging as a key variable for sampling functionally relevant conformations also observed in solution studies. Although concerns about radiation damage are warranted with higher temperature data collection, 'diffract and destroy' strategies on X-ray free electron lasers may permit radiation damage-free data collection. X-ray crystallography need not be confined to 'static unique snapshots'; these experimental and computational advances are revealing how the many conformations populated within a single crystal are used in biological mechanisms.