Conformational dynamics of chymotrypsin inhibitor 2 by coarse-grained simulations

A coarse-grained dynamic Monte Carlo (MC) simulation method is used to investigate the conformational dynamics of chymotrypsin inhibitor 2 (CI2). Each residue is represented therein by two interaction sites, one at the alpha-carbon and the other on the amino acid side-chain. The energy and geometry parameters extracted from databank structures are used. The calculated rms fluctuations of alpha-carbon atoms are in good agreement with crystallographic temperature factors. The two regions of the protein that pack against each other to form the main hydrophobic core exhibit negatively correlated fluctuations. The conformational dynamics could efficiently be probed by the time-delayed orientational and conformational correlation functions of the virtual bonds: the active site loop, excluding the active site bond, the turn region, and the N-terminal of the alpha-helix are relatively more mobile regions of the structure. A correlation is observed between the hydrogen/deuterium (H/D) exchange behavior and the long-time orientational and conformational autocorrelation function values for CI2. A cooperativity in the rotations of the bonds near in sequence is observed at all time windows, whereas the cooperative rotations of the bonds far along the sequence appear at long time windows; these correlations contribute to the stability of the secondary structures and the tertiary structure, respectively.     

Conformational dynamics of subtilisin-chymotrypsin inhibitor 2 complex by coarse-grained simulations

An off-lattice dynamic Monte Carlo (MC) method is used to investigate the conformational dynamics of chymotrypsin inhibitor 2 (CI2) and subtilisin in both free and complex forms over two time windows, referring to short and long time scales. The conformational dynamics of backbone bonds analysed from several independent trajectories reveal that: Both the inhibitor and the enzyme are restricted in their bond rotations, excluding a few bonds, upon binding; the effect being greatest for the loop regions, and for the inhibitor. A cooperativity in the near-neighbor bond rotations are observed on both time scales, whereas the cooperative rotations of the bonds far along the sequence appear only in the long time window, and the latter time window is where most of the interactions between the inhibitor and the enzyme are observed. Upon binding, the cooperatively rotating parts of the inhibitor and the enzyme are readjusted compared to their free forms, and new correlations appear. The binding loop, although it is the closest contact region, is not the only part of the inhibitor involved in the interactions with the enzyme. Loops 3 and 8 and the helices F and G in bound enzyme and the binding loop of the inhibitor contribute at the most to the collective motions of whole structure on the slow time scale and are apparently important for enzyme-inhibitor interactions and function. The results in general provide evidence for the contribution of the loops with cooperative motions to the extensive communication network of the complex.     

Conformational Dynamics of Subtilisin-Chymotrypsin Inhibitor 2 Complex by Coarse-Grained Simulations

Abstract

An off-lattice dynamic Monte Carlo (MC) method is used to investigate the conformational dynamics of chymotrypsin inhibitor 2 (CI2) and subtilisin in both free and complex forms over two time windows, referring to short and long time scales. The conformational dynamics of backbone bonds analysed from several independent trajectories reveal that: Both the inhibitor and the enzyme are restricted in their bond rotations, excluding a few bonds, upon binding; the effect being greatest for the loop regions, and for the inhibitor. A cooperativity in the near-neighbor bond rotations are observed on both time scales, whereas the cooperative rotations of the bonds far along the sequence appear only in the long time window, and the latter time window is where most of the interactions between the inhibitor and the enzyme are observed. Upon binding, the cooperatively rotating parts of the inhibitor and the enzyme are readjusted compared to their free forms, and new correlations appear. The binding loop, although it is the closest contact region, is not the only part of the inhibitor involved in the interactions with the enzyme. Loops 3 and 8 and the helices F and G in bound enzyme and the binding loop of the inhibitor contribute at the most to the collective motions of whole structure on the slow time scale and are apparently important for enzyme-inhibitor interactions and function. The results in general provide evidence for the contribution of the loops with cooperative motions to the extensive communication network of the complex.     

Hydrogen peroxide-induced structural alterations of RNAse A

Proteins exposed to oxidative stress are degraded via proteolytic pathways. In the present study, we undertook a series of in vitro experiments to establish a correlation between the structural changes induced by mild oxidation of the model protein RNase A and the proteolytic rate found upon exposure of the modified protein toward the isolated 20 S proteasome. Fourier transform infrared spectroscopy was used as a structure-sensitive probe. We report here strong experimental evidence for oxidation-induced conformational rearrangements of the model protein RNase A and, at the same time, for covalent modifications of amino acid side chains. Oxidation-related conformational changes, induced by H(2)O(2) exposure of the protein may be monitored in the amide I region, which is sensitive to changes in protein secondary structure. A comparison of the time- and H(2)O(2) concentration-dependent changes in the amide I region demonstrates a high degree of similarity to spectral alterations typical for temperature-induced unfolding of RNase A. In addition, spectral parameters of amino acid side chain marker bands (Tyr, Asp) revealed evidence for covalent modifications. Proteasome digestion measurements on oxidized RNase A revealed a specific time and H(2)O(2) concentration dependence; at low initial concentration of the oxidant, the RNase A turnover rate increases with incubation time and concentration. Based on these experimental findings, a correlation between structural alterations detected upon RNase A oxidation and proteolytic rates of RNase A is established, and possible mechanisms of the proteasome recognition process of oxidatively damaged proteins are discussed.     

Angular variances for internal bond rotations of side chains in GXG-based tripeptides derived from (13)C-NMR relaxation measurements: Implications to protein folding

The study of backbone and side-chain internal motions in proteins and peptides is crucial to having a better understanding of protein/peptide "structure" and to characterizing unfolded and partially folded states of proteins and peptides. To achieve this, however, requires establishing a baseline for internal motions and motional restrictions for all residues in the fully, solvent-exposed "unfolded state." GXG-based tripeptides are the simpliest peptides where residue X is fully solvent exposed in the context of an actual peptide. In this study, a series of GXG-based tripeptides has been synthesized with X being varied to include all twenty common amino acid residues. Proton-coupled and -decoupled (13)C-nmr relaxation measurements have been performed on these twenty tripeptides and various motional models (Lipari-Szabo model free approach, rotational anisotropic diffusion, rotational fluctuations within a potential well, rotational jump model) have been used to analyze relaxation data for derivation of angular variances and motional correlation times for backbone and side-chain chi(1) and chi(2) bonds and methyl group rotations. At 298 K, backbone motional correlation times range from about 50 to 85 ps, whereas side-chain motional correlation times show a much broader spread from about 18 to 80 ps. Angular variances for backbone phi,psi bond rotations range from 11 degrees to 23 degrees and those for side chains vary from 5 degrees to 24 degrees for chi(1) bond rotations and from 5 degrees to 27 degrees for chi(2) bond rotations. Even in these peptide models of the "unfolded state," side-chain angular variances can be as restricted as those for backbone and beta-branched (valine, threonine, and isoleucine) and aromatic side chains display the most restricted motions probably due to steric hinderence with backbone atoms. Comparison with motional data on residues in partially folded, beta-sheet-forming peptides indicates that side-chain motions of at least hydrophobic residues are less restricted in the partially folded state, suggesting that an increase in side-chain conformational entropy may help drive early-stage protein folding. Copyright 1999 John Wiley & Sons, Inc.     

Mechanism of action of Moloney murine leukemia virus RNA-directed DNA polymerase associated RNase H (RNase H I)

The mechanism of action of the ribonuclease H (RNase H) activity associated with Moloney murine leukemia virus RNA-directed DNA polymerase (RNase H I) and the two-subunit (alpha beta) form of avian myeloblastosis virus DNA polymerase were compared by utilizing the model substrate (A)n.(dT)n and polyacrylamide gel electrophoresis in 7 M urea to analyze digestion products. Examination on 25% polyacrylamide gels revealed that a larger proportion of the RNase H I oligonucleotide products generated by limited digestion of [3H](A)(1100).(dT)n were acid insoluble (15-26 nucleotides long) than acid soluble (less than 15 nucleotides long), while the opposite was true for products generated by alpha beta RNase H. RNase H I was capable of attacking RNA in RNA.DNA in the 5' to 3' and 3' to 5' directions, as demonstrated by the use of [3H,3'- or 5'-32P](A)(380).(dT)n and cellulose--[3H](A)n.(dT)n. Both RNase H I and alpha beta RNase H degraded [3H]-(A)n.(dT)n with a partially processive mechanism, based upon classical substrate competition experiments and analyses of the kinetics of degradation of [3H,3'- or 5'-32P](A)(380).(dT)n. That is, both enzymes remain bound to a RNA.DNA substrate through a finite number of hydrolytic events but dissociate before the RNA is completely degraded. Both RNase H I and alpha beta RNase H were capable of degrading [14C](A)n in [3H](C)n-[14C](A)n-[32P](dA)n.(dT)n, suggesting that retroviral RNase H is capable of removing the tRNA primer at the 5' terminus of minus strand DNA at the appropriate time during retroviral DNA synthesis in vitro.     

Two-phase unfolding pathway of ribonuclease A during denaturation induced by dithiothreitol

The dynamics of the unfolding process of bovine pancreatic ribonuclease A (RNase A) unfolded by dithiothreitol (DTT) at a low concentration of 1:30 were investigated in alkaline phosphate-buffered saline solutions at 303K and 313K by using proton nuclear magnetic resonance ((1)H NMR) spectra. Three NMR spectral parameters including Shannon entropy, mutual information, and correlation coefficient were introduced into the analysis. The results show that the unfolding process of RNase A was slowed to the order of many hours, and the kinetics of the unfolding pathway described by the three parameters is best fit by a model of two consecutive first-order reactions. Temperature greatly influences the rate constants of the unfolding kinetics with different temperature effects observed for the fast and the slow processes. The consistencies and the differences between the three sets of parameters show that they reflect the same relative denaturation pathway but different spectra windows of the unfolding process of RNase A. The results suggest that the unfolding process of RNase A induced by low concentrations of DTT is a two-phase pathway containing fast and slow first-order reactions.     

Molecular motions and dynamics of a diunsaturated acyl chain in a lipid bilayer: implications for the role of polyunsaturation in biological membranes.

The nature and dynamics of the motions of a diunsaturated fatty acyl chain in a lipid bilayer were examined using a comprehensive simulation program for 2H NMR line shapes developed by Wittebort et al. [Wittebort, R. J., Olejniczak, E. T., & Griffin, R. G. (1987) J. Chem. Phys. 36, 5411-5420]. A motional model in which the isolinoleoyl chain (18:2 delta 6,9) adopts two conformations consistent with the low energy structures proposed for 1,4-pentadiene [Applegate, K. R., & Glomset, J. A. (1986) J. Lipid Res. 27, 658-680], but undergoes a rapid jump between these states, is sufficient to account for the experimentally observed quadrupolar couplings, the 2H-2H and 1H-2H dipolar couplings, the longitudinal relaxation times, and the changes in the average conformation of the chain that occur with a variation in temperature. The jump motion originates via rotations about the C7-C8 and the C8-C9 carbon bonds and leads to the low order parameters assigned to the C8 methylene segment (0.18) and the C9-C10 double bond (0.11). In contrast, the C6-C7 double bond, which is not involved in the two-site jump, characterized by a relatively large order parameter (0.56). Fatty acyl chains containing three or more double bonds likely cannot undergo the same jump motion and consequently will be highly ordered structures. Correlation times for diffusion of the molecular long axis of the diunsaturated acyl chain about the bilayer normal (approximately 10(-10) s) and for the local jump motion (approximately 10(-10) s) were calculated.(ABSTRACT TRUNCATED AT 250 WORDS)     

ESR spectral analysis of the molecular motion of spin labels in lipid bilayers and membranes based on a model in terms of two angular motional parameters and rotational correlation times.

Electron spin resonance (ESR) spectral line shapes are calculated for a nitroxide spin-labeled molecule undergoing rapid restricted rotations (twisting) about its long molecular axis while simultaneously tumbling within a cone. Explicit expressions are derived for the hyperfine splittings and g-values, as well as for the secular contributions to the motionally modulated linewidths. The present model is useful for analyzing the restricted twisting and tumbling motions, and rotational correlation times, of spin-labeled molecules in bilayers. Simulated spectra compare well with experimental spectra of lecithin bilayers marked with cholestane spin label, over a wide temperature range.     

Effects of substrate binding on the dynamics of RNase A: Molecular dynamics simulations of UpA bound and native RNase A

RNase A has been extensively used as a model protein in several biophysical and biochemical studies. Using the available structural and biochemical results, RNase A-UpA interaction has been computationally modeled at an atomic level. In this study, the molecular dynamics (MD) simulations of native and UpA bound RNase A have been carried out. The gross dynamical behavior and atomic fluctuations of the free and UpA bound RNase A have been characterized. Principal component analysis is carried out to identify the important modes of collective motion and to analyze the changes brought out in these modes of RNase A upon UpA binding. The hydrogen bonds are monitored to study the atomic details of RNase A-UpA interactions and RNase A-water interactions. Based on these analysis, the stability of the free and UpA bound RNase A are discussed. © 1997 John Wiley & Sons, Inc. Biopoly 42: 505–520, 1997     

RNase H is an exo- and endoribonuclease with asymmetric directionality,depending on the binding mode to the structural variants of RNA:DNA hybrids

RNase H is involved in fundamental cellular processes and is responsible for removing the short stretch of RNA from Okazaki fragments and the long stretch of RNA from R-loops. Defects in RNase H lead to embryo lethality in mice and Aicardi-Goutieres syndrome in humans, suggesting the importance of RNase H. To date, RNase H is known to be a non-sequence-specific endonuclease, but it is not known whether it performs other functions on the structural variants of RNA:DNA hybrids. Here, we used Escherichia coli RNase H as a model, and examined its catalytic mechanism and its substrate recognition modes, using single-molecule FRET. We discovered that RNase H acts as a processive exoribonuclease on the 3′ DNA overhang side but as a distributive non-sequence-specific endonuclease on the 5′ DNA overhang side of RNA:DNA hybrids or on blunt-ended hybrids. The high affinity of previously unidentified double-stranded (ds) and single-stranded (ss) DNA junctions flanking RNA:DNA hybrids may help RNase H find the hybrid substrates in long genomic DNA. Our study provides new insights into the multifunctionality of RNase H, elucidating unprecedented roles of junctions and ssDNA overhang on RNA:DNA hybrids.     

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.     

ESR spectral analysis of the molecular motion of spin labels in lipid bilayers and membranes based on a model in terms of two angular motional parameters and rotational correlation times

Electron spin resonance (ESR) spectral line shapes are calculated for a nitroxide spin-labeled molecule undergoing rapid restricted rotations (twisting) about its long molecular axis while simultaneously tumbling within a cone. Explicit expressions are derived for the hyperfine splittings and $\text{g-}\text{values}$, as well as for the secular contributions to the motionally modulated linewidths. The present model is useful for analyzing the restricted twisting and tumbling motions, and rotational correlation times, of spin-labeled molecules in bilayers. Simulated spectra compare well with experimental spectra of lecithin bilayers marked with cholestane spin label, over a wide temperature range.     

The transmembrane 7-alpha-bundle of rhodopsin: distance geometry calculations with hydrogen bonding constraints.

A 3D model of the transmembrane 7-alpha-bundle of rhodopsin-like G-protein-coupled receptors (GPCRs) was calculated using an iterative distance geometry refinement with an evolving system of hydrogen bonds, formed by intramembrane polar side chains in various proteins of the family and collectively applied as distance constraints. The alpha-bundle structure thus obtained provides H bonding of nearly all buried polar side chains simultaneously in the 410 GPCRs considered. Forty evolutionarily conserved GPCR residues form a single continuous domain, with an aliphatic "core" surrounded by six clusters of polar and aromatic side chains. The 7-alpha-bundle of a specific GPCR can be calculated using its own set of H bonds as distance constraints and the common "average" model to restrain positions of the helices. The bovine rhodopsin model thus determined is closely packed, but has a few small polar cavities, presumably filled by water, and has a binding pocket that is complementary to 11-cis (6-s-cis, 12-s-trans, C = N anti)-retinal or to all-trans-retinal, depending on conformations of the Lys296 and Trp265 side chains. A suggested mechanism of rhodopsin photoactivation, triggered by the cis-trans isomerization of retinal, involves rotations of Glu134, Tyr223, Trp265, Lys296, and Tyr306 side chains and rearrangement of their H bonds. The model is in agreement with published electron cryomicroscopy, mutagenesis, chemical modification, cross-linking, Fourier transform infrared spectroscopy, Raman spectroscopy, electron paramagnetic resonance spectroscopy, NMR, and optical spectroscopy data. The rhodopsin model and the published structure of bacteriorhodopsin have very similar retinal-binding pockets.     

Ribonuclease A as a substrate of the protease from human immunodeficiency virus-1

Bovine pancreatic ribonuclease A (RNase) contains two bonds, Met29-Met30 and Tyr92-Pro93 which are representative of sites in the human immunodeficiency virus-1 (HIV-1) gag polyprotein precursors that are cleaved by the HIV-1 protease during viral maturation. Nevertheless, neither native nor performic acid-oxidized RNase is a substrate for the protease. However, RNase derivatives obtained by reduction and S-alkylation with iodoacetate or iodoacetamide undergo cleavage by the HIV-1 protease at a single site, Ala109-alkyl-Cys110, that is distinct from either of the two predicted bonds mentioned above. The neutral carboxyamido-methylcysteinyl derivative is cleaved 8 times faster than that containing the negatively charged carboxy-methyl substituent at P1'. Succinylation of these S-alkylated RNase derivatives creates a second site of cleavage by the protease between succinyl-Lys7 and Phe8. Thus, the pattern of cleavage of denatured RNase by the HIV-1 protease can be manipulated by chemical derivatization of the substrate, and the new sites of hydrolysis revealed by these studies add to our understanding of the specificity of this important enzyme.     

Adrenal Cortical Ribonuclease H (hybrid) (38494).

Ribonuclease Hybrid (RNaseH) from adrenal cortical tissue has been characterized. RNase H specifically degrades the RNA strand of purified RNA:DNA hybrids but is inactive on single- or double-stranded RNA or on DNA. The mode of clevage by RNase H is endonucleolytic, producing oligoribonucleoties with 3'-hydroxyl and 5'- phosphate termini. ACTH administration to guinea pigs results in no detectable change in adrenal cortical RNase H activity at times when changes in DNA synthesis are marked.     

Rescaled range analysis applied to the study delayed rectifier potassium channel kinetics

The gating of ion channels has widely been modeled by assuming that the transitions between open and closed states are a memoryless process. Nevertheless, analysis of records of unitary current events suggests that the kinetic process presents long lags (antipersistent correlation). Here, using the patch-voltage clamp technique and the rescaled range method, activity of single-channel delayed rectifier K(+) channels was studied. The experiment result showed that reversal potential was -73.3 mV in cell-attached mode. For the sequences of alternating open and shut time intervals, the Hurst coefficients were calculated for four different pipette potentials in rat dorsal root ganglion neurons. H=0.34169+/-0.00672 (n=4) for V=-30 mV; H=0.34632+/-0.0142 (n=3) for V=-40 mV; H=0.39237+/-0.0113 (n=4) for V=-50 mV; H=0.3954+/-0.0012 (n=4) for V=-60 mV. When the Hurst method was applied to the results from a simulated four-state Markovian model, it showed that it had different experimental data H coefficient, the distribution of the data values had no correlations between them, in particular, H=0.2531+/-0.00403 (n=50) for V=-40 mV. This indicates that open-dwell times and closed-dwell times are long lag (namely, antipersistent correlation) and do not change with the pipette potential applied to the patch.     

An inserted Gly residue fine tunes dynamics between mesophilic and thermophilic ribonucleases H

Dynamic processes are inherent properties of proteins and are crucial for a wide range of biological functions. To address how changes in protein sequence and structure affect dynamic processes, a quantitative comparison of microsecond-to-microsecond time scale conformational changes, measured by solution NMR spectroscopy, within homologous mesophilic and thermophilic ribonuclease H (RNase H) enzymes is presented. Kinetic transitions between the observed major state (high population) and alternate (low population) conformational state(s) of the substrate-binding handle region in RNase H from the mesophile Escherichia coli (ecRNH) and thermophile Thermus thermophilus (ttRNH) occur with similar kinetic exchange rate constants, but the difference in stability between exchanging conformers is smaller in ttRNH compared to ecRNH. The altered thermodynamic equilibrium between kinetically exchanging conformers in the thermophile is recapitulated in ecRNH by the insertion of a Gly residue within a putative hinge between alpha-helices B and C. This Gly insertion is conserved among thermophilic RNases H, and allows the formation of additional intrahelical hydrogen bonds. A Gly residue inserted between alpha-helices B and C appears to relieve unfavorable interactions in the transition state and alternate conformer(s) and represents an important adaptation to adjust conformational changes within RNase H for activity at high temperatures.