A newly introduced salt bridge cluster improves structural and biophysical properties of de novo TIM barrels

Protein stability can be fine‐tuned by modifying different structural features such as hydrogen‐bond networks, salt bridges, hydrophobic cores, or disulfide bridges. Among these, stabilization by salt bridges is a major challenge in protein design and engineering since their stabilizing effects show a high dependence on the structural environment in the protein, and therefore are difficult to predict and model. In this work, we explore the effects on structure and stability of an introduced salt bridge cluster in the context of three different de novo TIM barrels. The salt bridge variants exhibit similar thermostability in comparison with their parental designs but important differences in the conformational stability at 25°C can be observed such as a highly stabilizing effect for two of the proteins but a destabilizing effect to the third. Analysis of the formed geometries of the salt bridge cluster in the crystal structures show either highly ordered salt bridge clusters or only single salt bridges. Rosetta modeling of the salt bridge clusters results in a good prediction of the tendency on stability changes but not the geometries observed in the three‐dimensional structures. The results show that despite the similarities in protein fold, the salt bridge clusters differently influence the structural and stability properties of the de novo TIM barrel variants depending on the structural background where they are introduced.     

Frontotemporal dementia‐linked P112H mutation of TDP‐43 induces protein structural change and impairs its RNA binding function

TDP‐43 forms the primary constituents of the cytoplasmic inclusions contributing to various neurodegenerative diseases, including amyotrophic lateral sclerosis and frontotemporal dementia (FTD). Over 60 TDP‐43 mutations have been identified in patients suffering from these two diseases, but most variations are located in the protein''s disordered C‐terminal glycine‐rich region. P112H mutation of TDP‐43 has been uniquely linked to FTD, and is located in the first RNA recognition motif (RRM1). This mutation is thought to be pathogenic, but its impact on TDP‐43 at the protein level remains unclear. Here, we compare the biochemical and biophysical properties of TDP‐43 truncated proteins with or without P112H mutation. We show that P112H‐mutated TDP‐43 proteins exhibit higher thermal stability, impaired RNA‐binding activity, and a reduced tendency to aggregate relative to wild‐type proteins. Near‐UV CD, 2D‐nuclear‐magnetic resonance, and intrinsic fluorescence spectrometry further reveal that the P112H mutation in RRM1 generates local conformational changes surrounding the mutational site that disrupt the stacking interactions of the W113 side chain with nucleic acids. Together, these results support the notion that P112H mutation of TDP‐43 contributes to FTD through functional impairment of RNA metabolism and/or structural changes that curtail protein clearance.     

High‐resolution analysis of the conformational transition of pro‐apoptotic Bak at the lipid membrane

Permeabilization of the outer mitochondrial membrane by pore‐forming Bcl2 proteins is a crucial step for the induction of apoptosis. Despite a large set of data suggesting global conformational changes within pro‐apoptotic Bak during pore formation, high‐resolution structural details in a membrane environment remain sparse. Here, we used NMR and HDX‐MS (Hydrogen deuterium exchange mass spectrometry) in lipid nanodiscs to gain important high‐resolution structural insights into the conformational changes of Bak at the membrane that are dependent on a direct activation by BH3‐only proteins. Furthermore, we determined the first high‐resolution structure of the Bak transmembrane helix. Upon activation, α‐helix 1 in the soluble domain of Bak dissociates from the protein and adopts an unfolded and dynamic potentially membrane‐bound state. In line with this finding, comparative protein folding experiments with Bak and anti‐apoptotic BclxL suggest that α‐helix 1 in Bak is a metastable structural element contributing to its pro‐apoptotic features. Consequently, mutagenesis experiments aimed at stabilizing α‐helix 1 yielded Bak variants with delayed pore‐forming activity. These insights will contribute to a better mechanistic understanding of Bak‐mediated membrane permeabilization.     

Structure‐based analysis and evolution of a monomerized red‐colored chromoprotein from the Olindias formosa jellyfish

GFP‐like chromoproteins (CPs) with non‐fluorescence ability have been used as bioimaging probes. Existing CPs have voids in the optical absorption window which limits their extensibility. The development of new CP color is therefore ongoing. Here, we cloned CPs from the jellyfish, Olindias formosa, and developed a completely non‐fluorescent monomeric red CP, R‐Velour, with an absorption peak at 528 nm. To analyze the photophysical properties from a structural aspect, we determined the crystal structure of R‐Velour at a 2.1 Å resolution. R‐Velour has a trans‐chromophore similar to the green fluorescence protein, Gamillus, derived from the same jellyfish. However, in contrast to the two coplanar chromophoric rings in Gamillus, R‐Velour has a large torsion inducing non‐fluorescence property. Through site‐directed mutagenesis, we surveyed residues surrounding the chromophore and found a key residue, Ser155, which contributes to the generation of four‐color variants with the bathochromic and hypsochromic shift of the absorption peak, ranging from 506 to 554 nm. The recently proposed spectrum shift theory, based on the Marcus–Hush model, supports the spectrum shift of these mutants. These findings may support further development of R‐Velour variants with useful absorption characteristics for bioimaging, including fluorescence lifetime imaging and photoacoustic imaging.     

Complementing machine learning‐based structure predictions with native mass spectrometry

The advent of machine learning‐based structure prediction algorithms such as AlphaFold2 (AF2) and RoseTTa Fold have moved the generation of accurate structural models for the entire cellular protein machinery into the reach of the scientific community. However, structure predictions of protein complexes are based on user‐provided input and may require experimental validation. Mass spectrometry (MS) is a versatile, time‐effective tool that provides information on post‐translational modifications, ligand interactions, conformational changes, and higher‐order oligomerization. Using three protein systems, we show that native MS experiments can uncover structural features of ligand interactions, homology models, and point mutations that are undetectable by AF2 alone. We conclude that machine learning can be complemented with MS to yield more accurate structural models on a small and large scale.     

Enzyme catalysis prior to aromatic residues: Reverse engineering of a dephospho‐CoA kinase

The wide variety of protein structures and functions results from the diverse properties of the 20 canonical amino acids. The generally accepted hypothesis is that early protein evolution was associated with enrichment of a primordial alphabet, thereby enabling increased protein catalytic efficiencies and functional diversification. Aromatic amino acids were likely among the last additions to genetic code. The main objective of this study was to test whether enzyme catalysis can occur without the aromatic residues (aromatics) by studying the structure and function of dephospho‐CoA kinase (DPCK) following aromatic residue depletion. We designed two variants of a putative DPCK from Aquifex aeolicus by substituting (a) Tyr, Phe and Trp or (b) all aromatics (including His). Their structural characterization indicates that substituting the aromatics does not markedly alter their secondary structures but does significantly loosen their side chain packing and increase their sizes. Both variants still possess ATPase activity, although with 150–300 times lower efficiency in comparison with the wild‐type phosphotransferase activity. The transfer of the phosphate group to the dephospho‐CoA substrate becomes heavily uncoupled and only the His‐containing variant is still able to perform the phosphotransferase reaction. These data support the hypothesis that proteins in the early stages of life could support catalytic activities, albeit with low efficiencies. An observed significant contraction upon ligand binding is likely important for appropriate organization of the active site. Formation of firm hydrophobic cores, which enable the assembly of stably structured active sites, is suggested to provide a selective advantage for adding the aromatic residues.     

Cooperative hydrophobic core interactions in the β‐trefoil architecture

Symmetric protein architectures have a compelling aesthetic that suggests a plausible evolutionary process (i.e., gene duplication/fusion) yielding complex architecture from a simpler structural motif. Furthermore, symmetry inspires a practical approach to computational protein design that substantially reduces the combinatorial explosion problem, and may provide practical solutions for structure optimization. Despite such broad relevance, the role of structural symmetry in the key area of hydrophobic core‐packing cooperativity has not been adequately studied. In the present report, the threefold rotational symmetry intrinsic to the β‐trefoil architecture is shown to form a geometric basis for highly‐cooperative core‐packing interactions that both stabilize the local repeating motif and promote oligomerization/long‐range contacts in the folding process. Symmetry in the β‐trefoil structure also permits tolerance towards mutational drift that involves a structural quasi‐equivalence at several key core positions.     

Evidence for a synergistic effect of post‐translational modifications and genomic composition of eEF‐1α on the adaptation of Phytophthora infestans

Genetic variation plays a fundamental role in pathogen''s adaptation to environmental stresses. Pathogens with low genetic variation tend to survive and proliferate more poorly due to their lack of genotypic/phenotypic polymorphisms in responding to fluctuating environments. Evolutionary theory hypothesizes that the adaptive disadvantage of genes with low genomic variation can be compensated for structural diversity of proteins through post‐translation modification (PTM) but this theory is rarely tested experimentally and its implication to sustainable disease management is hardly discussed. In this study, we analyzed nucleotide characteristics of eukaryotic translation elongation factor‐1α (eEF‐lα) gene from 165 Phytophthora infestans isolates and the physical and chemical properties of its derived proteins. We found a low sequence variation of eEF‐lα protein, possibly attributable to purifying selection and a lack of intra‐genic recombination rather than reduced mutation. In the only two isoforms detected by the study, the major one accounted for >95% of the pathogen collection and displayed a significantly higher fitness than the minor one. High lysine representation enhances the opportunity of the eEF‐1α protein to be methylated and the absence of disulfide bonds is consistent with the structural prediction showing that many disordered regions are existed in the protein. Methylation, structural disordering, and possibly other PTMs ensure the ability of the protein to modify its functions during biological, cellular and biochemical processes, and compensate for its adaptive disadvantage caused by sequence conservation. Our results indicate that PTMs may function synergistically with nucleotide codes to regulate the adaptive landscape of eEF‐1α, possibly as well as other housekeeping genes, in P. infestans. Compensatory evolution between pre‐ and post‐translational phase in eEF‐1α could enable pathogens quickly adapting to disease management strategies while efficiently maintaining critical roles of the protein playing in biological, cellular, and biochemical activities. Implications of these results to sustainable plant disease management are discussed.     

Cyclic nucleotide selectivity of protein kinase G isozymes

The intrinsic activity of the C‐terminal catalytic (C) domain of cyclic guanosine monophosphate (cGMP)‐dependent protein kinases (PKG) is inhibited by interactions with the N‐terminal regulatory (R) domain. Selective binding of cGMP to cyclic nucleotide binding (CNB) domains within the R‐domain disrupts the inhibitory R–C interaction, leading to the release and activation of the C‐domain. Affinity measurements of mammalian and plasmodium PKG CNB domains reveal different degrees of cyclic nucleotide affinity and selectivity; the CNB domains adjacent to the C‐domain are more cGMP selective and therefore critical for cGMP‐dependent activation. Crystal structures of isolated CNB domains in the presence and absence of cyclic nucleotides reveal isozyme‐specific contacts that explain cyclic nucleotide selectivity and conformational changes that accompany CNB. Crystal structures of tandem CNB domains identify two types of CNB‐mediated dimeric contacts that indicate cGMP‐driven reorganization of domain–domain interfaces that include large conformational changes. Here, we review the available structural and functional information of PKG CNB domains that further advance our understanding of cGMP mediated regulation and activation of PKG isozymes.     

TSGIT: An N‐ and C‐terminal tandem tag system for purification of native and intein‐mediated ligation‐ready proteins

A large variety of fusion tags have been developed to improve protein expression, solubilization, and purification. Nevertheless, these tags have been combined in a rather limited number of composite tags and usually these composite tags have been dictated by traditional commercially‐available expression vectors. Moreover, most commercially‐available expression vectors include either N‐ or C‐terminal fusion tags but not both. Here, we introduce TSGIT, a fusion‐tag system composed of both N‐ and a C‐terminal composite fusion tags. The system includes two affinity tags, two solubilization tags and two cleavable tags distributed at both termini of the protein of interest. Therefore, the N‐ and the C‐terminal composite fusion tags in TSGIT are fully orthogonal in terms of both affinity selection and cleavage. For using TSGIT, we streamlined the cloning, expression, and purification procedures. Each component tag is selected to maximize its benefits toward the final construct. By expressing and partially purifying the protein of interest between the components of the TSGIT fusion, the full‐length protein is selected over truncated forms, which has been a long‐standing problem in protein purification. Moreover, due to the nature of the cleavable tags in TSGIT, the protein of interest is obtained in its native form without any additional undesired N‐ or C‐terminal amino acids. Finally, the resulting purified protein is ready for efficient ligation with other proteins or peptides for downstream applications. We demonstrate the use of this system by purifying a large amount of native fluorescent mRuby3 protein and bacteriophage T7 gp2.5 ssDNA‐binding protein.     

Early‐life conditions impact juvenile telomere length,but do not predict later life‐history strategies or fitness in a wild vertebrate

Environmental conditions experienced during early life may have long‐lasting effects on later‐life phenotypes and fitness. Individuals experiencing poor early‐life conditions may suffer subsequent fitness constraints. Alternatively, individuals may use a strategic “Predictive Adaptive Response” (PAR), whereby they respond—in terms of physiology or life‐history strategy—to the conditions experienced in early life to maximize later‐life fitness. Particularly, the Future Lifespan Expectation (FLE) PAR hypothesis predicts that when poor early‐life conditions negatively impact an individual''s physiological state, it will accelerate its reproductive schedule to maximize fitness during its shorter predicted life span. We aimed to measure the impact of early‐life conditions and resulting fitness across individual lifetimes to test predictions of the FLE hypothesis in a wild, long‐lived model species. Using a long‐term individual‐based dataset, we investigated how early‐life conditions are linked with subsequent fitness in an isolated population of the Seychelles warbler Acrocephalus sechellensis. How individuals experience early‐life environmental conditions may vary greatly, so we also tested whether telomere length—shorter telomers are a biomarker of an individual''s exposure to stress—can provide an effective measure of the individual‐specific impact of early‐life conditions. Specifically, under the FLE hypothesis, we would expect shorter telomeres to be associated with accelerated reproduction. Contrary to expectations, shorter juvenile telomere length was not associated with poor early‐life conditions, but instead with better conditions, probably as a result of faster juvenile growth. Furthermore, neither juvenile telomere length, nor other measures of early‐life conditions, were associated with age of first reproduction or the number of offspring produced during early life in either sex. We found no support for the FLE hypothesis. However, for males, poor early‐life body condition was associated with lower first‐year survival and reduced longevity, indicating that poor early‐life conditions pose subsequent fitness constraints. Our results also showed that using juvenile telomere length as a measure of early‐life conditions requires caution, as it is likely to not only reflect environmental stress but also other processes such as growth.     

Exploitation of the ribosomal protein L10 R98S mutation to enhance recombinant protein production in mammalian cells

Mammalian cells are commonly used to produce recombinant protein therapeutics, but suffer from a high cost per mg of protein produced. There is therefore great interest in improving protein yields to reduce production cost. We present an entirely novel approach to reach this goal through direct engineering of the cellular translation machinery by introducing the R98S point mutation in the catalytically essential ribosomal protein L10 (RPL10‐R98S). Our data support that RPL10‐R98S enhances translation levels and fidelity and reduces proteasomal activity in lymphoid Ba/F3 and Jurkat cell models. In HEK293T cells cultured in chemically defined medium, knock‐in of RPL10‐R98S was associated with a 1.7‐ to 2.5‐fold increased production of four transiently expressed recombinant proteins and 1.7‐fold for one out of two stably expressed proteins. In CHO‐S cells, eGFP reached a 2‐fold increased expression under stable but not transient conditions, but there was no production benefit for monoclonal antibodies. The RPL10‐R98S associated production gain thus depends on culture conditions, cell type, and the nature of the expressed protein. Our study demonstrates the potential for using a ribosomal protein mutation for pharmaceutical protein production gains, and further research on how various factors influence RPL10‐R98S phenotypes can maximize its exploitability for the mammalian protein production industry.     

Reproductive consequences of transient pathogen exposure across host genotypes and generations

To maximize fitness upon pathogenic infection, host organisms might reallocate energy and resources among life‐history traits, such as reproduction and defense. The fitness costs of infection can result from both immune upregulation and direct pathogen exploitation. The extent to which these costs, separately and together, vary by host genotype and across generations is unknown. We attempted to disentangle these costs by transiently exposing wild isolates and a lab‐domesticated strain of Caenorhabditis elegans nematodes to the pathogen Staphylococcus aureus, using exposure to heat‐killed pathogens to distinguish costs due to immune upregulation and pathogen exploitation. We found that host nematodes exhibit a short‐term delay in offspring production when exposed to live and heat‐killed pathogen, but their lifetime fecundity (total offspring produced) recovered to control levels. We also found genetic variation between host isolates for both cumulative offspring production and magnitude of fitness costs. We further investigated whether there were maternal pathogen exposure costs (or benefits) to offspring and revealed a positive correlation between the magnitude of the pathogen‐induced delay in the parent''s first day of reproduction and the cost to offspring population growth. Our findings highlight the capacity for hosts to recover fecundity after transient exposure to a pathogen.     

X‐ray scattering reveals disordered linkers and dynamic interfaces in complexes and mechanisms for DNA double‐strand break repair impacting cell and cancer biology

Evolutionary selection ensures specificity and efficiency in dynamic metastable macromolecular machines that repair DNA damage without releasing toxic and mutagenic intermediates. Here we examine non‐homologous end joining (NHEJ) as the primary conserved DNA double‐strand break (DSB) repair process in human cells. NHEJ has exemplary key roles in networks determining the development, outcome of cancer treatments by DSB‐inducing agents, generation of antibody and T‐cell receptor diversity, and innate immune response for RNA viruses. We determine mechanistic insights into NHEJ structural biochemistry focusing upon advanced small angle X‐ray scattering (SAXS) results combined with X‐ray crystallography (MX) and cryo‐electron microscopy (cryo‐EM). SAXS coupled to atomic structures enables integrated structural biology for objective quantitative assessment of conformational ensembles and assemblies in solution, intra‐molecular distances, structural similarity, functional disorder, conformational switching, and flexibility. Importantly, NHEJ complexes in solution undergo larger allosteric transitions than seen in their cryo‐EM or MX structures. In the long‐range synaptic complex, X‐ray repair cross‐complementing 4 (XRCC4) plus XRCC4‐like‐factor (XLF) form a flexible bridge and linchpin for DNA ends bound to KU heterodimer (Ku70/80) and DNA‐PKcs (DNA‐dependent protein kinase catalytic subunit). Upon binding two DNA ends, auto‐phosphorylation opens DNA‐PKcs dimer licensing NHEJ via concerted conformational transformations of XLF‐XRCC4, XLF–Ku80, and LigIV^BRCT–Ku70 interfaces. Integrated structures reveal multifunctional roles for disordered linkers and modular dynamic interfaces promoting DSB end processing and alignment into the short‐range complex for ligation by LigIV. Integrated findings define dynamic assemblies fundamental to designing separation‐of‐function mutants and allosteric inhibitors targeting conformational transitions in multifunctional complexes.     

The order of PDZ3 and TrpCage in fusion chimeras determines their properties—a biophysical characterization

Most of the structural proteins known today are composed of domains that carry their own functions while keeping their structural properties. It is supposed that such domains, when taken out of the context of the whole protein, can retain their original structure and function to a certain extent. Information on the specific functional and structural characteristics of individual domains in a new context of artificial fusion proteins may help to reveal the rules of internal and external domain communication. Moreover, this could also help explain the mechanism of such communication and address how the mutual allosteric effect plays a role in a such multi‐domain protein system. The simple model system of the two‐domain fusion protein investigated in this work consisted of a well‐folded PDZ3 domain and an artificially designed small protein domain called Tryptophan Cage (TrpCage). Two fusion proteins with swapped domain order were designed to study their structural and functional features as well as their biophysical properties. The proteins composed of PDZ3 and TrpCage, both identical in amino acid sequence but different in composition (PDZ3‐TrpCage, TrpCage‐PDZ3), were studied using circualr dichroism (CD) spectrometry, analytical ultracentrifugation, and molecular dynamic simulations. The biophysical analysis uncovered different structural and denaturation properties of both studied proteins, revealing their different unfolding pathways and dynamics.     

Enhancing protein perdeuteration by experimental evolution of Escherichia coli K‐12 for rapid growth in deuterium‐based media

Deuterium is a natural low abundance stable hydrogen isotope that in high concentrations negatively affects growth of cells. Here, we have studied growth of Escherichia coli MG1655, a wild‐type laboratory strain of E. coli K‐12, in deuterated glycerol minimal medium. The growth rate and final biomass in deuterated medium is substantially reduced compared to cells grown in ordinary medium. By using a multi‐generation adaptive laboratory evolution‐based approach, we have isolated strains that show increased fitness in deuterium‐based growth media. Whole‐genome sequencing identified the genomic changes in the obtained strains and show that there are multiple routes to genetic adaptation to growth in deuterium‐based media. By screening a collection of single‐gene knockouts of nonessential genes, no specific gene was found to be essential for growth in deuterated minimal medium. Deuteration of proteins is of importance for NMR spectroscopy, neutron protein crystallography, neutron reflectometry, and small angle neutron scattering. The laboratory evolved strains, with substantially improved growth rate, were adapted for recombinant protein production by T7 RNA polymerase overexpression systems and shown to be suitable for efficient production of perdeuterated soluble and membrane proteins for structural biology applications.     

Identification of 121 variants of honey bee Vitellogenin protein sequences with structural differences at functional sites

Proteins are under selection to maintain central functions and to accommodate needs that arise in ever‐changing environments. The positive selection and neutral drift that preserve functions result in a diversity of protein variants. The amount of diversity differs between proteins: multifunctional or disease‐related proteins tend to have fewer variants than proteins involved in some aspects of immunity. Our work focuses on the extensively studied protein Vitellogenin (Vg), which in honey bees (Apis mellifera) is multifunctional and highly expressed and plays roles in immunity. Yet, almost nothing is known about the natural variation in the coding sequences of this protein or how amino acid‐altering variants might impact structure–function relationships. Here, we map out allelic variation in honey bee Vg using biological samples from 15 countries. The successful barcoded amplicon Nanopore sequencing of 543 bees revealed 121 protein variants, indicating a high level of diversity in Vg. We find that the distribution of non‐synonymous single nucleotide polymorphisms (nsSNPs) differs between protein regions with different functions; domains involved in DNA and protein–protein interactions contain fewer nsSNPs than the protein''s lipid binding cavities. We outline how the central functions of the protein can be maintained in different variants and how the variation pattern may inform about selection from pathogens and nutrition.     

Improving integrative 3D modeling into low‐ to medium‐resolution electron microscopy structures with evolutionary couplings

Electron microscopy (EM) continues to provide near‐atomic resolution structures for well‐behaved proteins and protein complexes. Unfortunately, structures of some complexes are limited to low‐ to medium‐resolution due to biochemical or conformational heterogeneity. Thus, the application of unbiased systematic methods for fitting individual structures into EM maps is important. A method that employs co‐evolutionary information obtained solely from sequence data could prove invaluable for quick, confident localization of subunits within these structures. Here, we incorporate the co‐evolution of intermolecular amino acids as a new type of distance restraint in the integrative modeling platform in order to build three‐dimensional models of atomic structures into EM maps ranging from 10–14 Å in resolution. We validate this method using four complexes of known structure, where we highlight the conservation of intermolecular couplings despite dynamic conformational changes using the BAM complex. Finally, we use this method to assemble the subunits of the bacterial holo‐translocon into a model that agrees with previous biochemical data. The use of evolutionary couplings in integrative modeling improves systematic, unbiased fitting of atomic models into medium‐ to low‐resolution EM maps, providing additional information to integrative models lacking in spatial data.     

Structure,dynamics, and stability of the globular domain of human linker histone H1.0 and the role of positive charges

Linker histone H1 (H1) is an abundant chromatin‐binding protein that acts as an epigenetic regulator binding to nucleosomes and altering chromatin structures and dynamics. Nonetheless, the mechanistic details of its function remain poorly understood. Recent work suggest that the number and position of charged side chains on the globular domain (GD) of H1 influence chromatin structure and hence gene repression. Here, we solved the solution structure of the unbound GD of human H1.0, revealing that the structure is almost completely unperturbed by complex formation, except for a loop connecting two antiparallel β‐strands. We further quantified the role of the many positive charges of the GD for its structure and conformational stability through the analysis of 11 charge variants. We find that modulating the number of charges has little effect on the structure, but the stability is affected, resulting in a difference in melting temperature of 26 K between GD of net charge +5 versus +13. This result suggests that the large number of positive charges on H1‐GDs have evolved for function rather than structure and high stability. The stabilization of the GD upon binding to DNA can thus be expected to have a pronounced electrostatic component, a contribution that is amenable to modulation by posttranslational modifications, especially acetylation and phosphorylation.     

What the protein data bank tells us about the evolutionary conservation of protein conformational diversity

Proteins sample a multitude of different conformations by undergoing small‐ and large‐scale conformational changes that are often intrinsic to their functions. Information about these changes is often captured in the Protein Data Bank by the apparently redundant deposition of independent structural solutions of identical proteins. Here, we mine these data to examine the conservation of large‐scale conformational changes between homologous proteins. This is important for both practical reasons, such as predicting alternative conformations of a protein by comparative modeling, and conceptual reasons, such as understanding the extent of conservation of different features in evolution. To study this question, we introduce a novel approach to compare conformational changes between proteins by the comparison of their difference distance maps (DDMs). We found that proteins undergoing similar conformational changes have similar DDMs and that this similarity could be quantified by the correlation between the DDMs. By comparing the DDMs of homologous protein pairs, we found that large‐scale conformational changes show a high level of conservation across a broad range of sequence identities. This shows that conformational space is usually conserved between homologs, even relatively distant ones.