Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

Anisotropic refinement of the structure of Thermoascus aurantiacusxylanase I

The isotropic crystallographic model of the structure of xylanase I from Thermoascus aurantiacus (TAXI) has now been refined anisotropically at 1.14 Å resolution to a standard residual of R = 11.1% for all data. TAXI is amongst the five largest proteins deposited in the Protein Data Bank to have been refined with anisotropic displacement parameters (ADPs) at this level of resolution. The anisotropy analysis revealed a more isotropic distribution of anisotropy than usually observed previously. Adding ADPs resulted in high-quality electron-density maps which revealed discrepancies from the previously suggested primary sequences for this enzyme. Side-chain conformational disorder was modelled for 16 residues, including Trp275, a bulky residue at the active site. An unrestrained refinement was consistent with the protonation of the catalytic acid/base glutamate and the deprotonation of the nucleophile glutamate, as required for catalysis. The thermal stability of TAXI is reinterpreted in the light of the new refined model.

Estudo dos aspectos eletrostáticos da interação entre polieletrólitos e macroíons

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Conformational dynamics of thyroid hormones by variable temperature nuclear magnetic resonance: The role of side chain rotations and cisoid/transoid interconversions

H-1 NMR spectra of the thyroid hormone thyroxine recorded at low temperature and high field show splitting into two peaks of the resonance due to the H2,6 protons of the inner (tyrosyl) ring. A single resonance is observed in 600 MHz spectra at temperatures above 185 K. An analysis of the line shape as a function of temperature shows that the coalescence phenomenon is due to an exchange process with a barrier of 37 kJ mol(-1). This is identical to the barrier for coalescence of the H2',6' protons of the outer (phenolic) ring reported previously for the thyroid hormones and their analogues. It is proposed that the separate peaks at low temperature are due to resonances for H2,6 in cisoid and transoid conformers which are populated in approximately equal populations. These two peaks are averaged resonances for the individual H2 and H6 protons. Conversion of cisoid to transoid forms can occur via rotation of either the alanyl side chain or the outer ring, from one face of the inner ring to the other. It is proposed that the latter process is the one responsible for the observed coalescence phenomenon. The barrier to rotation of the alanyl side chain is greater than or equal to 37 kJ mol(-1), which is significantly larger than has previously been reported for Csp(2)-Csp(3) bonds in other Ph-CH2-X systems. The recent crystal structure of a hormone agonist bound to the ligand-binding domain of the rat thyroid hormone receptor (Wagner et al. Nature 1995, 378, 690-697) shows the transoid form to be the bound conformation. The significant energy barrier to cisoid/transoid interconversion determined in the current study combined with the tight fit of the hormone to its receptor suggests that interconversion between the forms cannot occur at the receptor site but that selection for the preferred bound form occurs from the 50% population of the transoid form in solution.

Kinetics of baculovirus replication and release using real-time quantitative polymerase chain reaction

The study of viral-based processes is hampered by (a) their complex, transient nature, (b) the instability of products, and (c) the lack of accurate diagnostic assays. Here, we describe the use of real-time quantitative polymerase chain reaction to characterize baculoviral infection. Baculovirus DNA content doubles every 1.7 h from 6 h post-infection until replication is halted at the onset of budding. No dynamic equilibrium exists between replication and release, and the kinetics are independent of the cell density at the time of infection. No more than 16% of the intracellular virus copies bud from the cell. (C) 2002 John Wiley & Sons, Inc. Biotechnol Bioeng 77: 476-480, 2002; DOI 10.1002/bit.10126.

Steady-state brain glucose transport kinetics re-evaluated with a four-state conformational model.

Glucose supply from blood to brain occurs through facilitative transporter proteins. A near linear relation between brain and plasma glucose has been experimentally determined and described by a reversible model of enzyme kinetics. A conformational four-state exchange model accounting for trans-acceleration and asymmetry of the carrier was included in a recently developed multi-compartmental model of glucose transport. Based on this model, we demonstrate that brain glucose (G(brain)) as function of plasma glucose (G(plasma)) can be described by a single analytical equation namely comprising three kinetic compartments: blood, endothelial cells and brain. Transport was described by four parameters: apparent half saturation constant K(t), apparent maximum rate constant T(max), glucose consumption rate CMR(glc), and the iso-inhibition constant K(ii) that suggests G(brain) as inhibitor of the isomerisation of the unloaded carrier. Previous published data, where G(brain) was quantified as a function of plasma glucose by either biochemical methods or NMR spectroscopy, were used to determine the aforementioned kinetic parameters. Glucose transport was characterized by K(t) ranging from 1.5 to 3.5 mM, T(max)/CMR(glc) from 4.6 to 5.6, and K(ii) from 51 to 149 mM. It was noteworthy that K(t) was on the order of a few mM, as previously determined from the reversible model. The conformational four-state exchange model of glucose transport into the brain includes both efflux and transport inhibition by G(brain), predicting that G(brain) eventually approaches a maximum concentration. However, since K(ii) largely exceeds G(plasma), iso-inhibition is unlikely to be of substantial importance for plasma glucose below 25 mM. As a consequence, the reversible model can account for most experimental observations under euglycaemia and moderate cases of hypo- and hyperglycaemia.

Hydration, conformational states and kinetics of yeast hexokinase PII /

The kinetic study of the coupled enzymatic reaction involving monomeric yeast hexokinase PII (HK) and yeast glucose-6-phosphate dehydrogenase (G-6-PDH) yields a Michaelis constant of 0.15 ± 0.01 mM for D-glucose. At pH 8.7 HK is present in monomeric form. The addition of polyethylene glycol (PEG), to the reaction mixture increased the affinity of HK for glucose, independent ofMW of the PEG from 2000 to 10000. The osmotic stress exerted by PEG can be used to measure the change in number of water molecules that accompany enzyme conformational changes (Rand, et al., 1993). Results indicate that the G-6-PDH is not osmotically sensitive and thus, the change in the number of PEG-inaccessible water molecules (ANw) measured in the coupled reaction is only the difference between the glucose-bound and glucosefree conformations of HK. ANw ~ 450 with PEGs of MW > 2000 under conditions for both binding (Reid and Rand, 1997) and kinetic assays. The contribution water may play in the binding of ATP (Km = 0.24 + 0.02 mM) has also been examined. It was found that in this case ANw = (for osmotic pressures < 2.8x10* dynes/cm^), suggesting no additional numbers of waters are displaced when ATP binds to HK. Osmotic pressure experiments were also performed with dimeric HK. It was determined that both the monomeric and dimeric forms of HK give the same ANw under low pressures. If this large ANw is due to conformational flexibility, it would appear that the flexibility is not reduced upon dimerization ofthe enzyme.

Effects of large anphiphilic ligands upon the spectra and kinetics of cytochrome C oxidase

Cytoch ro me c oxidase (ferrocytochrome c : 02 oxidoreductase ; EC 1.9. 3.1) is the terminal enzyme in the mitochondrial electron transport chain, catalyzing the transfer of electrons from ferrocytochrome c to molecular oxygen. The effects of two large amphiphilic molecules .. valinomycin and dibucaine upon the spectra of the isolated enzyme and upon the activity of both isolated enzyme and enzyme in membrane systems are investigated by using spectrophotometric and oxygen electrode techniques. The results show that both valinomycin and dibucaine change the Soret region of the spectrum and cause a partial inhibition in a concentration range higher than that in which they act as ionophores. It is concluded that both valinomycin and dibucain~ binding induce a conformational change of the protein structure which modifies the spectrum of the a3 CUB centre and diminishes the rate of electron transfer between cytochrome a and the binuclear centre.

Effects of large anphiphilic ligands upon the spectra and kinetics of cytochrome C oxidase

Cytoch ro me c oxidase (ferrocytochrome c : 02 oxidoreductase ; EC 1.9. 3.1) is the terminal enzyme in the mitochondrial electron transport chain, catalyzing the transfer of electrons from ferrocytochrome c to molecular oxygen. The effects of two large amphiphilic molecules - valinomycin and dibucaine upon the spectra of the isolated enzyme and upon the activity of both isolated enzyme and enzyme in membrane systems are investigated by using spectrophotometric and oxygen electrode techniques. The results show that both valinomycin and dibucaine change the Soret region of the speetrum and cause a partial inhibition in a concentration range higher than that in which they act as ionophores. It is concluded that both valinomycin and dibucaine binding induce a conformational change of the protein structure which modifies the spectrum of the a3 CUB centre and diminishes the rate of electron transfer between cytochrome a and the binuclear centre.

Thiol and Sulfenic Acid Oxidation of AhpE, the One-Cysteine Peroxiredoxin from Mycobacterium tuberculosis: Kinetics, Acidity Constants, and Conformational Dynamics

Drug resistance and virulence of Mycobacterium tuberculosis are partially related to the pathogen`s antioxidant systems. Peroxide detoxification in this bacterium is achieved by the heme-containing catalase peroxidase and different two-cysteine peroxiredoxins. M. tuberculosis genome also codifies for a putative one-cysteine peroxiredoxin, alkyl hydroperoxide reductase E (MtAhpE). Its expression was previously demonstrated at a transcriptional level, and the crystallographic structure of the recombinant protein was resolved under reduced and oxidized states. Herein, we report that the conformation of MtAhpE changed depending on its single cysteine redox state, as reflected by different tryptophan fluorescence properties and changes in quaternary structure. Dynamics of fluorescence changes, complemented by competition kinetic assays, were used to perform protein functional studies. MtAhE reduced peroxynitrite 2 orders of magnitude faster than hydrogen peroxide (1.9 x 10(7) M(-1) s(-1) vs 8.2 x 10(4) M(-1) s(-1) at pH 7.4 and 25 degrees C, respectively). The latter also caused cysteine overoxidation to sulfinic acid, but at much slower rate constant (40 M(-1) s(-1)). The pK(a) of the thiol in the reduced enzyme was 5.2, more than one unit lower than that of the sulfenic acid in the oxidized enzyme. The pH profile of hydrogen peroxide-mediated thiol and sulfenic acid oxidations indicated thiolate and sulfenate as the reacting species. The formation of sulfenic acid as well as the catalytic peroxidase activity of MtAhpE was demonstrated using the artificial reducing substrate thionitrobenzoate. Taken together, our results indicate that MtAhpE is a relevant component in the antioxidant repertoire of M. tuberculosis probably involved in peroxide and specially peroxynitrite detoxification.

A Single Mutation at the Sheet Switch Region Results in Conformational Changes Favoring lambda 6 Light-Chain Fibrillogenesis

Systemic amyloid light-chain (LC) amyloidosis is a disease process characterized by the pathological deposition of monoclonal LCs in tissue. All LC subtypes are capable of fibril formation although lambda chains, particularly those belonging to the lambda 6 type, are overrepresented. Here, we report the thermodynamic and in vitro fibrillogenic properties of several mutants of the lambda 6 protein 6aJL2 in which Pro7 and/or His8 was substituted by Ser or Pro. The H8P and H8S mutants were almost as stable as the wildtype protein and were poorly fibrillogenic. In contrast, the P7S mutation decreased the thermodynamic stability of 6aJL2 and greatly enhanced its capacity to form amyloid-like fibrils in vitro. The crystal structure of the P7S mutant showed that the substitution induced both local and long-distance effects, such as the rearrangement of the V(L) (variable region of the light chain)-V(L) interface. This mutant crystallized in two orthorhombic polymorphs, P2(1)2(1)2(1) and C222(1). In the latter, a monomer that was not arranged in the typical Bence-Jones dimer was observed for the first time. Crystal-packing analysis of the C222(1) lattice showed the establishment of intermolecular beta-beta interactions that involved the N-terminus and beta-strand B and that these could be relevant in the mechanism of LC fibril formation. Our results strongly suggest that Pro7 is a key residue in the conformation of the N-terminal sheet switch motif and, through long-distance interactions, is also critically involved in the contacts that stabilized the V(L) interface in lambda 6 LCs. (C) 2009 Elsevier Ltd. All rights reserved.

First report on the kinetics of the uncatalyzed esterification of cellulose under homogeneous reaction conditions: a rationale for the effect of carboxylic acid anhydride chain-length on the degree of biopolymer substitution

The kinetics of the homogeneous acylation of microcrystalline cellulose, MCC, with carboxylic acid anhydrides with different acyl chain-length (Nc; ethanoic to hexanoic) in LiCl/N,N-dimethylacetamide have been studied by conductivity measurements from 65 to 85 A degrees C. We have employed cyclohexylmethanol, CHM, and trans-1,2-cyclohexanediol, CHD, as model compounds for the hydroxyl groups of the anhydroglucose unit of cellulose. The ratios of rate constants of acylation of primary (CHM; Prim-OH) and secondary (CHD; Sec-OH) groups have been employed, after correction, in order to split the overall rate constants of the reaction of MCC into contributions from the discrete OH groups. For the model compounds, we have found that k((Prim-OH))/k((Sec-OH)) > 1, akin to reactions of cellulose under heterogeneous conditions; this ratio increases as a function of increasing Nc. The overall, and partial rate constants of the acylation of MCC decrease from ethanoic- to butanoic-anhydride and then increase for pentanoic- and hexanoic anhydride, due to subtle changes in- and compensations of the enthalpy and entropy of activation.

Conformational changes between the active-site and regulatory light chain of myosin as determined by luminescence resonance energy transfer: The effect of nucleotides and actin

Myosin is thought to generate movement of actin filaments via a conformational change between its light-chain domain and its catalytic domain that is driven by the binding of nucleotides and actin. To monitor this change, we have measured distances between a gizzard regulatory light chain (Cys 108) and the active site (near or at Trp 130) of skeletal myosin subfragment 1 (S1) by using luminescence resonance energy transfer and a photoaffinity ATP-lanthanide analog. The technique allows relatively long distances to be measured, and the label enables site-specific attachment at the active-site with only modest affect on myosin’s enzymology. The distance between these sites is 66.8 ± 2.3 Å when the nucleotide is ADP and is unchanged on binding to actin. The distance decreases slightly with ADP-BeF3, (−1.6 ± 0.3 Å) and more significantly with ADP-AlF4 (−4.6 ± 0.2 Å). During steady-state hydrolysis of ATP, the distance is temperature-dependent, becoming shorter as temperature increases and the complex with ADP⋅Pi is favored over that with ATP. We conclude that the distance between the active site and the light chain varies as Acto-S1-ADP ≈ S1-ADP > S1-ADP-BeF3 > S1-ADP-AlF4 ≈ S1-ADP-Pi and that S1-ATP > S1-ADP-Pi. The changes in distance are consistent with a substantial rotation of the light-chain binding domain of skeletal S1 between the prepowerstroke state, simulated by S1-ADP-AlF4, and the post-powerstroke state, simulated by acto-S1-ADP.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

Pseudoperiodic "Living" and/or controlled cationic ring-opening copolymerization of oxetane with tetrahydropyran:microstructure of polymers vs kinetics of chain growth

The "living" and/or controlled cationic ring-opening bulk copolymerization of oxetane (Ox) with tetrahydropyran (THP) (cyclic ether with no homopolymerizability) at 35°C was examined using ethoxymethyl-1 -oxoniacyclohexane hexafluoroantimonate (EMOA) and (BF3 · CH3OH)THP as fast and slow initiator, respectively, yielding living and nonliving polymers with pseudoperiodic sequences (i.e., each pentamethylene oxide fragment inserted into the polymer is flanked by two trimethylene oxide fragments). Good control over number-average molecular weight (Mn up to 150000 g mol-1) with molecular weight distribution (MWD ∼ 1.4-1, 5) broader than predicted by the Poison distribution (MWDs > 1 +1/DPn) was attained using EMOA as initiating system, i.e., C 2H5OCH2Cl with 1.1 equiv of AgSbF6 as a stable catalyst and 1.1 equiv of 2,6-di-tert-butylpyridine used as a non-nucleophilic proton trap. With (BF3 · CH 3OH)THP, a drift of the linear dependence M n(GPC) vs Mn(theory) to lower molecular weight was observed together with the production of cyclic oligomers, ∼3-5% of the Ox consumed in THP against ∼30% in dichloromethane. Structural and kinetics studies highlighted a mechanism of chains growth where the rate of mutual conversion between "strain ACE species" (chain terminated by a tertiary 1-oxoniacyclobutane ion, Al) and "strain-free ACE species" (chain terminated by a tertiary 1-oxoniacyclohexane ion, Tl) depends on the rate at which Ox converts the stable species T1 (kind of "dormant" species) into a living "propagating" center A1 (i.e., k aapp[Ox]). The role of the THP solvent associated with the suspension of irreversible and reversible transfer reactions to polymer, when the polymerization is initiated with EMOA, was predicted by our kinetic considerations. The activation -deactivation pseudoequilibrium coefficient (Qt) was then calculated in a pure theoretical basis. From the measured apparent rate constant of Ox (kOxapp) and THP (kTHPapp = ka(endo)app) consumption, Qt and reactivity ratio (kp/kd, k a(endo)/ka(exo), and ks/ka(endo) were calculated, which then allow the determination of the transition rate constant of elementary step reactions that governs the increase of Mu with conversion. © 2009 American Chemical Society.