Open "back door" in a molecular dynamics simulation of acetylcholinesterase.

The enzyme acetylcholinesterase generates a strong electrostatic field that can attract the cationic substrate acetylcholine to the active site. However, the long and narrow active site gorge seems inconsistent with the enzyme's high catalytic rate. A molecular dynamics simulation of acetylcholinesterase in water reveals the transient opening of a short channel, large enough to pass a water molecule, through a thin wall of the active site near tryptophan-84. This simulation suggests that substrate, products, or solvent could move through this "back door," in addition to the entrance revealed by the crystallographic structure. Electrostatic calculations show a strong field at the back door, oriented to attract the substrate and the reaction product choline and to repel the other reaction product, acetate. Analysis of the open back door conformation suggests a mutation that could seal the back door and thus test the hypothesis that thermal motion of this enzyme may open multiple routes of access to its active site.

Covalent modifications of the amyloid beta peptide by hydroxynonenal: Effects on metal ion binding by monomers and insights into the fibril topology.

Amyloid ß peptides (Aß) and metal ions are associated with oxidative stress in Alzheimer's disease (AD). Oxidative stress, acting on ω-6 polyunsaturated fatty acyl chains, produces diverse products, including 4-hydroxy-2-nonenal (HNE), which can covalently modify the Aß that helped to produce it. To examine possible feedback mechanisms involving Aß, metal ions and HNE production, the effects of HNE modification and fibril formation on metal ion binding was investigated. Results indicate that copper(II) generally inhibits the modification of His side chains in Aß by HNE, but that once modified, copper(II) still binds to Aß with high affinity. Fibril formation protects only one of the three His residues in Aß from HNE modification, and this protection is consistent with proposed models of fibril structure. These results provide insight into a network of biochemical reactions that may be operating as a consequence of oxidative stress in AD, or as part of the pathogenic process.

Orientational order determination by internal reflection infrared spectroscopy.

PATIR-FTIR spectroscopy is a powerful technique for the determination of molecular order in thin films such as supported lipid membranes, but it relies on electromagnetic theory which is incomplete and potentially misleading. A complete derivation of the current theory for two, three and four phase system has been reported. The two phase and thin film approximations most commonly used in practice have been shown to represent the thickness-dependent expressions from which they are derived with a high degree of accuracy. However, these expressions are based on the macroscopic behavior of dielectric materials, and may not be accurate when applied to microscopic circumstances. The potential error introduced is qualitatively and quantitatively significant. Further experimental and theoretical work is needed to verify the accuracy of this theory, or to refute and refine it. This effort to do this is warranted by the power and increasing popularity of the technique.

Folding of beta-sheet membrane proteins: a hydrophobic hexapeptide model.

Beta-sheets, in the form of the beta-barrel folding motif, are found in several constitutive membrane proteins (porins) and in several microbial toxins that assemble on membranes to form oligomeric transmembrane channels. We report here a first step towards understanding the principles of beta-sheet formation in membranes. In particular, we describe the properties of a simple hydrophobic hexapeptide, acetyl-Trp-Leu5 (AcWL5), that assembles cooperatively into beta-sheet aggregates upon partitioning into lipid bilayer membranes from the aqueous phase where the peptide is strictly monomeric and random coil. The aggregates, containing 10 to 20 monomers, undergo a relatively sharp and reversible thermal unfolding at approximately 60 degreesC. No pores are formed by the aggregates, but they do induce graded leakage of vesicle contents at very high peptide to lipid ratios. Because beta-sheet structure is not observed when the peptide is dissolved in n-octanol, trifluoroethanol or sodium dodecyl sulfate micelles, aggregation into beta-sheets appears to be an exclusive property of the peptide in the bilayer membrane interface. This is an expected consequence of the hypothesis that a reduction in the free energy of partitioning of peptide bonds caused by hydrogen bonding drives secondary structure formation in membrane interfaces. But, other features of interfacial partitioning, such as side-chain interactions and reduction of dimensionality, must also contribute. We estimate from our partitioning data that the free energy reduction per residue for aggregation is about 0.5 kcal mol-1. Although modest, its aggregate effect on the free energy of assembling beta-sheet proteins can be huge. This surprising finding, that a simple hydrophobic hexapeptide readily assembles into oligomeric beta-sheets in membranes, reveals the potent ability of membranes to promote secondary structure in peptides, and shows that the formation of beta-sheets in membranes is more facile than expected. Furthermore, it provides a basis for understanding the observation that membranes promote self-association of beta-amyloid peptides. AcWL5 and related peptides thus provide a good starting point for designing peptide models for exploring the principles of beta-sheet formation in membranes.

A ligand-mediated dimerization mode for vancomycin.

BACKGROUND: Vancomycin and related glycopeptide antibiotics exert their antimicrobial effect by binding to carboxy-terminal peptide targets in the bacterial cell wall and preventing the biosynthesis of peptidoglycan. Bacteria can resist the action of these agents by replacing the peptide targets with depsipeptides. Rational efforts to design new agents effective against resistant bacteria require a thorough understanding of the structural determinants of peptide recognition by vancomycin. RESULTS: The crystal structure of vancomycin in complex with N-acetyl-D-alanine has been determined at atomic resolution. Two different oligomeric interactions are seen in the structure: back-to-back dimers, as previously described for the vancomycin-acetate complex, and novel face-to-face dimers, mediated largely by the bound ligands. Models of longer, naturally occurring peptide ligands may be built by extension of N-acetyl-D-alanine. These larger ligands can form an extensive array of polar and nonpolar interactions with two vancomycin monomers in the face-to-face configuration. CONCLUSIONS: A new dimeric form of vancomycin has been found in which two monomers are related in a face-to-face configuration, and bound ligands comprise a large portion of the dimer interface. The relative importance of face-to-face and back-to-back dimers to the antimicrobial activity of vancomycin remains to be established, but face-to-face interactions appear to explain how increased antimicrobial activity may arise in covalent vancomycin dimers.

Lipoprotein A-I structure.

High density lipoproteins are produced by the liver as protein-lipid complexes with a characteristic discoidal shape. A crystal structure is available for the chief protein component of these complexes, apolipoprotein A-I, but controversy about how this protein is situated with respect to the lipid components has flourished for lack of experimental techniques that can characterize protein structure in a lipid environment. New spectroscopic techniques developed to address this problem now indicate that apolipoprotein A-I is arranged as a helical belt around a bilayer of phospholipids. This is an important step towards understanding how these lipoproteins regulate cholesterol transport.

Structure and dynamics of the active site gorge of acetylcholinesterase: synergistic use of molecular dynamics simulation and X-ray crystallography.

The active site of acetylcholinesterase (AChE) from Torpedo californica is located 20 A from the enzyme surface at the bottom of a narrow gorge. To understand the role of this gorge in the function of AChE, we have studied simulations of its molecular dynamics. When simulations were conducted with pure water filling the gorge, residues in the vicinity of the active site deviated quickly and markedly from the crystal structure. Further study of the original crystallographic data suggests that a bis-quaternary decamethonium (DECA) ion, acquired during enzyme purification, residues in the gorge. There is additional electron density within the gorge that may represent small bound cations. When DECA and 2 cations are placed within the gorge, the simulation and the crystal structure are dramatically reconciled. The small cations, more so than DECA, appear to stabilize part of the gorge wall through electrostatic interactions. This part of the gorge wall is relatively thin and may regulate substrate, product, and water movement through the active site.

Vancomycin binding to low-affinity ligands: delineating a minimum set of interactions necessary for high-affinity binding.

Bacterial resistance to vancomycin has been attributed to the loss of an intermolecular hydrogen bond between vancomycin and its peptidoglycan target when cell wall biosynthesis proceeds via depsipeptide intermediates rather than the usual polypeptide intermediates. To investigate the relative importance of this hydrogen bond to vancomycin binding, we have determined crystal structures at 1.0 A resolution for the vancomycin complexes with three ligands that mimic peptides and depsipeptides found in vancomycin-sensitive and vancomycin-resistant bacteria: N-acetylglycine, D-lactic acid, and 2-acetoxy-D-propanoic acid. These, in conjunction with structures that have been reported previously, indicate higher-affinity ligands elicit a structural change in the drug not seen with these low-affinity ligands. They also enable us to define a minimal set of drug-ligand interactions necessary to confer higher-affinity binding on a ligand. Most importantly, these structures point to factors in addition to the loss of an intermolecular hydrogen bond that must be invoked to explain the weaker affinity of vancomycin for depsipeptide ligands. These factors are important considerations for the design of vancomycin analogues to overcome vancomycin resistance.

The role of sugar residues in molecular recognition by vancomycin.

The sugar residues of the glycopeptide antibiotic vancomycin contribute to the cooperativity of ligand binding, thereby increasing ligand affinity and enhancing antimicrobial activity. To assess the structural basis for these effects, we determined a 0.98 A X-ray crystal structure of the vancomycin aglycon and compared it to structures of several intact vancomycin:ligand complexes. The crystal structure reveals that the aglycon binds acetate anions and forms back-to-back dimeric complexes in a manner similar to that of intact vancomycin. However, the four independent copies of the aglycon in each asymmetric unit of the crystal exhibit a high degree of conformational heterogeneity. These results suggest that the sugar residues, in addition to enlarging and strengthening the dimer interface, provide steric constraints that limit the vancomycin molecule to a relatively small number of productive conformations.

Toward a rational design of peptide inhibitors of ribonucleotide reductase: structure-function and modeling studies.

Mammalian ribonucleotide reductase, a chemotherapeutic target, has two subunits, mR1 and mR2, and is inhibited by AcF(1)TLDADF(7), denoted P7. P7 corresponds to the C-terminus of mR2 and competes with mR2 for binding to mR1. We report results of a structure-function analysis of P7, obtained using a new assay measuring peptide ligand binding to mR1, that demonstrate stringent specificity for Phe at F(7), high specificity for Phe at F(1), and little specificity for the N-acyl group. They support a structural model in which the dominant interactions of P7 occur at two mR1 sites, the F(1) and F(7) subsites. The model is constructed from the structure of Escherichia coli R1 (eR1) complexed with the C-terminal peptide from eR2, aligned sequences of mR1 and eR1, and the trNOE-derived structure of mR1-bound P7. Comparison of this model with similar models constructed for mR1 complexed with other inhibitory ligands indicates that increased F(1) subsite interaction can offset lower F(7) subsite interaction and suggests strategies for the design of new, higher affinity inhibitors.

Electronic transitions in molecules in static external fields. I. Indole and Trp-59 in ribonuclease T1.

Electronic transition properties of indole perturbed by its environment were calculated by use of quantum-mechanical semi-empirical numerical methods. The environment was represented by a discrete set of charges placed at different positions around the indole ring. Wavelength shifts and transition intensity changes in indole were evaluated for several, specifically modeled geometries of external charges. This methodology was employed to estimate the extent of spectroscopic changes induced by small nonprotein polar species on the Trp-59 residue in the anisotropic environment of the protein ribonuclease T1. The geometry of the residue environment was obtained from dynamically equilibrated X-ray crystallographic data of the protein.

Infrared spectroscopy of supported lipid monolayer, bilayer, and multibilayer membranes.

Cecropin A was examined in supported monolayer, bilayer, and multibilayer lipid membranes using attenuated total internal reflection Fourier-transform infrared spectroscopy. The spectral features provide an abundance of information about the conformation and orientation of the peptide, as well as about the effects of the peptide on lipid order. In this case, they serve to contrast results from the three preparations. The results of monolayer and bilayer studies are generally similar, although differences in the nature of the membranes appear to cause minor changes in the conformation and orientation of the peptide. The results of the multibilayer studies are different in many respects from those of the monolayer and bilayer studies, suggesting that fundamentally different peptide-lipid interactions occur in multibilayers.

A rational strategy for enhancing the affinity of vancomycin towards depsipeptide ligands.

Glycopeptide antibiotics with enhanced affinity for model depsipeptide ligands may also exhibit enhanced efficacy against bacteria exhibiting the vanA resistance phenotype. To design modified agents with enhanced affinity for these ligands, and better understand why traditional agents have low affinity for depsipeptide ligands, free energy perturbation studies were performed on vancomycin derivatives by means of molecular dynamics simulation. The results suggest that modifications of the asparagine side chain on residue 3 of the antibiotic which enhance its hydrophobicity will enhance the affinity of glycopeptide antibiotics for depsipeptide ligands, and act synergistically with other modifications that enhance the efficacy of these agents against vanA-positive bacteria.

Fourier transform infrared linked analysis of conformational changes in annexin V upon membrane binding.

Annexins are ubiquitous cellular proteins of unknown primary function that bind to anionic phospholipid membranes in a calcium-dependent manner. Correlative studies involving X-ray crystallography and electron microscopy suggest that annexins undergo a structural change upon binding to supported lipid monolayer membranes. In this investigation, novel spectroscopic and analytical techniques have been applied to verify and characterize this change. Soluble annexin V was examined with ordinary transmission infrared spectroscopy, while membrane-bound annexin V was examined with both transmission and internal reflection infrared spectroscopy. Spectra were processed by linked analysis, whereby multiple spectra are fit simultaneously with component bands that are constrained to share common fitting parameters. This approach is shown to enhance the sensitivity and accuracy of the bandfitting procedure. Our results are consistent with the general mode of membrane binding inferred from electron microscopy studies, and they provide independent support for the conclusion that annexin V undergoes a conformational change upon binding to lipid monolayer membranes. Most likely, this change involves the formation of new beta structure in which interstrand hydrogen bonds orient parallel to the membrane surface.

Determination of molecular order in supported lipid membranes by internal reflection Fourier transform infrared spectroscopy.

When polarized internal reflection infrared spectroscopy is used to determine molecular order in supported lipid membranes, the results are critically dependent on the accuracy of assumptions made about the evanescent electric field amplitudes in the membrane. In this work, we examine several expressions used for calculating evanescent electric field amplitudes in supported lipid monolayers and bilayers, and test their validity by measuring the infrared dichroism of poly-gamma-benzyl-L-glutamate and poly-beta-benzyl-L-aspartate under conditions in which their molecular order is known. Our results indicate that treating such systems as a simple single interface between two semi-infinite bulk phases is more accurate than the commonly employed thin-film approximation. This implies that earlier conclusions about molecular order in supported lipid membranes may require substantial revision.

Buried water in homologous serine proteases.

Buried water molecules in the structurally homologous family of eukaryotic serine proteases were examined to determine whether buried waters and their protein environments are conserved in these proteins. We found 16 equivalent water sites conserved in trypsin/ogen, chymotrypsin/ogen, elastase, kallikrein, thrombin, rat tonin and rat mast cell protease, and 5 additional water sites in enzymes which share the primary specificity of trypsin. Based on an alignment of 30 serine protease sequences, it appears that the protein environments of these 21 conserved buried waters are highly conserved. The protein environments of buried waters are comprised primarily of atoms from highly conserved residues or main chain atoms from nonconserved residues. In one instance, the protein environment of a water is conserved even in the presence of an unlikely Pro/Ala substitution. We also note 3 instances in which a histidine side chain substitutes for water, suggesting that the structural role of water at these sites is satisfied by the presence of an alternative hydrogen bonding partner. Buried waters appear to be integral structural components of these proteins and should be incorporated into protein structures predicted on the basis of sequence homology to this family, including the catalytic domains of coagulation proteases.

A constraint on possible stoichiometries of myocardial sodium-calcium exchange.

We have examined the sodium-calcium exchange stoichiometry in Langendorff-perfused rabbit hearts using gamma-emitting tracers under conditions of sodium pump inhibition. Following a 60-min perfusion with 10(-5) acetylstrophanthidin, and extracellular concentrations [Na]o = 70 mM and [Ca]o = 300 microM, intracellular sodium rose to 59.2 mM. At this point an increase in extracellular calcium [Ca]o = 1.52 mM) caused a net efflux of sodium, but an increase in sodium [Na]o = 105 mM) caused no measurable change. When sodium and calcium were simultaneously increased according to the ratio [Na]o)n/[Ca]o = [Na]'o)n/[Ca]'o, a sodium efflux is observed when n = 4, but not when n = 3. These results are consistent with an exchange stoichiometry of 3 Na+ for each Ca2+ ion, but not values of 4 or more.

Molecular dynamics of tryptophan in ribonuclease-T1. II. Correlations with fluorescence.

The interactions of tryptophan-59 (TRP-59) and its protein environment in ribonuclease-T1 (RNAse-T1) were examined in a 50-ps molecular dynamics simulation. The simulation used was previously shown to demonstrate a fluorescence anisotropy decay that closely agreed with the experimentally determined limiting anisotropy for RNAse-T1 (Axelsen, P. H., C. Haydock, and F. G. Prendergast. 1988. Biophys. J. 54:249-258). Further characterization of TRP-59 side chain dynamics and its protein environment has now been completed and correlated to other photophysical properties of this protein. Angular fluctuations of the side chain occur at rates of 1-10 cycles/ps and are limited to +/- 0.3 radians in all directions. Side chain motions are primarily limited by nonpolar collisions, although most side chain atoms have some collisional contact with polar atoms in the adjacent protein matrix or water. The steric relationship between PRO-39 and TRP-59 changes abruptly at 16 ps into the simulation. Two types of interaction with water are observed. First, a structural water appears to H-bond with the greater than N-H group of TRP-59. Second, water frequently contacts the six-atom ring. The electrostatic field experienced by the TRP-59 rings appears to be relatively constant and featureless regardless of ring orientation. We make the following interferences from our data: The fluorescent emission of TRP-59 may be red-shifted relative to TRP in nonpolar solvents either as a result of specific interactions with the structural water or relaxations of proximal bulk water and polar protein moieties. The quenching efficiency of polar interactions with TRP-59 must be extremely low given their frequency and the high quantum yield of RNAse-T1. This low efficiency may be due to restricted and unfavorable interaction geometries. PRO-39 is located near two titratable HIS residues in RNAse-T1 and may be involved in pH-dependent fluorescence phenomena by virtue of a metastable interaction with TRP-59. The interaction of bulk water with TRP-59 illustrates features of the gated transition state model for transient exposure to exogenously added collisional quenching agents. The restrictive environment of TRP-59 suggests that extrinsic quenching can only occur via interactions with the edge of the indole six-atom ring and that the efficiency of a quencher in a protein environment is likely to be a function of molecular symmetry.

Membrane-induced folding of cecropin A.

Lipid membranes manifest a diverse array of surface forces that can fold and orient an approaching protein. To better understand these forces and their ability to influence protein function, we have used infrared spectroscopy with isotopic editing to characterize the 37-residue membrane-active antimicrobial polypeptide cecropin A as it approached, adsorbed onto, and finally penetrated various lipid membranes. Intermediate stages in this process were isolated for study by the use of internal reflection and Langmuir trough techniques. Results indicate that this peptide adopts well-ordered secondary structure while superficially adsorbed to a membrane surface. Its conformation is predominantly alpha-helical, although some beta structure is likely to be present. The longitudinal axis of the helical structure, and the transverse axes of any beta structure, are preferentially oriented parallel to the membrane surface. The peptide expands the membrane against pressure when it penetrates the membrane surface, but its structure and orientation do not change. These observations indicate that interactions between the peptide and deeper hydrophobic regions of the membrane provide energy to perform thermodynamic work, but separate and distinct interactions between the peptide and superficial components of the membrane are responsible for peptide folding. These results have broad implications for our understanding of the mechanism of action and the specificity of these antimicrobial peptides.