Protein flexibility: coordinate uncertainties and interpretation of structural differences

AlphaFold2 Can Predict Single-Mutation Effects

AlphaFold2 (AF) is a promising tool, but is it accurate enough to predict single mutation effects? Here, we report that the localized structural deformation between protein pairs differing by only 1-3 mutations-as measured by the effective strain-is correlated across 3901 experimental and AF-predicted structures. Furthermore, analysis of ∼11 000 proteins shows that the local structural change correlates with various phenotypic changes. These findings suggest that AF can predict the range and magnitude of single-mutation effects on average, and we propose a method to improve precision of AF predictions and to indicate when predictions are unreliable.

Patterns in Protein Flexibility: A Comparison of NMR "Ensembles", MD Trajectories, and Crystallographic B-Factors

Proteins are molecular machines requiring flexibility to function. Crystallographic B-factors and Molecular Dynamics (MD) simulations both provide insights into protein flexibility on an atomic scale. Nuclear Magnetic Resonance (NMR) lacks a universally accepted analog of the B-factor. However, a lack of convergence in atomic coordinates in an NMR-based structure calculation also suggests atomic mobility. This paper describes a pattern in the coordinate uncertainties of backbone heavy atoms in NMR-derived structural “ensembles” first noted in the development of FindCore2 (previously called Expanded FindCore: DA Snyder, J Grullon, YJ Huang, R Tejero, GT Montelione, Proteins: Structure, Function, and Bioinformatics 82 (S2), 219–230) and demonstrates that this pattern exists in coordinate variances across MD trajectories but not in crystallographic B-factors. This either suggests that MD trajectories and NMR “ensembles” capture motional behavior of peptide bond units not captured by B-factors or indicates a deficiency common to force fields used in both NMR and MD calculations.

High-performance transformation of protein structure representation from internal to Cartesian coordinates

We present a highly parallel algorithm to convert internal coordinates of a polymeric molecule into Cartesian coordinates. Traditionally, converting the structures of polymers (e.g., proteins) from internal to Cartesian coordinates has been performed serially, due to an inherent linear dependency along the polymer chain. We show this dependency can be removed using a tree-based concatenation of coordinate transforms between segments, and then parallelized efficiently on graphics processing units (GPUs). The conversion algorithm is applicable to protein engineering and fitting protein structures to experimental data, and we observe an order of magnitude speedup using parallel processing on a GPU compared to serial execution on a CPU.

Structural and dynamical description of the enzymatic reaction of a phosphohexomutase

Enzymes are known to adopt various conformations at different points along their catalytic cycles. Here, we present a comprehensive analysis of 15 isomorphous, high resolution crystal structures of the enzyme phosphoglucomutase from the bacterium Xanthomonas citri. The protein was captured in distinct states critical to function, including enzyme-substrate, enzyme-product, and enzyme-intermediate complexes. Key residues in ligand recognition and regions undergoing conformational change are identified and correlated with the various steps of the catalytic reaction. In addition, we use principal component analysis to examine various subsets of these structures with two goals: (1) identifying sites of conformational heterogeneity through a comparison of room temperature and cryogenic structures of the apo-enzyme and (2) a priori clustering of the enzyme-ligand complexes into functionally related groups, showing sensitivity of this method to structural features difficult to detect by traditional methods. This study captures, in a single system, the structural basis of diverse substrate recognition, the subtle impact of covalent modification, and the role of ligand-induced conformational change in this representative enzyme of the α-D-phosphohexomutase superfamily.

Structural origins of misfolding propensity in the platelet adhesive von Willebrand factor A1 domain

The von Willebrand factor (VWF) A1 and A3 domains are structurally isomorphic yet exhibit distinct mechanisms of unfolding. The A1 domain, responsible for platelet adhesion to VWF in hemostasis, unfolds through a molten globule intermediate in an apparent three-state mechanism, while A3 unfolds by a classical two-state mechanism. Inspection of the sequences or structures alone does not elucidate the source of this thermodynamic conundrum; however, the three-state character of the A1 domain suggests that it has more than one cooperative substructure yielding two separate unfolding transitions not present in A3. We investigate the extent to which structural elements contributing to intermediate conformations can be identified using a residue-specific implementation of the structure-energy-equivalence-of-domains algorithm (SEED), which parses proteins of known structure into their constituent thermodynamically cooperative components using protein-group-specific, transfer free energies. The structural elements computed to contribute to the non-two-state character coincide with regions where Von Willebrand disease mutations induce misfolded molten globule conformations of the A1 domain. This suggests a mechanism for the regulation of rheological platelet adhesion to A1 based on cooperative flexibility of the α2 and α3 helices flanking the platelet GPIbα receptor binding interface.

Clusters of Structurally Similar MHC I HLA-A2 Molecules, Found with a New Method, Suggest Mechanisms of T-Cell Receptor Avidity

Only α1 and α2 domains of the α-chain of the major histocompatibility complex (MHC) directly bind peptide antigens (Ag-s) and the T-cell receptor (TCR). Significant plasticity was found in the TCR but only minor in (α1 + α2). The α3-domain position variation was noted only in connection to its binding the coreceptor CD8. We apply our methods for identifying functional conformational changes in proteins to a systematic study of similarities between 43 X-ray structures of the entire α chains of MHC-I HLA-A2. Out of 903 different αHLA-A2 pairs 203 show similarities within the earlier determined uncertainty threshold and unexpectedly form three similarity clusters (SCs) with all/most structures in a cluster similar within the uncertainty threshold. Pairs from different SCs always differ above the threshold, mainly due to variations in the α3 position/structure. All structures in SC3 cannot bind the CD8 coreceptor. Strong hydrogen bonds between (α1 + α2) and α3 differ between SC1 and SC2 but are nearly invariant within each SC. Small conformational changes in the (α1 + α2), caused by Ag-s differences, act as an α3 "allosteric switch" between SC2 and SC1. Binding of CD8 to SC2-HLA-A2 (Tax-type Ag-s) changes it to SC1-HLA-A2 (HuD-type Ag-s). HuD binding to HLA-A2 is much less stable than Tax binding. CD8-liganded HLA-A2 preference for binding HuD suggests that CD8-HLA-A2 may present a weakly binding peptide for TCR recognition, supporting the hypothesis that CD8 increases TCR avidity to weak Ag-s. Other HLA-A2 functions may involve α3. TCR-A6-liganded-Tax-type-HLA-A2s form two small clusters, similar to either A6-liganded-HuD or A6-liganded-native-Tax HLA-A2s.

Parsimony in Protein Conformational Change

Protein conformational change is analyzed by finding the minimalist backbone torsion angle rotations that superpose crystal structures within experimental error. Of several approaches for enforcing parsimony during flexible least-squares superposition, an ℓ(1)-norm restraint provided greatest consistency with independent indications of flexibility from nuclear magnetic resonance relaxation dispersion and chemical shift perturbation in arginine kinase and four previously studied systems. Crystallographic cross-validation shows that the dihedral parameterization describes conformational change more accurately than rigid-group approaches. The rotations that superpose the principal elements of structure constitute a small fraction of the raw (φ, ψ) differences that also reflect local conformation and experimental error. Substantial long-range displacements can be mediated by modest dihedral rotations, accommodated even within α helices and β sheets without disruption of hydrogen bonding at the hinges. Consistency between ligand-associated and intrinsic motions (in the unliganded state) implies that induced changes tend to follow low-barrier paths between conformational sub-states that are in intrinsic dynamic equilibrium.

Factors correlating with significant differences between X-ray structures of myoglobin

Validation of general ideas about the origins of conformational differences in proteins is critical in order to arrive at meaningful functional insights. Here, principal component analysis (PCA) and distance difference matrices are used to validate some such ideas about the conformational differences between 291 myoglobin structures from sperm whale, horse and pig. Almost all of the horse and pig structures form compact PCA clusters with only minor coordinate differences and outliers that are easily explained. The 222 whale structures form a few dense clusters with multiple outliers. A few whale outliers with a prominent distortion of the GH loop are very similar to the cluster of horse structures, which all have a similar GH-loop distortion apparently owing to intermolecular crystal lattice hydrogen bonds to the GH loop from residues near the distal histidine His64. The variations of the GH-loop coordinates in the whale structures are likely to be owing to the observed alternative intermolecular crystal lattice bond, with the change to the GH loop distorting bonds correlated with the binding of specific `unusual' ligands. Such an alternative intermolecular bond is not observed in horse myoglobins, obliterating any correlation with the ligands. Intermolecular bonds do not usually cause significant coordinate differences and cannot be validated as their universal cause. Most of the native-like whale myoglobin structure outliers can be correlated with a few specific factors. However, these factors do not always lead to coordinate differences beyond the previously determined uncertainty thresholds. The binding of unusual ligands by myoglobin, leading to crystal-induced distortions, suggests that some of the conformational differences between the apo and holo structures might not be `functionally important' but rather artifacts caused by the binding of `unusual' substrate analogs. The causes of P6 symmetry in myoglobin crystals and the relationship between crystal and solution structures are also discussed.

Comparison of tertiary structures of proteins in protein-protein complexes with unbound forms suggests prevalence of allostery in signalling proteins

BACKGROUND

Most signalling and regulatory proteins participate in transient protein-protein interactions during biological processes. They usually serve as key regulators of various cellular processes and are often stable in both protein-bound and unbound forms. Availability of high-resolution structures of their unbound and bound forms provides an opportunity to understand the molecular mechanisms involved. In this work, we have addressed the question "What is the nature, extent, location and functional significance of structural changes which are associated with formation of protein-protein complexes?"

RESULTS

A database of 76 non-redundant sets of high resolution 3-D structures of protein-protein complexes, representing diverse functions, and corresponding unbound forms, has been used in this analysis. Structural changes associated with protein-protein complexation have been investigated using structural measures and Protein Blocks description. Our study highlights that significant structural rearrangement occurs on binding at the interface as well as at regions away from the interface to form a highly specific, stable and functional complex. Notably, predominantly unaltered interfaces interact mainly with interfaces undergoing substantial structural alterations, revealing the presence of at least one structural regulatory component in every complex.Interestingly, about one-half of the number of complexes, comprising largely of signalling proteins, show substantial localized structural change at surfaces away from the interface. Normal mode analysis and available information on functions on some of these complexes suggests that many of these changes are allosteric. This change is largely manifest in the proteins whose interfaces are altered upon binding, implicating structural change as the possible trigger of allosteric effect. Although large-scale studies of allostery induced by small-molecule effectors are available in literature, this is, to our knowledge, the first study indicating the prevalence of allostery induced by protein effectors.

CONCLUSIONS

The enrichment of allosteric sites in signalling proteins, whose mutations commonly lead to diseases such as cancer, provides support for the usage of allosteric modulators in combating these diseases.

Molecular determinants for antibody binding on group 1 house dust mite allergens

House dust mites produce potent allergens, Der p 1 and Der f 1, that cause allergic sensitization and asthma. Der p 1 and Der f 1 are cysteine proteases that elicit IgE responses in 80% of mite-allergic subjects and have proinflammatory properties. Their antigenic structure is unknown. Here, we present crystal structures of natural Der p 1 and Der f 1 in complex with a monoclonal antibody, 4C1, which binds to a unique cross-reactive epitope on both allergens associated with IgE recognition. The 4C1 epitope is formed by almost identical amino acid sequences and contact residues. Mutations of the contact residues abrogate mAb 4C1 binding and reduce IgE antibody binding. These surface-exposed residues are molecular targets that can be exploited for development of recombinant allergen vaccines.

Percentile-based spread: a more accurate way to compare crystallographic models

The comparison of biomacromolecular crystal structures is traditionally based on the root-mean-square distance between corresponding atoms. This measure is sensitive to the presence of outliers, which inflate it disproportionately to their fraction. An alternative measure, the percentile-based spread (p.b.s.), is proposed and is shown to represent the average variation in atomic positions more adequately. It is discussed in the context of isomorphous crystal structures, conformational changes and model ensembles generated by repetitive automated rebuilding.

Diversity of function-related conformational changes in proteins: coordinate uncertainty, fragment rigidity, and stability

It was found that the variety of function-related conformational changes ("movements") in proteins is beyond the earlier simple classifications. Here we offer biochemists a more comprehensive, transparent, and easy-to-use approach allowing a detailed and accurate interpretation of such conformational changes. It makes possible a more multifaceted characterization of protein flexibility via identification of rigidly and nonrigidly repositioned fragments, stable and nonstable fragments, and domain and nondomain repositioning. "Coordinate uncertainty thresholds" derived from computed differences between independently determined coordinates of the same molecules are used as the criteria for conformational identity. "Identical" rigid substructures are localized in the distance difference matrices (DDMs). A sequence of simple transformations determines whether a structural change occurs by rigid-body movements of fragments or largely through non-rigid-body deformations. We estimate the stability of protein fragments and compare stable and rigidly moving fragments. The motions computed with the coarse-grained elastic networks are also compared to those of their DDM analogues. We study and suggest a classification for 17 structural pairs, differing in their functional states. For five of the 17 proteins, conformational change cannot be accomplished by rigid-body transformations and requires significant non-rigid-body deformations. Stable fragments rarely coincide with rigidly moving fragments and often disagree with the CATH identifications of domains. Almost all monomeric apo chains, containing stable fragments and/or domains, indicate instability of the entire molecule, suggesting the importance of fragments and domains motions prior to stabilization by substrate binding or crystallization. Notably, kinases exhibit the greatest extent of nonrigidity among the proteins investigated.