The structure of bovine F1-ATPase complexed with the antibiotic inhibitor aurovertin B.

In the structure of bovine mitochondrial F1-ATPase that was previously determined with crystals grown in the presence of adenylyl-imidodiphosphate (AMP-PNP) and ADP, the three catalytic beta-subunits have different conformations and nucleotide occupancies. Adenylyl-imidodiphosphate is bound to one beta-subunit (betaTP), ADP is bound to the second (betaDP), and no nucleotide is bound to the third (betaE). Here we show that the uncompetitive inhibitor aurovertin B binds to bovine F1 at two equivalent sites in betaTP and betaE, in a cleft between the nucleotide binding and C-terminal domains. In betaDP, the aurovertin B pocket is incomplete and is inaccessible to the inhibitor. The aurovertin B bound to betaTP interacts with alpha-Glu399 in the adjacent alphaTP subunit, whereas the aurovertin B bound to betaE is too distant from alphaE to make an equivalent interaction. Both sites encompass betaArg-412, which was shown by mutational studies to be involved in binding aurovertin. Except for minor changes around the aurovertin pockets, the structure of bovine F1-ATPase is the same as determined previously. Aurovertin B appears to act by preventing closure of the catalytic interfaces, which is essential for a catalytic mechanism involving cyclic interconversion of catalytic sites.

Protonatable residues at the cytoplasmic end of transmembrane helix-2 in the signal transducer HtrI control photochemistry and function of sensory rhodopsin I.

Neutral residue replacements were made of 21 acidic and basic residues within the N-terminal half of the Halobacterium salinarium signal transducer HtrI [the halobacterial transducer for sensory rhodopsin I (SRI)] by site-specific mutagenesis. The replacements are all within the region of HtrI that we previously concluded from deletion analysis to contain sites of interaction with the phototaxis receptor SRI. Immunoblotting shows plasmid expression of the htrI-sopI operon containing the mutations produces SRI and mutant HtrI in cells at near wild-type levels. Six of the HtrI mutations perturb photochemical kinetics of SRI and one reverses the phototaxis response. Substitution with neutral amino acids of Asp-86, Glu-87, and Glu-108 accelerate, and of Arg-70, Arg-84, and Arg-99 retard, the SRI photocycle. Opposite effects on photocycle rate cancel in double mutants containing one replaced acidic and one replaced basic residue. Laser flash spectroscopy shows the kinetic perturbations are due to alteration of the rate of reprotonation of the retinylidene Schiff base. All of these mutations permit normal attractant and repellent signaling. On the other hand, the substitution of Glu-56 with the isosteric glutamine converts the normally attractant effect of orange light to a repellent signal in vivo at neutral pH (inverted signaling). Low pH corrects the inversion due to Glu-56 -> Gln and the apparent pK of the inversion is increased when arginine is substituted at position 56. The results indicate that the cytoplasmic end of transmembrane helix-2 and the initial part of the cytoplasmic domain contain interaction sites with SRI. To explain these and previous results, we propose a model in which (i) the HtrI region identified here forms part of an electrostatic bonding network that extends through the SRI protein and includes its photoactive site; (ii) alteration of this network by photoisomerization-induced Schiff base deprotonation and reprotonation shifts HtrI between attractant and repellent conformations; and (iii) HtrI mutations and extracellular pH alter the equilibrium ratios of these conformations.

Alpha-helical, but not beta-sheet, propensity of proline is determined by peptide environment.

Proline is established as a potent breaker of both alpha-helical and beta-sheet structures in soluble (globular) proteins. Thus, the frequent occurrence of the Pro residue in the putative transmembrane helices of integral membrane proteins, particularly transport proteins, presents a structural dilemma. We propose that this phenomenon results from the fact that the structural propensity of a given amino acid may be altered to conform to changes imposed by molecular environment. To test this hypothesis on proline, we synthesized model peptides of generic sequence H2N-(Ser-LyS)2-Ala- Leu-Z-Ala-Leu-Z-Trp-Ala-Leu-Z-(Lys-Ser)3-OH (Z = Ala and/or Pro). Peptide conformations were analyzed by circular dichroism spectroscopy in aqueous buffer, SDS, lysophosphatidylglycerol micelles, and organic solvents (methanol, trifluoroethanol, and 2-propanol). The helical propensity of Pro was found to be greatly enhanced in the membrane-mimetic environments of both lipid micelles and organic solvents. Proline was found to stabilize the alpha-helical conformation relative to Ala at elevated temperatures in 2-propanol, an observation that argues against the doctrine that Pro is the most potent alpha-helix breaker as established in aqueous media. Parallel studies in deoxycholate micelles of the temperature-induced conformational transitions of the single-spanning membrane bacteriophage IKe major coat protein, in which the Pro-containing wild type was compared with Pro30 --> Ala mutant, Pro was found to protect the helix, but disrupt the beta-sheet structure as effectively as it does to model peptides in water. The intrinsic capacity of Pro to disrupt beta-sheets was further reflected in a survey of porins where Pro was found to be selectively excluded from the core of membrane-spanning beta-sheet barrels. The overall data provide a rationale for predicting and understanding the structural consequences when Pro occurs in the context of a membrane.

The hard-soft acid-base principle in enzymatic catalysis: dual reactivity of phosphoenolpyruvate.

In this paper, the chemical reactivity of C3 of phosphoenolpyruvate (PEP) has been analyzed in terms of density functional theory quantified through quantum chemistry calculations. PEP is involved in a number of important enzymatic reactions, in which its C3 atom behaves like a base. In three different enzymatic reactions analyzed here, C3 sometimes behaves like a soft base and sometimes behaves like a hard base in terms of the hard-soft acid-base principle. This dual nature of C3 of PEP was found to be related to the conformational change of the molecule. This leads to a testable hypothesis: that PEP adopts particular conformations in the enzyme-substrate complexes of different PEP-using enzymes, and that the enzymes control the reactivity through controlling the dihedral angle between the carboxylate and the C==C double bond of PEP.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

The multiple phenotypes of allosteric receptor mutants.

Channel-linked neurotransmitter receptors are membrane-bound heterooligomers made up of distinct, although homologous, subunits. They mediate chemo-electrical signal transduction and its regulation via interconversion between multiple conformations that exhibit distinct pharmacological properties and biological activities. The large diversity of functional properties and the widely pleiotropic phenotypes, which arise from point mutations in their subunits (or from subunit substitutions), are interpreted in terms of an allosteric model that incorporates multiple discrete conformational states. The model predicts that three main categories of phenotypes may result from point mutations, altering selectively one (or more) of the following features: (i) the properties of individual binding sites (K phenotype), (ii) the biological activity of the ion channel (gamma phenotype) of individual conformations, or (iii) the isomerization constants between receptor conformations (L phenotype). Several nicotinic acetylcholine and glycine receptor mutants with complex phenotypes are quantitatively analyzed in terms of the model, and the analogies among phenotypes are discussed.

Structure of the N intermediate of bacteriorhodopsin revealed by x-ray diffraction.

X-ray diffraction experiments revealed the structure of the N photointermediate of bacteriorhodopsin. Since the retinal Schiff base is reprotonated from Asp-96 during the M to N transition in the photocycle, and Asp-96 is reprotonated during the lifetime of the N intermediate, or immediately after, N is a key intermediate for understanding the light-driven proton pump. The N intermediate accumulates in large amounts during continuous illumination of the F171C mutant at pH 7 and 5 degrees Celsius. Small but significant changes of the structure were detected in the x-ray diffraction profile under these conditions. The changes were reversible and reproducible. The difference Fourier map indicates that the major change occurs near helix F. The observed diffraction changes between N and the original state were essentially identical to the diffraction changes reported for the M intermediate of the D96N mutant of bacteriorhodopsin. Thus, we find that the protein conformations of the M and N intermediates of the photocycle are essentially the same, in spite of the fact that in M the Schiff base is unprotonated and in N it is protonated. The observed structural change near helix F will increase access of the Schiff base and Asp-96 to the cytoplasmic surface and facilitate the proton transfer events that begin with the decay of the M state.

Crystallographic analysis of endogenous peptides associated with HLA-DR1 suggests a common, polyproline II-like conformation for bound peptides.

The structure of the human major histocompatibility complex (MHC) class II molecule HLA-DR1 derived from the human lymphoblastoid cell line LG-2 has been determined in a complex with the Staphylococcus aureus enterotoxin B superantigen. The HLA-DR1 molecule contains a mixture of endogenous peptides derived from cellular or serum proteins bound in the antigen-binding site, which copurify with the class II molecule. Continuous electron density for 13 amino acid residues is observed in the MHC peptide-binding site, suggesting that this is the core length of peptide that forms common interactions with the MHC molecule. Electron density is also observed for side chains of the endogenous peptides. The electron density corresponding to peptide side chains that interact with the DR1-binding site is more clearly defined than the electron density that extends out of the binding site. The regions of the endogenous peptides that interact with DRI are therefore either more restricted in conformation or sequence than the peptide side chains or amino acids that project out of the peptide-binding site. The hydrogen-bond interactions and conformation of a peptide model built into the electron density are similar to other HLA-DR-peptide structures. The bound peptides assume a regular conformation that is similar to a polyproline type II helix. The side-chain pockets and conserved asparagine residues of the DR1 molecule are well-positioned to interact with peptides in the polyproline type II conformation and may restrict the range of acceptable peptide conformations.

Thermodynamic parameters for the helix-coil transition of oligopeptides: molecular dynamics simulation with the peptide growth method.

The helix-coil transition equilibrium of polypeptides in aqueous solution was studied by molecular dynamics simulation. The peptide growth simulation method was introduced to generate dynamic models of polypeptide chains in a statistical (random) coil or an alpha-helical conformation. The key element of this method is to build up a polypeptide chain during the course of a molecular transformation simulation, successively adding whole amino acid residues to the chain in a predefined conformation state (e.g., alpha-helical or statistical coil). Thus, oligopeptides of the same length and composition, but having different conformations, can be incrementally grown from a common precursor, and their relative conformational free energies can be calculated as the difference between the free energies for growing the individual peptides. This affords a straightforward calculation of the Zimm-Bragg sigma and s parameters for helix initiation and helix growth. The calculated sigma and s parameters for the polyalanine alpha-helix are in reasonable agreement with the experimental measurements. The peptide growth simulation method is an effective way to study quantitatively the thermodynamics of local protein folding.

Growth hormone secretagogues: characterization, efficacy, and minimal bioactive conformation.

Another class of growth hormone (GH) secretagogues has been discovered by altering the backbone structure of a flexible linear GH-releasing peptide (GHRP). In vitro and in vivo characterization confirms these GH secretagogues as the most potent and smallest (M(r) < 500) reported. Anabolic efficacy is demonstrated in rodents with intermittent delivery. A convergent model of the bioactive conformation of GHRPs is developed and is supported by the NMR structure of a highly potent cyclic analog of GHRP-2. The model and functional data provide a logical framework for the further design of low-molecular weight secretagogues and illustrate the utility of an interdisciplinary approach to elucidating potential bound-state conformations of flexible peptide ligands.

In vitro integration of human immunodeficiency virus type 1 cDNA into targets containing protein-induced bends.

Integration of human immunodeficiency virus type 1 cDNA into a target DNA can be strongly influenced by the conformation of the target. For example, integration in vitro is sometimes favored in target DNAs containing sequence-directed bends or DNA distortions caused by bound proteins. We have analyzed the effect of DNA bending by studying integration into two well-characterized protein-DNA complexes: Escherichia coli integration host factor (IHF) protein bound to a phage IHF site, and the DNA binding domain of human lymphoid enhancer factor (LEF) bound to a LEF site. Both of these proteins have previously been reported to bend DNA by approximately 140 degrees. Binding of IHF greatly increases the efficiency of in vitro integration at hotspots within the IHF site. We analyzed a series of mutants in which the IHF site was modified at the most prominent hotspot. We found that each variant still displayed enhanced integration upon IHF binding. Evidently the local sequence is not critical for formation of an IHF hotspot. LEF binding did not create preferred sites for integration. The different effects of IHF and LEF binding can be rationalized in terms of the different proposed conformations of the two protein-DNA complexes.

Lattice model for rapidly folding protein-like heteropolymers.

Protein folding is a relatively fast process considering the astronomical number of conformations in which a protein could find itself. Within the framework of a lattice model, we show that one can design rapidly folding sequences by assigning the strongest attractive couplings to the contacts present in a target native state. Our protein design can be extended to situations with both attractive and repulsive contacts. Frustration is minimized by ensuring that all the native contacts are again strongly attractive. Strikingly, this ensures the inevitability of folding and accelerates the folding process by an order of magnitude. The evolutionary implications of our findings are discussed.

On constructing folding heteropolymers.

Simplified models of the protein-folding process have led to valuable insights into the generic properties of the folding of heteropolymers. On the basis of theoretical arguments, Shakhnovich and Gutin [(1993) Proc. Natl. Acad. Sci. USA 90, 7195-7199] have proposed a specific method to generate folding sequences for one of these. Here we present a model of folding in heteropolymers that is comparable in simplicity but different in spirit to the one studied by Shakhnovich and Gutin. In our model, the proposed recipe for constructing folding sequence fails. We find that, as a rule, the construction of folding sequences is impossible to achieve by looking at the native conformation only. Rather, competing conformations have to be taken into account too. An evolutionary algorithm that generates folding sequences by optimizing both stability of the native state and folding time is described. Remarkably, this algorithm produces, among others, sequences that fold reproducibly to metastable states.

Trapping of branched DNA in microfabricated structures.

We have observed electrostatic trapping of tribranched DNA molecules undergoing electrophoresis in a microfabricated pseudo-two-dimensional array of posts. Trapping occurs in a unique transport regimen in which the electrophoretic mobility is extremely sensitive to polymer topology. The arrest of branched polymers is explained by considering their center-of-mass motion; in certain conformations, owing to the constraints imposed by the obstacles a molecule cannot advance without the center of mass first moving a short distance backwards. The depth of the resulting local potential well can be much greater than the thermal energy so that escape of an immobilized molecule can be extremely slow. We summarize the expected behavior of the mobility as a function of field strength and topology and point out that the microfabricated arrays are highly suitable for detecting an extremely small number of branched molecules in a very large population of linear molecules.

Crystal and molecular structure of paclitaxel (taxol).

Paclitaxel (formerly called taxol), an important anticancer drug, inhibits cell replication by binding to and stabilizing microtubule polymers. As drug-receptor interactions are governed by the three-dimensional stereochemistries of both participants, we have determined the crystal structure of paclitaxel to identify its conformational preferences that may be related to biological activity. The monoclinic crystals contain two independent paclitaxel molecules in the asymmetric unit plus several water and dioxane solvent molecules. Taxane ring conformation is very similar in both paclitaxel molecules and is similar to the taxane ring conformation found in the crystal structure of the paclitaxel analogue docetaxel (formerly called taxotere). The two paclitaxel molecules have carbon-13 side-chain conformations that differ from each other and from that of the corresponding side chain in the docetaxel crystal structure. The carbon-13 side-chain conformation of one paclitaxel molecule is similar to what was proposed from NMR studies done in polar solvents, while that of the other paclitaxel molecule is different and hitherto unobserved. The paclitaxel molecules interact with each other and with solvent atoms through an extensive network of hydrogen bonds. Analysis of the hydrogen-bonding network together with structure-activity studies may suggest which atoms of paclitaxel are important for binding to microtubule receptors.