Protein Dynamics, Function, and Design

Much of our understanding of protein structure and mechanistic function has been derived from static high-resolution structures. As structural biology has continued to evolve it has become clear that high-resolution structures alone are unable to fully capture the mechanistic basis for protein structure and function in solution. Recently Hydrogen/Deuterium-exchange Mass Spectrometry (HDX-MS) has developed into a powerful and versatile tool for structural biologists that provides novel insights into protein structure and function. HDX-MS enables direct monitoring of a protein's structural fluctuations and conformational changes under native conditions in solution even as it is carrying out its functions. In this review, we focus on the use of HDX-MS to monitor these dynamic changes in proteins. We examine how HDX-MS has been applied to study protein structure and function in systems ranging from large, complex assemblies to intrinsically disordered proteins, and we discuss its use in probing conformational changes during protein folding and catalytic function. STATEMENT FOR A BROAD AUDIENCE: The biophysical and structural characterization of proteins provides novel insight into their functionalities. Protein motions, ranging from small scale local fluctuations to larger concerted structural rearrangements, often determine protein function. Hydrogen/Deuterium-exchange Mass Spectrometry (HDX-MS) has proven a powerful biophysical tool capable of probing changes in protein structure and dynamic protein motions that are often invisible to most other techniques.

Protein engineering, the process by which novel proteins with desired properties are developed, has grown by leaps and bounds since the first examples of protein mutagenesis were described over three decades ago. In the last 30 years, protein engineers have successfully created a wide range of proteins tailored to specific industrial, medical, and research applications. Yet challenges remain that prevent the engineering of complex protein functions on demand. With the aim of creating a forum for the protein engineering community, the inaugural Protein Engineering Canada Conference was held on June 20–22, 2014, in Ottawa, Canada. Over the course of the 2-day conference, 115 protein scientists from over 30 institutions representing 5 countries shared ideas, networked with colleagues, and learned about exciting new research covering a broad array of topics such as enzyme engineering, computational protein design, X-ray crystallography, protein NMR, molecular modeling, and protein evolution. This special issue highlights the breadth of research that was covered during the Protein Engineering Canada Conference by including previously unpublished research articles submitted by speakers as well as some selections that were submitted directly to Protein Science.

Apart from playing key roles in drug metabolism and adverse drug-drug interactions, CYPs are potential drug targets to treat a variety of diseases. The intervention of over expression of P450 1A1 (CYP1A1) in tumor cells is identified as a novel strategy for anticancer therapy. We investigated three isoforms of CYP1 family (CYP1A1, CYP1A2, and CYP1B1) for their substrate specificity. The understanding of macromolecular features that govern substrate specificity is required to understand the interplay between the protein function and dynamics. This can help in design of new antitumor molecule specifically metabolized by CYP1A1 to mediate their antitumor activity. In the present study, we carried out the comparative protein structure analysis of the three isoforms. Sequence alignment, root mean square deviation (RMSD) analysis, B-factor analysis was performed to give a better understanding of the macromolecular features involved in substrate specificity and to understand the interplay between protein dynamics and functions which will have important implications on rational design of anticancer drugs. We identified the differences in amino acid residues among the three isoforms of CYP1 family, which may account for differential substrate specificity. Six putative substrate recognition sequences are characterized along with the regions they form in the protein structure. Further the RMSD and B-factor analysis provides the information about the identified residues having the maximum RMSD and B-factor deviations.

An integrated view of protein structure, dynamics, and function is emerging, where proteins are considered as dynamically active assemblies and internal motions are closely linked to function such as enzyme catalysis. Further, the motion of solvent bound to external regions of protein impacts internal motions and, therefore, protein function. Recently, we discovered a network of protein vibrations in enzyme cyclophilin A, coupled to its catalytic activity of peptidyl-prolyl cis-trans isomerization. Detailed studies suggest that this network, extending from surface regions to active site, is a conserved part of enzyme structure and has a role in promoting catalysis. In this report, theoretical investigations of concerted conformational fluctuations occurring on microsecond and longer time scales within the discovered network are presented. Using a new technique, kinetic energy was added to protein vibrational modes corresponding to conformational fluctuations in the network. The results reveal that protein dynamics promotes catalysis by altering transition state barrier crossing behavior of reaction trajectories. An increase in transmission coefficient and number of productive trajectories with increasing amounts of kinetic energy in vibrational modes is observed. Variations in active site enzyme-substrate interactions near transition state are found to be correlated with barrier recrossings. Simulations also showed that energy transferred from first solvation shell to surface residues impacts catalysis through network fluctuations. The detailed characterization of network presented here indicates that protein dynamics plays a role in rate enhancement by enzymes. Therefore, coupled networks in enzymes have wide implications in understanding allostericity and cooperative effects, as well as protein engineering and drug design.

Proteins are inherently plastic molecules, whose function often critically depends on excursions between different molecular conformations (conformers). However, a rigorous understanding of the relation between a protein's structure, dynamics and function remains elusive. This is because many of the conformers on its energy landscape are only transiently formed and marginally populated (less than a few per cent of the total number of molecules), so that they cannot be individually characterized by most biophysical tools. Here we study a lysozyme mutant from phage T4 that binds hydrophobic molecules and populates an excited state transiently (about 1 ms) to about 3% at 25 °C (ref. 5). We show that such binding occurs only via the ground state, and present the atomic-level model of the 'invisible', excited state obtained using a combined strategy of relaxation-dispersion NMR (ref. 6) and CS-Rosetta model building that rationalizes this observation. The model was tested using structure-based design calculations identifying point mutants predicted to stabilize the excited state relative to the ground state. In this way a pair of mutations were introduced, inverting the relative populations of the ground and excited states and altering function. Our results suggest a mechanism for the evolution of a protein's function by changing the delicate balance between the states on its energy landscape. More generally, they show that our approach can generate and validate models of excited protein states.

Chemical cross-linking coupled with mass spectroscopy (CXMS) provides proximity information for the cross-linked residues and is used increasingly for modeling protein structures. However, experimentally identified cross-links are sometimes incompatible with the known structure of a protein, as the distance calculated between the cross-linked residues far exceeds the maximum length of the cross-linker. The discrepancies may persist even after eliminating potentially false cross-links and excluding intermolecular ones. Thus the "over-length" cross-links may arise from alternative excited-state conformation of the protein. Here we present a method and associated software DynaXL for visualizing the ensemble structures of multidomain proteins based on intramolecular cross-links identified by mass spectrometry with high confidence. Representing the cross-linkers and cross-linking reactions explicitly, we show that the protein excited-state structure can be modeled with as few as two over-length cross-links. We demonstrate the generality of our method with three systems: calmodulin, enzyme I, and glutamine-binding protein, and we show that these proteins alternate between different conformations for interacting with other proteins and ligands. Taken together, the over-length chemical cross-links contain valuable information about protein dynamics, and our findings here illustrate the relationship between dynamic domain movement and protein function.

A protein's structure is most often used to explain its function, but function also depends on dynamics. To date, protein dynamics have been studied only in vitro under dilute solution conditions where solute concentrations are typically less than 10 g/L, yet proteins function in a crowded environment where the solute concentration can exceed 400 g/L. Does the intracellular environment affect protein dynamics? The answer will help in assessing the biological significance of the NMR-derived dynamics data collected to date. We investigated fast protein dynamics inside living Escherichia coli by using in-cell NMR. The backbone dynamics of apocytochrome b5 were quantified using {1H}-15N nuclear Overhauser effect (nOe) measurements, which characterize motions on the pico- to nanosecond time scale. The overall trend of backbone dynamics remains the same in cells. Some of the nOe values differ, but most of the differences track the increased intracellular viscosity rather than a change in dynamics. Therefore, it appears that dilute solution steady-state {1H}-15N nOe measurements provide biologically relevant information about pico- to nanosecond backbone motion in proteins.

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using “Backrub” motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.

Proteins have evolved to carry out very specific functions within the cell by interacting with a diverse set of biomolecules. Understanding how a protein's higher order structure relates to its function is important for defining the molecular basis of these interactions. In recent years, mass spectrometry has become an important tool for dissecting protein structure and function. Using electrospray ionization (ESI)- and matrix-assisted laser desorption/ionization (MALDI)-based approaches, it has been possible to monitor protein folding, characterize noncovalent protein complexes, and assess the contribution of individual amino acid residues to a protein's function. Here, it is our goal to summarize these approaches and highlight recent, biologically relevant applications where mass spectrometry has provided unique insight into the mysteries of protein structure and function.

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

The unique biophysical properties of tryptophan residues have been exploited for decades to monitor protein structure and dynamics using a variety of spectroscopic techniques, such as fluorescence and nuclear magnetic resonance (NMR). We recently designed a tryptophan mutant in the regulatory N-domain of cardiac troponin C (F77W-cNTnC) to study the domain orientation of troponin C in muscle fibers using solid-state NMR. In our previous study, we determined the NMR structure of calcium-saturated mutant F77W-V82A-cNTnC in the presence of 19% 2,2,2-trifluoroethanol (TFE). TFE is a widely used cosolvent in the biophysical characterization of the solution structures of peptides and proteins. It is generally assumed that the structures are unchanged in the presence of cosolvents at relatively low concentrations, and this has been verified for TFE at the level of the overall secondary and tertiary structure for several calcium regulatory proteins. Here, we present the NMR solution structure of the calcium saturated F77W-cNTnC in presence of its biological binding partner troponin I peptide (cTnI(144-163)) and in the absence of TFE. We have also characterized a panel of six F77W-cNTnC structures in the presence and absence TFE, cTnI(144-163), and the extra mutation V82A, and used (19)F NMR to characterize the effect of TFE on the F77(5fW) analog. Our results show that although TFE did not perturb the overall protein structure, TFE did induce a change in the orientation of the indole ring of the buried tryptophan side chain from the anticipated position based upon homology with other proteins, highlighting the potential dangers of the use of cosolvents.

Protein dynamics play a crucial role in their diverse functions. The intracellular environment significantly influences protein dynamics, particularly for intrinsically disordered proteins (IDPs). To comprehensively capture structural information from various proteins wihtin cells and characterize protein dynamics, chemical cross-linking mass spectrometry was employed. In this study, we introduce a hierarchical decoding strategy that enables the investigation of protein dynamics in vivo. Computational analysis based on distance restraints derived from cross-links is used to infer protein dynamics in cells. To facilitate this analysis, we leverage the prior structure obtained from AlphaFold2. By employing this strategy, we can characterize the full-length structure of multi-domain proteins taking into account their distinct dynamic features. Furthermore, by combining restraint sampling with an unbiased sampling and evaluation approach, we can provide a comprehensive description of the intrinsic motion of IDPs. Consequently, the hierarchical strategy we propose holds significant potential in advancing our understanding of the molecular mechanisms that undelie protein functions in cells.

Measuring the positions and dynamics of proteins in intact tissues or whole animals is key to understanding protein function. However, to date, this is challenging, as the accessibility of large antibodies to dense tissues is often limited, and fluorescent proteins inserted close to a domain of interest may affect protein function. These complications apply in particular to muscle sarcomeres, arguably one of the most protein-dense assemblies in nature, which complicates studying sarcomere morphogenesis at molecular resolution. Here, we introduce a toolbox of nanobodies recognising various domains of the two Drosophila titin homologs, Sallimus and Projectin, as well as the key sarcomeric proteins Obscurin, α-Actinin, and Zasp52. We verified the superior labelling qualities of our nanobodies in muscle tissue as compared to antibodies. By applying our toolbox to larval muscles, we found a gigantic Sallimus isoform stretching more than 2 µm to bridge the sarcomeric I-band, while Projectin covers almost the entire myosin filaments in a polar orientation. Transgenic expression of tagged nanobodies confirmed their high affinity-binding without affecting target protein function. Finally, adding a degradation signal to anti-Sallimus nanobodies suggested that it is difficult to fully degrade Sallimus in mature sarcomeres; however, expression of these nanobodies caused developmental lethality. These results may inspire the generation of similar toolboxes for other large protein complexes in Drosophila or mammals.

During the last decade, network approaches became a powerful tool to describe protein structure and dynamics. Here we review the links between disordered proteins and the associated networks, and describe the consequences of local, mesoscopic and global network disorder on changes in protein structure and dynamics. We introduce a new classification of protein networks into 'cumulus-type', i.e., those similar to puffy (white) clouds, and 'stratus-type', i.e., those similar to flat, dense (dark) low-lying clouds, and relate these network types to protein disorder dynamics and to differences in energy transmission processes. In the first class, there is limited overlap between the modules, which implies higher rigidity of the individual units; there the conformational changes can be described by an 'energy transfer' mechanism. In the second class, the topology presents a compact structure with significant overlap between the modules; there the conformational changes can be described by 'multi-trajectories'; that is, multiple highly populated pathways. We further propose that disordered protein regions evolved to help other protein segments reach 'rarely visited' but functionally-related states. We also show the role of disorder in 'spatial games' of amino acids; highlight the effects of intrinsically disordered proteins (IDPs) on cellular networks and list some possible studies linking protein disorder and protein structure networks.

Ensemble-function relationships: From qualitative to quantitative relationships between protein structure and function