Nonlinear optical studies of heme protein dynamics: implications for proteins as hybrid states of matter

Protein structure is fundamentally related to function. However, static structures alone are insufficient to understand how a protein works. Dynamics play an equally important role. Given that proteins are highly associated aperiodic systems, it may be expected that protein dynamics would follow glass-like dynamics. However, protein functions occur on time scales orders of magnitude faster than the time scales typically associated with glassy systems. It is becoming clear that the reaction forces driving functions do not sample entirely the large number of configurations available to a protein but are highly directed along an optimized pathway. Could there be any correlation between specific topological features in protein structures and dynamics that leads to strongly correlated atomic displacements in the dynamical response to a perturbation? This review will try to provide an answer by focusing upon recent nonlinear optical studies with the aim of directly observing functionally important protein motions over the entire dynamic range of the protein response function. The specific system chosen is photoinduced dynamics of ligand dissociation at the active site in heme proteins, with myoglobin serving as the simplest model system. The energetics and nuclear motions from the very earliest events involved in bond breaking on the femtosecond time scale all the way out to ligand escape and bimolecular rebinding on the microsecond and millisecond time scale have been mapped out. The picture that is emerging is that the system consists of strongly coupled motions from the very instant the bond breaks at the active site that cascade into low frequency collective modes specific to the protein structure. It is this coupling that imparts the ability of a protein to function on time scales more commensurate with liquids while simultaneously conserving structural integrity akin to solids.     

Hierarchical Description and Extensive Classification of Protein Structural Changes by Motion Tree

The structures of the same protein, determined under different conditions, provide clues toward understanding the role of structural changes in the protein's function. Structural changes are usually identified as rigid-body motions, which are defined using a particular threshold of rigidity, such as domain motions. However, each protein actually undergoes motions with various size and magnitude ranges. In this study, to describe protein structural changes more comprehensively, we propose a method based on hierarchical clustering. This method enables the illustration of a wide range of protein motions in a single tree diagram, named the “Motion Tree”. We applied the method to 432 proteins exhibiting large structural changes and classified their Motion Trees in terms of the characteristic indices of the trees. This classification of the Motion Trees revealed clear relationships to their protein functions. Especially, complex structural changes are significantly correlated with multi-step protein functions.     

Close correspondence between the motions from principal component analysis of multiple HIV-1 protease structures and elastic network modes

The large number of available HIV-1 protease structures provides a remarkable sampling of conformations of the different conformational states, which can be viewed as direct structural information about the dynamics of the HIV-1 protease. After structure matching, we apply principal component analysis (PCA) to obtain the important apparent motions for both bound and unbound structures. There are significant similarities between the first few key motions and the first few low-frequency normal modes calculated from a static representative structure with an elastic network model (ENM), strongly suggesting that the variations among the observed structures and the corresponding conformational changes are facilitated by the low-frequency, global motions intrinsic to the structure. Similarities are also found when the approach is applied to an NMR ensemble, as well as to molecular dynamics (MD) trajectories. Thus, a sufficiently large number of experimental structures can directly provide important information about protein dynamics, but ENM can also provide similar sampling of conformations.     

Hierarchical domain‐motion analysis of conformational changes in sarcoplasmic reticulum Ca2+‐ATPase

Sarco(endo)plasmic reticulum Ca²⁺‐ATPase transports two Ca²⁺ per ATP‐hydrolyzed across biological membranes against a large concentration gradient by undergoing large conformational changes. Structural studies with X‐ray crystallography revealed functional roles of coupled motions between the cytoplasmic domains and the transmembrane helices in individual reaction steps. Here, we employed “Motion Tree (MT),” a tree diagram that describes a conformational change between two structures, and applied it to representative Ca²⁺‐ATPase structures. MT provides information of coupled rigid‐body motions of the ATPase in individual reaction steps. Fourteen rigid structural units, “common rigid domains (CRDs)” are identified from seven MTs throughout the whole enzymatic reaction cycle. CRDs likely act as not only the structural units, but also the functional units. Some of the functional importance has been newly revealed by the analysis. Stability of each CRD is examined on the morphing trajectories that cover seven conformational transitions. We confirmed that the large conformational changes are realized by the motions only in the flexible regions that connect CRDs. The Ca²⁺‐ATPase efficiently utilizes its intrinsic flexibility and rigidity to response different switches like ligand binding/dissociation or ATP hydrolysis. The analysis detects functional motions without extensive biological knowledge of experts, suggesting its general applicability to domain movements in other membrane proteins to deepen the understanding of protein structure and function. Proteins 2015; 83:746–756. © 2015 Wiley Periodicals, Inc.     

Molecular motions within the pore of voltage-dependent sodium channels.

The pores of ion channel proteins are often modeled as static structures. In this view, selectivity reflects rigidly constrained backbone orientations. Such a picture is at variance with the generalization that biological proteins are flexible, capable of major internal motions on biologically relevant time scales. We tested for motions in the sodium channel pore by systematically introducing pairs of cysteine residues throughout the pore-lining segments. Two distinct pairs of residues spontaneously formed disulfide bonds bridging domains I and II. Nine other permutations, involving all four domains, were capable of disulfide bonding in the presence of a redox catalyst. The results are inconsistent with a single fixed backbone structure for the pore; instead, the segments that line the permeation pathway appear capable of sizable motions.     

Envisioning the Loop Movements and Rotation of the Two Subdomains of Dihydrofolate Reductase by Elastic Normal Mode Analysis

Abstract

Normal mode analysis, using the elastic network model, has been employed to envision the low frequency normal mode motion trends in the structures of five intermediates and a transition state in the kinetic pathway of E. coli dihydrofolate reductase (DHFR). Five of the reaction pathway analog structures and a crystal structure resembling the transition state, using X-ray analyses determined by Kraut et al., have been adapted as structural models. The motions that poise pathways of the M20 loop transitions from closed to occluded conformations and sub domain rotation to close the substrate cleft, have been predicted and envisioned for the first time by this study. Pathway entries to the movement of the substrate binding cleft helices are also envisioned. These motions play roles in transition structure stabilization and in regulating the release of the product tetrahydrofolate (THF). The motions observed push the ground state conformation of each intermediate towards a higher energy sub state conformation. A set of conserved residues involved in the catalytic reactions and conformational changes, previously studied by kinetic, theoretical and NMR, have been analyzed. The importance of these motions in terms of protein dynamics are revealed and envisioned by the normal mode analysis. Additional residues are proposed as candidates for further study of their potential promotional function.     

Classification and annotation of the relationship between protein structural change and ligand binding

The causal relationship between protein structural change and ligand binding was classified and annotated for 839 nonredundant pairs of crystal structures in the Protein Data Bank—one with and the other without a bound low-molecular-weight ligand molecule. Protein structural changes were first classified into either domain or local motions depending on the size of the moving protein segments. Whether the protein motion was coupled with ligand binding was then evaluated based on the location of the ligand binding site and by application of the linear response theory of protein structural change. Protein motions coupled with ligand binding were further classified into either closure or opening motions. This classification revealed the following: (i) domain motions coupled with ligand binding are dominated by closure motions, which can be described by the linear response theory; (ii) local motions frequently accompany order-disorder or α-helix-coil conformational transitions; and (iii) transferase activity (Enzyme Commission   number 2) is the predominant function among coupled domain closure motions. This could be explained by the closure motion acting to insulate the reaction site of these enzymes from environmental water.     

Protein motions visualized by femtosecond time-resolved crystallography: The case of photosensory vs photosynthetic proteins

Proteins are dynamic objects and undergo conformational changes when functioning. These changes range from interconversion between states in equilibrium to ultrafast and coherent structural motions within one perturbed state. Time-resolved serial femtosecond crystallography at free-electron X-ray lasers can unravel structural changes with atomic resolution and down to femtosecond time scales. In this review, we summarize recent advances on detecting structural changes for phytochrome photosensor proteins and a bacterial photosynthetic reaction center. In the phytochrome structural changes are extensive and involve major rearrangements of many amino acids and water molecules, accompanying the regulation of its biochemical activity, whereas in the photosynthetic reaction center protein the structural changes are smaller, more localized, and are optimized to facilitate electron transfer along the chromophores. The detected structural motions underpin the proteins’ function, providing a showcase for the importance of detecting ultrafast protein structural dynamics.     

Coupling of structural fluctuations to deamidation reaction in triosephosphate isomerase by Gaussian network model

We study the structural fluctuations of triosephosphate isomerase (TIM) by an elastic model, namely, the Gaussian network model (GNM), to identify a network of coupled motions in the allosteric communication between its deamidation and catalytic sites, and the promoting motions for the deamidation activity. For this, three TIM structures have been studied: one crystal structure and two model structures designed to describe different putative models for the deamidation reaction taking place at the subunit interface. The structural fluctuations have been mapped on the functional properties; then the differences in the fluctuations between the two models in relation to the deamidation reaction have been considered. The results demonstrate that the qualitative picture of the mean-square fluctuations and the correlations between the fluctuations are similar in both, but the differences may affect the observed barrier height of the deamidation reaction. The higher packing density at regions close to deamidation sites, reflected by the high-frequency fluctuating residues in the respective regions, the stronger positive correlation between the fluctuations of the deamidation sites, and enhanced positive correlation of the primary deamidation site with the extended vicinity of the catalytic region on the juxtaposed unit promote the probability of the deamidation reaction. The results in general emphasize the importance of structural fluctuations in enzyme reactions, as well as proposing the present methodology as a plausible approach for studies on the network of coupled promoting motions in protein functions.     

High-throughput modeling and analysis of protein structural dynamics

Protein function is a dynamic property closely related to the conformational mechanisms of protein structure in its physiological environment. To understand and control the function of target proteins, it becomes increasingly important to develop methods and tools for predicting collective motions at the molecular level. In this article, we review computational methods for predicting conformational dynamics and discuss software tools for data analysis. In particular, we discuss a high-throughput, web-based system called iGNM for protein structural dynamics. iGNM contains a database of protein motions for more than 20 000 PDB structures and supports online calculations for newly deposited PDB structures or user-modified structures. iGNM allows dynamics analysis of protein structures ranging from enzymes to large complexes and assemblies, and enables the exploration of protein sequence-structure-dynamics-function relations.     

Concerted motions in copper plastocyanin and azurin: an essential dynamics study

Essential dynamics analysis of molecular dynamics simulation trajectories (1.1 ns) of two copper containing electron transfer proteins, plastocyanin and azurin, has been performed. The protein essential modes have been analysed in order to identify large concerted motions which could be relevant for the electron transfer function exerted by these proteins. The analysis, conducted for temporal windows of different lengths along the protein trajectories, shows a rapid convergence and indicates that for both the proteins the predominant internal motions occur in a subspace of only a few degrees of freedom. Moreover, it is found that for both the proteins the likely binding sites (i.e. the hydrophobic and negative patches) with the reaction partners move in a concerted fashion with a few structural regions far from the active site. Such results are discussed in connection with the possible involvement of large concerted motions in the recognition and binding interaction with physiological electron transfer partners.     

Insect lipoprotein biogenesis depends on an amphipathic beta cluster in apolipophorin II/I and is stimulated by microsomal triglyceride transfer protein

Lipoproteins transport lipids in the circulation of an evolutionally wide diversity of animals. The pathway for lipoprotein biogenesis has been revealed to a large extent in mammals only, in which apolipoprotein B (apoB) acquires lipids via the assistance of microsomal triglyceride transfer protein (MTP) and binds them by means of amphipathic protein structures. To investigate whether this is a common mechanism for lipoprotein biogenesis in animals, we studied the structural elements involved in the assembly of the insect lipoprotein, lipophorin. LOCATE sequence analysis predicted that the insect lipoprotein precursor, apolipophorin II/I (apoLp-II/I), contains clusters of amphipathic alpha-helices and beta-strands, organized along the protein as N-alpha(1)-beta-alpha(2)-C, reminiscent of a truncated form of apoB. Recombinant expression of a series of C-terminal truncation variants of Locusta migratoria apoLp-II/I in an insect cell (Sf9) expression system revealed that the formation of a buoyant high density lipoprotein requires the amphipathic beta cluster. Coexpression of apoLp-II/I with the MTP homolog of Drosophila melanogaster affected insect lipoprotein biogenesis quantitatively as well as qualitatively, as the secretion of apoLp-II/I proteins was increased several-fold and the buoyant density of the secreted lipoprotein decreased concomitantly, indicative of augmented lipidation. Based on these findings, we propose that, despite specific modifications, the assembly of lipoproteins involves MTP as well as amphipathic structures in the apolipoprotein carrier, both in mammals and insects. Thus, lipoprotein biogenesis in animals appears to rely on structural elements that are of early metazoan origin.     

Protein structural change upon ligand binding correlates with enzymatic reaction mechanism

Overall structural changes of enzymes in response to ligand binding were investigated by database analysis of 62 non-redundant enzymes whose ligand-unbound and ligand-bound forms were available in the Protein Data Bank. The results of analysis indicate that transferases often undergo large rigid-body domain motions upon ligand binding, while other enzymes, most typically, hydrolases, change their structures to a small extent. It was also found that the solvent accessibility of the substrate molecule was low in transferases but high in hydrolases. These differences are explained by the enzymatic reaction mechanisms. The transferase reaction requires the catalytic groups to be insulated from the water environment, and thus transferases bury the ligand molecule inside the protein by closing the cleft. On the other hand, the hydrolase reaction involves the surrounding water molecules and occurs at the protein surface, requiring only a small structural change.     

Using cross-links to study ribosomal dynamics

After publications of 3-D models of a static ribosome and its large and small subunits, one of the next tasks is to recognize movable ribosomal elements responsible for mechanical shifts during protein synthesis. Statistic analysis of available cross-linking data allowed us to reveal three well separated groups of motions in the ribosome: I, mean magnitude of 10 A; II, most abundant, centered at 20 A and of wide dispersion, and III, sparsely populated, with large distances up to 95 A. The last group, III, comprises elements, like the L7/12-stalk and the L1- protuberance, that adopt different positions in crystallographic or electron micrographic structures, and neighboring hairpins 88 and 89, indicating mobility. We demonstrate that the cross-linking method can be applied to study ribosomal dynamics, including large-scale functional movements and, in particular, to estimate which structures participate in molecular switches.     

Exploring the structural dynamics of the E.coli chaperonin GroEL using translation-libration-screw crystallographic refinement of intermediate states

Large rigid-body domain movements are critical to GroEL-mediated protein folding, especially apical domain elevation and twist associated with the formation of a folding chamber upon binding ATP and co-chaperonin GroES. Here, we have modeled the anisotropic displacements of GroEL domains from various crystallized states, unliganded GroEL, ATPgammaS-bound, ADP-AlFx/GroES-bound, and ADP/GroES bound, using translation-libration-screw (TLS) analysis. Remarkably, the TLS results show that the inherent motions of unliganded GroEL, a polypeptide-accepting state, are biased along the transition pathway that leads to the folding-active state. In the ADP-AlFx/GroES-bound folding-active state the dynamic modes of the apical domains become reoriented and coupled to the motions of bound GroES. The ADP/GroES complex exhibits these same motions, but they are increased in magnitude, potentially reflecting the decreased stability of the complex after nucleotide hydrolysis. Our results have allowed the visualization of the anisotropic molecular motions that link the static conformations previously observed by X-ray crystallography. Application of the same analyses to other macromolecules where rigid body motions occur may give insight into the large scale dynamics critical for function and thus has the potential to extend our fundamental understanding of molecular machines.     

A comparative study of dynamic structures between phage 434 Cro and repressor proteins by normal mode analysis

Two DNA binding proteins, Cro and the amino-terminal domain of the repressor of bacteriophage 434 (434 Cro and 434 repressor) that regulate gene expression and contain a helix-turn-helix (HTH) motif responsible for their site-specific DNA recognition adopt very similar three-dimensional structures when compared to each other. To reveal structural differences between these two similar proteins, their dynamic structures, as examined by normal mode analysis, are compared in this paper. Two kinds of structural data, one for the monomer and the other for a complex with DNA, for each protein, are used in the analyses. From a comparison between the monomers it is found that the interactions of Ala-24 in 434 Cro or Val-24 in 434 repressor, both located in the HTH motif, with residues 44, 47, 48, and 51 located in the domain facing the motif, and the interactions between residues 17, 18, 28, and 32, located in the HTH motif, cause significant differences in the correlative motions of these residues. From the comparison between the monomer and the complex with DNA for each protein, it was found that the first helix in the HTH motif is distorted in the complex form. While the residues in the HTH motif in 434 Cro have relatively larger positive correlation coefficients of motions with other residues within the HTH motif, such correlations are not large in the HTH motif of 434 repressor. It is suggestive to their specificity because the 434 repressor is less specific than 434 Cro. Although a structural comparison of proteins has been performed mainly from a static or geometrical point of view, this study demonstrates that the comparison from a dynamic point of view, using the normal mode analysis, is useful and convenient to explore a difference that is difficult to find only from a geometrical point of view, especially for proteins very similar in structure. © 1996 Wiley-Liss, Inc.     

Using Cross-links to Study Ribosomal Dynamics

Abstract After publications of 3-D models of a static ribosome and its large and small subunits, one of the next tasks is to recognize movable ribosomal elements responsible for mechanical shifts during protein synthesis. Statistic analysis of available cross-linking data allowed us to reveal three well separated groups of motions in the ribosome: I, mean magnitude of 10?; II, most abundant, centered at 20? and of wide dispersion, and III, sparsely populated, with large distances up to 95?. The last group, III, comprises elements, like the L7/12-stalk and the L1- protuberance, that adopt different positions in crystallographic or electron micrographic structures, and neighboring hairpins 88 and 89, indicating mobility. We demonstrate that the cross- linking method can be applied to study ribosomal dynamics, including large-scale functional movements and, in particular, to estimate which structures participate in molecular switches.     

Energy Landscapes and Catalysis in Nitric-oxide Synthase

Nitric oxide (NO) plays diverse roles in mammalian physiology. It is involved in blood pressure regulation, neurotransmission, and immune response, and is generated through complex electron transfer reactions catalyzed by NO synthases (NOS). In neuronal NOS (nNOS), protein domain dynamics and calmodulin binding are implicated in regulating electron flow from NADPH, through the FAD and FMN cofactors, to the heme oxygenase domain, the site of NO generation. Simple models based on crystal structures of nNOS reductase have invoked a role for large scale motions of the FMN-binding domain in shuttling electrons from the FAD-binding domain to the heme oxygenase domain. However, molecular level insight of the dynamic structural transitions in NOS enzymes during enzyme catalysis is lacking. We use pulsed electron-electron double resonance spectroscopy to derive inter-domain distance relationships in multiple conformational states of nNOS. These distance relationships are correlated with enzymatic activity through variable pressure kinetic studies of electron transfer and turnover. The binding of NADPH and calmodulin are shown to influence interdomain distance relationships as well as reaction chemistry. An important effect of calmodulin binding is to suppress adventitious electron transfer from nNOS to molecular oxygen and thereby preventing accumulation of reactive oxygen species. A complex landscape of conformations is required for nNOS catalysis beyond the simple models derived from static crystal structures of nNOS reductase. Detailed understanding of this landscape advances our understanding of nNOS catalysis/electron transfer, and could provide new opportunities for the discovery of small molecule inhibitors that bind at dynamic protein interfaces of this multidimensional energy landscape.     

Crystal structure of the plant epigenetic protein arginine methyltransferase 10

Protein arginine methyltransferase 10 (PRMT10) is a type I arginine methyltransferase that is essential for regulating flowering time in Arabidopsis thaliana. We present a 2.6 Å resolution crystal structure of A. thaliana PRMT 10 (AtPRMT10) in complex with a reaction product, S-adenosylhomocysteine. The structure reveals a dimerization arm that is 12-20 residues longer than PRMT structures elucidated previously; as a result, the essential AtPRMT10 dimer exhibits a large central cavity and a distinctly accessible active site. We employ molecular dynamics to examine how dimerization facilitates AtPRMT10 motions necessary for activity, and we show that these motions are conserved in other PRMT enzymes. Finally, functional data reveal that the 10 N-terminal residues of AtPRMT10 influence substrate specificity, and that enzyme activity is dependent on substrate protein sequences distal from the methylation site. Taken together, these data provide insights into the molecular mechanism of AtPRMT10, as well as other members of the PRMT family of enzymes. They highlight differences between AtPRMT10 and other PRMTs but also indicate that motions are a conserved element of PRMT function.