Proton magnetic resonance studies of the tyrosine residues of hen lysozyme-assignment and detection of conformational mobility

Formerly degenerate seventh zinc finger domain from transcription factor ZNF711 rehabilitated by experimental NMR structure

Domain Z7 of nuclear transcription factor ZNF711 has the consensus last metal-ligand H23 found in odd-numbered zinc-fingers of this protein replaced by a phenylalanine. Ever since the discovery of ZNF711 it has been thought that Z7 is probably non-functional because of the H23F substitution. The presence of H26 three positions downstream prompted us to examine if this histidine could substitute as the last metal ligand. The Z7 domain adopts a stable tertiary structure upon metal binding. The NMR structure of Zn2+-bound Z7 shows the classical ββα-fold of CCHH zinc fingers. Mutagenesis and pH titration experiments indicate that H26 is not involved in metal binding and that Z7 has a tridentate metal-binding site comprised of only residues C3, C6, and H19. By contrast, an F23H mutation that introduces a histidine in the consensus position forms a tetradentate ligand. The structure of the WT Z7 is stable causing restricted ring-flipping of phenyalanines 10 and 23. Dynamics are increased with either the H26A or F23H substitutions and aromatic ring rotation is no longer hindered in the two mutants. The mutations have only small effects on the Kd values for Zn2+ and Co2+ and retain the high thermal stability of the WT domain above 80 °C. Like two previously reported designed zinc fingers with the last ligand replaced by water, the WT Z7 domain is catalytically active, hydrolyzing 4-nitophenyl acetate. We discuss the implications of naturally occurring tridentate zinc fingers for cancer mutations and drug targeting of notoriously undruggable transcription factors. Our findings that Z7 can fold with only a subset of three metal ligands suggests the recent view that most everything about protein structure can be predicted through homology modeling might be premature for at least the resilient and versatile zinc-finger motif.

Bioactives from medicinal herb against bedaquiline resistant tuberculosis: removing the dark clouds from the horizon

Tuberculosis is a contagious bacterial ailment that primarily affects the lungs and is brought on by the bacterium Mycobacterium tuberculosis (MTB). An antimycobacterial medication called bedaquiline (BQ) is specified to treat multidrug-resistant tuberculosis (MDR-TB). Despite its contemporary use in clinical practice, the mutations (D32 A/G/N/V/P) constrain the potential of BQ by causing transitions in the structural conformation of the atpE subunit-c after binding. In this study, we have taken the benzylisoquinoline alkaloids from thalictrum foliolosum due to its antimicrobial activity reported in prior literature. We used an efficient and optimized structure-based strategy to examine the wild type (WT) and mutated protein upon molecule binding. Our results emphasize the drastic decline in BQ binding affinity of mutant and WT atpE subunit-c complexes compared to thalirugidine (top hit) from thalictrum foliolosum. The decrease in BQ binding free energy is due to electrostatic energy because nearly every atom in a macromolecule harbors a partial charge, and molecules taking part in molecular recognition will interact electrostatically. Similarly, the high potential mean force of thalirugidine than BQ in WT and mutant complexes demonstrated the remarkable ability to eradicate mycobacteria efficiently. Furthermore, the Alamar blue cell viability and ATP determination assay were performed to validate the computational outcomes in search of novel antimycobacterial. Upon closer examination of the ATP determination assay, it became apparent that both BQ and thalirugidine showed similar reductions in ATP levels at their respective MICs, presenting a potential common mechanism of action.

NMR Studies of Aromatic Ring Flips to Probe Conformational Fluctuations in Proteins

Aromatic residues form a significant part of the protein core, where they make tight interactions with multiple surrounding side chains. Despite the dense packing of internal side chains, the aromatic rings of phenylalanine and tyrosine residues undergo 180° rotations, or flips, which are mediated by transient and large-scale "breathing" motions that generate sufficient void volume around the aromatic ring. Forty years after the seminal work by Wagner and Wüthrich, NMR studies of aromatic ring flips are now undergoing a renaissance as a powerful means of probing fundamental dynamic properties of proteins. Recent developments of improved NMR methods and isotope labeling schemes have enabled a number of advances in addressing the mechanisms and energetics of aromatic ring flips. The nature of the transition states associated with ring flips can be described by thermodynamic activation parameters, including the activation enthalpy, activation entropy, activation volume, and also the isothermal volume compressibility of activation. Consequently, it is of great interest to study how ring flip rate constants and activation parameters might vary with protein structure and external conditions like temperature and pressure. The field is beginning to gather such data for aromatic residues in a variety of environments, ranging from surface exposed to buried. In the future, the combination of solution and solid-state NMR spectroscopy together with molecular dynamics simulations and other computational approaches is likely to provide detailed information about the coupled dynamics of aromatic rings and neighboring residues. In this Perspective, we highlight recent developments and provide an outlook toward the future.

Global insights into the fine tuning of human A2AAR conformational dynamics in a ternary complex with an engineered G protein viewed by NMR

G protein-coupled receptor (GPCR) conformational plasticity enables formation of ternary signaling complexes with intracellular proteins in response to binding extracellular ligands. We investigate the dynamic process of GPCR complex formation in solution with the human A_2A adenosine receptor (A_2AAR) and an engineered G_s protein, mini-G_s. 2D nuclear magnetic resonance (NMR) data with uniform stable isotope-labeled A_2AAR enabled a global comparison of A_2AAR conformations between complexes with an agonist and mini-G_s and with an agonist alone. The two conformations are similar and show subtle differences at the receptor intracellular surface, supporting a model whereby agonist binding alone is sufficient to populate a conformation resembling the active state. However, an A_2AAR "hot spot" connecting the extracellular ligand-binding pocket to the intracellular surface is observed to be highly dynamic in the ternary complex, suggesting a mechanism for allosteric connection between the bound G protein and the drug-binding pocket involving structural plasticity of the "toggle switch" tryptophan.

Slow ring flips in aromatic cluster of GB1 studied by aromatic 13C relaxation dispersion methods

Ring flips of phenylalanine and tyrosine are a hallmark of protein dynamics. They report on transient breathing motions of proteins. In addition, flip rates also depend on stabilizing interactions in the ground state, like aromatic stacking or cation-π interaction. So far, experimental studies of ring flips have almost exclusively been performed on aromatic rings without stabilizing interactions. Here we investigate ring flip dynamics of Phe and Tyr in the aromatic cluster in GB1. We found that all four residues of the cluster, Y3, F30, Y45 and F52, display slow ring flips. Interestingly, F52, the central residue of the cluster, which makes aromatic contacts with all three others, is flipping significantly faster, while the other rings are flipping with the same rates within margin of error. Determined activation enthalpies and activation volumes of these processes are in the same range of other reported ring flips of single aromatic rings. There is no correlation of the number of aromatic stacking interactions to the activation enthalpy, and no correlation of the ring's extent of burying to the activation volume. Because of these findings, we speculate that F52 is undergoing concerted ring flips with each of the other rings.

Aromatic Ring Dynamics, Thermal Activation, and Transient Conformations of a 468 kDa Enzyme by Specific 1H-13C Labeling and Fast Magic-Angle Spinning NMR

Aromatic residues are located at structurally important sites of many proteins. Probing their interactions and dynamics can provide important functional insight but is challenging in large proteins. Here, we introduce approaches to characterize the dynamics of phenylalanine residues using ¹H-detected fast magic-angle spinning (MAS) NMR combined with a tailored isotope-labeling scheme. Our approach yields isolated two-spin systems that are ideally suited for artifact-free dynamics measurements, and allows probing motions effectively without molecular weight limitations. The application to the TET2 enzyme assembly of ∼0.5 MDa size, the currently largest protein assigned by MAS NMR, provides insights into motions occurring on a wide range of time scales (picoseconds to milliseconds). We quantitatively probe ring-flip motions and show the temperature dependence by MAS NMR measurements down to 100 K. Interestingly, favorable line widths are observed down to 100 K, with potential implications for DNP NMR. Furthermore, we report the first ¹³C R_1ρ MAS NMR relaxation-dispersion measurements and detect structural excursions occurring on a microsecond time scale in the entry pore to the catalytic chamber and at a trimer interface that was proposed as the exit pore. We show that the labeling scheme with deuteration at ca. 50 kHz MAS provides superior resolution compared to 100 kHz MAS experiments with protonated, uniformly ¹³C-labeled samples.

An NMR View of Protein Dynamics in Health and Disease

Biological molecules are often highly dynamic, and this flexibility can be critical for function. The large range of sampled timescales and the fact that many of the conformers that are continually explored are only transiently formed and sparsely populated challenge current biophysical approaches. Solution nuclear magnetic resonance (NMR) spectroscopy has emerged as a powerful method for characterizing biomolecular dynamics in detail, even in cases where excursions involve short-lived states. Here, we briefly review a number of NMR experiments for studies of biomolecular dynamics on the microsecond-to-second timescale and focus on applications to protein and nucleic acid systems that clearly illustrate the functional relevance of motion in both health and disease.

General trends of dihedral conformational transitions in a globular protein

Dihedral conformational transitions are analyzed systematically in a model globular protein, cytochrome P450cam, to examine their structural and chemical dependences through combined conventional molecular dynamics (cMD), accelerated molecular dynamics (aMD) and adaptive biasing force (ABF) simulations. The aMD simulations are performed at two acceleration levels, using dihedral and dual boost, respectively. In comparison with cMD, aMD samples protein dihedral transitions approximately two times faster on average using dihedral boost, and ∼ 3.5 times faster using dual boost. In the protein backbone, significantly higher dihedral transition rates are observed in the bend, coil, and turn flexible regions, followed by the β bridge and β sheet, and then the helices. Moreover, protein side chains of greater length exhibit higher transition rates on average in the aMD-enhanced sampling. Side chains of the same length (particularly Nχ = 2) exhibit decreasing transition rates with residues when going from hydrophobic to polar, then charged and aromatic chemical types. The reduction of dihedral transition rates is found to be correlated with increasing energy barriers as identified through ABF free energy calculations. These general trends of dihedral conformational transitions provide important insights into the hierarchical dynamics and complex free energy landscapes of functional proteins.

Cavities and Excited States in Proteins

Protein cavities or voids are observed as defects in atomic packing. Cavities have long been suggested to play important roles in protein dynamics and function, but the underlying origin and mechanism remains elusive. Here, recent studies about the cavities characterized by high-pressure NMR spectroscopy have been reviewed. Analysis of the pressure-dependent chemical shifts showed both linear and nonlinear response of proteins to pressure. The linear response corresponded to compression within the native ensemble, while the nonlinear response indicated the involvement of low-lying excited states that were different from the native state. The finding of non-linear pressure shifts in various proteins suggested that the existence of the low-lying excited states was common for globular proteins. However, the absolute nonlinear coefficient values varied significantly from protein to protein, and showed a good correlation with the density of cavities. Extensive studies on hen lysozyme as a model system showed that the cavity hydration and water penetration into the interior of proteins was an origin of the conformational transition to the excited states. The importance of cavities for protein function and evolution has also been explained. In addition to these "equilibrium" cavities, there are also "transient" cavities formed in the interior of the protein structure, as manifested by the ring flip motions of aromatic rings. The significance of transient cavities, reflecting an intrinsic dynamic nature within the native state, has also been discussed.

Off-resonance rotating-frame relaxation dispersion experiment for 13C in aromatic side chains using L-optimized TROSY-selection

Protein dynamics on the microsecond-millisecond time scales often play a critical role in biological function. NMR relaxation dispersion experiments are powerful approaches for investigating biologically relevant dynamics with site-specific resolution, as shown by a growing number of publications on enzyme catalysis, protein folding, ligand binding, and allostery. To date, the majority of studies has probed the backbone amides or side-chain methyl groups, while experiments targeting other sites have been used more sparingly. Aromatic side chains are useful probes of protein dynamics, because they are over-represented in protein binding interfaces, have important catalytic roles in enzymes, and form a sizable part of the protein interior. Here we present an off-resonance R 1ρ experiment for measuring microsecond to millisecond conformational exchange of aromatic side chains in selectively (13)C labeled proteins by means of longitudinal- and transverse-relaxation optimization. Using selective excitation and inversion of the narrow component of the (13)C doublet, the experiment achieves significant sensitivity enhancement in terms of both signal intensity and the fractional contribution from exchange to transverse relaxation; additional signal enhancement is achieved by optimizing the longitudinal relaxation recovery of the covalently attached (1)H spins. We validated the L-TROSY-selected R 1ρ experiment by measuring exchange parameters for Y23 in bovine pancreatic trypsin inhibitor at a temperature of 328 K, where the ring flip is in the fast exchange regime with a mean waiting time between flips of 320 μs. The determined chemical shift difference matches perfectly with that measured from the NMR spectrum at lower temperatures, where separate peaks are observed for the two sites. We further show that potentially complicating effects of strong scalar coupling between protons (Weininger et al. in J Phys Chem B 117: 9241-9247, 2013b) can be accounted for using a simple expression, and provide recommendations for data acquisition when the studied system exhibits this behavior. The present method extends the repertoire of relaxation methods tailored for aromatic side chains by enabling studies of faster processes and improved control over artifacts due to strong coupling.

Protein conformational space populated in solution probed with aromatic residual dipolar (13) C-(1) H couplings

The use of aromatic (13) C-(1) H residual dipolar couplings (RDCs) to probe the conformational space populated in solution is demonstrated for the protein BPTI. RDCs allow one to assess accuracy of atomic resolution structures and potentially to characterize low-populated subspaces corresponding to "excited states" in conformationally labile systems. They also allow one to assess sampling accuracy of molecular dynamics simulations.

NMR spectroscopy brings invisible protein states into focus

Molecular dynamics are essential for protein function. In some cases these dynamics involve the interconversion between ground state, highly populated conformers and less populated higher energy structures ('excited states') that play critical roles in biochemical processes. Here we describe recent advances in NMR spectroscopy methods that enable studies of these otherwise invisible excited states at an atomic level and that help elucidate their important relation to function. We discuss a range of examples from molecular recognition, ligand binding, enzyme catalysis and protein folding that illustrate the role that motion plays in 'funneling' conformers along preferred pathways that facilitate their biological function.

Selective characterization of microsecond motions in proteins by NMR relaxation

The three-dimensional structures of macromolecules fluctuate over a wide range of time-scales. Separating the individual dynamic processes according to frequency is of importance in relating protein motions to biological function and stability. We present here a general NMR method for the specific characterization of microsecond motions at backbone positions in proteins even in the presence of other dynamics such as large-amplitude nanosecond motions and millisecond chemical exchange processes. The method is based on measurement of relaxation rates of four bilinear coherences and relies on the ability of strong continuous radio frequency fields to quench millisecond chemical exchange. The utility of the methodology is demonstrated and validated through two specific examples focusing on the thermo-stable proteins, ubiquitin and protein L, where it is found that small-amplitude microsecond dynamics are more pervasive than previously thought. Specifically, these motions are localized to alpha helices, loop regions, and regions along the rim of beta sheets in both of the proteins examined. A third example focuses on a 28 kDa ternary complex of the chaperone Chz1 and the histones H2A.Z/H2B, where it is established that pervasive microsecond motions are localized to a region of the chaperone that is important for stabilizing the complex. It is further shown that these motions can be well separated from extensive millisecond dynamics that are also present and that derive from exchange of Chz1 between bound and free states. The methodology is straightforward to implement, and data recorded at only a single static magnetic field are required.

Observing biological dynamics at atomic resolution using NMR

Biological macromolecules are highly flexible and continually undergo conformational fluctuations on a broad spectrum of timescales. It has long been recognized that dynamics have an important role in the action of these molecules. However, the relationship between molecular function and motion is extremely challenging to delineate, because the conformational space available to macromolecules is vast and the relevant excursions can be infrequent and short-lived. Recent advances in solution nuclear magnetic resonance (NMR) spectroscopy permit biomolecular dynamics to be observed with unprecedented detail. Applications of these new NMR techniques to the study of fundamental processes such as binding and catalysis have provided new insights into how living systems operate at an atomic level.

What Contributions to Protein Side-chain Dynamics are Probed by NMR Experiments? A Molecular Dynamics Simulation Analysis

Molecular dynamics simulations of the structurally homologous proteins TNfn3 and FNfn10 have been used to investigate the contributions to side-chain dynamics measured by NMR relaxation experiments. The results reproduce the variation in core side-chain dynamics observed by NMR and highlight the relevance of anharmonic motion and transitions between local minima for explaining NMR side-chain order parameters. A method is described for calculating converged order parameters by use of replica exchange molecular dynamics in conjunction with an implicit solvent model. These simulations allow the influence of various factors, such as the flexibility of side-chains and their free volume, on the mobility to be tested by perturbing the system. Deletion mutations are found to have the largest effect on the more densely packed FNfn10. Some counterintuitive effects are seen, such as an increase in order parameters close to deletion mutation sites, but these can be rationalized in terms of direct interactions with the modified side-chains. A statistical analysis of published order parameters supports the conclusions drawn from the simulations.

Infrequent cavity-forming fluctuations in HPr from Staphylococcus carnosus revealed by pressure- and temperature-dependent tyrosine ring flips

Infrequent structural fluctuations of a globular protein is seldom detected and studied in detail. One tyrosine ring of HPr from Staphylococcus carnosus, an 88-residue phosphocarrier protein with no disulfide bonds, undergoes a very slow ring flip, the pressure and temperature dependence of which is studied in detail using the on-line cell high-pressure nuclear magnetic resonance technique in the pressure range from 3 MPa to 200 MPa and in the temperature range from 257 K to 313 K. The ring of Tyr6 is buried sandwiched between a beta-sheet and alpha-helices (the water-accessible area is less than 0.26 nm2), its hydroxyl proton being involved in an internal hydrogen bond. The ring flip rates 10(1)-10(5) s(-1) were determined from the line shape analysis of H(delta1, delta2) and H(epsilon1,epsilon2) of Tyr6, giving an activation volume DeltaV++ of 0.044 +/- 0.008 nm3 (27 mL mol(-1)), an activation enthalpy DeltaH++ of 89 +/- 10 kJ mol(-1), and an activation entropy DeltaS++ of 16 +/- 2 JK(-1) mol(-1). The DeltaV++) and DeltaH++ values for HPr found previously for Tyr and Phe ring flips of BPTI and cytochrome c fall within the range of DeltaV(double dagger) of 28 to 51 mL mol(-1) and DeltaH++ of 71 to 155 kJ mol(-1). The fairly common DeltaV++ and DeltaH++ values are considered to represent the extra space or cavity required for the ring flip and the extra energy required to create a cavity, respectively, in the core part of a globular protein. Nearly complete cold denaturation was found to take place at 200 MPa and 257 K independently from the ring reorientation process.

Steric effects of isoleucine 107 on heme reorientation reaction in human myoglobin

Structural factors to regulate the heme reorientation reaction in myoglobin were examined and we found that the side chain at position 107 (Ile107), which is located between the 2-vinyl and 3-methyl groups of heme, forms a kinetic barrier for the heme rotation about the alpha-gamma axis. The phenylalanine-substituted mutant showed an extremely slow heme reorientation rate, compared to that of the wild-type protein, while replacement by the decreased side chain, valine, at position 107 accelerated the reorientation reaction. Considering that the spectroscopic data show only minor structural changes in the heme environments of the Ile107 mutants, the side chain at position 107 sterically interacts with the heme peripheral groups in the activation state for the heme reorientation, which supports the intramolecular mechanism that the heme rotates about the alpha-gamma axis without leaving the "protein cage."

Aromatic Ring Currents at a Protein Surface: Use of (1)H-NMR Chemical Shifts to Refine the Structure of a Naked β Sheet

Abstract The naked β sheet, a newly recognized motif of protein structure, exhibits ordered surfaces in the absence of a conventional hydrophobic core. A model is provided by an archaeal Zn ribbon homologous to eukaryotic RNA polymerase II subunit 9 (RPB9). This subunit, which regulates transcriptional start-site selection and downstream pausing, contains Zn(2+)-binding motifs similar to those of general transcription factors TFIIB and TFIIS. Interestingly, distance-geometry yields two models of the archaeal Zn ribbon differing in the orientation of a conserved tyrosine side chain on the well-ordered surface of the naked β-sheet. The models are equally consistent with conventional restraints and otherwise contain indistinguishable structural features, including a tetrahedral Cys(4) Zn(2+)-binding sites, four antiparallel β-strands, and disordered loop. Due to the change in tyrosine orientation and correlated changes in the configuration of neighboring side chains, the two models predict inequivalent patterns of aromatic ring-current shifts. The observed secondary shifts of adjoining resonances are shown to be consistent with one model but not the other. In the consistent model the surface of the β-sheet contains successive aromatic edge-to-face contacts in accord with semi-classical and ab initio potentials. We speculate that the aromatic-rich surface of the hyperthermophilic RPB9 domain contributes its thermodynamic stability and provides a nucleic-acid-binding site in the eukaryotic and archaeal transcriptional machinery. The present study demonstrates how the reduced dimensionality of a surface can lead to ambiguities in the interpretation of nuclear Overhauser enhancements. The results illustrate the utility of chemical shifts at such a surface in the cross-validation of a high-resolution solution structure.

Structure-based analysis of protein dynamics: comparison of theoretical results for hen lysozyme with X-ray diffraction and NMR relaxation data

An analytical approach based on Gaussian network model (GNM) is proposed for predicting the rotational dynamics of proteins. The method, previously shown to successfully reproduce X-ray crystallographic temperature factors for a series of proteins is extended here to predict bond torsional mobilities and reorientation of main chain amide groups probed by 15N-H nuclear magnetic resonance (NMR) relaxation. The dynamics of hen egg-white lysozyme (HEWL) in the folded state is investigated using the proposed approach. Excellent agreement is observed between theoretical results and experimental (X-ray diffraction and NMR relaxation) data. The analysis reveals the important role of coupled rotations, or cross-correlations between dihedral angle librations, in defining the relaxation mechanism on a local scale. The crystal and solution structures exhibit some differences in their local motions, but their global motions are identical. Hinge residues mediating the cooperative movements of the alpha- and beta-domains are identified, which comprise residues in helix C, Glu35 and Ser36 on the loop succeeding helix B, Ile55 and Leu56 at the turn between strands II and III. The central part of the beta-domain long loop and the turn between strands I and II display an enhanced mobility. Finally, kinetically hot residues and key interactions are identified, which point at helix B and beta-strand III as the structural elements underlying the stability of the tertiary structure.

Access of ligands to cavities within the core of a protein is rapid

We have investigated the magnitude and timescale of fluctuations within the core of a protein using the exchange kinetics of indole and benzene binding to engineered hydrophobic cavities in T4 lysozyme. The crystal structures of variant-benzene complexes suggest that relatively large scale fluctuations (1-2 angstrom) of backbone atoms are required for entry of these ligands into the core. Nonetheless, these ligands enter the cavities rapidly, with bimolecular rate constants of approximately 10(6)-10(7) M(-1) s(-1) and a low activation barrier, 2-5 kcal mol(-1). These results suggest that protein cores undergo substantial fluctuations on the millisecond to microsecond timescale and that entry of small molecules into protein interiors is not strongly limited by steric occlusion.

Structural determinants of protein dynamics: analysis of 15N NMR relaxation measurements for main-chain and side-chain nuclei of hen egg white lysozyme

15N-labeled hen lysozyme has been studied by 2D and 3D NMR in order to characterize its dynamic behavior. The resonances of all main-chain amide nitrogen atoms were assigned, as were resonances of nitrogen atoms in 28 side chains. Relaxation measurements for the main-chain and arginine and tryptophan side-chain 15N nuclei used standard methods, and those for the 15N nuclei of asparagine and glutamine side chains used pulse sequences designed to remove unwanted relaxation pathways in the NH2 groups. The calculated order parameters (S2) show that the majority of main-chain amides undergo only small amplitude librational motions on a fast time scale (S2 > or = 0.8). Increased main-chain motion (0.5 < S2 < 0.8) is observed for a total of 19 residues located at the C-terminus, in loop and turn regions, and in the first strand of the main beta-sheet. Order parameters derived for the side chains range from 0.05 to 0.9; five of the six tryptophan residues have high order parameters (S2 > or = 0.8), consistent with their location in the closely packed core of the protein, whereas the order parameters between 0.05 and 0.3 for arginine residues confirm increased side-chain mobility at the protein surface. Order parameters for the side chains of asparagine and glutamine residues range from 0.2 to 0.8; high values are found for side chains that have low solvent accessible surfaces and well-defined chi 1 values, as measured by 3J alpha beta coupling constants. Many of the main-chain and side-chain groups with low order parameters have higher than average temperature factors in X-ray crystal structures and increased positional uncertainty in NMR solution structures. They also tend to lack persistent hydrogen bond interactions and protection against amide hydrogen exchange. The most significant correlations are found between residues with low order parameters and high surface accessibility in both crystal and solution structures. The results suggest that a lack of van der Waals contacts is a major determinant of side-chain and main-chain mobility in proteins.

Probing weakly polar interactions in cytochrome c

Theoretical, statistical, and model studies suggest that proteins are stabilized by weakly polar attractions between sulfur atoms and properly oriented aromatic rings. The two sulfur-containing amino acids, methionine and cysteine, occur frequently among functional alleles in random mutant libraries of Saccharomyces cerevisiae iso-1-cytochrome c genes at positions that form a weakly polar aromatic-aromatic interaction, the wild-type protein. To determine if a weakly polar sulfur-aromatic interaction replaced the aromatic-aromatic interaction, the structure and stability of two variants were examined. Phenylalanine 10, which interacts with tyrosine 97, was replaced by methionine and cysteine. The cysteine was modified to form the methionine and cysteine analog, S-methyl cysteine (CysSMe). Proton NMR studies indicate that changing Phe 10 to Met or CysSMe affects only local structure and that the structures of sulfur-containing variants are nearly identical. Analysis of chemical shifts and nuclear Overhauser effect data indicates that both sulfur-containing side chains are in position to form a weakly polar interaction with Tyr 97. The F10M and F10CSMe variants are 2-3 kcal mol-1 less stable than iso-1-cytochrome c at 300 K. Comparison of the stabilities of the F10M and F10CSMe variants allows evaluation of the potential weakly polar interaction between the additional sulfur atom of F10CSMe and the aromatic moiety of Tyr 97. The F10CSMe;C102T variant is 0.7 +/- 0.3 kcal mol-1 more stable than the F10M;C102T protein. The increased stability is explained by the difference in hydrophobicity of the sulfur-containing side chains. We conclude that any weakly polar interaction between the additional sulfur and the aromatic ring is too weak to detect or is masked by destabilizing contributions to the free energy of denaturation.

Changes in global stability and local structure of cytochrome c upon substituting phenylalanine-82 with tyrosine

We have examined the F82Y;C102T variant of Saccharomyces cerevisiae iso-1-cytochrome c using high-resolution proton nuclear magnetic resonance spectroscopy, chemical denaturation, and differential scanning calorimetry. Comparison of proton chemical shifts, paramagnetic shifts, and nuclear Overhauser effects indicates structural changes are localized to the vicinity of position 82. One alteration involves the rearrangement of the side chain of leucine-85. Using many more proton assignments than were available in the initial report [G. J. Pielak, R. A. Atkinson, J. Boyd, and R. J. P. Williams, Eur. J. Biochem. 177, 179-185 (1988)], a second alteration involving an interaction between arginine-13 and tyrosine-82 is observed. The interaction appears to involve a hydrogen bond with the eta-protons of arginine's guanido group acting as donor and tyrosine's phenolic eta-oxygen as acceptor. In spite of this potentially-stabilizing interaction, the free energy of denaturation decreases by approximately 2.4 kcal mol-1. Results are discussed with respect to alterations in the native and denatured states.

19F n.m.r. studies of conformational changes accompanying cyclic AMP binding to 3-fluorophenylalanine-containing cyclic AMP receptor protein from Escherichia coli

A fluorine-containing analogue of the cyclic AMP (cAMP) receptor protein (CRP) from Escherichia coli was prepared by biosynthetic incorporation of 3-fluorophenylalanine (3-F-Phe). 19F n.m.r. studies on this protein have provided direct evidence for cAMP-induced conformational changes not only within the cAMP-binding domain but also within the hinge region connecting the cAMP-binding domain to the DNA-binding headpiece. At 313 K, the 19F n.m.r. spectrum of [3-F-Phe]CRP showed five signals corresponding to the five phenylalanine residues as expected for a symmetrical dimer. Proteolysis of [3-F-Phe]CRP with subtilisin produced a fragment (the alpha-fragment) containing the cAMP-binding domain. The alpha-fragment contains all the phenylalanines except for Phe-136, a residue located in the hinge region. By comparing the 19F spectra of [3-F-Phe]CRP and its alpha-fragment, the signal for Phe-136 was assigned. The chemical shifts of the corresponding signals in the two spectra are similar, indicating that the alpha-fragment retains the structure it has in the intact protein. The largest cAMP-induced shift was observed for the signal from Phe-136 providing direct evidence for a conformational change in the hinge region. However, whereas binding of a single cAMP molecule to a CRP dimer is known to be sufficient to activate the DNA binding, the n.m.r. data indicate that the hinge region does not have the same conformation in both subunits when only one cAMP molecule is bound.

Comparison of reduced and oxidized yeast iso-1-cytochrome c using proton paramagnetic shifts

Dipolar paramagnetic shifts for protons of yeast iso-1-cytochrome c have been calculated by using an optimized g-tensor and the X-ray crystallographic coordinates of the reduced form of yeast iso-1-cytochrome c [Louie, G. V., & Brayer, G. D. (1990) J. Mol. Biol. 214, 527-555]. The calculated values are compared with the observed paramagnetic shift determined from over 450 nonequivalent protons that have been assigned in both oxidation states [Gao, Y., Boyd, J., Williams, R. J. P., & Pielak, G. J. (1990) Biochemistry 29, 6994-7003]. There is good agreement between the calculated and the experimental data with a few exceptions. This indicates that, overall, the solution structures must be very similar in both the reduced and oxidized states in solution as is the case in crystals. The differences between observed and calculated shift values for the molecule in solution are most readily explained by slight movement of the heme and certain changes in diamagnetic shift due to small rearrangements of a few residues and some considerable changes in a few hydrogen bonds. It is also known that small differences exist between the structures of the two oxidation states in crystals but the hydrogen-bond changes are not so easily observed there. Structural changes from nuclear magnetic resonance data are in reasonable agreement with those deduced from crystallography, but additional information is clearly available concerning changes in hydrogen bonding.

NMR comparison of prokaryotic and eukaryotic cytochromes c

1H NMR spectroscopy has been used to examine ferrocytochrome c-551 from Pseudomonas aeruginosa (ATCC 19429) over the pH range 3.5-10.6 and the temperature range 4-60 degrees C. Resonance assignments are proposed for main-chain and side-chain protons. Comparison of results for cytochrome c-551 to recently assigned spectra for horse cytochrome c (Wand et al. (1989) Biochemistry 28, 186-194) and mutants of yeast iso-1 cytochrome (Pielak et al. (1988) Eur. J. Biochem. 177, 167-177) reveals some unique resonances with unusual chemical shifts in all cytochromes that may serve as markers for the heme region. Results for cytochrome c-551 indicate that in the smaller prokaryotic cytochrome, all benzoid side chains are rapidly flipping on the NMR time scale. In contrast, in eukaryotic cytochromes there are some rings flipping slowly on the NMR time scale. The ferrocytochrome c-551 undergoes a transition linked to pH with a pK around 7. The pH behavior of assigned resonances provides evidence that the site of protonation is the inner or buried 17-propionic acid heme substituent (IUPAC-IUB porphyrin nomenclature). Conformational heterogeneity has been observed for segments near the inner heme propionate substituent.

13C and 15N NMR and time-resolved fluorescence depolarization study of bovine carbonic anhydrase--4-methylbenzenesulfonamide complex

The influence of the binding of the high-affinity inhibitor, 4-methylbenzenesulfonamide, to the active site of bovine carbonic anhydrase B was studied by 15N- and 13C-NMR spectroscopy. The rotational correlation time dependence on temperature and concentration of the complex was determined by time-resolved fluorescence depolarization measurements. Our experiment provides evidence that the stoichiometry of the interaction of 4-methylbenzenesulfonamide with carbonic anhydrase B is 1:1 and the inhibitor is bound in anionic form. The 15N-NMR relaxation parameters confirm our previous conclusions about the presence of librational motions in the active site of carbonic anhydrase and indicate that the internal motion in the enzyme-inhibitor complex is more restricted than the backbone motion in the uncomplexed native enzyme.

Proton-NMR studies show that the Thr-102 mutant of yeast iso-1-cytochrome c is a typical member of the eukaryotic cytochrome c family

The Thr-102 mutant of yeast iso-1-cytochrome c is a useful system for the study of structure-function relationships in this important class of electron transfer proteins, but little is known about its structure. Furthermore, few assignments of individual amino acid residues in yeast iso-1-cytochrome c have been made by proton NMR. Here we report assignments for nearly half of the amino acids in the reduced Thr-102 mutant of yeast iso-1-cytochrome c. We also report assignments for the oxidized Thr-102 mutant. While the crystal structure of the reduced iso-1-cytochrome c (N. B. not the Thr-102 mutant) has been reported, there is currently little structural information concerning its solution structure and none concerning the oxidized protein. There is also no information concerning the structure of either oxidation state of the Thr-102 mutant. Comparison of the chemical shift and NOE data for the reduced Thr-102 mutant and comparison of paramagnetic shifts for analogous residues between this mutant and horse-heart and tuna cytochromes c reveal that both the basic fold of Thr-102 yeast iso-1-cytochrome c and the region around the site of the mutation are the same as those found in the latter two proteins. It is concluded that the results from structure function studies using the Thr-102 mutant will be applicable to eukaryotic cytochrome c in general. This knowledge allows us to proceed to a description of some mutants of yeast iso-1-cytochrome c in the next paper.

Two-dimensional NMR as a probe of structural similarity applied to mutants of cytochrome c

Using site-directed mutagenesis, it is possible to prepare many mutants of a protein in a short time, and to uncover differences in function. To understand the changes in function, it is essential to understand the effect(s) of the mutation in terms of structural and dynamic changes. It is particularly important to establish a rapid method for comparing the structure of the mutants with that of the wild-type protein. We propose that a combination of overlayed and difference two-dimensional NOE spectra between the wild-type and mutant protein provide a rapid method for determination of structural similarity. The observation of differences other than those due directly to the field effects of the exchanged side chain allow both local and distant conformational changes to be assessed. Here we compare NOESY spectra from a mutant of yeast iso-1-ferrocytochrome c in which the invariant residue Phe-82 has been changed to a Tyr. We conclude that NMR can show subtle changes in protein structure. Specifically, we show the change must involve the reorientation of the side chain of Leu-85 which is proximal to the mutation. The dynamics of the aromatic side chain at position 82 are shown not to give rise to measurable differences between the wild-type and mutant protein. Structural changes are not propagated to a measurable degree in other parts of the protein.

Analysis of the aromatic 1H-NMR spectrum of the kringle 5 domain from human plasminogen. Evidence for a conserved kringle fold

A kringle 5 domain fragment from human plasminogen has been investigated by 1H-NMR spectroscopy at 300 MHz and 620 MHz. The study focuses on the kringle 5 aromatic spectrum as aromatic side chains appear to mediate the binding of benzamidine. Spin-echo experiments and acid/base-titration studies in conjunction with two-dimensional double-quantum and chemical-shift-correlated spectroscopies were used to identify individual spin systems. Sequence-specific assignments of aromatic resonances are derived from direct comparison of the kringle 5 spectrum with spectra of the homologous kringle 1 and kringle 4 domains of plasminogen. As previously observed for kringles 1 and 4, the pattern we detect for Tyr9 in kringle 5 reflects a slow conformational exchange between two states in equilibrium, one in which the Tyr9 ring is freely mobile and one in which its flip dynamics are constrained. Proton Overhauser experiments in 1H2O and in 2H2O have been used to probe aromatic ring interactions and to identify residues which are part of the hydrophobic core centered at the Leu46 side chain. Overall, the data indicate a strong structural homology among the three plasminogen kringles.

Dynamics of trimethoprim bound to dihydrofolate reductase

The conformation of a small molecule in its binding site on a protein is a major factor in the specificity of the interaction between them. In this paper, we report the use of 1H and 13C NMR spectroscopy to study the fluctuations in conformation of the anti-bacterial drug trimethoprim when it is bound to its "target," dihydrofolate reductase. 13C relaxation measurements reveal dihedral angle changes of +/- 25 degrees to +/- 35 degrees on the subnanosecond time scale, while 13C line-shape analysis demonstrates dihedral angle changes of at least +/- 65 degrees on the millisecond time scale. 1H NMR shows that a specific hydrogen bond between the inhibitor and enzyme, which is believed to make an important contribution to binding, makes and breaks rapidly at room temperature.

A unique pair of zinc binding sites in the human alpha 2-macroglobulin tetramer. A 35Cl and 37Cl NMR study

35Cl NMR has been used to demonstrate that human alpha 2-macroglobulin tetramer possesses a unique pair of zinc binding sites. Zinc bound at these sites does not affect the 35Cl NMR line width of free Cl-. Additional lower affinity zinc sites exist that bind chloride weakly and cause broadening of the free chloride resonance through fast exchange with bound chloride. Using both 35Cl and 37Cl relaxation measurements it has been shown that chloride bound at these sites has an internal correlation time of 5.1 ns and a quadrupolar interaction, chi, of 4.2 MHz with zinc. Manganese binds to apo-alpha 2-macroglobulin analogously to zinc. alpha 2-macroglobulin that has been reacted with methylamine still possesses two classes of zinc sites per tetramer, but their relative affinities differ more than for unreacted alpha 2-macroglobulin. These data are discussed with respect to possible models for the subunit arrangement in the tetramer.

Influence of inhibitor binding on the internal motions of lysozyme

Time-resolved laser-induced fluorescence depolarization measurements of internal motions in lysozyme are presented. The fluorescent dye eosin binds in a one-to-one complex with the enzyme, and is used both to measure the overall tumbling time constants and to probe the motions of residues in the region of binding. The precision and accuracy of the present method for determining the overall tumbling time constants compare favorably with those from other methods used in the literature. The extent of the internal motions, as described by a model independent order parameter, S2, is temperature dependent, and changes when the inhibitor N,N',N"-triacetylchitotriose, (GlcNAc)3, is bound to the active site of the enzyme. The observed temperature dependence and changes in S2 upon binding of (GlcNAc)3 are interpreted in terms of a nonharmonic model of the effective potential that is consistent with the picture of concerted motions in the protein. The values of the parameters of the potential that reproduce the data with and without the bound inhibitor imply that (GlcNAc)3 binding causes an increase in the rigidity of the protein, which agree qualitatively with other results on the lysozyme-(GlcNAc)3 system.

Surface accessibility of aromatic residues in human lysozyme using photochemically induced dynamic nuclear polarization NMR spectroscopy

Resonances in the photo-CIDNP spectrum of human lysozyme have been assigned to specific spin systems despite extensive spectral overlap using the two-dimensional photo-CIDNP COSY experiment. Five of the 12 tyrosine, tryptophan and histidine residues of human lysozyme are found to be accessible to flavin dye in solution. This result is in good agreement with surface accessibility calculations carried out on the human lysozyme crystal structure. When amino acid differences are considered the photo-CIDNP results obtained for human lysozyme are in good agreement with results obtained previously for hen lysozyme.

19F-n.m.r. studies of 3',5'-difluoromethotrexate binding to Lactobacillus casei dihydrofolate reductase. Molecular motion and coenzyme-induced conformational changes

19F-n.m.r. spectroscopy was used to study the binding of 3',5'-difluoromethotrexate to dihydrofolate reductase (tetrahydrofolate dehydrogenase) from Lactobacillus casei. The benzoyl ring of the bound difluoromethotrexate was found to 'flip' about its symmetry axis, and the rate (7.3 X 10(3) s-1 at 298 K) and activation parameters for this process were determined by lineshape analysis of the 19F-n.m.r. spectrum at a series of temperatures in the range 273-308 K. The contributions to the barrier for this process are discussed. Addition of NADP+ or NADPH to form the enzyme-difluoromethotrexate-coenzyme ternary complex led to an increase in the rate of benzoyl ring flipping by a factor of 2.6-2.8-fold, and to substantial changes in the 19F-n.m.r. chemical shifts. The possible nature of the coenzyme-induced conformational changes responsible for these effects is discussed.

Pyrrolidine ring conformations in prolyl peptides from 13C spin-lattice relaxation times

A method was developed in the framework of a bistable jump model to obtain the pyrrolidine ring conformations in proline peptides from 13C spin-lattice relaxation times. Equations are presented expressing the ring torsions in terms of the 13C spin-lattice relaxation times of the ring carbons. This method was applied to 26 pyrrolidine ring systems and acceptable conformations were obtained.

The assignment of the 1H nuclear magnetic resonance spectrum of azurin. An investigation of the 1H NMR spectrum of the blue copper protein, azurin, from Pseudomonas aeruginosa, with reference to the previously determined crystal structure

A detailed assignment of the 1H nuclear magnetic resonance spectrum of azurin has been made. Resonances associated with the single tryptophan residue, all six phenylalanine residues, one of the two tyrosine residues and all four histidine residues, as well as most of the resonances from the ring-current shifted methyl groups have been assigned. These assignments have been used to study the pH dependence of the structure of the protein and binding of analogues of redox-active reagents to the protein.

Structural and functional aspects of domain motions in proteins

Three distinct categories of large-scale flexibility in proteins have been documented by single-crystal X-ray diffraction studies: the relatively free movement of essentially rigid globular domains that are connected by a flexible segment of polypeptide, the reorientation of essentially rigid domains among a few distinct conformations, and the concerted transition of a contiguous region of the surface of a protein from a disordered state to an ordered state. In a number of examples, well-defined functions can be assigned to these large-scale structural changes. The occurrence of such motions in proteins of known structure is reviewed, and the best-studied examples are discussed in detail to allow a critical evaluation of the methods used to identify and study these motions.

Nuclear magnetic resonance studies on calmodulin: spectral assignments in the calcium-free state

The 400-MHz proton magnetic resonance spectra of calcium-free scallop testis calmodulin (CaM) and pig brain CaM were observed. Detailed spectral assignments were made by resolution enhancement, spin decoupling, and nuclear Overhauser enhancement (NOE) experiments as well as pH titration. Comparison between spectra of scallop testis CaM and pig brain CaM were also utilized for the assignment. Previous assignments for tyrosine-99, histidine-107, epsilon-trimethyllysine-115, and tyrosine-138 [Seamon, K. B. (1980) Biochemistry 19, 207; Krebs, J., & Carafoli, E. (1982) Eur. J. Biochem. 124, 619] were confirmed. Phenylalanine-99 and threonine-143 of scallop testis CaM were identified. Sixteen methyl resonances from one isoleucine, two valines, nine methionines, and the amino-terminal acetyl group were identified. First-stage assignments were made of resonances arising from seven phenylalanines. The uniquely high field shifted phenylalanine resonance previously reported by Seamon was found to consist of two doublets from the two pairs of delta protons of two phenylalanines. The NOE experiments showed that the two phenylalanines are located closely to each other. The large high-field shifts of these phenylalanines were accounted for the ring-current effects due to their proximity. An isoleucine and a valine of which methyl resonances appear at high fields were found to be situated closely to each other. It was found that two delta protons and two epsilon protons of almost all aromatic residues are magnetically equivalent, suggesting that the local structure of aromatic residues is so flexible as to permit the rapid flipping motion of the ring.

Complete assignment of the 1H NMR spectrum of the aromatic residues of lysozyme

Assignment of the nuclear magnetic resonance (NMR) spectrum of the 59 non-exchangeable protons of the 13 aromatic amino acid residues of lysozyme has been completed using a combination of selective spin decoupling and nuclear Overhauser experiments. The experimental chemical shift data were used to simulate the aromatic region of the NMR spectrum at 300 MHz and at 470 MHz. Excellent agreement with the experimental spectra was obtained and the simulations permitted more accurate chemical shift values and resonance linewidths to be obtained. The current set of assignments is compared with those of the assignments made previously. Evidence for the motional behaviour of the aromatic amino acid side chains is presented.

Hydrogen exchange and the dynamic structure of proteins

In native proteins, buried, labile protons undergo isotope exchange with solvent hydrogens, but the kinetics of exchange are markedly slower than in unfolded polypeptides. This indicates that, whereas buried protein atoms are shielded from solvent, the protein fluctuates around the time average structure and occasionally exposes buried sites to solvent. Generally, hydrogen exchange studies are designed to characterize the nature of the fluctuations between conformational substates, to monitor the shift in conformational equilibria among protein substates due to ligand binding or other factors, or to monitor the major cooperative denaturation transition. In this article, we review the recent reports of hydrogen exchange in proteins, focusing on recent advances in methodology, especially with regard to the implications of the results for the mechanism of hydrogen exchange in folded proteins.

Structural role of the tyrosine residues of cytochrome c

The tertiary structures of horse, tuna, Neurospora crassa, horse [Hse65,Leu67]- and horse [Hse65,Leu74]-cytochromes c were studied with high-resolution 1H n.m.r. spectroscopy. The amino acid sequences of these proteins differ at position 46, which is occupied by phenylalanine in the horse proteins but by tyrosine in the remaining two, and at positions 67, 74 and 97, which are all occupied by tyrosine residues in horse and tuna cytochrome c but in the other proteins are substituted by phenylalanine or leucine, though there is only one such substitution per protein. The various aromatic-amino-acid substitutions do not seriously affect the protein structure.