Structural heterogeneity in protein crystals

Extensive conformational heterogeneity is reported in highly refined crystallographic models for the proteins crambin, erabutoxin, myohemerythrin, and lamprey hemoglobin. From 6% to 13% of the amino acid side chains of these four proteins are seen in multiple, discrete conformations. Most common are flexible side chains on the molecular surface, but structural heterogeneity occasionally extends to buried side chains or to the polypeptide backbone. A few instances of sequence heterogeneity are also very clear. Numerous solvent sites are multiplets, and at high resolution, multiple, mutually exclusive solvent networks are observed. The proteins have been studied with X-ray diffraction data extending to spacings of from 0.945 to 2.0 A. The extensive heterogeneity observed here provides detailed, accurate structures for conformational substates of these molecules and sets a lower bound on the number of substates accessible to each protein molecule in solution. Electron density is missing or very weak for only a few side chains in these protein crystals, revealing a strong preference for discrete over continuous conformational perturbations. The results at very high resolution further suggest that even rather small conformational fluctuations produce discrete substates and that unresolved conformers are accommodated in increased atomic thermal parameters.     

Advances in all atom sampling methods for modeling protein-ligand binding affinities

Conformational dynamics plays a fundamental role in the regulation of molecular recognition processes. Conformational heterogeneity and entropy variations upon binding, although not always evident from the analysis of structural data, can substantially affect affinity and specificity. Computer modeling is able to provide some of the most direct insights into these aspects of molecular recognition. We review recent physics-based computational studies that employ advanced conformational sampling algorithms and effective potentials to model the three main classes of degrees of freedom relevant to the binding process: ligand positioning relative to the receptor, ligand and receptor internal reorganization, and hydration. Collectively these studies show that all of these elements are important for proper modeling of protein-ligand interactions.     

Integrating Molecular Models Into CryoEM Heterogeneity Analysis Using Scalable High-resolution Deep Gaussian Mixture Models

Resolving the structural variability of proteins is often key to understanding the structure–function relationship of those macromolecular machines. Single particle analysis using Cryogenic electron microscopy (CryoEM), combined with machine learning algorithms, provides a way to reveal the dynamics within the protein system from noisy micrographs. Here, we introduce an improved computational method that uses Gaussian mixture models for protein structure representation and deep neural networks for conformation space embedding. By integrating information from molecular models into the heterogeneity analysis, we can analyze continuous protein conformational changes using structural information at the frequency of 1/3 Å⁻¹, and present the results in a more interpretable form.     

Conformational control through translocational regulation: a new view of secretory and membrane protein folding

We suggest a new view of secretory and membrane protein folding that emphasizes the role of pathways of biogenesis in generating functional and conformational heterogeneity. In this view, heterogeneity results from action of accessory factors either directly binding specific sequences of the nascent chain, or indirectly, changing the environment in which a particular domain is synthesized. Entrained by signaling pathways, these variables create a combinatorial set of necessary-but-not-sufficient conditions that enhance synthesis and folding of particular alternate, functional, conformational forms. We therefore propose that protein conformation is productively regulated by the cell during translocation across the endoplasmic reticulum (ER), a concept that may account for currently poorly understood aspects of physiological function, natural selection, and disease pathogenesis.     

Effects of temperature and Y21M mutation on conformational heterogeneity of the major coat protein (pVIII) of filamentous bacteriophage fd

Solid-state NMR spectroscopy was used to analyze the conformational heterogeneity of the major coat protein (pVIII) of filamentous bacteriophage fd. Both one and two-dimensional solid-state NMR spectra of magnetically aligned samples of fd bacteriophage reveal that an increase in temperature and a single site substitution (Tyr21 to Met, Y21M) reduce the conformational heterogeneity observed throughout wild-type pVIII. The NMR results are consistent with previous studies indicating that conformational flexibility in the hinge-bend segment that links the amphipathic and hydrophobic helices in the membrane-bound form of the protein plays an essential role during phage assembly, which involves a major change in the tertiary, but not secondary, structure of the coat protein.     

Hemocyanin from E. californicum encapsulated in silica gels: oxygen binding and conformational states

Cooperativity depends on the existence of equilibria among functionally distinct conformational states that are affected by homo and heterotropic effectors. In order to isolate the quaternary conformations of hemocyanin from E. californicum, the 24-meric giant protein was encapsulated in wet, nanoporous silica gels, either in the absence or presence of oxygen. The deoxy- and oxy-hemocyanin gels exhibit a p50 for oxygen of 11 and 2.5 torr, respectively, values in close agreement with those for hemocyanin in solution. The observed Hill coefficients are lower than unity, indicating a conformational heterogeneity within each locked conformational state, a finding in agreement with the assumption that at least four conformational states are required to explain the oxygen binding properties of E. californicum hemocyanin in solution.     

Effects of pH and temperature on cardioactive polypeptides from sea anemones: a 1H-NMR study

The effects of pH and temperature on the 300-MHz ¹H-nmr spectra of three cardioactive polypeptides from sea anemones, anthopleurin-A from Anthopleura xanthogrammica (AP-A) and Anemonia sulcata toxins I and II (ATX I and II), are described. AP-A and ATX II exhibit major spectral heterogeneity. Evidence from the pH and temperature studies and from high performance liquid chromatography indicates that this heterogeneity is conformational rather than chemical in origin. By contrast, purified isotoxins of ATX I show no evidence of conformational heterogeneity. The pKa values of most of the ionizable groups in these polypeptides are not strongly perturbed by interactions in the tertiary structure, with the exception of one of the Asp carboxylates, which has a pKa of ? 2 in AP-A and ATX II and 3.0 in ATX I. Protonation of this carboxylate, suggested to be Asp-9, leads to a conformational change in all three molecules. All three polypeptides are thermally stable, showing some conformational changes but not major unfolding at elevated temperatures.     

Probing the conformation of the prion protein within a single amyloid fibril using a novel immunoconformational assay

The coexistence of multiple strains or subtypes of the disease-related isoform of prion protein (PrP) in natural isolates, together with the observed conformational heterogeneity of PrP amyloid fibrils generated in vitro, indicates the importance of probing the conformation of single particles within heterogeneous samples. Using an array of PrP-specific antibodies, we report the development of a novel immunoconformational assay. Uniquely, application of this new technology allows the conformation of multimeric PrP within a single fibril or particle to be probed without pretreatment of the sample with proteinase K. Using amyloid fibrils prepared from full-length recombinant PrP, we demonstrated the utility of this assay to define (i) PrP regions that are surface-exposed or buried, (ii) the susceptibility of defined PrP regions to GdnHCl-induced denaturation, and (iii) the conformational heterogeneity of PrP fibrils as measured for either the entire fibrillar population or for individual fibrils. Specifically, PrP regions 159-174 and 224-230 were shown to be buried and were the most resistant to denaturation. The 132-156 segment of PrP was found to be cryptic under native conditions and solvent-exposed under partially denaturing conditions, whereas the region 95-105 was solvent-accessible regardless of the solvent conditions. Remarkably, a subfraction of fibrils showed immunoreactivity to PrPSc-specific antibodies designated as IgGs 89-112 and 136-158. The immunoreactivity of the conformational epitopes was reduced upon exposure to partially denaturing conditions. Unexpectedly, PrPSc -specific antibodies revealed conformational polymorphisms even within individual fibrils. Our studies provide valuable new insight into fibrillar substructure and offer a new tool for probing the conformation of single PrP fibrils.     

Structure and dynamics of an RNA tetraloop: a joint molecular dynamics and NMR study

Molecular dynamics simulations of the RNA tetraloop 5'-CGCUUUUGCG-3' with high melting temperature and significant conformational heterogeneity in explicit water solvent are presented and compared to NMR studies. The NMR data allow for a detailed test of the theoretical model, including the quality of the force field and the conformational sampling. Due to the conformational heterogeneity of the tetraloop, high temperature (350 K) and locally enhanced sampling simulations need to be invoked. The Amber98 force field leads to a good overall agreement with experimental data. Based on NMR data and a principal component analysis of the 350 K trajectory, the dynamic structure of the tetraloop is revealed. The principal component free energy surface exhibits four minima, which correspond to well-defined conformational structures that differ mainly by their base stacking in the loop region. No correlation between the motion of the sugar rings and the stacking dynamics of the loop bases is found.     

Modeling conformational states of proteins with AlphaFold

Many proteins exert their function by switching among different structures. Knowing the conformational ensembles affiliated with these states is critical to elucidate key mechanistic aspects that govern protein function. While experimental determination efforts are still bottlenecked by cost, time, and technical challenges, the machine-learning technology AlphaFold showed near experimental accuracy in predicting the three-dimensional structure of monomeric proteins. However, an AlphaFold ensemble of models usually represents a single conformational state with minimal structural heterogeneity. Consequently, several pipelines have been proposed to either expand the structural breadth of an ensemble or bias the prediction toward a desired conformational state. Here, we analyze how those pipelines work, what they can and cannot predict, and future directions.     

1D2DSimScore: A novel method for comparing contacts in biomacromolecules and their complexes

The biologically relevant structures of proteins and nucleic acids and their complexes are dynamic. They include a combination of regions ranging from rigid structural segments to structural switches to regions that are almost always disordered, which interact with each other in various ways. Comparing conformational changes and variation in contacts between different conformational states is essential to understand the biological functions of proteins, nucleic acids, and their complexes. Here, we describe a new computational tool, 1D2DSimScore, for comparing contacts and contact interfaces in all kinds of macromolecules and macromolecular complexes, including proteins, nucleic acids, and other molecules. 1D2DSimScore can be used to compare structural features of macromolecular models between alternative structures obtained in a particular experiment or to score various predictions against a defined “ideal” reference structure. Comparisons at the level of contacts are particularly useful for flexible molecules, for which comparisons in 3D that require rigid-body superpositions are difficult, and in biological systems where the formation of specific inter-residue contacts is more relevant for the biological function than the maintenance of a specific global 3D structure. Similarity/dissimilarity scores calculated by 1D2DSimScore can be used to complement scores describing 3D structural similarity measures calculated by the existing tools.     

Fuzziness in Protein Interactions—A Historical Perspective

The proposal that coupled folding to binding is not an obligatory mechanism for intrinsically disordered (ID) proteins was put forward 10 years ago. The notion of fuzziness implies that conformational heterogeneity can be maintained upon interactions of ID proteins, which has a functional impact either on regulated assembly or activity of the corresponding complexes. Here I review how the concept has evolved in the past decade, via increasing experimental data providing insights into the mechanisms, pathways and regulatory modes. The effects of structural diversity and transient contacts on protein assemblies have been collected and systematically analyzed (Fuzzy Complexes Database, http://protdyn-database.org). Fuzziness has also been exploited as a framework to decipher molecular organization of higher-order protein structures. Quantification of conformational heterogeneity opens exciting future perspectives for drug discovery from small molecule–ID protein interactions to supramolecular assemblies.     

Antigenically distinct conformations of CXCR4

The major human immunodeficiency virus type 1 (HIV-1) coreceptors are the chemokine receptors CCR5 and CXCR4. The patterns of expression of the major coreceptors and their use by HIV-1 strains largely explain viral tropism at the level of entry. However, while virus infection is dependent upon the presence of CD4 and an appropriate coreceptor, it can be influenced by a number of factors, including receptor concentration, affinity between envelope gp120 and receptors, and potentially receptor conformation. Indeed, seven-transmembrane domain receptors, such as CCR5, can exhibit conformational heterogeneity, although the significance for virus infection is uncertain. Using a panel of monoclonal antibodies (MAbs) to CXCR4, we found that CXCR4 on both primary and transformed T cells as well as on primary B cells exhibited considerable conformational heterogeneity. The conformational heterogeneity of CXCR4 explains the cell-type-dependent ability of CXCR4 antibodies to block chemotaxis to stromal cell-derived factor 1 alpha and to inhibit HIV-1 infection. In addition, the MAb most commonly used to study CXCR4 expression, 12G5, recognizes only a subpopulation of CXCR4 molecules on all primary cell types analyzed. As a result, CXCR4 concentrations on these important cell types have been underestimated to date. Finally, while the factors responsible for altering CXCR4 conformation are not known, we found that they do not involve CXCR4 glycosylation, sulfation of the N-terminal domain of CXCR4, or pertussis toxin-sensitive G-protein coupling. The fact that this important HIV-1 coreceptor exists in multiple conformations could have implications for viral entry and for the development of receptor antagonists.     

Heterogeneity of Protein Substates Visualized by Spin-label EPR

The energy landscape of proteins is characterized by a hierarchy of substates, which give rise to conformational heterogeneity at low temperatures. In multiply spin-labeled membranous Na,K-ATPase, this heterogeneous population of conformations is manifest by strong inhomogeneous broadening of the electron paramagnetic resonance (EPR) line shapes and nonexponential spin-echo decays, which undergo a transition to homogeneous broadening and exponential relaxation at higher temperatures (previous study). In this study, we apply these EPR methods to small water-soluble proteins, of the type for which the existence of conformational substates is well established. Both α-helical and β-sheet aqueous proteins that are spin-labeled on a single cysteine residue display spin-echo decays with a single phase-memory time T_2M and conventional EPR line shapes with predominantly homogeneous broadening, over a broad range of temperatures from 77 K to ∼ 250 K or higher. Above ∼ 200 K, the residual inhomogeneous broadening is reduced almost to zero. In contrast, both the proteins and the spin label alone, when in a glycerol-water mixture below the glass transition, display heterogeneity in spin-echo phase-memory time and a stronger inhomogeneous broadening of the conventional line shapes, similar to multiply spin-labeled membranous Na,K-ATPase below 200 K. Above 200 K (or the glass transition), a single phase-memory time and predominantly homogeneous broadening are found in both spin-label systems. The results are discussed in terms of solvent-mediated protein transitions, the ability of single spin-label sites to detect conformational heterogeneity, and the desirability of exploring multiple sites for proteins with the size and complexity of the Na,K-ATPase.     

Heterogeneity of Protein Substates Visualized by Spin-label EPR

The energy landscape of proteins is characterized by a hierarchy of substates, which give rise to conformational heterogeneity at low temperatures. In multiply spin-labeled membranous Na,K-ATPase, this heterogeneous population of conformations is manifest by strong inhomogeneous broadening of the electron paramagnetic resonance (EPR) line shapes and nonexponential spin-echo decays, which undergo a transition to homogeneous broadening and exponential relaxation at higher temperatures (previous study). In this study, we apply these EPR methods to small water-soluble proteins, of the type for which the existence of conformational substates is well established. Both α-helical and β-sheet aqueous proteins that are spin-labeled on a single cysteine residue display spin-echo decays with a single phase-memory time T_2M and conventional EPR line shapes with predominantly homogeneous broadening, over a broad range of temperatures from 77 K to ∼ 250 K or higher. Above ∼ 200 K, the residual inhomogeneous broadening is reduced almost to zero. In contrast, both the proteins and the spin label alone, when in a glycerol-water mixture below the glass transition, display heterogeneity in spin-echo phase-memory time and a stronger inhomogeneous broadening of the conventional line shapes, similar to multiply spin-labeled membranous Na,K-ATPase below 200 K. Above 200 K (or the glass transition), a single phase-memory time and predominantly homogeneous broadening are found in both spin-label systems. The results are discussed in terms of solvent-mediated protein transitions, the ability of single spin-label sites to detect conformational heterogeneity, and the desirability of exploring multiple sites for proteins with the size and complexity of the Na,K-ATPase.     

Conformational transitions of adenylate kinase: switching by cracking

Conformational heterogeneity in proteins is known to often be the key to their function. We present a coarse grained model to explore the interplay between protein structure, folding and function which is applicable to allosteric or non-allosteric proteins. We employ the model to study the detailed mechanism of the reversible conformational transition of Adenylate Kinase (AKE) between the open to the closed conformation, a reaction that is crucial to the protein's catalytic function. We directly observe high strain energy which appears to be correlated with localized unfolding during the functional transition. This work also demonstrates that competing native interactions from the open and closed form can account for the large conformational transitions in AKE. We further characterize the conformational transitions with a new measure Phi(Func), and demonstrate that local unfolding may be due, in part, to competing intra-protein interactions.     

Persistent conformational heterogeneity of triosephosphate isomerase: separation and characterization of conformational isomers in solution

Subunit dissociation of dimeric rabbit muscle triosephosphate isomerase (TIM) by hydrostatic pressure has previously been shown not to follow the expected dependence on protein concentration [Rietveld and Ferreira (1996) Biochemistry 35, 7743-7751]. This anomalous behavior was attributed to persistent conformational heterogeneity (i.e., the coexistence of long-lived conformational isomers) in the ensemble of TIM dimers. Here, we initially show that subunit dissociation/unfolding of TIM by guanidine hydrochloride (GdnHCl) also exhibits an anomalous dependence on protein concentration. Dissociation/unfolding of TIM by GdnHCl was investigated by intrinsic fluorescence and circular dichroism spectroscopies and was found to be a highly cooperative transition in which the tertiary and secondary structures of the protein were concomitantly lost. A procedure based on size-exclusion chromatography in the presence of intermediate (0.6 M) GdnHCl concentrations was developed to isolate two conformational isomers of TIM that exhibit significantly different stabilities and kinetics of unfolding by GdnHCl. Complete unfolding of the two isolated conformers at a high GdnHCl concentration (1.5 M), followed by refolding by removal of the denaturant, completely abolished the differences in their unfolding kinetics. These results indicate that such differences stem from conformational heterogeneity of TIM and are not related to any chemical modification of the protein. Furthermore, they add support to the notion that long-lived conformational isomers of TIM coexist in solution and provide a basis for the interpretation of the persistent heterogeneity of this protein.