1
|
Kulkarni M, Söderhjelm P. Free-Energy Landscape and Rate Estimation of the Aromatic Ring Flips in Basic Pancreatic Trypsin Inhibitors Using Metadynamics. J Chem Theory Comput 2023; 19:6605-6618. [PMID: 37698852 PMCID: PMC10569046 DOI: 10.1021/acs.jctc.3c00460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2023] [Indexed: 09/13/2023]
Abstract
Aromatic side chains (phenylalanine and tyrosine) of a protein flip by 180° around the Cβ-Cγ axis (χ2 dihedral of the side chain), producing two symmetry-equivalent states. The study of ring flip dynamics with nuclear magnetic resonance (NMR) experiments helps to understand local conformational fluctuations. Ring flips are categorized as slow (milliseconds and onward) or fast (nanoseconds to near milliseconds) based on timescales accessible to NMR experiments. In this study, we investigated the ability of the infrequent metadynamics approach to estimate the flip rate and discriminate between slow and fast ring flips for eight individual aromatic side chains (F4, Y10, Y21, F22, Y23, F33, Y35, and F45) of the basic pancreatic trypsin inhibitor. Well-tempered metadynamics simulations were performed to estimate the ring-flipping free-energy surfaces for all eight aromatic residues. The results indicate that χ2 as a standalone collective variable (CV) is not sufficient to obtain computationally consistent results. Inclusion of a complementary CV, such as χ1(Cα-Cβ), solved the problem for most residues and enabled us to classify fast and slow ring flips. This indicates the importance of librational motions in ring flips. Multiple pathways and mechanisms were observed for residues F4, Y10, and F22. Recrossing events were observed for residues F22 and F33, indicating a possible role of friction effects in ring flipping. The results demonstrate the successful application of infrequent metadynamics to estimate ring flip rates and identify certain limitations of the approach.
Collapse
Affiliation(s)
- Mandar Kulkarni
- Division of Biophysical Chemistry, Lund University, Chemical Center, 22100 Lund, Sweden
| | - Pär Söderhjelm
- Division of Biophysical Chemistry, Lund University, Chemical Center, 22100 Lund, Sweden
| |
Collapse
|
2
|
Weininger U, Brath U, Modig K, Teilum K, Akke M. Off-resonance rotating-frame relaxation dispersion experiment for 13C in aromatic side chains using L-optimized TROSY-selection. JOURNAL OF BIOMOLECULAR NMR 2014; 59:23-9. [PMID: 24706175 PMCID: PMC4003406 DOI: 10.1007/s10858-014-9826-2] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/26/2013] [Accepted: 03/25/2014] [Indexed: 05/04/2023]
Abstract
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.
Collapse
Affiliation(s)
- Ulrich Weininger
- Department of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, 22100 Lund, Sweden
| | - Ulrika Brath
- Department of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, 22100 Lund, Sweden
- Present Address: Department of Chemistry and Molecular Biology, University of Gothenburg, 41296 Göteborg, Sweden
| | - Kristofer Modig
- Department of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, 22100 Lund, Sweden
| | - Kaare Teilum
- Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, 2200 Copenhagen, Denmark
| | - Mikael Akke
- Department of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, 22100 Lund, Sweden
| |
Collapse
|
3
|
Hattori M, Li H, Yamada H, Akasaka K, Hengstenberg W, Gronwald W, Kalbitzer HR. Infrequent cavity-forming fluctuations in HPr from Staphylococcus carnosus revealed by pressure- and temperature-dependent tyrosine ring flips. Protein Sci 2005; 13:3104-14. [PMID: 15557257 PMCID: PMC2287304 DOI: 10.1110/ps.04877104] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
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.
Collapse
Affiliation(s)
- Mineyuki Hattori
- Department of Molecular Science, Graduate School of Science and Technology, Kobe University, Japan
| | | | | | | | | | | | | |
Collapse
|
4
|
Mummert ME, Voss EW. Transition-state theory and secondary forces in antigen--antibody complexes. Biochemistry 1996; 35:8187-92. [PMID: 8679572 DOI: 10.1021/bi9604791] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Secondary forces, defined as those interactions between the antigen (epitope including the surrounding environment) and areas immediately adjacent to the antibody active site, were investigated using monofluorescein-derivatized synthetic peptides of varying electrostatic properties. Secondary forces were quantitated by measuring the unimolecular rate constants at two different temperatures using the high-affinity anti-fluorescein monoclonal antibody 4-4-20 complexed with fluorescein-derivatized synthetic peptides. Unimolecular rate constants were correlated with transition-state theory to explain secondary effects. An acidic peptide produced a large temperature-dependent effect upon binding including a significant enthalpic factor (+33.28 kcal/mol) relative to the binding of fluorescein ligand (+23.96 kcal/mol). Binding of a basic peptide produced both a relatively smaller temperature effect and enthalpy factor than fluorescein ligand. The antibody-ligand binding results were interpreted invoking the concepts of thermally averaged metatypic (liganded) states of the antibody as well as potential biochemical interactions between the antigen and accessible surface regions of the antibody's complementarity determining regions.
Collapse
Affiliation(s)
- M E Mummert
- Department of Microbiology, University of Illinois, Urbana 61801-3704, USA
| | | |
Collapse
|
5
|
Brunne RM, van Gunsteren WF. Dynamical properties of bovine pancreatic trypsin inhibitor from a molecular dynamics simulation at 5000 atm. FEBS Lett 1993; 323:215-7. [PMID: 7684708 DOI: 10.1016/0014-5793(93)81342-w] [Citation(s) in RCA: 40] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Molecular dynamics simulations of bovine pancreatic trypsin inhibitor in water have been performed with coupling to pressure baths at 1 atm and at 5000 atm. The positional fluctuations of atoms in the alpha- and gamma-positions are slightly decreased at 5000 atm. The mobility of the backbone phi- and psi-angles is not affected with respect to the root mean square fluctuations and the rate of torsional angle transitions. The amplitude of libration of sidechain chi-angles remains nearly the same for both pressures, but the rate of torsional angle transitions decreases on average by 30% when increasing the hydrostatic pressure to 5000 atm.
Collapse
Affiliation(s)
- R M Brunne
- Laboratorium für Physikalische Chemie, Eidgenössische Technische Hochschule, ETH Zentrum, Zürich, Switzerland
| | | |
Collapse
|
6
|
Jona I, Martonosi A. The effect of high pressure on the conformation, interactions and activity of the Ca(2+)-ATPase of sarcoplasmic reticulum. BIOCHIMICA ET BIOPHYSICA ACTA 1991; 1070:355-73. [PMID: 1837234 DOI: 10.1016/0005-2736(91)90077-l] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
High pressure (100-150 MPa) increases the intensity and polarization of fluorescence of FITC-labeled Ca(2+)-ATPase in a medium containing 0.1 mM Ca2+, suggesting a reversible pressure-induced transition from the E1 into an E2-like state with dissociation of ATPase oligomers. Under similar conditions but using unlabeled sarcoplasmic reticulum vesicles, high pressure caused the reversible release of Ca2+ from the high-affinity Ca2+ sites of Ca(2+)-ATPase, as indicated by changes in the fluorescence of the Ca2+ indicator, Fluo-3; this was accompanied by reversible inhibition of the Ca(2+)-stimulated ATPase activity measured in a coupled enzyme system of pyruvate kinase and lactate dehydrogenase, and by redistribution of Prodan in the lipid phase of the membrane, as shown by marked changes in its fluorescence emission characteristics. In a Ca(2+)-free medium where the equilibrium favors the E2 conformation of Ca(2+)-ATPase the fluorescence intensity of FITC-ATPase was not affected or only slightly reduced by high pressure. The enhancement of TNP-AMP fluorescence by 100 mM inorganic phosphate in the presence of EGTA and 20% dimethylsulfoxide was essentially unaffected by 150 MPa pressure at pH 6.0 and was only slightly reduced at pH 8.0. As the enhancement of TNP-AMP fluorescence by Pi is associated with the Mg(2+)-dependent phosphorylation of the enzyme and the formation of Mg.E2-P intermediate, it appears that the reactions of Ca(2+)-ATPase associated with the E2 state are relatively insensitive to high pressure. These observations suggest that high pressure stabilizes the enzyme in an E2-like state characterized by low reactivity with ATP and Ca2+ and high reactivity with Pi. The transition from the E1 to the E2-like state involves a decrease in the effective volume of Ca(2+)-ATPase.
Collapse
Affiliation(s)
- I Jona
- Department of Biochemistry and Molecular Biology, State University of New York, Syracuse 13210
| | | |
Collapse
|
7
|
Abstract
The effect of hydrostatic pressure (0-2.6 kbar) on the acrylamide quenching of the fluorescence of indole derivatives and several single-tryptophan-containing proteins has been studied using phase fluorometry at 25 degrees C. For the model system, N-acetyl-L-tryptophanamide in water, there is essentially no pressure dependence of the quenching rate constant, kappa q. For the internal Trp residue of ribonuclease T1 and cod parvalbumin, there also is essentially no pressure dependence of the apparent kappa q at low pressure. Thus, the activation volume, delta V not equal to, for these quenching processes is approximately zero. Such small delta V not equal to values are expected for diffusion-limited reactions in water at this temperature. The low, apparent delta V not equal to values for the globular proteins characterize these quenching processes as involving very small amplitude fluctuations in the protein structures. Only for the poised tetramer in equilibrium monomer equilibrium of melittin were we able to observe a significant effect of pressure on kappa q and this is due to the pressure-induced shift in the equilibrium position.
Collapse
Affiliation(s)
- M R Eftink
- Department of Chemistry, University of Mississippi, University 38655
| | | |
Collapse
|
8
|
Barman T, Travers F, Balny C, Hui Bon Hoa G, Douzou P. New trends in cryoenzymology: probing the functional role of protein dynamics by single-step kinetics. Biochimie 1986; 68:1041-51. [PMID: 3096382 DOI: 10.1016/s0300-9084(86)80178-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Cryoenzymology was initially used to slow down enzyme-catalyzed reactions so as to stabilize intermediates for further study. During the course of this early work, it became clear that cryoenzymology could be extended to other ends and some of these are described. First, the use of a cryosolvent on its own (or together with temperature) as a perturbant has allowed a resolution of the substrate binding steps of certain enzymes (myosin, D-amino acid oxidase, peroxidase and cytochrome P450). Second, by the use of cryosolvent and temperature, coupled with the classical physico-chemical perturbants, one can selectively modulate the various steps of an enzyme pathway. This approach can lead to an understanding of the mechanism of enzyme regulation. Finally, by carrying out experiments over a wide range of temperatures (-30 degrees C- +30 degrees C) and pressure (up to several kbars) in specially constructed fast reaction equipment, one can study the thermodynamic properties of the individual rate constants describing the interconversions of reaction intermediates. Experiments with creatine kinase, cytochrome P450 and peroxidase are described. The thermodynamic parameters delta H, delta G, delta S and delta V are thus measured and when this is done under different solvent conditions one can, at least within the theories available, attempt an approach to the problem of protein dynamics.
Collapse
|
9
|
Wagner G, Wüthrich K. Observation of internal motility of proteins by nuclear magnetic resonance in solution. Methods Enzymol 1986; 131:307-26. [PMID: 3773764 DOI: 10.1016/0076-6879(86)31047-4] [Citation(s) in RCA: 39] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
|
10
|
Wagner G. Characterization of the distribution of internal motions in the basic pancreatic trypsin inhibitor using a large number of internal NMR probes. Q Rev Biophys 1983; 16:1-57. [PMID: 6878622 DOI: 10.1017/s0033583500004911] [Citation(s) in RCA: 191] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
The experimental observations described in this article indicated that a distribution of many different fluctuations is present in a globular protein. These fluctuations were characterized by observation of many natural internal probes such as the labile peptide protons and the aromatic side chains. The conditions which are necessary to get reactions of the internal probes have been discussed in detail. The structural interpretation of the data was facilitated by the development and the use of new NMR techniques which provided the identification of the resonances of all the labile peptide protons. With NOE measurements a distinction between correlated and uncorrelated exchange events was obtained. This enabled us to elucidate the exchange mechanism over a wide range of p2H and temperature and to classify different subsets of fluctuations with respect to their lifetimes. It was further demonstrated that a change of external conditions such as temperature, p2H or pressure can change the distribution of fluctuations in the protein. The mechanisms responsible for rotation of internal aromatic side chains were also found to change with temperature, and mechanistic aspects of these fluctuations were discussed. This demonstration of a manifold of spatial fluctuations in a small protein provides an impression on the kind of fluctuations which have to be expected for larger proteins. When studying protein reactions one should therefore consider the presence of a large number of different, transiently formed, spatial structures available for the partner in the reaction, which may pick out only that structure which will optimally perform a particular reaction with the highest efficiency.
Collapse
|