Slow Dynamics of Tryptophan-Water Networks in Proteins

Dynamics underlie the drug recognition mechanism by the efflux transporter EmrE

The multidrug efflux transporter EmrE from Escherichia coli requires anionic residues in the substrate binding pocket for coupling drug transport with the proton motive force. Here, we show how protonation of a single membrane embedded glutamate residue (Glu14) within the homodimer of EmrE modulates the structure and dynamics in an allosteric manner using NMR spectroscopy. The structure of EmrE in the Glu14 protonated state displays a partially occluded conformation that is inaccessible for drug binding by the presence of aromatic residues in the binding pocket. Deprotonation of a single Glu14 residue in one monomer induces an equilibrium shift toward the open state by altering its side chain position and that of a nearby tryptophan residue. This structural change promotes an open conformation that facilitates drug binding through a conformational selection mechanism and increases the binding affinity by approximately 2000-fold. The prevalence of proton-coupled exchange in efflux systems suggests a mechanism that may be shared in other antiporters where acid/base chemistry modulates access of drugs to the substrate binding pocket.

Heterogeneous and Allosteric Role of Surface Hydration for Protein-Ligand Binding

Atomistic-level understanding of surface hydration mediating protein-protein interactions and ligand binding has been a challenge due to the dynamic nature of water molecules near the surface. We develop a computational method to evaluate the solvation free energy based on the density map of the first hydration shell constructed from all-atom molecular dynamics simulation and use it to examine the binding of two intrinsically disordered ligands to their target protein domain. One ligand is from the human protein, and the other is from the 1918 Spanish flu virus. We find that the viral ligand incurs a 6.9 kcal/mol lower desolvation penalty upon binding to the target, which is consistent with its stronger binding affinity. The difference arises from the spatially fragmented and nonuniform water density profiles of the first hydration shell. In particular, residues that are distal from the ligand-binding site contribute to a varying extent to the desolvation penalty, among which the "entropy hotspot" residues contribute significantly. Thus, ligand binding alters hydration on remote sites in addition to affecting the binding interface. The nonlocal effect disappears when the conformational motion of the protein is suppressed. The present results elucidate the interplay between protein conformational dynamics and surface hydration. Our approach of measuring the solvation free energy based on the water density of the first hydration shell is tolerant of the conformational fluctuation of protein, and we expect it to be applicable to investigating a broad range of biomolecular interfaces.

Slow conformational changes in the rigid and highly stable chymotrypsin inhibitor 2

Slow conformational changes are often directly linked to protein function. It is however less clear how such processes may perturb the overall folding stability of a protein. We previously found that the stabilizing double mutant L49I/I57V in the small protein chymotrypsin inhibitor 2 from barley led to distributed increased nanosecond and faster dynamics. Here we asked what effects the L49I and I57V substitutions, either individually or together, have on the slow conformational dynamics of CI2. We used ¹⁵ N CPMG spin relaxation dispersion experiments to measure the kinetics, thermodynamics and structural changes associated with slow conformational change in CI2. These changes result in an excited state that is populated to 4.3% at 1 °C. As the temperature is increased the population of the excited state decreases. Structural changes in the excited state are associated with residues that interact with water molecules that have well defined positions and are found at these positions in all crystal structures of CI2. The substitutions in CI2 have only little effect on the structure of the excited state whereas the stability of the excited state to some extent follows the stability of the main state. The minor state is thus most populated for the most stable CI2 variant and least populated for the least stable variant. We hypothesize that the interactions between the substituted residues and the well-ordered water molecules links subtle structural changes around the substituted residues to the region in the protein that experience slow conformational changes. This article is protected by copyright. All rights reserved.

Understanding Conformational Properties and Role of Hydrogen Bonds in Glycosaminoglycans-Interleukin8 Complexes in Aqueous Medium by Molecular Dynamics Simulation

Atomistic molecular dynamics simulations were performed under ambient conditions to explore the conformational features and binding affinities of hexameric glycosaminoglycans (GAGs) with chemokine Interleukin8 (IL8) in an aqueous medium. We tried to understand the role of hydrogen bonds (HBs) involving conserved water in mediating the interactions. The Luzar-Chandler model was adopted to study the kinetics of HB breaking and formation concerning different water-mediated HBs. The conformational flexibilities of bound GAGs are due to the flexible glycosidic linkages than the occasional/rare ring pucker conformation. The free energy landscape constructed with ϕ, and ψ, depicted that different conformational minima associated with the glycosidic linkage flexibility of the GAGs in bound states are separated by energy barriers. The binding affinities of IL8 towards GAGs are favored through the electrostatic and non-polar solvation interactions. 4-different types of conserved water were explored in the solvent-mediated binding of GAGs with IL8. The average lifetime of the IL8-GAG direct HB pairs was ∼ten times less than the IL8-GAG-shared water HBs. This is due to the rapid establishment of HB breaking and reformation kinetics involving water of a shared layer. We find that despite the highly negatively charged surface of GAGs, the IL8 surface populated by non-cationic amino acids could serve as a promising binding site in addition to the cationic surface of the protein.

Environment-Driven Coherent Population Transfer Governs the Ultrafast Photophysics of Tryptophan

By combining UV transient absorption spectroscopy with sub-30-fs temporal resolution and CASPT2/MM calculations, we present a complete description of the primary photoinduced processes in solvated tryptophan. Our results shed new light on the role of the solvent in the relaxation dynamics of tryptophan. We unveil two consecutive coherent population transfer events involving the lowest two singlet excited states: a sub-50-fs nonadiabatic L_a → L_b transfer through a conical intersection and a subsequent 220 fs reverse L_b → L_a transfer due to solvent-assisted adiabatic stabilization of the L_a state. Vibrational fingerprints in the transient absorption spectra provide compelling evidence of a vibronic coherence established between the two excited states from the earliest times after photoexcitation and lasting until the back-transfer to L_a is complete. The demonstration of response to the environment as a driver of coherent population dynamics among the excited states of tryptophan closes the long debate on its solvent-assisted relaxation mechanisms and extends its application as a local probe of protein dynamics to the ultrafast time scales.

Roles of interfacial water states on advanced biomedical material design

When biomedical materials come into contact with body fluids, the first reaction that occurs on the material surface is hydration; proteins are then adsorbed and denatured on the hydrated material surface. The amount and degree of denaturation of adsorbed proteins affect subsequent cell behavior, including cell adhesion, migration, proliferation, and differentiation. Biomolecules are important for understanding the interactions and biological reactions of biomedical materials to elucidate the role of hydration in biomedical materials and their interaction partners. Analysis of the water states of hydrated materials is complicated and remains controversial; however, knowledge about interfacial water is useful for the design and development of advanced biomaterials. Herein, we summarize recent findings on the hydration of synthetic polymers, supramolecular materials, inorganic materials, proteins, and lipid membranes. Furthermore, we present recent advances in our understanding of the classification of interfacial water and advanced polymer biomaterials, based on the intermediate water concept.

Ultraslow Biological Water-Like Dynamics in Waterless Liquid Protein

Fluorescence correlation spectroscopy and time-dependent fluorescence Stokes shift have been employed to elucidate dynamics in different time scales, ranging from picoseconds to nanoseconds, for human serum albumin, in its native and cationized forms as well as in the self-assembled complex of the cationized protein with the polymer surfactant (PS) glycolic acid ethoxylate lauryl ether. The effect of crowding in this complex, especially in the waterless condition, is of prime importance in this context. Excellent correlation of the dynamics with the structures, obtained by circular dichroism and Fourier transform infrared spectroscopy, has been observed. Slow solvation, associated classically with biological water, has been observed in these systems, even in the waterless condition. This apparently intriguing observation has been rationalized by the relaxation of segments of the protein and the PS in the microenvironment of the fluorescent probe.

Similarities and Differences among Protein Dynamics Studied by Variable Temperature Nuclear Magnetic Resonance Relaxation

Understanding and describing the dynamics of proteins is one of the major challenges in biology. Here, we use multifield variable-temperature NMR longitudinal relaxation (R₁) measurements to determine the hierarchical activation energies of motions of four different proteins: two small globular proteins (GB1 and the SH3 domain of α-spectrin), an intrinsically disordered protein (the C-terminus of the nucleoprotein of the Sendai virus, Sendai Ntail), and an outer membrane protein (OmpG). The activation energies map the motions occurring in the side chains, in the backbone, and in the hydration shells of the proteins. We were able to identify similarities and differences in the average motions of the proteins. We find that the NMR relaxation properties of the four proteins do share similar features. The data characterizing average backbone motions are found to be very similar, the same for methyl group rotations, and similar activation energies are measured. The main observed difference occurs for the intrinsically disordered Sendai Ntail, where we observe much lower energy of activation for motions of protons associated with the protein-solvent interface as compared to the others. We also observe variability between the proteins regarding side chain ¹⁵N relaxation of lysine residues, with a higher activation energy observed in OmpG. This hints at strong interactions with negatively charged lipids in the bilayer and provides a possible mechanistic clue for the "positive-inside" rule for helical membrane proteins. Overall, these observations refine the understanding of the similarities and differences between hierarchical dynamics in proteins.

Structure-based mechanisms: On the way to apply alcohol dehydrogenases/reductases to organic-aqueous systems

Alcohol dehydrogenases/reductases catalyze enantioselective syntheses of versatile chiral compounds relying on direct hydride transfer from cofactor to substrates, or to an intermediate and then to substrates. Since most of the substrates catalyzed by alcohol dehydrogenases/reductases are insoluble in aqueous solutions, increasing interest has been turning to organic-aqueous systems. However, alcohol dehydrogenases/reductases are normally instable in organic solvents, leading to the unsatisfied enantioselective synthesis efficiency. The behaviors of these enzymes in organic solvents at an atomic level are unclear, thus it is of great importance to understand its structure-based mechanisms in organic-aqueous systems to improve their relative stability. Here, we summarized the accessible structures of alcohol dehydrogenases/reductases in Protein Data Bank crystallized in organic-aqueous systems, and compared the structures of alcohol dehydrogenases/reductases which have different tolerance towards organic solvents. By understanding the catalytic behaviors and mechanisms of these enzymes in organic-aqueous systems, the efficient enantioselective syntheses mediated by alcohol dehydrogenases/reductases and further challenges are also discussed through solvent engineering and enzyme-immobilization in the last decade.

Differential impact of BTK active site inhibitors on the conformational state of full-length BTK

Bruton’s tyrosine kinase (BTK) is targeted in the treatment of B-cell disorders including leukemias and lymphomas. Currently approved BTK inhibitors, including Ibrutinib, a first-in-class covalent inhibitor of BTK, bind directly to the kinase active site. While effective at blocking the catalytic activity of BTK, consequences of drug binding on the global conformation of full-length BTK are unknown. Here, we uncover a range of conformational effects in full-length BTK induced by a panel of active site inhibitors, including large-scale shifts in the conformational equilibria of the regulatory domains. Additionally, we find that a remote Ibrutinib resistance mutation, T316A in the BTK SH2 domain, drives spurious BTK activity by destabilizing the compact autoinhibitory conformation of full-length BTK, shifting the conformational ensemble away from the autoinhibited form. Future development of BTK inhibitors will need to consider long-range allosteric consequences of inhibitor binding, including the emerging application of these BTK inhibitors in treating COVID-19.

Treatments for blood cancers, such as leukemia and lymphoma, rely heavily on chemotherapy, using drugs that target a vulnerable aspect of the cancer cells. B-cells, a type of white blood cell that produces antibodies, require a protein called Bruton’s tyrosine kinase, or BTK for short, to survive. The drug ibrutinib (Imbruvica) is used to treat B-cell cancers by blocking BTK. The BTK protein consists of several regions. One of them, known as the kinase domain, is responsible for its activity as an enzyme (which allows it to modify other proteins by adding a ‘tag’ known as a phosphate group). The other regions of BTK, known as regulatory modules, control this activity. In BTK’s inactive form, the regulatory modules attach to the kinase domain, blocking the regulatory modules from interacting with other proteins. When BTK is activated, it changes its conformation so the regulatory regions detach and become available for interactions with other proteins, at the same time exposing the active kinase domain.

Ibrutinib and other BTK drugs in development bind to the kinase domain to block its activity. However, it is not known how this binding affects the regulatory modules. Previous efforts to study how drugs bind to BTK have used a version of the protein that only had the kinase domain, instead of the full-length protein.

Now, Joseph et al. have studied full-length BTK and how it binds to five different drugs. The results reveal that ibrutinib and another drug called dasatinib both indirectly disrupt the normal position of the regulatory domains pushing BTK toward a conformation that resembles the activated state. By contrast, the three other compounds studied do not affect the inactive structure. Joseph et al. also examined a mutation in BTK that confers resistance against ibrutinib. This mutation increases the activity of BTK by disrupting the inactive structure, leading to B cells surviving better.

Understanding how drug resistance mechanisms can work will lead to better drug treatment strategies for cancer. BTK is also a target in other diseases such as allergies or asthma and even COVID-19. If interactions between partner proteins and the regulatory domain are important in these diseases, then they may be better treated with drugs that maintain the regulatory modules in their inactive state. This research will help to design drugs that are better able to control BTK activity.

Phonon-mediated lipid raft formation in biological membranes

Short-wavelength collective molecular motions, also known as phonons, have recently attracted much interest in revealing dynamic properties of biological membranes through the use of neutron and X-ray scattering, infrared and Raman spectroscopies, and molecular dynamics simulations. Experimentally detecting unique vibrational patterns such as, shear phonon excitations, viscoelastic crossovers, transverse acoustic phonon gaps, and continuous and truncated optical phonon modes in cellular membranes, to name a few, has proven non-trivial. Here, we review recent advances in liquid thermodynamics that have resulted in the development of the phonon theory of liquids. The theory has important predictions regarding the shear vibrational spectra of fluids, namely the emergence of viscoelastic crossovers and transverse acoustic phonon gaps. Furthermore, we show that these vibrational patterns are common in soft (non-crystalline) materials, including, but not limited to liquids, colloids, liquid crystals (mesogens), block copolymers, and biological membranes. The existence of viscoelastic crossovers and acoustic phonon gaps define the self-diffusion properties of cellular membranes and provide a molecular picture of the transient nature of lipid rafts (Bolmatov et al., 2020). Importantly, the timescales (picoseconds) for the formation and dissolution of transient lipid rafts match the lifetime of the formation and breakdown of interfacial water hydrogen bonds. Apart from acoustic propagating phonon modes, biological membranes can also support more energetic non-propagating optical phonon excitations, also known as standing waves or breathing modes. Importantly, optical phonons can be truncated due to the existence of finite size nanodomains made up of strongly correlated lipid-cholesterol molecular pairs. These strongly coupled molecular pairs can serve as nucleation centers for the formation of stable rafts at larger length scales, due to correlations of spontaneous fluctuations (Onsager's regression hypothesis). Finally and importantly, molecular level viscoelastic crossovers, acoustic phonon gaps, and continuous and truncated optical phonon modes may offer insights as to how lipid-lipid and lipid-protein interactions enable biological function.

Migrations and Catalytic Action of Water Molecules in the Ionized Formamide-(H2O)2 Cluster

Isomerization dynamics involving the migrations, proton transfer reaction, and catalytic actions of water molecules upon vertical ionization of the formamide (FA)-(H₂O)₂ cluster is investigated by the infrared spectroscopy and theoretical reaction path search calculation. The infrared spectroscopic result indicates the [FA-(H₂O)₂]⁺ cation has the hydrogen-bonded structure of the enol isomer cation of formamide and the water dimer. This structure is formed by proton transfer from the CH bond to the carbonyl group through the catalytic action of the water molecules. The isomerization paths involving this enolization in ionized FA-(H₂O)₂ are explored by using the anharmonic downward distortion following method. We found multiple enolization paths which accompany proton exchanges among the formamide moiety and water molecules through the catalytic actions of the water molecules.

Self-Assembly of Aromatic Amino Acid Enantiomers into Supramolecular Materials of High Rigidity

Most natural biomolecules may exist in either of two enantiomeric forms. Although in nature, amino acid biopolymers are characterized by l-type homochirality, incorporation of d-amino acids in the design of self-assembling peptide motifs has been shown to significantly alter enzyme stability, conformation, self-assembly behavior, cytotoxicity, and even therapeutic activity. However, while functional metabolite assemblies are ubiquitous throughout nature and play numerous important roles including physiological, structural, or catalytic functions, the effect of chirality on the self-assembly nature and function of single amino acids is not yet explored. Herein, we investigated the self-assembly mechanism of amyloid-like structure formation by two aromatic amino acids, phenylalanine (Phe) and tryptophan (Trp), both previously found as extremely important for the nucleation and self-assembly of aggregation-prone peptide regions into functional structures. Employing d-enantiomers, we demonstrate the critical role that amino acid chirality plays in their self-assembly process. The kinetics and morphology of pure enantiomers is completely altered upon their coassembly, allowing to fabricate different nanostructures that are mechanically more robust. Using diverse experimental techniques, we reveal the different molecular arrangement and self-assembly mechanism of the dl-racemic mixtures that resulted in the formation of advanced supramolecular materials. This study provides a simple yet sophisticated engineering model for the fabrication of attractive materials with bionanotechnological applications.

Solid-State NMR Reveals the Structural Transformation of the TDP-43 Amyloidogenic Region upon Fibrillation

TDP-43 is a primary pathological hallmark protein of amyotrophic lateral sclerosis and frontotemporal lobar degeneration, which may exist in the form of amyloid inclusions in the cells of patients. In addition to serving as a biomarker for these diseases, TDP-43 can also directly trigger neurodegeneration. We previously determined the amyloidogenic core region of TDP-43 (residues 311-360) and showed by solution NMR that this region includes two α-helices [(321-330) and (335-343)] in solution. We suggested that the helix-to-sheet structural transformation initiates TDP-43 aggregation. In the present study, X-ray diffraction shows that TDP-43 (311-360) aggregates adopt a cross-β structure. Thioredoxin (Trx)-fused TDP-43 (311-360) can undergo liquid-liquid phase separation (LLPS) before fibrillation, suggesting that phase separation is an intermediate step before amyloid formation. Solid-state NMR (SSNMR), carried out to elucidate the structural changes of TDP-43 (311-360) at the atomic level, indicates five β-strands of the amyloids formed, with the major two β-strands contributed by the first helical region in the solution structure. The NMR evidence is also in support of the fibril having a parallel in-register conformation, implying a mechanism in which the helix-helix interactions in LLPS are converted into β-strand parallel lateral association upon fibrillation. Our studies have assigned many key interresidue interactions that contribute to the stability of the fibril, including F316 with I318 and Q327 and W334 with A325, A326, A329, and S332. SSNMR with ¹H detection reveals a unique close interaction between the indole Nε1-Hε1 of W334 and the side-chain carbonyl of Q343. This interaction could be a very important factor in initiating TDP-43 (311-360) folding/misfolding in LLPS.

Temperature-jump solution X-ray scattering reveals distinct motions in a dynamic enzyme

Correlated motions of proteins are critical to function, but these features are difficult to resolve using traditional structure determination techniques. Time-resolved X-ray methods hold promise for addressing this challenge, but have relied on the exploitation of exotic protein photoactivity, and are therefore not generalizable. Temperature jumps, through thermal excitation of the solvent, have been utilized to study protein dynamics using spectroscopic techniques, but their implementation in X-ray scattering experiments has been limited. Here, we perform temperature-jump small- and wide-angle X-ray scattering measurements on a dynamic enzyme, cyclophilin A, demonstrating that these experiments are able to capture functional intramolecular protein dynamics on the microsecond timescale. We show that cyclophilin A displays rich dynamics following a temperature jump, and use the resulting time-resolved signal to assess the kinetics of conformational changes. Two relaxation processes are resolved: a fast process is related to surface loop motions, and a slower process is related to motions in the core of the protein that are critical for catalytic turnover.

Time-resolved crystallography reveals allosteric communication aligned with molecular breathing

A comprehensive understanding of protein function demands correlating structure and dynamic changes. Using time-resolved serial synchrotron crystallography, we visualized half-of-the-sites reactivity and correlated molecular-breathing motions in the enzyme fluoroacetate dehalogenase. Eighteen time points from 30 milliseconds to 30 seconds cover four turnover cycles of the irreversible reaction. They reveal sequential substrate binding, covalent-intermediate formation, setup of a hydrolytic water molecule, and product release. Small structural changes of the protein mold and variations in the number and placement of water molecules accompany the various chemical steps of catalysis. Triggered by enzyme-ligand interactions, these repetitive changes in the protein framework’s dynamics and entropy constitute crucial components of the catalytic machinery.

Water-Mediated Selenium Hydrogen-Bonding in Proteins: PDB Analysis and Gas-Phase Spectroscopy of Model Complexes

High-resolution X-ray crystallography and two-dimensional NMR studies demonstrate that water-mediated conventional hydrogen-bonding interactions (N-H···N, O-H···N, etc.) bridging two or more amino acid residues contribute to the stability of proteins and protein-ligand complexes. In this work, we have investigated single water-mediated selenium hydrogen-bonding interactions (unconventional hydrogen-bonding) between amino acid residues in proteins through extensive protein data bank (PDB) analysis coupled with gas-phase spectroscopy and quantum chemical calculation of a model complex consisting of indole, dimethyl selenide, and water. Here, indole and dimethyl selenide represent the amino acid residues tryptophan and selenomethionine, respectively. The current investigation demonstrates that the most stable structure of the model complex observed in the IR spectroscopy mimics single water-mediated selenium hydrogen-bonded structural motifs present in the crystal structures of proteins. The present work establishes that water-mediated Se hydrogen-bonding interactions are ubiquitous in proteins and the number of these interactions observed in the PDB is more than that of direct Se hydrogen-bonds present there.

Ultrafast energy relaxation dynamics of amide I vibrations coupled with protein-bound water molecules

The influence of hydration water on the vibrational energy relaxation in a protein holds the key to understand ultrafast protein dynamics, but its detection is a major challenge. Here, we report measurements on the ultrafast vibrational dynamics of amide I vibrations of proteins at the lipid membrane/H₂O interface using femtosecond time-resolved sum frequency generation vibrational spectroscopy. We find that the relaxation time of the amide I mode shows a very strong dependence on the H₂O exposure, but not on the D₂O exposure. This observation indicates that the exposure of amide I bond to H₂O opens up a resonant relaxation channel and facilitates direct resonant vibrational energy transfer from the amide I mode to the H₂O bending mode. The protein backbone motions can thus be energetically coupled with protein-bound water molecules. Our findings highlight the influence of H₂O on the ultrafast structure dynamics of proteins.

Functional Hydration Behavior: Interrelation between Hydration and Molecular Properties at Lipid Membrane Interfaces

Water is an abundant commodity and has various important functions. It stabilizes the structure of biological macromolecules, controls biochemical activities, and regulates interfacial/intermolecular interactions. Common aspects of interfacial water can be obtained by overviewing fundamental functions and properties at different temporal and spatial scales. It is important to understand the hydrogen bonding and structural properties of water and to evaluate the individual molecular species having different hydration properties. Water molecules form hydrogen bonds with biomolecules and contribute to the adjustment of their properties, such as surface charge, hydrophilicity, and structural flexibility. In this review, the fundamental properties of water molecules and the methods used for the analyses of water dynamics are summarized. In particular, the interrelation between the hydration properties, determined by molecules, and the properties of molecules, determined by their hydration properties, are discussed using the lipid membrane as an example. Accordingly, interesting water functions are introduced that provide beneficial information in the fields of biochemistry, medicine, and food chemistry.