Neighboring aliphatic/aromatic side chain interactions between residues 9 and 10 in gramicidin channels

The interactions between an aliphatic or phenyl side chain and an indole ring in a phospholipid environment were investigated by synthesizing and characterizing gramicidins in which Trp(9) was ring-labeled and D-Leu(10) was replaced by D-Val, D-Ala, or D-Phe. All three analogues form conducting channels, with conductances that are lower than that of gramicidin A (gA) channels. The channel lifetimes vary by less than 50% from that of gA channels. Circular dichroism spectra and size-exclusion chromatography show that the conformation of each analogue in dimyristoylphosphatidylcholine (DMPC) vesicles is similar to the right-handed beta(6.3)-helical conformation that is observed for gA. (2)H NMR spectra of oriented samples in DMPC show large changes for the Trp(9) ring when residue 10 is modified, suggesting a steric interaction between D-Leu(10) and Trp(9), in agreement with previous acylation studies (R. E. Koeppe II et al. (1995) Biochemistry 34, 9299-9307). The outer quadrupolar splitting for Trp(9) is unchanged with D-Phe(10), at approximately 153 kHz, but increases by approximately 25 kHz with D-Val(10) and decreases by approximately 10 kHz with D-Ala(10). With D-Ala(10) or D-Val(10), the outer resonance splits into two in a temperature-dependent manner. The NMR spectra indicate that the side chain torsion angles chi1 and chi2 for Trp(9) change when residue 10 is substituted. The changes in chi1 are small, in all cases less than 10 degrees, as is Deltachi2 when D-Ala(10) is introduced, but with D-Val(10) and D-Phe(10) Deltachi2 is at least 25 degrees. We conclude that D-Leu(10) helps to stabilize an optimal orientation of Trp(9) in gA channels in lipid bilayers and that changes in Trp orientation alter channel conductance and lifetime without affecting the basic channel fold.     

The membrane interface dictates different anchor roles for "inner pair" and "outer pair" tryptophan indole rings in gramicidin A channels

We investigated the effects of substituting two of the four tryptophans (the "inner pair" Trp(9) and Trp(11) or the "outer pair" Trp(13) and Trp(15)) in gramicidin A (gA) channels. The conformational preferences of the doubly substituted gA analogues were assessed using circular dichroism spectroscopy and size-exclusion chromatography, which show that the inner tryptophans 9 and 11 are critical for the gA's conformational preference in lipid bilayer membranes. [Phe(13,15)]gA largely retains the single-stranded helical channel structure, whereas [Phe(9,11)]gA exists primarily as double-stranded conformers. Within this context, the (2)H NMR spectra from labeled tryptophans were used to examine the changes in average indole ring orientations, induced by the Phe substitutions and by the shift in conformational preference. Using a method for deuterium labeling of already synthesized gAs, we introduced deuterium selectively onto positions C2 and C5 of the remaining tryptophan indole rings in the substituted gA analogues for solid-state (2)H NMR spectroscopy. The (least possible) changes in orientation and overall motion of each indole ring were estimated from the experimental spectra. Regardless of the mixture of backbone folds, the indole ring orientations observed in the analogues are similar to those found previously for gA channels. Both Phe-substituted analogues form single-stranded channels, as judged from the formation of heterodimeric channels with the native gA. [Phe(13,15)]gA channels have Na(+) currents that are ~50% and lifetimes that are ~80% of those of native gA channels. The double-stranded conformer(s) of [Phe(9,11)]gA do not form detectable channels. The minor single-stranded population of [Phe(9,11)]gA forms channels with Na(+) currents that are ~25% and single-channel lifetimes that are ~300% of those of native gA channels. Our results suggest that Trp(9) and Trp(11), when "reaching" for the interface, tend to drive both monomer folding (to "open" a channel) and dimer dissociation (to "close" a channel). Furthermore, the dipoles of Trp(9) and Trp(11) are relatively more important for the single-channel conductance than are the dipoles of Trp(13) and Trp(15).     

Noncontact dipole effects on channel permeation. II. Trp conformations and dipole potentials in gramicidin A

The four Trp dipoles in the gramicidin A (gA) channel modulate channel conductance, and their side chain conformations should therefore be important, but the energies of different conformations are unknown. A conformational search for the right-handed helix based on molecular mechanics in vacuo yielded 46 conformations within 20 kcal/mol of the lowest energy conformation. The two lowest energy conformations correspond to the solid-state and solution-state NMR conformations, suggesting that interactions within the peptide determine the conformation. For representative conformations, the electrostatic potential of the Trp side chains on the channel axis was computed. A novel application of the image-series method of. Biophys. J. 9:1160-1170) was introduced to simulate the polarization of bulk water by the Trp side chains. For the experimentally observed structures, the CHARm toph19 potential energy (PE) of a cation in the channel center is -1.65 kcal/mol without images. With images, the PE is -1.9 kcal/mol, demonstrating that the images further enhance the direct dipole effect. Nonstandard conformations yielded less favorable PEs by 0.4-1.1 kcal/mol.     

Orientations of the tryptophan 9 and 11 side chains of the gramicidin channel based on deuterium nuclear magnetic resonance spectroscopy.

Deuterium nuclear magnetic resonance spectroscopy was used to investigate the orientations of the indole rings of Trp9 and Trp11 in specific indole-d5-labeled samples of gramicidin A incorporated into dimyristoyl phosphatidylcholine bilayers in the beta 6.3 channel conformation. The magnitudes and signs of the deuterium quadrupolar splittings were fit to the rings and assigned to specific ring bonds, using a full rotation search of the chi 1 and chi 2 angles of each Trp and a least-squares method. Unique assignments were obtained. The data and assignments are in close agreement with four sets of (chi 1, chi 2) angles for each Trp in which the indole N-H is oriented toward the membrane's exterior surface. (Four additional sets of (chi 1, chi 2) angles with the N-H's pointing toward the membrane interior are inconsistent with previous observations.) One of the sets of (chi 1, chi 2) angles for each Trp is consistent with the corresponding Trp orientation found by Arsen'ev et al. (1986. Biol. Membr. 3:1077-1104) for gramicidin in sodium dodecyl sulfate micelles. Together, the 1H and 2H nuclear magnetic resonance methods suggest that the Trp9 and Trp11 side chain orientations could be very similar in dimyristoyl phosphatidylcholine membranes and in sodium dodecyl sulfate micelles. The data for Trp11 could be fit using a static quadrupolar coupling constant of 180 kHz under the assumption that the ring is essentially immobile. By contrast, Trp9 could be fit only under the assumption that the quadrupolar splittings for ring 9 are reduced by approximately 14% due to motional averaging. Such a difference in motional averaging between rings 11 and 9 is also consistent with the 15N data of Hu et al. (1993. Biochemistry. 32:7035-7047).     

The preference of tryptophan for membrane interfaces: insights from N-methylation of tryptophans in gramicidin channels

To better understand the structural and functional roles of tryptophan at the membrane/water interface in membrane proteins, we examined the structural and functional consequences of Trp --> 1-methyl-tryptophan substitutions in membrane-spanning gramicidin A channels. Gramicidin A channels are miniproteins that are anchored to the interface by four Trps near the C terminus of each subunit in a membrane-spanning dimer. We masked the hydrogen bonding ability of individual or multiple Trps by 1-methylation of the indole ring and examined the structural and functional changes using circular dichroism spectroscopy, size exclusion chromatography, solid state (2)H NMR spectroscopy, and single channel analysis. N-Methylation causes distinct changes in the subunit conformational preference, channel-forming propensity, single channel conductance and lifetime, and average indole ring orientations within the membrane-spanning channels. The extent of the local ring dynamic wobble does not increase, and may decrease slightly, when the indole NH is replaced by the non-hydrogen-bonding and more bulky and hydrophobic N-CH(3) group. The changes in conformational preference, which are associated with a shift in the distribution of the aromatic residues across the bilayer, are similar to those observed previously with Trp --> Phe substitutions. We conclude that indole N-H hydrogen bonding is of major importance for the folding of gramicidin channels. The changes in ion permeability, however, are quite different for Trp --> Phe and Trp --> 1-methyl-tryptophan substitutions, indicating that the indole dipole moment and perhaps also ring size and are important for ion permeation through these channels.     

Mechanism of the efficient tryptophan fluorescence quenching in human gammaD-crystallin studied by time-resolved fluorescence

Human gammaD-crystallin (HgammaD-Crys) is a two-domain, beta-sheet eye lens protein found in the lens nucleus. Its long-term solubility and stability are important to maintain lens transparency throughout life. HgammaD-Crys has four highly conserved buried tryptophans (Trps), with two in each of the homologous beta-sheet domains. In situ, these Trps will be absorbing ambient UV radiation that reaches the lens. The dispersal of the excited-state energy to avoid covalent damage is likely to be physiologically relevant for the lens crystallins. Trp fluorescence is efficiently quenched in native HgammaD-Crys. Previous steady-state fluorescence measurements provide strong evidence for energy transfer from Trp42 to Trp68 in the N-terminal domain and from Trp130 to Trp156 in the C-terminal domain [Chen, J., et al. (2006) Biochemistry 45, 11552-11563]. Hybrid quantum mechanical-molecular mechanical (QM-MM) simulations indicated that the fluorescence of Trp68 and Trp156 is quenched by fast electron transfer to the amide backbone. Here we report additional information obtained using time-resolved fluorescence spectroscopy. In the single-Trp-containing proteins (Trp42-only, Trp68-only, Trp130-only, and Trp156-only), the highly quenched Trp68 and Trp156 have very short lifetimes, tau approximately 0.1 ns, whereas the moderately fluorescent Trp42 and Trp130 have longer lifetimes, tau approximately 3 ns. In the presence of the energy acceptor (Trp68 or Trp156), the lifetime of the energy donor (Trp42 or Trp130) decreased from approximately 3 to approximately 1 ns. The intradomain energy transfer efficiency is 56% in the N-terminal domain and is 71% in the C-terminal domain. The experimental values of energy transfer efficiency are in good agreement with those calculated theoretically. The absence of a time-dependent red shift in the time-resolved emission spectra of Trp130 proves that its local environment is very rigid. Time-resolved fluorescence anisotropy measurements with the single-Trp-containing proteins, Trp42-only and Trp130-only, indicate that the protein rotates as a rigid body and no segmental motion is detected. A combination of energy transfer with electron transfer results in short excited-state lifetimes of all Trps, which, together with the high rigidity of the protein matrix around Trps, could protect HgammaD-Crys from excited-state reactions causing permanent covalent damage.     

Human glycolipid transfer protein: probing conformation using fluorescence spectroscopy

Glycolipid transfer protein (GLTP) is a soluble 24 kDa protein that selectively accelerates the intermembrane transfer of glycolipids in vitro. Little is known about the GLTP structure and dynamics. Here, we report the cloning of human GLTP and characterize the environment of the three tryptophans (Trps) of the protein using fluorescence spectroscopy. Excitation at 295 nm yielded an emission maximum (lambda(max)) near 347 nm, indicating a relatively polar average environment for emitting Trps. Quenching with acrylamide at physiological ionic strength or with potassium iodide resulted in linear Stern-Volmer plots, suggesting accessibility of emitting Trps to soluble quenchers. Insights into reversible conformational changes accompanying changes in GLTP activity were provided by addition and rapid dilution of urea while monitoring changes in Trp or 1-anilinonaphthalene-8-sulfonic acid fluorescence. Incubation of GLTP with glycolipid liposomes caused a blue shift in the Trp emission maximum but diminished the fluorescence intensity. The blue-shifted emission maximum, centered near 335 nm, persisted after separation of glycolipid liposomes from GLTP, consistent with formation of a GLTP-glycolipid complex at a glycolipid-liganding site containing Trp. The results provide the first insights into human GLTP structural dynamics by fluorescence spectroscopy, including global conformational changes that accompany GLTP folding into an active conformational state as well as more subtle conformational changes that play a role in GLTP-mediated transfer of glycolipids between membranes, and establish a foundation for future studies of membrane rafts using GLTP.     

Effect of avidin on channel kinetics of biotinylated gramicidin

Membrane protein functioning basically depends on the supramolecular structure of the proteins which can be modulated by specific interactions with external ligands. The effect of a water-soluble protein bearing specific binding sites on the kinetics of ionic channels formed by gramicidin A (gA) in planar bilayer lipid membranes (BLM) has been studied using three independent approaches: (1) sensitized photoinactivation, (2) single-channel, and (3) autocorrelation measurements of current fluctuations. As shown previously [Rokitskaya, T. I., et al. (1996) Biochim. Biophys. Acta 1275, 221], the time course of the flash-induced current decrease in most cases follows a single-exponential decay with an exponential factor (tau) that corresponds to the gA single-channel lifetime. Addition of avidin does not affect tau for gA channels, but causes a dramatic increase in tau for channels formed by gA5XB, a biotinylated analogue of gA. This effect is reversed by addition of an excess of biotin to the bathing solution. The average single-channel duration of gA5XB was about 3.6 s as revealed by single-channel recording of the BLM current. After prolonged incubation with avidin, a long-lasting open state of the gA5XB channel appeared which did not close for more than 10 min. The data on gA5XB photoinactivation kinetics and single-channel measurements were confirmed by analysis of the corresponding power spectra of the current fluctuations obtained in the control, in the presence of avidin, and after the addition of biotin. We infer that avidin produces a deceleration of gA5XB channel kinetics by motional restriction of gA5XB monomers and dimers upon the formation of avidin and gA5XB complexes, which would stabilize the channel state and thus increase the single-channel lifetime.     

Luminescence studies of perturbation of tryptophan residues of tubulin in the complexes of tubulin with colchicine and colchicine analogues

Tubulin, a heterodimeric (alphabeta) protein, the main constituent of microtubules, binds efficiently with colchicine (consisting of a trimethoxybenzene ring, a seven-member ring and methoxy tropone moiety) and its analogues, viz., demecolcine and AC [2-methoxy-5-(2',3',4'-trimethoxyphenyl)tropone]. Tubulin contains eight tryptophan (Trp) residues at A21, A346, A388, A407, B21, B103, B346, and B407 in the two subunits. The role of these eight Trp residues in this interaction and also their perturbation due to binding have been explored via time-resolved fluorescence at room temperature and low-temperature (77 K) phosphorescence in a suitable cryosolvent. Both the time-resolved fluorescence data and 77 K phosphorescence spectra indicate that the emitting residues move toward a more hydrophobic and less polar environment after complex formation. The environment of emitting Trps in the complex also becomes slightly more heterogeneous. Our analysis using the experimental results, the calculation of the accessible surface area (ASA) of all the Trps in the wild type and tubulin-colchicine complex [Ravelli, R. B. G., et al. (2004) Nature 428, 198-202], the distance of the Trp residues from the different moieties of the colchicine molecule, the knowledge of the nature of the immediate residues (<5 A) present near each Trp residue, and the calculation of the intramolecular Trp-Trp energy transfer efficiencies indicate that Trp A346, Trp A407, Trp B21, and Trp B407 are the major contributors to the emission in the free protein, while Trp B21 and Trp B103 are mainly responsible for the emission of the complexes. A comparative account of the photophysical aspects of the drug molecules bound to protein in aqueous buffer and in buffer containing 40% ethylene glycol has been presented. The quantum yield and average lifetime of fluorescence in tubulin and its complexes with colchicine are used to predict the possible donors and the energy transfer (ET) efficiency in the ET process from Trps to colchicine in the complex. This study is a unique attempt to identify the Trp residues contributing to the emission in the free protein and in a complex of a multi-Trp protein with a drug molecule without performing the mutation of the protein.     

Involvement of tryptophans at the catalytic and subunit-binding domains of transcarboxylase

Transcarboxylase from Propionibacterium shermanii is a multisubunit enzyme. It consists of one central hexameric subunit to which six outer dimeric subunits are attached through twelve biotinyl subunits. Both the central and the outer subunits are multi-tryptophan (Trp) proteins, and each contains 5 Trps per monomer. The roles of the Trps during catalysis and assembly of the enzyme have been studied by using N-bromosuccinimide (NBS) oxidation as a probe. Modification of approximately 10 Trps of the total 90 Trps of the intact enzyme results in loss of activity. Both the substrates, viz., methylmalonyl-CoA and pyruvate, afford protection (approximately 50%) against inactivation caused by NBS. Analyses of tryptic peptide maps and intrinsic fluorescence studies have indicated that modification of 10 Trps of the whole enzyme does not cause extensive conformational changes. Therefore, the Trps appear to be essential for catalytic activity. NBS modification of the individual subunits at pH 6.5 has demonstrated differential reactivity of their Trps. Modification of the exposed/reactive Trps of either one of the subunits significantly affects the subunit assembly with the complementary unmodified subunits to form active enzyme. It is proposed that Trps are involved at the subunit-binding domains of either the central or the outer subunit of transcarboxylase, in addition to those critical for catalysis.     

Agonist-mediated conformational changes in acetylcholine-binding protein revealed by simulation and intrinsic tryptophan fluorescence

We delineated acetylcholine (ACh)-dependent conformational changes in a prototype of the nicotinic receptor ligand binding domain by molecular dynamics simulation and changes in intrinsic tryptophan (Trp) fluorescence. Prolonged molecular dynamics simulation of ACh-binding protein showed that binding of ACh establishes close register of Trps from adjacent subunits, Trp(143) and Trp(53), and draws the peripheral C-loop inward to occlude the entrance to the binding cavity. Close register of Trp(143) and Trp(53) was demonstrated by ACh-mediated quenching of intrinsic Trp fluorescence, elimination of quenching by mutation of one or both Trps to Phe, and decreased lifetime of Trp fluorescence by bound ACh. Occlusion of the binding cavity by the C-loop was demonstrated by restricted access of an extrinsic quencher of binding site Trp fluorescence by ACh. The collective findings showed that ACh initially establishes close register of conserved Trps from adjacent subunits and then draws the C-loop inward to occlude the entrance to the binding cavity.     

K channel subconductance levels result from heteromeric pore conformations

Voltage-gated K channels assemble from four identical subunits symmetrically arranged around a central permeation pathway. Each subunit harbors a voltage-sensing domain. The sigmoidal nature of the activation kinetics suggests that multiple sensors need to undergo a conformational change before the channel can open. Following activation, individual K channels alternate stochastically between two main permeation states, open and closed. This binary character of single channel behavior suggests the presence of a structure in the permeation pathway that can exist in only two conformations. However, single channel analysis of drk1 (K(v)2.1) K channels demonstrated the existence of four additional, intermediate conductance levels. These short-lived subconductance levels are visited when the channel gate moves between the closed and fully open state. We have proposed that these sublevels arise from transient heteromeric pore conformations, in which some, but not all, subunits are in the "open" state. A minimal model based on this hypothesis relates specific subconductance states with the number of activated subunits (Chapman et al., 1997). To stringently test this hypothesis, we constructed a tandem dimer that links two K channel subunits with different activation thresholds. Activation of this dimer by strong depolarizations resulted in the characteristic binary open-close behavior. However, depolarizations to membrane potentials in between the activation thresholds of the two parents elicited highly unusual single channel gating, displaying frequent visits to two subconductance levels. The voltage dependence and kinetics of the small and large sublevels associate them with the activation of one and two subunits, respectively. The data therefore support the hypothesis that subconductance levels result from heteromeric pore conformations. In this model, both sensor movement and channel opening have a subunit basis and these processes are allosterically coupled.     

Excluded volume in protein side-chain packing

The excluded volume occupied by protein side-chains and the requirement of high packing density in the protein interior should severely limit the number of side-chain conformations compatible with a given native backbone. To examine the relationship between side-chain geometry and side-chain packing, we use an all-atom Monte Carlo simulation to sample the large space of side-chain conformations. We study three models of excluded volume and use umbrella sampling to effectively explore the entire space. We find that while excluded volume constraints reduce the size of conformational space by many orders of magnitude, the number of allowed conformations is still large. An average repacked conformation has 20 % of its chi angles in a non-native state, a marked reduction from the expected 67 % in the absence of excluded volume. Interestingly, well-packed conformations with up to 50 % non-native chi angles exist. The repacked conformations have native packing density as measured by a standard Voronoi procedure. Entropy is distributed non-uniformly over positions, and we partially explain the observed distribution using rotamer probabilities derived from the Protein Data Bank database. In several cases, native rotamers that occur infrequently in the database are seen with high probability in our simulation, indicating that sequence-specific excluded volume interactions can stabilize rotamers that are rare for a given backbone. In spite of our finding that 65 % of the native rotamers and 85 % of chi(1) angles can be predicted correctly on the basis of excluded volume only, 95 % of positions can accommodate more than one rotamer in simulation. We estimate that, in order to quench the side-chain entropy observed in the presence of excluded volume interactions, other interactions (hydrophobic, polar, electrostatic) must provide an additional stabilization of at least 0.6 kT per residue in order to single out the native state.     

Mechanism of the highly efficient quenching of tryptophan fluorescence in human gammaD-crystallin

Quenching of the fluorescence of buried tryptophans (Trps) is an important reporter of protein conformation. Human gammaD-crystallin (HgammaD-Crys) is a very stable eye lens protein that must remain soluble and folded throughout the human lifetime. Aggregation of non-native or covalently damaged HgammaD-Crys is associated with the prevalent eye disease mature-onset cataract. HgammaD-Crys has two homologous beta-sheet domains, each containing a pair of highly conserved buried tryptophans. The overall fluorescence of the Trps is quenched in the native state despite the absence of the metal ligands or cofactors. We report the results of detailed quantitative measurements of the fluorescence emission spectra and the quantum yields of numerous site-directed mutants of HgammaD-Crys. From fluorescence of triple Trp to Phe mutants, the homologous pair Trp68 and Trp156 were found to be extremely quenched, with quantum yields close to 0.01. The homologous pair Trp42 and Trp130 were moderately fluorescent, with quantum yields of 0.13 and 0.17, respectively. In an attempt to identify quenching and/or electrostatically perturbing residues, a set of 17 candidate amino acids around Trp68 and Trp156 were substituted with neutral or hydrophobic residues. None of these mutants showed significant changes in the fluorescence intensity compared to their own background. Hybrid quantum mechanical-molecular mechanical (QM-MM) simulations with the four different excited Trps as electron donors strongly indicate that electron transfer rates to the amide backbone of Trp68 and Trp156 are extremely fast relative to those for Trp42 and Trp130. This is in agreement with the quantum yields measured experimentally and consistent with the absence of a quenching side chain. Efficient electron transfer to the backbone is possible for Trp68 and Trp156 because of the net favorable location of several charged residues and the orientation of nearby waters, which collectively stabilize electron transfer electrostatically. The fluorescence emission spectra of single and double Trp to Phe mutants provide strong evidence for energy transfer from Trp42 to Trp68 in the N-terminal domain and from Trp130 to Trp156 in the C-terminal domain. The backbone conformation of tryptophans in HgammaD-Crys may have evolved in part to enable the lens to become a very effective UV filter, while the efficient quenching provides an in situ mechanism to protect the tryptophans of the crystallins from photochemical degradation.     

Trp42 rotamers report reduced flexibility when the inhibitor acetyl-pepstatin is bound to HIV-1 protease

The Q7K/L331/L631 HIV-1 protease mutant was expressed in Escherichia coli and the effect of binding a substrate-analog inhibitor, acetyl-pepstatin, was investigated by fluorescence spectroscopy and molecular dynamics. The dimeric enzyme has four intrinsic tryptophans, located at positions 6 and 42 in each monomer. Fluorescence spectra and acrylamide quenching experiments show two differently accessible Trp populations in the apoenzyme with k(q1) = 6.85 x 10(9) M(-1) s(-1) and k(q2) = 1.88 x 10(9) M(-1) s(-1), that merge into one in the complex with k(q) = 1.78 x 10(9) M(-1) s(-1). 500 ps trajectory analysis of Trp X1/X2 rotameric interconversions suggest a model to account for the observed Trp fluorescence. In the simulations, Trp6/Trp6B rotameric interconversions do not occur on this timescale for both HIV forms. In the apoenzyme simulations, however, both Trp42s and Trp42Bs are flipping between X1/X2 states; in the complexed form, no such interconverions occur. A detailed investigation of the local Trp environments sampled during the molecular dynamics simulation suggests that one of the apoenzyme Trp42B rotameric interconversions would allow indole-quencher contact, such as with nearby Tyr59. This could account for the short lifetime component. The model thus interprets the experimental data on the basis of the conformational fluctuations of Trp42s alone. It suggests that the rotameric interconversions of these Trps, located relatively far from the active site and at the very start of the flap region, becomes restrained when the apoenzyme binds the inhibitor. The model is thus consistent with associating components of the fluorescence decay in HIV-1 protease to ground state conformational heterogeneity.     

Using robotics to fold proteins and dock ligands

The problems of protein folding and ligand docking have been explored largely using molecular dynamics or Monte Carlo methods. These methods are very compute intensive because they often explore a much wider range of energies, conformations and time than necessary. In addition, Monte Carlo methods often get trapped in local minima. We initially showed that robotic motion planning permitted one to determine the energy of binding and dissociation of ligands from protein binding sites (Singh et al., 1999). The robotic motion planning method maps complicated three-dimensional conformational states into a much simpler, but higher dimensional space in which conformational rearrangements can be represented as linear paths. The dimensionality of the conformation space is of the same order as the number of degrees of conformational freedom in three-dimensional space. We were able to determine the relative energy of association and dissociation of a ligand to a protein by calculating the energetics of interaction for a few thousand conformational states in the vicinity of the protein and choosing the best path from the roadmap. More recently, we have applied roadmap planning to the problem of protein folding (Apaydin et al., 2002a). We represented multiple conformations of a protein as nodes in a compact graph with the edges representing the probability of moving between neighboring states. Instead of using Monte Carlo simulation to simulate thousands of possible paths through various conformational states, we were able to use Markov methods to calculate the steady state occupancy of each conformation, needing to calculate the energy of each conformation only once. We referred to this Markov method of representing multiple conformations and transitions as stochastic roadmap simulation or SRS. We demonstrated that the distribution of conformational states calculated with exhaustive Monte Carlo simulations asymptotically approached the Markov steady state if the same Boltzman energy distribution was used in both methods. SRS permits one to calculate contributions from all possible paths simultaneously with far fewer energy calculations than Monte Carlo or molecular dynamics methods. The SRS method also permits one to represent multiple unfolded starting states and multiple, near-native, folded states and all possible paths between them simultaneously. The SRS method is also independent of the function used to calculate the energy of the various conformational states. In a paper to be presented at this conference (Apaydin et al., 2002b) we have also applied SRS to ligand docking in which we calculate the dynamics of ligand-protein association and dissociation in the region of various binding sites on a number of proteins. SRS permits us to determine the relative times of association to and dissociation from various catalytic and non-catalytic binding sites on protein surfaces. Instead of just following the best path in a roadmap, we can calculate the contribution of all the possible binding or dissociation paths and their relative probabilities and energies simultaneously.     

Modification of myosin subfragment 1 tryptophans by dimethyl(2-hydroxy-5-nitrobenzyl)sulfonium bromide

Modification of tryptophanyl residues (Trps) of myosin subfragments 1 (S-1) was performed with dimethyl(2-hydroxy-5-nitrobenzyl)sulfonium bromide (DHNBS). Under controlled conditions, pH 6 at 0 degrees C and 10-min reaction with 10-100-fold molar excess, K+(EDTA) activity was reduced down to less than half, whereas Ca2+-ATPase activity increased and acto-S-1-ATPase was not affected. The number of modified Trps (up to 2.5) agreed well with the number of 2-hydroxy-5-nitrobenzyl moieties incorporated in S-1. The thiol groups of S-1 were not affected up to 50-fold molar excess of DHNBS, thus indicating that the modification was selective for Trps. The modification of as few as one Trp caused a blue shift of the emission spectrum, accompanied by a reduction in the fluorescence quantum yield. The accessibility of Trps to the fluorescence quencher acrylamide is drastically reduced upon modification, indicating that DHNBS-reactive Trps are more "exposed" than the DHNBS-refractive ones. DHNBS modification did not seem to affect the ATP-induced tryptophan fluorescence enhancement of S-1. The effect of DHNBS modification of the intrinsic fluorescence of S-1 indicates that the modified Trps are located in a polar environment and that they may be identical with the long-lifetime Trps of Torgerson [Torgerson, P. (1984) Biochemistry 23, 3002-3007]. The most reactive Trp is located in the N-terminal 27-kDa fragment of the S-1 heavy chain. It might also be inferred from the above data that the nonexposed and ATP-perturbed Trp(s) is (are) located in the 50-kDa fragment.     

The role of Trp side chains in tuning single proton conduction through gramicidin channels

We present an extensive set of measurements of proton conduction through gramicidin A (gA), B (gB), and M (gM) homodimer channels which have 4, 3, or 0 Trp residues at each end of the channel, respectively. In gA we find a shoulder separating two domains of conductance increasing with concentration, confirming the results of Eisenman, G., B. Enos, J. Hagglund, and J. Sandblom. 1980. Ann. NY. Acad. Sci. 339:8-20. In gB, the shoulder is shifted by approximately 1/2 pH unit to higher H(+) concentrations and is very sharply defined. No shoulder appears in the gM data, but an associated transition from sublinear to superlinear I-V values occurs at a 100-fold higher [H(+)] in gM than in gA. The data in the low concentration domain are analyzed using a configuration space model of single-proton conduction, assuming that the difference in the proton potential of mean force (PMF) between gA and its analogs is constant, similar to the results of Anderson, D., R. B. Shirts, T. A. Cross, and D. D. Busath. 2001. Biophys. J. 81:1255-1264. Our results suggest that the average amplitudes of the calculated proton PMFs are nearly correct, but that the water reorientation barrier calculated for gA by molecular dynamics using the PM6 water model (Pomès, R., and B. Roux. 1997. Biophys. J. 72:246a) must be reduced in amplitude by 1.5 kcal/mol or more, and is not rate-limiting for gA.     

Comparison of conformational features of staphylococcal nuclease in ternary complexes with pdTp, pdGp, and nitrophenyl-pdTp.

The conformations of wild-type staphylococcal nuclease (SNase) in the ternary complexes with thymidine 3',5'-bisphosphate (pdTp), 2'-deoxyguanine 3',5'-bisphosphate (pdGp), and thymidine 3'-phosphate 5'-(p-nitrophenylphosphate) (NpdTp) with Ca2+ were examined by two-dimensional NMR NOESY and ROESY experiments. The results of these experiments indicate that the conformational features of the SNase are quite similar in the three ternary complexes. This suggests that the conformational features of SNase, in these ternary complexes, are not strongly dependent on whether the 5'-phosphate is a mono- or diester. This is in contrast to our prior studies on substitutions of active site charged amino acids which indicated that the conformational features of SNase in the ternary complex are quite sensitive to substitutions for active site charged amino acids (Hibler et al., 1987; Wilde et al., 1988; Pourmotabbed et al., 1990). The similarity of the SNase conformational features in the ternary complexes with pdTp and pdGp indicates that the features of the nucleotide bound at the active site are not strong determinants of the enzyme conformation in the ternary complexes. These conclusions are in general agreement with the results on pdApdT ternary complexes with SNase which suggested that it is the conformational features of the bound nucleic acid which determine the differences in catalysis observed for SNase with different substrates (Weber et al., 1991), more so than the conformational features of the enzyme.     

Protein's native state stability in a chemically induced denaturation mechanism

In this work, we present a generalization of Zwanzig's protein unfolding analysis [Zwanzig, R., 1997. Two-state models of protein folding kinetics. Proc. Natl Acad. Sci. USA 94, 148-150; Zwanzig, R., 1995. Simple model of protein folding kinetics. Proc. Natl Acad. Sci. USA 92, 9801], in order to calculate the free energy change Delta(N)(D)F between the protein's native state N and its unfolded state D in a chemically induced denaturation. This Extended Zwanzig Model (EZM) is both based on an equilibrium statistical mechanics approach and the inclusion of experimental denaturation curves. It enables us to construct a suitable partition function Z and to derive an analytical formula for Delta(N)(D)F in terms of the number K of residues of the macromolecule, the average number nu of accessible states for each single amino acid and the concentration C(1/2) where the midpoint of the N<==>D transition occurs. The results of the EZM for proteins where chemical denaturation follows a sigmoidal-type profile, as it occurs for the case of the T70N human variant of lysozyme (PDB code: T70N) [Esposito, G., et al., 2003. J. Biol. Chem. 278, 25910-25918], can be splitted into two lines. First, EZM shows that for sigmoidal denaturation profiles, the internal degrees of freedom of the chain play an outstanding role in the stability of the native state. On the other hand, that under certain conditions DeltaF can be written as a quadratic polynomial on concentration C(1/2), i.e., DeltaF approximately aC(1/2)(2)+bC(1/2)+c, where a,b,c are constant coefficients directly linked to protein's size K and the averaged number of non-native conformations nu. Such functional form for DeltaF has been widely known to fit experimental measures in chemically induced protein denaturation [Yagi, M., et al., 2003. J. Biol. Chem. 278, 47009-47015; Asgeirsson, B., Guojonsdottir, K., 2006. Biochim. Biophys. Acta 1764, 190-198; Sharma, S., et al., 2006. Protein Pept. Lett. 13(4), 323-329; Salem, M., et al., 2006. Biochim. Biophys. Acta 1764(5), 903-912] so EZM can shed some light into the physical meaning of the experimental values for the a,b,c coefficients.