Segment swapping aided the evolution of enzyme function: The case of uroporphyrinogen III synthase

In an earlier study, we showed that two‐domain segment‐swapped proteins can evolve by domain swapping and fusion, resulting in a protein with two linkers connecting its domains. We proposed that a potential evolutionary advantage of this topology may be the restriction of interdomain motions, which may facilitate domain closure by a hinge‐like movement, crucial for the function of many enzymes. Here, we test this hypothesis computationally on uroporphyrinogen III synthase, a two‐domain segment‐swapped enzyme essential in porphyrin metabolism. To compare the interdomain flexibility between the wild‐type, segment‐swapped enzyme (having two interdomain linkers) and circular permutants of the same enzyme having only one interdomain linker, we performed geometric and molecular dynamics simulations for these species in their ligand‐free and ligand‐bound forms. We find that in the ligand‐free form, interdomain motions in the wild‐type enzyme are significantly more restricted than they would be with only one interdomain linker, while the flexibility difference is negligible in the ligand‐bound form. We also estimated the entropy costs of ligand binding associated with the interdomain motions, and find that the change in domain connectivity due to segment swapping results in a reduction of this entropy cost, corresponding to ～20% of the total ligand binding free energy. In addition, the restriction of interdomain motions may also help the functional domain‐closure motion required for catalysis. This suggests that the evolution of the segment‐swapped topology facilitated the evolution of enzyme function for this protein by influencing its dynamic properties. Proteins 2016; 85:46–53. © 2016 Wiley Periodicals, Inc.     

Molecular dynamics simulations elucidate the mode of protein recognition by Skp1 and the F‐box domain in the SCF complex

Polyubiquitination of the target protein by a ubiquitin transferring machinery is key to various cellular processes. E3 ligase Skp1‐Cul1‐F‐box (SCF) is one such complex which plays crucial role in substrate recognition and transfer of the ubiquitin molecule. Previous computational studies have focused on S‐phase kinase‐associated protein 2 (Skp2), cullin, and RING‐finger proteins of this complex, but the roles of the adapter protein Skp1 and F‐box domain of Skp2 have not been determined. Using sub‐microsecond molecular dynamics simulations of full‐length Skp1, unbound Skp2, Skp2‐Cks1 (Cks1: Cyclin‐dependent kinases regulatory subunit 1), Skp1‐Skp2, and Skp1‐Skp2‐Cks1 complexes, we have elucidated the function of Skp1 and the F‐box domain of Skp2. We found that the L₁₆ loop of Skp1_, which was deleted in previous X‐ray crystallography studies, can offer additional stability to the ternary complex via its interactions with the C‐terminal tail of Skp2. Moreover, Skp1 helices H6, H7, and H8 display vivid conformational flexibility when not bound to Skp2, suggesting that these helices can recognize and lock the F‐box proteins. Furthermore, we observed that the F‐box domain could rotate (5°–129°), and that the binding partner determined the degree of conformational flexibility. Finally, Skp1 and Skp2 were found to execute a domain motion in Skp1‐Skp2 and Skp1‐Skp2‐Cks1 complexes that could decrease the distance between ubiquitination site of the substrate and the ubiquitin molecule by 3 nm. Thus, we propose that both the F‐box domain of Skp2 and Skp1‐Skp2 domain motions displaying preferential conformational control can together facilitate polyubiquitination of a wide variety of substrates. Proteins 2016; 84:159–171. © 2015 Wiley Periodicals, Inc.     

Common functionally important motions of the nucleotide‐binding domain of Hsp70

The 70 kDa heat shock proteins (Hsp70) are a family of molecular chaperones involved in protein folding, aggregate prevention, and protein disaggregation. They consist of the substrate‐binding domain (SBD) that binds client substrates, and the nucleotide‐binding domain (NBD), whose cycles of nucleotide hydrolysis and exchange underpin the activity of the chaperone. To characterize the structure–function relationships that link the binding state of the NBD to its conformational behavior, we analyzed the dynamics of the NBD of the Hsp70 chaperone from Bos taurus (PDB 3C7N:B) by all‐atom canonical molecular dynamics simulations. It was found that essential motions within the NBD fall into three major classes: the mutual class, reflecting tendencies common to all binding states, and the ADP‐ and ATP‐unique classes, which reflect conformational trends that are unique to either the ADP‐ or ATP‐bound states, respectively. “Mutual” class motions generally describe “in‐plane” and/or “out‐of‐plane” (scissor‐like) rotation of the subdomains within the NBD. This result is consistent with experimental nuclear magnetic resonance data on the NBD. The “unique” class motions target specific regions on the NBD, usually surface loops or sites involved in nucleotide binding and are, therefore, expected to be involved in allostery and signal transmission. For all classes, and especially for those of the “unique” type, regions of enhanced mobility can be identified; these are termed “hot spots,” and their locations generally parallel those found by NMR spectroscopy. The presence of magnesium and potassium cations in the nucleotide‐binding pocket was also found to influence the dynamics of the NBD significantly. Proteins 2015; 83:282–299. © 2014 Wiley Periodicals, Inc.     

NMR characterization of the conformational fluctuations of the human lymphocyte function‐associated antigen‐1 I‐domain

Lymphocyte function‐associated antigen‐1 (LFA‐1) is an integrin protein that transmits information across the plasma membrane through the so‐called inside‐out and outside‐in signaling mechanisms. To investigate these mechanisms, we carried out an NMR analysis of the dynamics of the LFA‐1 I‐domain, which has enabled us to characterize the motions of this domain on a broad range of timescales. We studied first the internal motions on the nanosecond timescale by spin relaxation measurements and model‐free analysis. We then extended this analysis to the millisecond timescale motions by measuring ¹⁵N‐¹H residual dipolar couplings of the backbone amide groups. We analyzed these results in the context of the three major conformational states of the I‐domain using their corresponding X‐ray crystallographic structures. Our results highlight the importance of the low‐frequency motions of the LFA‐1 I‐domain in the inactive apo‐state. We found in particular that α‐helix 7 is in a position in the apo‐closed state that cannot be fully described by any of the existing X‐ray structures, as it appears to be in dynamic exchange between different conformations. This type of motion seems to represent an inherent property of the LFA‐1 I‐domain and might be relevant for controlling the access to the allosteric binding pocket, as well as for the downward displacement of α‐helix 7 that is required for the activation of LFA‐1.     

Nanoscale protein domain motion and long-range allostery in signaling proteins—a view from neutron spin echo spectroscopy

Many cellular proteins are multi-domain proteins. Coupled domain–domain interactions in these multidomain proteins are important for the allosteric relay of signals in the cellular signaling networks. We have initiated the application of neutron spin echo spectroscopy to the study of nanoscale protein domain motions on submicrosecond time scales and on nanometer length scale. Our NSE experiments reveal the activation of protein domain motions over a long distance of over more than 100 Å in a multidomain scaffolding protein NHERF1 upon binding to another protein, Ezrin. Such activation of nanoscale protein domain motions is correlated with the allosteric assembly of multi-protein complexes by NHERF1 and Ezrin. Here, we summarize the theoretical framework that we have developed, which uses simple concepts from nonequilibrium statistical mechanics to interpret the NSE data, and employs a mobility tensor to describe nanoscale protein domain motion. Extracting nanoscale protein domain motion from the NSE does not require elaborate molecular dynamics simulations, nor complex fits to rotational motion, nor elastic network models. The approach is thus more robust than multiparameter techniques that require untestable assumptions. We also demonstrate that an experimental scheme of selective deuteration of a protein subunit in a complex can highlight and amplify specific domain dynamics from the abundant global translational and rotational motions in a protein. We expect NSE to provide a unique tool to determine nanoscale protein dynamics for the understanding of protein functions, such as how signals are propagated in a protein over a long distance to a distal domain.     

The red‐edge effects: 30 years of exploration

In 1970, three laboratories independently made a discovery that, for aromatic fluorophores embedded into different rigid and highly viscous media, the spectroscopic properties do not conform to classical rules. The fluorescence spectra can depend on excitation wavelength, and the excited‐state energy transfer, if present, fails at the ‘red’ excitation edge. These red‐edge effects were related to the existence of excited‐state distribution of fluorophores on their interaction energy with the environment and the slow rate of dielectric relaxation of this environment. In these conditions the site‐selection can be provided by variation of the energy of illuminating light quanta, and the behaviour of selected species can be followed as a function of time and other variables. These observations found extensive application in different areas of research: colloid and polymer science, molecular biophysics, photochemistry and photobiology. In particular, they led to the development of very productive methods of studying the dynamics of dielectric relaxations in protein and membranes, using the tryptophan emission and the emission of a variety of probes. These studies were extended to the time domain with the observation of new site‐selective effects in emission intensity and anisotropy decays. They stimulated the emergence and development of cryogenic energy‐selective and single‐molecular techniques that became valuable tools in their own right in chemistry and biophysics research. Site‐selection effects were discovered for electron‐transfer and proton‐transfer reactions if they depended on the dynamics of the environment. This review is focused on the progress in the field of red‐edge effects, their applications and prospects. Copyright © 2002 John Wiley & Sons, Ltd.     

Method for identification of rigid domains and hinge residues in proteins based on exhaustive enumeration

Many proteins undergo large‐scale motions where relatively rigid domains move against each other. The identification of rigid domains, as well as the hinge residues important for their relative movements, is important for various applications including flexible docking simulations. In this work, we develop a method for protein rigid domain identification based on an exhaustive enumeration of maximal rigid domains, the rigid domains not fully contained within other domains. The computation is performed by mapping the problem to that of finding maximal cliques in a graph. A minimal set of rigid domains are then selected, which cover most of the protein with minimal overlap. In contrast to the results of existing methods that partition a protein into non‐overlapping domains using approximate algorithms, the rigid domains obtained from exact enumeration naturally contain overlapping regions, which correspond to the hinges of the inter‐domain bending motion. The performance of the algorithm is demonstrated on several proteins. Proteins 2015; 83:1054–1067. © 2015 Wiley Periodicals, Inc.     

Thermal‐ and urea‐induced unfolding processes of glutathione S‐transferase by molecular dynamics simulation

The Schistosoma juponicum 26 kDa glutathione S‐transferase (sj26GST) consists of the N‐terminal domain (N‐domain), containing three alpha‐helices (named H1‐H3) and four anti‐parallel beta‐strands (S1‐S4), and the C‐terminal domain (C‐domain), comprising five alpha‐helices (named H4‐H8). In present work, molecular dynamics simulations and fluorescence spectroscopic were used to gain insights into the unfolding process of sj26GST. The molecular dynamics simulations on sj26GST subunit both in water and in 8 M urea were carried out at 300 K, 400 K and 500 K, respectively. Spectroscopic measurements were employed to monitor structural changes. Molecular dynamics simulations of sj26GST subunit induced by urea and temperature showed that the initial unfolding step of sj26GST both in water and urea occurred on N‐domain, involving the disruption of helices H2, H3 and strands S3 and S4, whereas H6 was the last region exposed to solution and was the last helix to unfold. Moreover, simulations analyses combining with fluorescence and circular dichroism spectra indicated that N‐domain could not fold independent, suggesting that correct folding of N‐domain depended on its interactions with C‐domain. We further proposed that the folding of GSTs could begin with the hydrophobic collapse of C‐domain whose H4, H5, H6 and H7 could move close to each other and form a hydrophobic core, especially H6 wrapped in the hydrophobic center and beginning spontaneous formation of the helix. S3, S4, H3, and H2 could form in the wake of the interaction between C‐domain and N‐domain. The paper can offer insights into the molecular mechanism of GSTs unfolding. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 247–259, 2015.     

Molecular dynamics study of time-correlated protein domain motions and molecular flexibility: cytochrome P450BM-3.

Time-correlated atomic motions were used to characterize protein domain boundaries from atomic coordinates generated by molecular dynamics simulations. A novel application of the dynamical cross-correlation matrix (DCCM) analysis tool was used to help identify putative protein domains. In implementing this new approach, several DCCM maps were calculated, each using a different coordinate reference frame from which protein domain boundaries and protein domain residue constituents could be identified. Cytochrome P450BM-3, from Bacillus megaterium, was used as the model protein in this study. The analyses indicated that the simulated protein comprises three distinct domain regions; in contrast, only two protein domains were identified in the original crystal structure report. Specifically, the DCCM analyses showed that the F-G helix region was a separate domain entity and not a part of the alpha domain, as previously designated. The simulations demonstrated that the domain motions of the F-G helix region effected both the size and shape of the enzyme active site, and that the dynamics of the F-G helix domain could possibly control access of substrate to the binding pocket.     

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.     

Toward decrypting the allosteric mechanism of the ryanodine receptor based on coarse‐grained structural and dynamic modeling

The ryanodine receptors (RyRs) are a family of calcium (Ca) channels that regulate Ca release by undergoing a closed‐to‐open gating transition in response to action potential or Ca binding. The allosteric mechanism of RyRs gating, which is activated/regulated by ligand/protein binding >200 Å away from the channel gate, remains elusive for the lack of high‐resolution structures. Recent solution of the closed‐form structures of the RyR1 isoform by cryo‐electron microscopy has paved the way for detailed structure‐driven studies of RyRs functions. Toward elucidating the allosteric mechanism of RyRs gating, we performed coarse‐grained modeling based on the newly solved closed‐form structures of RyR1. Our normal mode analysis captured a key mode of collective motions dominating the observed structural variations in RyR1, which features large outward and downward movements of the peripheral domains with the channel remaining closed, and involves hotspot residues that overlap well with key functional sites and disease mutations. In particular, we found a key interaction between a peripheral domain and the Ca‐binding EF hand domain, which may allow for direct coupling of Ca binding to the collective motions as captured by the above mode. This key mode was robustly reproduced by the normal mode analysis of the other two closed‐form structures of RyR1 solved independently. To elucidate the closed‐to‐open conformational changes in RyR1 with amino‐acid level of details, we flexibly fitted the closed‐form structures of RyR1 into a 10‐Å cryo‐electron microscopy map of the open state. We observed extensive structural changes involving the peripheral domains and the central domains, resulting in the channel pore opening. In sum, our findings have offered unprecedented structural and dynamic insights to the allosteric mechanism of RyR1 via modulation of the key collective motions involved in RyR1 gating. The predicted hotspot residues and open‐form conformation of RyR1 will guide future mutational and functional studies. Proteins 2015; 83:2307–2318. © 2015 Wiley Periodicals, Inc.     

Conformational frustration in calmodulin–target recognition

Calmodulin (CaM) is a primary calcium (Ca²⁺)‐signaling protein that specifically recognizes and activates highly diverse target proteins. We explored the molecular basis of target recognition of CaM with peptides representing the CaM‐binding domains from two Ca²⁺‐CaM‐dependent kinases, CaMKI and CaMKII, by employing experimentally constrained molecular simulations. Detailed binding route analysis revealed that the two CaM target peptides, although similar in length and net charge, follow distinct routes that lead to a higher binding frustration in the CaM–CaMKII complex than in the CaM–CaMKI complex. We discovered that the molecular origin of the binding frustration is caused by intermolecular contacts formed with the C‐domain of CaM that need to be broken before the formation of intermolecular contacts with the N‐domain of CaM. We argue that the binding frustration is important for determining the kinetics of the recognition process of proteins involving large structural fluctuations. Copyright © 2015 John Wiley & Sons, Ltd.     

Investigation of Domain Motions in Bacteriophage T4 Lysozyme

Abstract

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge loci were found in the simulation. One locus comprises residues 8–14 near the C-terminal of the A helix; the other site, residues 77–83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19–24) acted as a mobile ‘flap’ and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.     

Comprehensive analysis of motions in molecular dynamics trajectories of the actin capping protein and its inhibitor complexes

The actin capping protein (CP) binds to actin filaments to block further elongation. The capping activity is inhibited by proteins V‐1 and CARMIL interacting with CP via steric and allosteric mechanisms, respectively. The crystal structures of free CP, CP/V‐1, and CP/CARMIL complexes suggest that the binding of CARMIL alters the flexibility of CP rather than the overall structure of CP, and this is an allosteric inhibition mechanism. Here, we performed molecular dynamics (MD) simulations of CP in the free form, and in complex with CARMIL or V‐1. The resulting trajectories were analyzed exhaustively using Motion Tree, which identifies various rigid‐body motions ranging from small local motions to large domain motions. After enumerating all the motions, CP flexibilities with different ligands were characterized by a list of frequencies for 20 dominant rigid‐body motions, some of which were not identified in previous studies. The comparative analysis highlights the influence of the binding of the CARMIL peptide to CP flexibility. In free CP and the CP/V‐1 complex, domain motions around a large crevice between the N‐stalk and the CP‐S domain occur frequently. The CARMIL peptide binds the crevice and suppresses the motions effectively. In addition, the binding of the CARMIL peptide enhances and alters local motions around the pocket that participates in V‐1 binding. These newly identified motions are likely to suppress the binding of V‐1 to CP. The observed changes in CP motion provide insights that describe the mechanism of allosteric regulation by CARMIL through modulating CP flexibility. Proteins 2016; 84:948–956. © 2016 Wiley Periodicals, Inc.     

Differential backbone dynamics of companion helices in the extended helical coiled‐coil domain of a bacterial chemoreceptor

Cytoplasmic domains of transmembrane bacterial chemoreceptors are largely extended four‐helix coiled coils. Previous observations suggested the domain was structurally dynamic. We probed directly backbone dynamics of this domain of the transmembrane chemoreceptor Tar from Escherichia coli using site‐directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy. Spin labels were positioned on solvent‐exposed helical faces because EPR spectra for such positions reflect primarily polypeptide backbone movements. We acquired spectra for spin‐labeled, intact receptor homodimers solubilized in detergent or inserted into native E. coli lipid bilayers in Nanodiscs, characterizing 16 positions distributed throughout the cytoplasmic domain and on both helices of its helical hairpins, one amino terminal to the membrane‐distal tight turn (N‐helix), and the other carboxyl terminal (C‐helix). Detergent solubilization increased backbone dynamics for much of the domain, suggesting that loss of receptor activities upon solubilization reflects wide‐spread destabilization. For receptors in either condition, we observed an unanticipated difference between the N‐ and C‐helices. For bilayer‐inserted receptors, EPR spectra from sites in the membrane‐distal protein‐interaction region and throughout the C‐helix were typical of well‐structured helices. In contrast, for approximately two‐thirds of the N‐helix, from its origin as the AS‐2 helix of the membrane‐proximal HAMP domain to the beginning of the membrane‐distal protein‐interaction region, spectra had a significantly mobile component, estimated by spectral deconvolution to average approximately 15%. Differential helical dynamics suggests a four‐helix bundle organization with a pair of core scaffold helices and two more dynamic partner helices. This newly observed feature of chemoreceptor structure could be involved in receptor function.     

Dynamics of a truncated prion protein,PrP(113–231), from 15N NMR relaxation: Order parameters calculated and slow conformational fluctuations localized to a distinct region

Prion diseases are associated with the misfolding of the prion protein (PrP^C) from a largely α‐helical isoform to a β‐sheet rich oligomer (PrP^Sc). Flexibility of the polypeptide could contribute to the ability of PrP^C to undergo the conformational rearrangement during PrP^C–PrP^Sc interactions, which then leads to the misfolded isoform. We have therefore examined the molecular motions of mouse PrP^C, residues 113–231, in solution, using ¹⁵N NMR relaxation measurements. A truncated fragment has been used to eliminate the effect of the 90‐residue unstructured tail of PrP^C so the dynamics of the structured domain can be studied in isolation. ¹⁵N longitudinal (T₁) and transverse relaxation (T₂) times as well as the proton‐nitrogen nuclear Overhauser effects have been used to calculate the spectral density at three frequencies, 0, ω_N, and 0.87ω_H. Spectral densities at each residue indicate various time‐scale motions of the main‐chain. Even within the structured domain of PrP^C, a diverse range of motions are observed. We find that removal of the tail increases T₂ relaxation times significantly indicating that the tail is responsible for shortening of T₂ times in full‐length PrP^C. The truncated fragment of PrP has facilitated the determination of meaningful order parameters (S²) from the relaxation data and shows for the first time that all three helices in PrP^C have similar rigidity. Slow conformational fluctuations of mouse PrP^C are localized to a distinct region that involves residues 171 and 172. Interestingly, residues 170–175 have been identified as a segment within PrP that will form a steric zipper, believed to be the fundamental amyloid unit. The flexibility within these residues could facilitate the PrP^C–PrP^Sc recognition process during fibril elongation.     

Domain motions of hyaluronan lyase underlying processive hyaluronan translocation

Hyaluronan lyase (Hyal) is a surface enzyme occurring in many bacterial organisms including members of Streptococcus species. Streptococcal Hyal primarily degrades hyaluronan‐substrate (HA) of the extracellular matrix. This degradation appears to facilitate the spread of this bacterium throughout host tissues. Unlike purely endolytic degradation of its other substrates, unsulfated chondroitin or some chondroitin sulfates, the degradation of HA by Hyal proceeds by processive exolytic cleavage of one disaccharide at a time following an initial endolytic cut. Molecular dynamics (MD) studies of Hyal from Streptococcus pneumoniae are presented that address the enzyme's molecular mechanism of action and the role of domain motions for processive functionality. The analysis of extensive sub‐microsecond MD simulations of this enzyme action on HA‐substrates of different lengths and the connection between the domain dynamics of Hyal and the translocation of the HA‐substrate reveals that opening/closing and twisting domain motions of the Hyal are intimately linked to processive HA degradation. Enforced simulations confirmed this finding as the domain motions in SpnHyal were found to be induced by enforced substrate translocation. These results establish the dynamic interplay between Hyal flexibility and substrate translocation and provide insight into the processive mechanism of Hyal. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Elbow‐ and hinge‐bending motions of IgG: Dielectric response and dynamic feature

Immunoglobulin G (IgG) is a Y‐shaped globular protein consisting of two Fab segments connecting to an Fc segment with a flexible hinge region, in which the Fab segments show secondary flexibility at an “elbow” region. In the present work, the hinge‐bending and elbow‐bending motions of aqueous solutions of IgG by microwave dielectric measurements below the freezing point of bulk water was observed. The presence of unfreezable water around the macromolecules reduced the effects of steric hindrance normally generated by ice and enabled the intramolecular motions of IgG. At the same time, the overall IgG molecule rotation was restricted by ice. Papain digestion and reduction of the disulfide linkage at the hinge region was used to generate Fab and Fc fragments. In solutions of these fragments, the dielectric relaxation process of the hinge‐bending motion was absent, although the elbow‐bending motion remained. Three relaxation processes were observed for papain‐digested IgG. The high, middle, and low frequency processes were attributed to unfrozen water, local peptide motions cooperating with bound water, and the elbow‐bending motion, respectively. In the case of the intact IgG, an additional relaxation process due to the hinge‐bending motion was observed at frequencies lower than that of the elbow‐bending motion. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 626–632, 2016.     

On the dielectrically observable consequences of the diffusional motions of lipids and proteins in membranes

A system consisting of any array of cylindrical, polytopic membrane proteins (or protein complexes) possessed of a permanent dipole moment and immersed in a closed, spherical phospholipid bilayer sheet is considered. It is assumed that rotation of the protein (complex) in a plane normal to the membrane, if occurring, is restricted by viscous drag alone. Lateral diffusion is assumed either to be free and random or to be partially constrained by barriers of an unspecified nature. The dielectric relaxation times calculated for membrane protein rotation in a suspension of vesicles of the above type are much longer than those observed with globular proteins in aqueous solution, and fall in the mid-to-high audio frequency range. If the long range lateral diffusion of (charged) membrane protein complexes is essentially unrestricted, as in the "fluid mosaic" membrane model, dielectric relaxation times for lateral motions will lie, except in the case of the very smallest vesicles, in the sub-audio (ELF) range. If, in contrast, the lateral diffusion of membrane protein complexes is partially restricted by "barriers" or "long-range" interactions (of unspecified nature), significant dielectric dispersions may be expected in both audio- and radio-frequency ranges, the critical (characteristic) frequencies depending upon the average distance moved before a barrier is encountered. Similar analyses are given for rotational and translational motions of phospholipids. At very low frequencies, a dispersion due to vesicle orientation might in principle also be observed; the dielectrically observable extent of this rotation will depend, inter alia, upon the charge mobility and disposition of the membrane protein complexes, as well as, of course, on the viscosity of the aqueous phase. The role of electroosmotic interactions between double layer ions (and water dipoles) and proteins raised above the membrane surface is considered. In some cases, it seems likely that such interactions serve to raise the dielectric increment, relative to that which might otherwise have been expected, of dispersions due to protein motions in membranes. Depending upon the tortuosity of the ion-relaxation pathways, such a relaxation mechanism might lead to almost any characteristic frequency, and, even in the absence of protein/lipid motions, would cause dielectric spectra to be much broader than one might expect from a simple, macroscopic treatment.     

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.