Terahertz underdamped vibrational motion governs protein-ligand binding in solution

Low-frequency collective vibrational modes in proteins have been proposed as being responsible for efficiently directing biochemical reactions and biological energy transport. However, evidence of the existence of delocalized vibrational modes is scarce and proof of their involvement in biological function absent. Here we apply extremely sensitive femtosecond optical Kerr-effect spectroscopy to study the depolarized Raman spectra of lysozyme and its complex with the inhibitor triacetylchitotriose in solution. Underdamped delocalized vibrational modes in the terahertz frequency domain are identified and shown to blue-shift and strengthen upon inhibitor binding. This demonstrates that the ligand-binding coordinate in proteins is underdamped and not simply solvent-controlled as previously assumed. The presence of such underdamped delocalized modes in proteins may have significant implications for the understanding of the efficiency of ligand binding and protein–molecule interactions, and has wider implications for biochemical reactivity and biological function

Predominant Functional Expression of Kv1.3 by Activated Microglia of the Hippocampus after Status epilepticus

BACKGROUND:Growing evidence indicates that the functional state of microglial cells differs according to the pathological conditions that trigger their activation. In particular, activated microglial cells can express sets of Kv subunits which sustain delayed rectifying potassium currents (Kdr) and modulate differently microglia proliferation and ability to release mediators. We recently reported that hippocampal microglia is in a particular activation state after a status epilepticus (SE) and the present study aimed at identifying which of the Kv channels are functionally expressed by microglia in this model. METHODOLOGY/PRINCIPAL FINDINGS:SE was induced by systemic injection of kainate in CX3CR1(eGFP/+) mice and whole cell recordings of fluorescent microglia were performed in acute hippocampal slices prepared 48 h after SE. Microglia expressed Kdr currents which were characterized by a potential of half-maximal activation near -25 mV, prominent steady-state and cumulative inactivations. Kdr currents were almost abolished by the broad spectrum antagonist 4-Aminopyridine (1 mM). In contrast, tetraethylammonium (TEA) at a concentration of 1 mM, known to block Kv3.1, Kv1.1 and 1.2 subunits, only weakly reduced Kdr currents. However, at a concentration of 5 mM which should also affect Kv1.3 and 1.6, TEA inhibited about 30% of the Kdr conductance. Alpha-dendrotoxin, which selectively inhibits Kv1.1, 1.2 and 1.6, reduced only weakly Kdr currents, indicating that channels formed by homomeric assemblies of these subunits are not important contributors of Kdr currents. Finally, agitoxin-2 and margatoxin strongly inhibited the current. CONCLUSIONS/SIGNIFICANCE:These results indicate that Kv1.3 containing channels predominantly determined Kdr currents in activated microglia after SE

Coupled Motions Direct Electrons along Human Microsomal P450 Chains

Directional electron transfer through biological redox chains can be achieved by coupling reaction chemistry to conformational changes in individual redox enzymes

Enzyme Promiscuity in Enolase Superfamily. Theoretical Study of o-Succinylbenzoate Synthase Using QM/MM Methods

The promiscuous activity of the enzyme o-succinylbenzoate synthase (OSBS) from the actinobacteria Amycolatopsis is investigated by means of QM/MM methods, using both density functional theory and semiempirical Hamiltonians. This enzyme catalyzes not only the dehydration of 2-succinyl-6R-hydroxy-2,4-cyclohexadiene-1R-carboxylate but also catalyzes racemization of different acylamino acids, with N-succinyl-R-phenylglycine being the best substrate. We investigated the molecular mechanisms for both reactions exploring the potential energy surface. Then, molecular dynamics simulations were performed to obtain the free energy profiles and the averaged interaction energies of enzymatic residues with the reacting system. Our results confirm the plausibility of the reaction mechanisms proposed in the literature, with a good agreement between theoretical and experimentally derived activation free energies. Our simulations unravel the role played by the different residues in each of the two possible reactions. The presence of flexible loops in the active site and the selection of structural modifications in the substrate seem to be key elements to promote the promiscuity of this enzyme.This work was supported by the Spanish Ministerio de Economia y Competitividad project CTQ2012-36253-C03-03 ́ and FEDER funds. K.S. thanks the Polish National Science Center (NCN) for Grant 2011/02/A/ST4/00246. The authors acknowledge computational facilities of the Servei d’Informatica ̀ de la Universitat de Valencia in the ̀ “Tirant” supercomputer, which is part of the Spanish Supercomputing Network

Global Conformational Dynamics of a Y-Family DNA Polymerase during Catalysis

High-resolution analysis of protein, and DNA conformational changes during DNA polymerization, established relationships between the enzymatic function and conformational dynamics of individual domains for a DNA polymerase

Evolutionarily Conserved Linkage between Enzyme Fold, Flexibility, and Catalysis

Proteins are intrinsically flexible molecules. The role of internal motions in a protein's designated function is widely debated. The role of protein structure in enzyme catalysis is well established, and conservation of structural features provides vital clues to their role in function. Recently, it has been proposed that the protein function may involve multiple conformations: the observed deviations are not random thermodynamic fluctuations; rather, flexibility may be closely linked to protein function, including enzyme catalysis. We hypothesize that the argument of conservation of important structural features can also be extended to identification of protein flexibility in interconnection with enzyme function. Three classes of enzymes (prolyl-peptidyl isomerase, oxidoreductase, and nuclease) that catalyze diverse chemical reactions have been examined using detailed computational modeling. For each class, the identification and characterization of the internal protein motions coupled to the chemical step in enzyme mechanisms in multiple species show identical enzyme conformational fluctuations. In addition to the active-site residues, motions of protein surface loop regions (>10 Å away) are observed to be identical across species, and networks of conserved interactions/residues connect these highly flexible surface regions to the active-site residues that make direct contact with substrates. More interestingly, examination of reaction-coupled motions in non-homologous enzyme systems (with no structural or sequence similarity) that catalyze the same biochemical reaction shows motions that induce remarkably similar changes in the enzyme–substrate interactions during catalysis. The results indicate that the reaction-coupled flexibility is a conserved aspect of the enzyme molecular architecture. Protein motions in distal areas of homologous and non-homologous enzyme systems mediate similar changes in the active-site enzyme–substrate interactions, thereby impacting the mechanism of catalyzed chemistry. These results have implications for understanding the mechanism of allostery, and for protein engineering and drug design

A Rigidifying Salt-Bridge Favors the Activity of Thermophilic Enzyme at High Temperatures at the Expense of Low-Temperature Activity

Although enzymes from thermophiles thriving in hot habitats are more stable than their mesophilic homologs, they are often less active at low temperatures. One theory suggests that extra stabilizing interactions found in thermophilic enzymes may increase their rigidity and decrease enzymatic activity at lower temperatures. We used acylphosphatase as a model to study how flexibility affects enzymatic activity. This enzyme has a unique structural feature in that an invariant arginine residue, which takes part in catalysis, is restrained by a salt-bridge in the thermophilic homologs but not in its mesophilic homologs. Here, we demonstrate the trade-offs between flexibility and enzymatic activity by disrupting the salt-bridge in a thermophilic acylphosphatase and introducing it in the mesophilic human homolog. Our results suggest that the salt-bridge is a structural adaptation for thermophilic acylphosphatases as it entropically favors enzymatic activity at high temperatures by restricting the flexibility of the active-site residue. However, at low temperatures the salt-bridge reduces the enzymatic activity because of a steeper temperature-dependency of activity

Investigating Homology between Proteins using Energetic Profiles

Accumulated experimental observations demonstrate that protein stability is often preserved upon conservative point mutation. In contrast, less is known about the effects of large sequence or structure changes on the stability of a particular fold. Almost completely unknown is the degree to which stability of different regions of a protein is generally preserved throughout evolution. In this work, these questions are addressed through thermodynamic analysis of a large representative sample of protein fold space based on remote, yet accepted, homology. More than 3,000 proteins were computationally analyzed using the structural-thermodynamic algorithm COREX/BEST. Estimated position-specific stability (i.e., local Gibbs free energy of folding) and its component enthalpy and entropy were quantitatively compared between all proteins in the sample according to all-vs.-all pairwise structural alignment. It was discovered that the local stabilities of homologous pairs were significantly more correlated than those of non-homologous pairs, indicating that local stability was indeed generally conserved throughout evolution. However, the position-specific enthalpy and entropy underlying stability were less correlated, suggesting that the overall regional stability of a protein was more important than the thermodynamic mechanism utilized to achieve that stability. Finally, two different types of statistically exceptional evolutionary structure-thermodynamic relationships were noted. First, many homologous proteins contained regions of similar thermodynamics despite localized structure change, suggesting a thermodynamic mechanism enabling evolutionary fold change. Second, some homologous proteins with extremely similar structures nonetheless exhibited different local stabilities, a phenomenon previously observed experimentally in this laboratory. These two observations, in conjunction with the principal conclusion that homologous proteins generally conserved local stability, may provide guidance for a future thermodynamically informed classification of protein homology

Corresponding Functional Dynamics across the Hsp90 Chaperone Family: Insights from a Multiscale Analysis of MD Simulations

Understanding how local protein modifications, such as binding small-molecule ligands, can trigger and regulate large-scale motions of large protein domains is a major open issue in molecular biology. We address various aspects of this problem by analyzing and comparing atomistic simulations of Hsp90 family representatives for which crystal structures of the full length protein are available: mammalian Grp94, yeast Hsp90 and E.coli HtpG. These chaperones are studied in complex with the natural ligands ATP, ADP and in the Apo state. Common key aspects of their functional dynamics are elucidated with a novel multi-scale comparison of their internal dynamics. Starting from the atomic resolution investigation of internal fluctuations and geometric strain patterns, a novel analysis of domain dynamics is developed. The results reveal that the ligand-dependent structural modulations mostly consist of relative rigid-like movements of a limited number of quasi-rigid domains, shared by the three proteins. Two common primary hinges for such movements are identified. The first hinge, whose functional role has been demonstrated by several experimental approaches, is located at the boundary between the N-terminal and Middle-domains. The second hinge is located at the end of a three-helix bundle in the Middle-domain and unfolds/unpacks going from the ATP- to the ADP-state. This latter site could represent a promising novel druggable allosteric site common to all chaperones