Protein dynamics investigated by the neutron diffraction-hydrogen exchange technique

Interpretation of Hydrogen/Deuterium Exchange Mass Spectrometry

This paper sheds light on the meaning of hydrogen/deuterium exchange-mass spectrometry (HDX-MS) data. HDX-MS data provide not structural information but dynamic information on an analyte protein. First, the reaction mechanism of backbone amide HDX reaction is considered and the correlation between the parameters from an X-ray crystal structure and the protection factors of HDX reactions of cytochrome c is evaluated. The presence of H-bonds in a protein structure has a strong influence on HDX rates which represent protein dynamics, while the solvent accessibility only weakly affects the HDX rates. Second, the energy diagrams of the HDX reaction at each residue in the presence and absence of perturbation are described. Whereas the free energy change upon mutation can be directly measured by the HDX rates, the free energy change upon ligand binding may be complicated due to the presence of unbound analyte protein in the protein-ligand mixture. Third, the meanings of HDX and other biophysical techniques are explained using a hypothetical protein folding well. The shape of the protein folding well describes the protein dynamics and provides Boltzmann distribution of open and closed states which yield HDX protection factors, while a protein's crystal structure represents a snapshot near the bottom of the well. All biophysical data should be consistent yet provide different information because they monitor different parts of the same protein folding well.

Unveiling the complexity: assessing models describing the structure and function of the nuclear pore complex

The nuclear pore complex (NPC) serves as a pivotal subcellular structure, acting as a gateway that orchestrates nucleocytoplasmic transport through a selectively permeable barrier. Nucleoporins (Nups), particularly those containing phenylalanine-glycine (FG) motifs, play indispensable roles within this barrier. Recent advancements in technology have significantly deepened our understanding of the NPC's architecture and operational intricacies, owing to comprehensive investigations. Nevertheless, the conspicuous presence of intrinsically disordered regions within FG-Nups continues to present a formidable challenge to conventional static characterization techniques. Historically, a multitude of strategies have been employed to unravel the intricate organization and behavior of FG-Nups within the NPC. These endeavors have given rise to multiple models that strive to elucidate the structural layout and functional significance of FG-Nups. Within this exhaustive review, we present a comprehensive overview of these prominent models, underscoring their proposed dynamic and structural attributes, supported by pertinent research. Through a comparative analysis, we endeavor to shed light on the distinct characteristics and contributions inherent in each model. Simultaneously, it remains crucial to acknowledge the scarcity of unequivocal validation for any of these models, as substantiated by empirical evidence.

Description of peptide bond planarity from high-resolution neutron crystallography

Neutron crystallography is a highly effective method for visualizing hydrogen atoms in proteins. In our recent study, we successfully determined the high-resolution (1.2 Å) neutron structure of high-potential iron-sulfur protein, refining the coordinates of some amide protons without any geometric restraints. Interestingly, we observed that amide protons are deviated from the peptide plane due to electrostatic interactions. Moreover, the difference in the position of the amide proton of Cys75 between reduced and oxidized states is possibly attributed to the electron storage capacity of the iron-sulfur cluster. Additionally, we have discussed about the rigidity of the iron-sulfur cluster based on the results of the hydrogen-deuterium exchange. Our research underscores the significance of neutron crystallography in protein structure elucidation, enriching our understanding of protein functions at an atomic resolution.

Hydrogen Deuterium Exchange and other Mass Spectrometry-based Approaches for Epitope Mapping

Antigen-antibody interactions are a fundamental subset of protein-protein interactions responsible for the "survival of the fittest". Determining the interacting interface of the antigen, called an epitope, and that on the antibody, called a paratope, is crucial to antibody development. Because each antigen presents multiple epitopes (unique footprints), sophisticated approaches are required to determine the target region for a given antibody. Although X-ray crystallography, Cryo-EM, and nuclear magnetic resonance can provide atomic details of an epitope, they are often laborious, poor in throughput, and insensitive. Mass spectrometry-based approaches offer rapid turnaround, intermediate structural resolution, and virtually no size limit for the antigen, making them a vital approach for epitope mapping. In this review, we describe in detail the principles of hydrogen deuterium exchange mass spectrometry in application to epitope mapping. We also show that a combination of MS-based approaches can assist or complement epitope mapping and push the limit of structural resolution to the residue level. We describe in detail the MS methods used in epitope mapping, provide our perspective about the approaches, and focus on elucidating the role that HDX-MS is playing now and in the future by organizing a discussion centered around several improvements in prototype instrument/applications used for epitope mapping. At the end, we provide a tabular summary of the current literature on HDX-MS-based epitope mapping.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and x-ray diffraction at room temperature

The nonstructural protein 3 (NSP3) macrodomain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Mac1) removes adenosine diphosphate (ADP) ribosylation posttranslational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the coronavirus disease 2019 pandemic. Here, we determined neutron and x-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase the potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a reevaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and X-ray diffraction at room temperature

The NSP3 macrodomain of SARS CoV 2 (Mac1) removes ADP-ribosylation post-translational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the COVID-19 pandemic. Here, we determined neutron and X-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site, and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a re-evaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

Hydrogen/Deuterium Exchange Behavior During Denaturing/Refolding Processes Determined in Tetragonal Hen Egg-White Lysozyme Crystals

The hydrogen/deuterium (H/D) exchange of main-chain amide hydrogens in the protein that denatured and refolded in deuterated solvent is considered to contain the traces of hydrogen bond cleavages or the exposure to solvent of the buried part of the protein during the denaturing and refolding (denaturing/refolding) processes. Here, we report the H/D exchange behaviors in hen egg-white lysozymes denatured under acidic conditions, basic conditions, and thermal conditions and then refolded in deuterated solvents, using crystallographic methods. The results indicate that the space containing the Trp28 side chain was hardly exposed to the solvent in acidic conditions, but exposed under basic or heated conditions. Moreover, the β-bridges between Tyr53 and Ile58 in strands β2 and β3, which are in a highly conserved region, show some tolerance to changes in pD. The results indicate that crystallographic method is one of the powerful tools to analyze the denaturing/refolding processes of proteins.

Capturing the Catalytic Proton of Dihydrofolate Reductase: Implications for General Acid-Base Catalysis

Acid–base catalysis, which involves one or more proton transfer reactions, is a chemical mechanism commonly employed by many enzymes. The molecular basis for catalysis is often derived from structures determined at the optimal pH for enzyme activity. However, direct observation of protons from experimental structures is quite difficult; thus, a complete mechanistic description for most enzymes remains lacking. Dihydrofolate reductase (DHFR) exemplifies general acid–base catalysis, requiring hydride transfer and protonation of its substrate, DHF, to form the product, tetrahydrofolate (THF). Previous X-ray and neutron crystal structures coupled with theoretical calculations have proposed that solvent mediates the protonation step. However, visualization of a proton transfer has been elusive. Based on a 2.1 Å resolution neutron structure of a pseudo-Michaelis complex of E. coli DHFR determined at acidic pH, we report the direct observation of the catalytic proton and its parent solvent molecule. Comparison of X-ray and neutron structures elucidated at acidic and neutral pH reveals dampened dynamics at acidic pH, even for the regulatory Met20 loop. Guided by the structures and calculations, we propose a mechanism where dynamics are crucial for solvent entry and protonation of substrate. This mechanism invokes the release of a sole proton from a hydronium (H₃O⁺) ion, its pathway through a narrow channel that sterically hinders the passage of water, and the ultimate protonation of DHF at the N5 atom.

Hydrogen/deuterium exchange behavior in tetragonal hen egg-white lysozyme crystals affected by solution state

Neutron diffraction studies of hydrogen/deuterium-exchanged hen egg-white lysozyme were performed by a joint X-ray and neutron refinement to elucidate the hydrogen/deuterium exchange behavior. Large crystals for neutron work, consisting of molecules that were exchanged before crystallization, were obtained by repeatedly adding protein solution to the crystal batch using deuterated precipitant reagent. There are differences in hydrogen/deuterium exchange behavior compared with previous crystallographic or NMR studies, which could be due to intermolecular interactions in the crystal or to different lengths of exchange period.

Mass Spectrometry-Based Protein Footprinting for Higher-Order Structure Analysis: Fundamentals and Applications

Proteins adopt different higher-order structures (HOS) to enable their unique biological functions. Understanding the complexities of protein higher-order structures and dynamics requires integrated approaches, where mass spectrometry (MS) is now positioned to play a key role. One of those approaches is protein footprinting. Although the initial demonstration of footprinting was for the HOS determination of protein/nucleic acid binding, the concept was later adapted to MS-based protein HOS analysis, through which different covalent labeling approaches "mark" the solvent accessible surface area (SASA) of proteins to reflect protein HOS. Hydrogen-deuterium exchange (HDX), where deuterium in D2O replaces hydrogen of the backbone amides, is the most common example of footprinting. Its advantage is that the footprint reflects SASA and hydrogen bonding, whereas one drawback is the labeling is reversible. Another example of footprinting is slow irreversible labeling of functional groups on amino acid side chains by targeted reagents with high specificity, probing structural changes at selected sites. A third footprinting approach is by reactions with fast, irreversible labeling species that are highly reactive and footprint broadly several amino acid residue side chains on the time scale of submilliseconds. All of these covalent labeling approaches combine to constitute a problem-solving toolbox that enables mass spectrometry as a valuable tool for HOS elucidation. As there has been a growing need for MS-based protein footprinting in both academia and industry owing to its high throughput capability, prompt availability, and high spatial resolution, we present a summary of the history, descriptions, principles, mechanisms, and applications of these covalent labeling approaches. Moreover, their applications are highlighted according to the biological questions they can answer. This review is intended as a tutorial for MS-based protein HOS elucidation and as a reference for investigators seeking a MS-based tool to address structural questions in protein science.

High-resolution neutron crystallography visualizes an OH-bound resting state of a copper-containing nitrite reductase

X-ray crystallography often fails to determine the positions of hydrogen atoms, which play crucial roles in enzymatic reactions. Despite many X-ray crystallographic studies, the reaction mechanism of copper-containing nitrite reductases (CuNIRs), which reduce nitrite using two protons, has been controversial. The high-resolution neutron structure of a CuNIR reveals the protonation states of catalytic residues and key water molecules, thus providing insights into the catalytic mechanism. The catalytic Cu is shown to be coordinated by a hydroxide ion and not water. Furthermore, the hydrogen-deuterium exchange ratio suggests that intramolecular electron transfer is involved in a hydrogen-bond jump. These observations are consistent with previous computational chemistry; therefore, our study forms a bridge between the structural biology and quantum chemistry of CuNIRs.

Copper-containing nitrite reductases (CuNIRs) transform nitrite to gaseous nitric oxide, which is a key process in the global nitrogen cycle. The catalytic mechanism has been extensively studied to ultimately achieve rational control of this important geobiochemical reaction. However, accumulated structural biology data show discrepancies with spectroscopic and computational studies; hence, the reaction mechanism is still controversial. In particular, the details of the proton transfer involved in it are largely unknown. This situation arises from the failure of determining positions of hydrogen atoms and protons, which play essential roles at the catalytic site of CuNIRs, even with atomic resolution X-ray crystallography. Here, we determined the 1.50 Å resolution neutron structure of a CuNIR from Geobacillus thermodenitrificans (trimer molecular mass of ∼106 kDa) in its resting state at low pH. Our neutron structure reveals the protonation states of catalytic residues (deprotonated aspartate and protonated histidine), thus providing insights into the catalytic mechanism. We found that a hydroxide ion can exist as a ligand to the catalytic Cu atom in the resting state even at a low pH. This OH-bound Cu site is unexpected from previously given X-ray structures but consistent with a reaction intermediate suggested by computational chemistry. Furthermore, the hydrogen-deuterium exchange ratio in our neutron structure suggests that the intramolecular electron transfer pathway has a hydrogen-bond jump, which is proposed by quantum chemistry. Our study can seamlessly link the structural biology to the computational chemistry of CuNIRs, boosting our understanding of the enzymes at the atomic and electronic levels.

Visualization of H atoms in the X-ray crystal structure of photoactive yellow protein: Does it contain low-barrier hydrogen bonds?

The hydrogen bond (HB) between 4-hydroxycinnamic acid (HC4) and glutamic acid E46 of photoactive yellow protein is exceptionally strong. In the 0.82-å resolution X-ray structure for this protein (PDB ID: 1NWZ), the OH…O distance is only 2.57 å. The position of the H atom between these two O atoms has not been determined in that structure, and in the absence of that information, it is impossible to determine whether or not this HB is a low-barrier HB (LBHB), as was proposed recently based on neutron structures of this protein (Yamaguchi et al., Proceedings of the National Academy of Sciences of the United States of America, 2009, 106: 440-444). Residual electron density maps computed using the 1NWZ data reveal that this H atom is 0.92 å from the O_ε2 atom of E46 and 1.67 å from the O₄ ' of HC4, and that the OH…O bond angle is 167°. These observations indicate that E46 is protonated, and HC4 is deprotonated, as was originally suggested, and that the HB in question is not an LBHB.

Computational methods and challenges in hydrogen/deuterium exchange mass spectrometry

Hydrogen/Deuterium exchange (HDX) has been applied, since the 1930s, as an analytical tool to study the structure and dynamics of (small) biomolecules. The popularity of using HDX to study proteins increased drastically in the last two decades due to the successful combination with mass spectrometry (MS). Together with this growth in popularity, several technological advances have been made, such as improved quenching and fragmentation. As a consequence of these experimental improvements and the increased use of protein-HDXMS, large amounts of complex data are generated, which require appropriate analysis. Computational analysis of HDXMS requires several steps. A typical workflow for proteins consists of identification of (non-)deuterated peptides or fragments of the protein under study (local analysis), or identification of the deuterated protein as a whole (global analysis); determination of the deuteration level; estimation of the protection extent or exchange rates of the labile backbone amide hydrogen atoms; and a statistically sound interpretation of the estimated protection extent or exchange rates. Several algorithms, specifically designed for HDX analysis, have been proposed. They range from procedures that focus on one specific step in the analysis of HDX data to complete HDX workflow analysis tools. In this review, we provide an overview of the computational methods and discuss outstanding challenges. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 36:649-667, 2017.

Solid-State NMR Provides Evidence for Small-Amplitude Slow Domain Motions in a Multispanning Transmembrane α-Helical Protein

Proteins are dynamic entities and populate ensembles of conformations. Transitions between states within a conformational ensemble occur over a broad spectrum of amplitude and time scales, and are often related to biological function. Whereas solid-state NMR (SSNMR) spectroscopy has recently been used to characterize conformational ensembles of proteins in the microcrystalline states, its applications to membrane proteins remain limited. Here we use SSNMR to study conformational dynamics of a seven-helical transmembrane (TM) protein, Anabaena Sensory Rhodopsin (ASR) reconstituted in lipids. We report on site-specific measurements of the 15N longitudinal R1 and rotating frame R1ρ relaxation rates at two fields of 600 and 800 MHz and at two temperatures of 7 and 30 °C. Quantitative analysis of the R1 and R1ρ values and of their field and temperature dependencies provides evidence of motions on at least two time scales. We modeled these motions as fast local motions and slower collective motions of TM helices and of structured loops, and used the simple model-free and extended model-free analyses to fit the data and estimate the amplitudes, time scales and activation energies. Faster picosecond (tens to hundreds of picoseconds) local motions occur throughout the protein and are dominant in the middle portions of the TM helices. In contrast, the amplitudes of the slower collective motions occurring on the nanosecond (tens to hundreds of nanoseconds) time scales, are smaller in the central parts of helices, but increase toward their cytoplasmic sides as well as in the interhelical loops. ASR interacts with a soluble transducer protein on its cytoplasmic surface, and its binding affinity is modulated by light. The larger amplitude of motions on the cytoplasmic side of the TM helices correlates with the ability of ASR to undergo large conformational changes in the process of binding/unbinding the transducer.

Fifteen years of the Protein Crystallography Station: the coming of age of macromolecular neutron crystallography

The Protein Crystallography Station (PCS), located at the Los Alamos Neutron Scattering Center (LANSCE), was the first macromolecular crystallography beamline to be built at a spallation neutron source. Following testing and commissioning, the PCS user program was funded by the Biology and Environmental Research program of the Department of Energy Office of Science (DOE-OBER) for 13 years (2002-2014). The PCS remained the only dedicated macromolecular neutron crystallography station in North America until the construction and commissioning of the MaNDi and IMAGINE instruments at Oak Ridge National Laboratory, which started in 2012. The instrument produced a number of research and technical outcomes that have contributed to the field, clearly demonstrating the power of neutron crystallo-graphy in helping scientists to understand enzyme reaction mechanisms, hydrogen bonding and visualization of H-atom positions, which are critical to nearly all chemical reactions. During this period, neutron crystallography became a technique that increasingly gained traction, and became more integrated into macromolecular crystallography through software developments led by investigators at the PCS. This review highlights the contributions of the PCS to macromolecular neutron crystallography, and gives an overview of the history of neutron crystallography and the development of macromolecular neutron crystallography from the 1960s to the 1990s and onwards through the 2000s.

Chemical degradation of proteins in the solid state with a focus on photochemical reactions

Protein pharmaceuticals comprise an increasing fraction of marketed products but the limited solution stability of proteins requires considerable research effort to prepare stable formulations. An alternative is solid formulation, as proteins in the solid state are thermodynamically less susceptible to degradation. Nevertheless, within the time of storage a large panel of kinetically controlled degradation reactions can occur such as, e.g., hydrolysis reactions, the formation of diketopiperazine, condensation and aggregation reactions. These mechanisms of degradation in protein solids are relatively well covered by the literature. Considerably less is known about oxidative and photochemical reactions of solid proteins. This review will provide an overview over photolytic and non-photolytic degradation reactions, and specially emphasize mechanistic details on how solid structure may affect the interaction of protein solids with light.

How amide hydrogens exchange in native proteins

Amide hydrogen exchange (HX) is widely used in protein biophysics even though our ignorance about the HX mechanism makes data interpretation imprecise. Notably, the open exchange-competent conformational state has not been identified. Based on analysis of an ultralong molecular dynamics trajectory of the protein BPTI, we propose that the open (O) states for amides that exchange by subglobal fluctuations are locally distorted conformations with two water molecules directly coordinated to the N-H group. The HX protection factors computed from the relative O-state populations agree well with experiment. The O states of different amides show little or no temporal correlation, even if adjacent residues unfold cooperatively. The mean residence time of the O state is ∼100 ps for all examined amides, so the large variation in measured HX rate must be attributed to the opening frequency. A few amides gain solvent access via tunnels or pores penetrated by water chains including native internal water molecules, but most amides access solvent by more local structural distortions. In either case, we argue that an overcoordinated N-H group is necessary for efficient proton transfer by Grotthuss-type structural diffusion.

Analytical Aspects of Hydrogen Exchange Mass Spectrometry

This article reviews the analytical aspects of measuring hydrogen exchange by mass spectrometry (HX MS). We describe the nature of analytical selectivity in hydrogen exchange, then review the analytical tools required to accomplish fragmentation, separation, and the mass spectrometry measurements under restrictive exchange quench conditions. In contrast to analytical quantitation that relies on measurements of peak intensity or area, quantitation in HX MS depends on measuring a mass change with respect to an undeuterated or deuterated control, resulting in a value between zero and the maximum amount of deuterium that can be incorporated. Reliable quantitation is a function of experimental fidelity and to achieve high measurement reproducibility, a large number of experimental variables must be controlled during sample preparation and analysis. The method also reports on important qualitative aspects of the sample, including conformational heterogeneity and population dynamics.

Toward resolving the catalytic mechanism of dihydrofolate reductase using neutron and ultrahigh-resolution X-ray crystallography

Dihydrofolate reductase (DHFR) catalyzes the NADPH-dependent reduction of dihydrofolate (DHF) to tetrahydrofolate (THF). An important step in the mechanism involves proton donation to the N5 atom of DHF. The inability to determine the protonation states of active site residues and substrate has led to a lack of consensus regarding the catalytic mechanism involved. To resolve this ambiguity, we conducted neutron and ultrahigh-resolution X-ray crystallographic studies of the pseudo-Michaelis ternary complex of Escherichia coli DHFR with folate and NADP(+). The neutron data were collected to 2.0-Å resolution using a 3.6-mm(3) crystal with the quasi-Laue technique. The structure reveals that the N3 atom of folate is protonated, whereas Asp27 is negatively charged. Previous mechanisms have proposed a keto-to-enol tautomerization of the substrate to facilitate protonation of the N5 atom. The structure supports the existence of the keto tautomer owing to protonation of the N3 atom, suggesting that tautomerization is unnecessary for catalysis. In the 1.05-Å resolution X-ray structure of the ternary complex, conformational disorder of the Met20 side chain is coupled to electron density for a partially occupied water within hydrogen-bonding distance of the N5 atom of folate; this suggests direct protonation of substrate by solvent. We propose a catalytic mechanism for DHFR that involves stabilization of the keto tautomer of the substrate, elevation of the pKa value of the N5 atom of DHF by Asp27, and protonation of N5 by water that gains access to the active site through fluctuation of the Met20 side chain even though the Met20 loop is closed.

Neutron scattering techniques and applications in structural biology

Neutron scattering is exquisitely sensitive to the position, concentration, and dynamics of hydrogen atoms in materials and is a powerful tool for the characterization of structure-function and interfacial relationships in biological systems. Modern neutron scattering facilities offer access to a sophisticated, nondestructive suite of instruments for biophysical characterization that provides spatial and dynamic information spanning from Ångstroms to microns and from picoseconds to microseconds, respectively. Applications in structural biology range from the atomic-resolution analysis of individual hydrogen atoms in enzymes through to meso- and macro-scale analysis of complex biological structures, membranes, and assemblies. The large difference in neutron scattering length between hydrogen and deuterium allows contrast variation experiments to be performed and enables H/D isotopic labeling to be used for selective and systematic analysis of the local structure, dynamics, and interactions of multi-component systems. This overview describes the available techniques and summarizes their practical application to the study of biomolecular systems.

Direct observation of hydrogen atom dynamics and interactions by ultrahigh resolution neutron protein crystallography

The 1.1 Å, ultrahigh resolution neutron structure of hydrogen/deuterium (H/D) exchanged crambin is reported. Two hundred ninety-nine out of 315, or 94.9%, of the hydrogen atom positions in the protein have been experimentally derived and resolved through nuclear density maps. A number of unconventional interactions are clearly defined, including a potential O─H…π interaction between a water molecule and the aromatic ring of residue Y44, as well as a number of potential C─H…O hydrogen bonds. Hydrogen bonding networks that are ambiguous in the 0.85 Å ultrahigh resolution X-ray structure can be resolved by accurate orientation of water molecules. Furthermore, the high resolution of the reported structure has allowed for the anisotropic description of 36 deuterium atoms in the protein. The visibility of hydrogen and deuterium atoms in the nuclear density maps is discussed in relation to the resolution of the neutron data.

Room-temperature ultrahigh-resolution time-of-flight neutron and X-ray diffraction studies of H/D-exchanged crambin

The room-temperature (RT) X-ray structure of H/D-exchanged crambin is reported at 0.85 Å resolution. As one of the very few proteins refined with anisotropic atomic displacement parameters at two temperatures, the dynamics of atoms in the RT and 100 K structures are compared. Neutron diffraction data from an H/D-exchanged crambin crystal collected at the Protein Crystallography Station (PCS) showed diffraction beyond 1.1 Å resolution. This is the highest resolution neutron diffraction reported to date for a protein crystal and will reveal important details of the anisotropic motions of H and D atoms in protein structures.

The lighter side of a sweet reaction

The unique advantage of neutrons as biological probes is the ability to visualize hydrogen atoms in macromolecules. In this issue, Kovalevsky et al. (2010) solved an ensemble of xylose isomerase structures by neutron crystallography, and the determination of hydrogen atom rearrangements during the catalytic cycle provides insight into the enzyme's mechanism.

Normal mode paths for hydrogen exchange in the peptide ferrichrome

Possible paths for exposure to solvent and hydrogen exchange of the amide protons of ferrichrome, a cyclic hexapeptide, are examined. The paths are obtained from calculations of the vibrational normal modes of ferrichrome and correspond to low energy atomic displacements away from the local minimum in the multidimensional conformational space of the molecule. Exposure of exchangeable groups along the normal modes was determined by using the solvent accessible surface area algorithm of Lee and Richards. Three of the exchangeable protons (Gly(1,2,3,)) are largely exposed to solvent in the x-ray structure while the remaining three exchangeable protons of the ornithines are totally shielded from solvent. A very small number of normal mode displacements are found to expose the Orn(2) and Orn(3) amide groups while the Orn(1) amide proton remains shielded from solvent for all the paths studied. The effective paths for exposure of Orn(2) and Orn(3) correspond to the lowest frequency ( approximately 18 cm(-1)) motions. The paths are characterized in terms of the magnitude and energy of atomic displacements, correlated changes in dihedral angles, and the resulting changes in exposure and hydrogen bonding of exchangeable groups.

A joint x-ray and neutron study on amicyanin reveals the role of protein dynamics in electron transfer

The joint x-ray/neutron diffraction model of the Type I copper protein, amicyanin from Paracoccus denitrificans was determined at 1.8 A resolution. The protein was crystallized using reagents prepared in D(2)O. About 86% of the amide hydrogen atoms are either partially or fully exchanged, which correlates well with the atomic depth of the amide nitrogen atom and the secondary structure type, but with notable exceptions. Each of the four residues that provide copper ligands is partially deuterated. The model reveals the dynamic nature of the protein, especially around the copper-binding site. A detailed analysis of the presence of deuterated water molecules near the exchange sites indicates that amide hydrogen exchange is primarily due to the flexibility of the protein. Analysis of the electron transfer path through the protein shows that residues in that region are highly dynamic, as judged by hydrogen/deuterium exchange. This could increase the rate of electron transfer by transiently shortening through-space jumps in pathways or by increasing the atomic packing density. Analysis of C-HX bonding reveals previously undefined roles of these relatively weak H bonds, which, when present in sufficient number can collectively influence the structure, redox, and electron transfer properties of amicyanin.

Structure and dynamics of NBD1 from CFTR characterized using crystallography and hydrogen/deuterium exchange mass spectrometry

The DeltaF508 mutation in nucleotide-binding domain 1 (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR) is the predominant cause of cystic fibrosis. Previous biophysical studies on human F508 and DeltaF508 domains showed only local structural changes restricted to residues 509-511 and only minor differences in folding rate and stability. These results were remarkable because DeltaF508 was widely assumed to perturb domain folding based on the fact that it prevents trafficking of CFTR out of the endoplasmic reticulum. However, the previously reported crystal structures did not come from matched F508 and DeltaF508 constructs, and the DeltaF508 structure contained additional mutations that were required to obtain sufficient protein solubility. In this article, we present additional biophysical studies of NBD1 designed to address these ambiguities. Mass spectral measurements of backbone amide (1)H/(2)H exchange rates in matched F508 and DeltaF508 constructs reveal that DeltaF508 increases backbone dynamics at residues 509-511 and the adjacent protein segments but not elsewhere in NBD1. These measurements also confirm a high level of flexibility in the protein segments exhibiting variable conformations in the crystal structures. We additionally present crystal structures of a broader set of human NBD1 constructs, including one harboring the native F508 residue and others harboring the DeltaF508 mutation in the presence of fewer and different solubilizing mutations. The only consistent conformational difference is observed at residues 509-511. The side chain of residue V510 in this loop is mostly buried in all non-DeltaF508 structures but completely solvent exposed in all DeltaF508 structures. These results reinforce the importance of the perturbation DeltaF508 causes in the surface topography of NBD1 in a region likely to mediate contact with the transmembrane domains of CFTR. However, they also suggest that increased exposure of the 509-511 loop and increased dynamics in its vicinity could promote aggregation in vitro and aberrant intermolecular interactions that impede trafficking in vivo.

HAAD: A quick algorithm for accurate prediction of hydrogen atoms in protein structures

Hydrogen constitutes nearly half of all atoms in proteins and their positions are essential for analyzing hydrogen-bonding interactions and refining atomic-level structures. However, most protein structures determined by experiments or computer prediction lack hydrogen coordinates. We present a new algorithm, HAAD, to predict the positions of hydrogen atoms based on the positions of heavy atoms. The algorithm is built on the basic rules of orbital hybridization followed by the optimization of steric repulsion and electrostatic interactions. We tested the algorithm using three independent data sets: ultra-high-resolution X-ray structures, structures determined by neutron diffraction, and NOE proton-proton distances. Compared with the widely used programs CHARMM and REDUCE, HAAD has a significantly higher accuracy, with the average RMSD of the predicted hydrogen atoms to the X-ray and neutron diffraction structures decreased by 26% and 11%, respectively. Furthermore, hydrogen atoms placed by HAAD have more matches with the NOE restraints and fewer clashes with heavy atoms. The average CPU cost by HAAD is 18 and 8 times lower than that of CHARMM and REDUCE, respectively. The significant advantage of HAAD in both the accuracy and the speed of the hydrogen additions should make HAAD a useful tool for the detailed study of protein structure and function. Both an executable and the source code of HAAD are freely available at http://zhang.bioinformatics.ku.edu/HAAD.

Rapid determination of hydrogen positions and protonation states of diisopropyl fluorophosphatase by joint neutron and X-ray diffraction refinement

Hydrogen atoms constitute about half of all atoms in proteins and play a critical role in enzyme mechanisms and macromolecular and solvent structure. Hydrogen atom positions can readily be determined by neutron diffraction, and as such, neutron diffraction is an invaluable tool for elucidating molecular mechanisms. Joint refinement of neutron and X-ray diffraction data can lead to improved models compared with the use of neutron data alone and has now been incorporated into modern, maximum-likelihood based crystallographic refinement programs like CNS. Joint refinement has been applied to neutron and X-ray diffraction data collected on crystals of diisopropyl fluorophosphatase (DFPase), a calcium-dependent phosphotriesterase capable of detoxifying organophosphorus nerve agents. Neutron omit maps reveal a number of important features pertaining to the mechanism of DFPase. Solvent molecule W33, coordinating the catalytic calcium, is a water molecule in a strained coordination environment, and not a hydroxide. The smallest Ca-O-H angle is 53 degrees, well beyond the smallest angles previously observed. Residue Asp-229, is deprotonated, supporting a mechanism involving nucleophilic attack by Asp-229, and excluding water activation by the catalytic calcium. The extended network of hydrogen bonding interactions in the central water filled tunnel of DFPase is revealed, showing that internal solvent molecules form an important, integrated part of the overall structure.

Anion modulation of the 1H/2H exchange rates in backbone amide protons monitored by NMR spectroscopy

The ability of three anionic cosolutes (sulfate, thiocyanate, and chloride) in modulating the (1)H/(2)H exchange rates for backbone amide protons has been investigated using nuclear magnetic resonance (NMR) for two different proteins: the IGg-binding domain of protein L (ProtL) and the glucose-galactose-binding protein (GGBP). Our results show that moderate anion concentrations (0.2 M-1 M) regulate the exchange rate following the Hofmeister series: Addition of thiocyanate increases the exchange rates for both proteins, while sulfate and chloride (to a less extent) slow down the exchange reaction. In the presence of the salt, no alteration of the protein structure and minimal variations in the number of measurable peaks are observed. Experiments with model compounds revealed that the unfolded state is modulated in an equivalent way by these cosolutes. For ProtL, the estimated values for the local free energy change upon salt addition (m (3,DeltaG )) are consistent with the previously reported free energy contribution from the cosolute's preferential interaction/exclusion term indicating that nonspecific weak interactions between the anion and the amide groups constitute the dominant mechanism for the exchange-rate modulation. The same trend is also found for GGBP in the presence of thiocyanate, underlining the generality of the exchange-rate modulation mechanism, complementary to more investigated effects like the electrostatic interactions or specific anion binding to protein sites.

Neutron diffraction studies of Escherichia coli dihydrofolate reductase complexed with methotrexate

Hydrogen atoms play a central role in many biochemical processes yet are difficult to visualize by x-ray crystallography. Spallation neutron sources provide a new arena for protein crystallography with TOF measurements enhancing data collection efficiency and allowing hydrogen atoms to be located in smaller crystals of larger biological macromolecules. Here we report a 2.2-A resolution neutron structure of Escherichia coli dihydrofolate reductase (DHFR) in complex with methotrexate (MTX). Neutron data were collected on a 0.3-mm(3) D(2)O-soaked crystal at the Los Alamos Neutron Scattering Center. This study provides an example of using spallation neutrons to study protein dynamics, to identify protonation states directly from nuclear density maps, and to analyze solvent structure. Our structure reveals that the occluded loop conformation [monomer (mon.) A] of the DHFR.MTX complex undergoes greater H/D exchange compared with the closed-loop conformer (mon. B), partly because the Met-20 and beta(F-G) loops readily exchange in mon. A. The eight-stranded beta sheet of both DHFR molecules resists H/D exchange more than the helices and loops. However, the C-terminal strand, betaH, in mon. A is almost fully exchanged. Several D(2)Os form hydrogen bonds with exchanged amides. At the active site, the N1 atom of MTX is protonated and thus charged when bound to DHFR. Several D(2)Os are observed at hydrophobic surfaces, including two pockets near the MTX-binding site. A previously unidentified D(2)O hydrogen bonds with the catalytic D27 in mon. B, stabilizing its negative charge.

Hydrogen/deuterium-exchange (H/D-Ex) of PPARgamma LBD in the presence of various modulators

A nuclear receptor, peroxisome proliferator-activated receptor gamma (PPARgamma), is a ligand-dependent transcription factor involved in glucose homeostasis and adipocyte differentiation. PPARgamma is the molecular target of various natural and synthetic molecules, including anti-diabetic agents such as rosiglitazone. Amide hydrogen/deuterium-exchange (H/D-Ex), coupled with proteolysis and mass spectrometry, was applied to study the dynamics of the PPARgamma ligand binding domain (LBD) with or without molecules that modulate PPARgamma activity. The H/D-Ex patterns of ligand-free PPARgamma LBD show that the ligand binding pocket of LBD is significantly more dynamic than the rest of the LBD. Presumably, the binding pocket is intrinsically disordered in order to accommodate different ligands. The presence of two full agonists (rosiglitazone and GW1929), a partial agonist (nTZDpa), and a covalent antagonist (GW9662), changed the dynamics/conformation of PPARgamma LBD and slowed the H/D exchange rate in various regions of the protein. The full agonists slowed the H/D exchange more globally and to a greater extent than the partial agonist or the antagonist, indicating that the full agonist stabilizes the PPARgamma LBD more than the partial agonist or the antagonist. One interesting observation is that the two full agonists significantly stabilized helix 12 while the partial agonist and the antagonist did not perturb the H/D exchange of this region. The results showed that the change in protein dynamics induced by ligand binding may be an important factor for the activation of genes and that H/D-Ex is a useful method for analyzing the biological activity of drug leads.

New accessory for studies of isotopic 1H/2H exchange and biomolecular interactions using transmission infrared spectroscopy

We present here a new accessory for IR transmission measurements of 1H/2H exchange, as an ancillary tool for monitoring structural features of biomolecules in aqueous solution. This new accessory results from the combination of two dialysis membranes and a conventional liquid cell having two cylinders containing 2H2O buffer. When compared with conventional transmission measurements, carried out either after dissolving lyophilized biomolecules in 2H2O or after dialyzing the aqueous solution considered against 2H2O buffer, this accessory shows the following advantages: (1) controlled measurements over the initial steps of this isotopic exchange and absence of molecular aggregation, and (2) smaller sample amounts. This new Fourier transform IR cell can also be used to analyze ligand-biomolecule and drug-cell interactions.

Carbon-13 NMR method for the detection of correlated hydrogen exchange at adjacent backbone peptide amides and its application to hydrogen exchange in five antiparallel beta strands within the hydrophobic core of Streptomyces subtilisin inhibitor (SSI)

A novel method for monitoring proton-deuteron (H/D) exchange at backbone amides is based on the observation of H/D isotope effects on the (13)C NMR signals from peptide carbonyls. The line shape of the carbonyl (13)C(i) signal is influenced by differential H/D occupancy at the two adjacent amides: the H(N)(i)(+1) (beta site) and the H(N)(i) (gamma site). At a carbon frequency of 75.4 MHz, the H --> D isotope shifts on the (13)C signal are about 5-7 Hz for exchange at the beta site and 2 Hz or less for exchange at the gamma site. Because the effects at the two sites are additive, the time dependence of the line shape of a particular carbonyl resonance can report not only the exchange rates at the individual sites but also the level of dual exchange. Therefore, the data can be analyzed to determine the rate (k(c)) and degree of correlated exchange (X(betagamma)) at the two sites. We have applied this approach to the investigation of the pH dependence of hydrogen exchange at several adjacent residues in Streptomyces subtilisin inhibitor (SSI). Two selectively labeled SSI proteins were produced: one with selective (13)C' labeling at all valyl residues and one with selective (13)C' labeling at all leucyl residues. This permitted the direct observation by one-dimensional (13)C NMR of selected carbonyl signals from residues with slowly exchanging amides at the i and i + 1 positions. The residues investigated were located in an alpha helix and in a five-stranded antiparallel beta sheet. Samples of the two labeled proteins were prepared at various pH values, and (13)C NMR spectra were collected at 50 degrees C prior to and at various times after transferring the sample from H(2)O to (2)H(2)O. Most of the slowly exchanging amides studied were intramolecular hydrogen-bond donors. In agreement with prior studies, the results indicated that the exchange rates of the amide hydrogens in proteins are governed not only by hydrogen bonding but also by other factors. For example, the amide hydrogen of Thr34 exchanges rapidly even though it is an intramolecular hydrogen-bond donor. Over nearly the whole pH range studied, the apparent rates of uncorrelated exchange (k(beta) and k(gamma)) were proportional to [OH(-)] and the apparent rates of correlated exchange at two adjacent sites (k(c)) were roughly proportional to [OH(-)](2). This enabled us to extract the pH-independent exchange rates (k(beta) degrees , k(gamma) degrees , and k(c) degrees ). In all cases in which correlated exchange could be measured, the observed sigmoidal pH dependence of X(betagamma) could be replicated roughly from the derived pH-independent rates.

Mapping intersubunit interactions of the regulatory subunit (RIalpha) in the type I holoenzyme of protein kinase A by amide hydrogen/deuterium exchange mass spectrometry (DXMS)

Protein kinase A (PKA), a central locus for cAMP signaling in the cell, is composed of regulatory (R) and catalytic (C) subunits. The C-subunits are maintained in an inactive state by binding to the R-subunit dimer in a tetrameric holoenzyme complex (R(2)C(2)). PKA is activated by cAMP binding to the R-subunits which induces a conformational change leading to release of the active C-subunit. Enzymatic activity of the C-subunit is thus regulated by cAMP via the R-subunit, which toggles between cAMP and C-subunit bound states. The R-subunit is composed of a dimerization/docking (D/D) domain connected to two cAMP-binding domains (cAMP:A and cAMP:B). While crystal structures of the free C-subunit and cAMP-bound states of a deletion mutant of the R-subunit are known, there is no structure of the holoenzyme complex or of the cAMP-free state of the R-subunit. An important step in understanding the cAMP-dependent activation of PKA is to map the R-C interface and characterize the mutually exclusive interactions of the R-subunit with cAMP and C-subunit. Amide hydrogen/deuterium exchange mass spectrometry is a suitable method that has provided insights into the different states of the R-subunit in solution, thereby allowing mapping of the effects of cAMP and C-subunit on different regions of the R-subunit. Our study has localized interactions with the C-subunit to a small contiguous surface on the cAMP:A domain and the linker region. In addition, C-subunit binding causes increased amide hydrogen exchange within both cAMP-domains, suggesting that these regions become more flexible in the holoenzyme and are primed to bind cAMP. Furthermore, the difference in the protection patterns between RIalpha and the previously studied RIIbeta upon cAMP-binding suggests isoform-specific differences in cAMP-dependent regulation of PKA activity.

Conformational changes in chemically modified Escherichia coli thioredoxin monitored by H/D exchange and electrospray ionization mass spectrometry

Hydrogen/deuterium (H/D) exchange in combination with electrospray ionization mass spectrometry and near-ultraviolet (UV) circular dichroism (CD) was used to study the conformational properties and thermal unfolding of Escherichia coli thioredoxin and its Cys32-alkylated derivatives in 1% acetic acid (pH 2.7). Thermal unfolding of oxidized (Oxi) and reduced (Red) -thioredoxin (TRX) and Cys-32-ethylglutathionyl (GS-ethyl-TRX) and Cys-32-ethylcysteinyl (Cys-ethyl-TRX), which are derivatives of Red-TRX, follow apparent EX1 kinetics as charge-state envelopes, H/D mass spectral exchange profiles, and near-UV CD appear to support a two-state folding/unfolding model. Minor mass peaks in the H/D exchange profiles and nonsuperimposable MS- and CD-derived melting curves, however, suggest the participation of unfolding intermediates leading to the conclusion that the two-state model is an oversimplification of the process. The relative stabilities as measured by melting temperatures by both CD and mass spectral charge states are, Oxi-TRX, GS-ethyl-TRX, Cys-ethyl-TRX, and Red-TRX. The introduction of the Cys-32-ethylglutathionyl group provides extra stabilization that results from additional hydrogen bonding interactions between the ethylglutathionyl group and the protein. Near-UV CD data show that the local environment near the active site is perturbed to almost an identical degree regardless of whether alkylation at Cys-32 is by the ethylglutathionyl group, or the smaller, nonhydrogen-bonding ethylcysteinyl group. Mass spectral data, however, indicate a tighter structure for GS-ethyl-TRX.

Affinity labeling of cysteine-mutants evidences contact residues in modeled receptor binding sites

To investigate the topology of binding sites in two ionotropic receptors, we have initiated a strategy combining affinity labeling with cysteine-scanning mutagenesis. For the GABAA receptor we have used reactive derivatives of non-competitive blockers (NCBs) to explore interacting positions in its channel. The polypeptide positions of the M2 segment of the alpha1 subunit which we mutated into cysteine were selected for their established accessibility, as determined by the substituted-cysteine accessibility method (SCAM). Using the Xenopus oocyte expression system, we show that receptors containing mutations V257C and S272C are inactivated by several reactive NCBs. These position-selective inactivations lead to an analysis of NCB binding in the channel. For the NMDA receptor glycine-binding site, the prototype antagonist L-701,324 was derivatized at different positions with different reactive groups. The receptor positions to mutate into cysteine were selected after a 3-D homology model. The observed receptor inactivations are mutant- and probe-selective, leading to an unambiguous chemical docking of the antagonist pharmacophore and supporting the model. The site-specificity of the inactivating reactions is assessed by protection experiments and by mutant to wild-type (WT) comparisons. The scope and limitations of the method are briefly discussed.

A synthetic picture of intramolecular dynamics of proteins. Towards a contemporary statistical theory of biochemical processes

An increasing body of experimental evidence indicates the slow character of internal dynamics of native proteins. The important consequence of this is that theories of chemical reactions, used hitherto, appear inadequate for description of most biochemical reactions. Construction of a contemporary, truly advanced statistical theory of biochemical processes will need simple but realistic models of microscopic dynamics of biomolecules. In this review, intended to be a contribution towards this direction, three topics are considered. First, an intentionally simplified picture of dynamics of native proteins which emerges from recent investigations is presented. Fast vibrational modes of motion, of periods varying from 10(-14) to 10(-11) s, are contrasted with purely stochastic conformational transitions. Significant evidence is adduced that the relaxation time spectrum of the latter spreads in the whole range from 10(-11) to 10(5) s or longer, and up to 10(-7) s it is practically quasi-continuous. Next, the essential ideas of the theory of reaction rates based on stochastic models of intramolecular dynamics are outlined. Special attention is paid to reactions involving molecules in the initial conformational substrates confirmed to the transition state, which is realized in actual experimental situations. And finally, the two best experimentally justified classes of models of conformational transition dynamics, symbolically referred to as "protein glass" and "protein machine", are described and applied to the interpretation of a few simple biochemical processes, perhaps the most important result reported is the demonstration of the possibility of predominance of the short initial condition-dependent stage of protein involved reactions over the main stage described by the standard kinetics. This initial stage, and not the latter, is expected to be responsible for the coupling of component reactions in the complete enzymatic cycles as well as more complex processes of biological free energy transduction.

Ligand binding to proteins: the binding landscape model

Models of ligand binding are often based on four assumptions: (1) steric fit: that binding is determined mainly by shape complementarity; (2) native binding: that ligands mainly bind to native states; (3) locality: that ligands perturb protein structures mainly at the binding site; and (4) continuity: that small changes in ligand or protein structure lead to small changes in binding affinity. Using a generalization of the 2D HP lattice model, we study ligand binding and explore these assumptions. We first validate the model by showing that it reproduces typical binding behaviors. We observe ligand-induced denaturation, ANS and heme-like binding, and "lock-and-key" and "induced-fit" specific binding behaviors characterized by Michaelis-Menten or more cooperative types of binding isotherms. We then explore cases where the model predicts violations of the standard assumptions. For example, very different binding modes can result from two ligands of identical shape. Ligands can sometimes bind highly denatured states more tightly than native states and yet have Michaelis-Menten isotherms. Even low-population binding to denatured states can cause changes in global stability, hydrogen-exchange rates, and thermal B-factors, contrary to expectations, but in agreement with experiments. We conclude that ligand binding, similar to protein folding, may be better described in terms of energy landscapes than in terms of simpler mass-action models.

Hydrogen exchange: the modern legacy of Linderstrøm-Lang

This discussion, prepared for the Protein Society's symposium honoring the 100th anniversary of Kaj Linderstrøm-Lang, shows how hydrogen exchange approaches initially conceived and implemented by Lang and his colleagues some 50 years ago are contributing to current progress in structural biology. Examples are chosen from the active protein folding field. Hydrogen exchange methods now make it possible to define the structure of protein folding intermediates in various contexts: as tenuous molten globule forms at equilibrium under destabilizing conditions, in kinetic intermediates that exist for less than one second, and as infinitesimally populated excited state forms under native conditions. More generally, similar methods now find broad application in studies of protein structure, energetics, and interactions. This article considers the rise of these capabilities from their inception at the Carlsberg Labs to their contemporary role as a significant tool of modern structural biology.