Peptides and proteins in membranes: what can we learn via computer simulations?

Membrane and membrane-active peptides and proteins play a crucial role in numerous cell processes, such as signaling, ion conductance, fusion, and others. Many of them act as highly specific and efficient drugs or drug targets, and, therefore, attract growing interest of medicinal chemists. Because of experimental difficulties with characterization of their spatial structure and mode of membrane binding, essential attention is given now to molecular modeling techniques. During the last years an important progress has been achieved in molecular dynamics (MD) and Monte Carlo (MC) simulations of peptides and proteins with explicit and/or implicit theoretical models of membranes. The first ones allow atomic-resolution studies of peptides behavior on the membrane-water interfaces. Models with implicit consideration of membrane are of a special interest because of their computational efficiency and ability to account for principal trends in protein-lipid interactions. In this approximation, the bilayer is usually treated as continuum whose properties vary along the membrane thickness, and membrane insertion is simulated using either MC or MD methods. This review surveys recent applications of both types of lipid bilayer models in computer simulations of a wide variety of peptides and proteins with different biological activities. Theoretical background of the membrane models is considered with examples of their applications to biologically relevant problems. The emphasis of the review is made on recent MC and MD computations, on structural and/or functional information, which may be obtained via molecular modeling. The approximations and shortcomings of the models, along with their perspectives in design of new membrane active drugs, are discussed.

Factors important for fusogenic activity of peptides: molecular modeling study of analogs of fusion peptide of influenza virus hemagglutinin.

Nine analogs of fusion peptide of influenza virus hemagglutinin whose membrane perturbation activity has been thoroughly tested [Murata et al. (1992) Biochemistry 31, 1986-1992; Murata et al. (1993) Biophys. J. 64, 724-734] were characterized by molecular modeling techniques with the aim of delineating any specific structural and/or hydrophobic properties inherent in peptides with fusogenic activity. It was shown that, regardless of characteristics common to all analogs (peripheral disposition at the water-lipid interface, amphiphilic nature, alpha-helical structure, etc.), only fusion active peptides reveal a specific 'tilted oblique-oriented' pattern of hydrophobicity on their surfaces and a certain depth of penetration to the non-polar membrane core. The conclusion was reached that these factors are among the most important for the specific destabilization of a bilayer, which is followed by membrane fusion.

A solvent model for simulations of peptides in bilayers. I. Membrane-promoting alpha-helix formation.

We describe an efficient solvation model for proteins. In this model atomic solvation parameters imitating the hydrocarbon core of a membrane, water, and weak polar solvent (octanol) were developed. An optimal number of solvation parameters was chosen based on analysis of atomic hydrophobicities and fitting experimental free energies of gas-cyclohexane, gas-water, and octanol-water transfer for amino acids. The solvation energy term incorporated into the ECEPP/2 potential energy function was tested in Monte Carlo simulations of a number of small peptides with known energies of bilayer-water and octanol-water transfer. The calculated properties were shown to agree reasonably well with the experimental data. Furthermore, the solvation model was used to assess membrane-promoting alpha-helix formation. To accomplish this, all-atom models of 20-residue homopolypeptides-poly-Leu, poly-Val, poly-Ile, and poly-Gly in initial random coil conformation-were subjected to nonrestrained Monte Carlo conformational search in vacuo and with the solvation terms mimicking the water and hydrophobic parts of the bilayer. All the peptides demonstrated their largest helix-forming tendencies in a nonpolar environment, where the lowest-energy conformers of poly-Leu, Val, Ile revealed 100, 95, and 80% of alpha-helical content, respectively. Energetic and conformational properties of Gly in all environments were shown to be different from those observed for residues with hydrophobic side chains. Applications of the solvation model to simulations of peptides and proteins in the presence of membrane, along with limitations of the approach, are discussed.

A solvent model for simulations of peptides in bilayers. II. Membrane-spanning alpha-helices.

We describe application of the implicit solvation model (see the first paper of this series), to Monte Carlo simulations of several peptides in bilayer- and water-mimetic environments, and in vacuum. The membrane-bound peptides chosen were transmembrane segments A and B of bacteriorhodopsin, the hydrophobic segment of surfactant lipoprotein, and magainin2. Their conformations in membrane-like media are known from the experiments. Also, molecular dynamics study of surfactant lipoprotein with different explicit solvents has been reported (Kovacs, H., A. E. Mark, J. Johansson, and W. F. van Gunsteren. 1995. J. Mol. Biol. 247:808-822). The principal goal of this work is to compare the results obtained in the framework of our solvation model with available experimental and computational data. The findings could be summarized as follows: 1) structural and energetic properties of studied molecules strongly depend on the solvent; membrane-mimetic media significantly promote formation of alpha-helices capable of traversing the bilayer, whereas a polar environment destabilizes alpha-helical conformation via reduction of solvent-exposed surface area and packing; 2) the structures calculated in a membrane-like environment agree with the experimental ones; 3) noticeable differences in conformation of surfactant lipoprotein assessed via Monte Carlo simulation with implicit solvent (this work) and molecular dynamics in explicit solvent were observed; 4) in vacuo simulations do not correctly reproduce protein-membrane interactions, and hence should be avoided in modeling membrane proteins.

Three-dimensional structure of binase in solution.

We present the spatial structure of binase, a small extracellular ribonuclease, derived from 1H-NMR* data in aqueous solution. The total of 20 structures were obtained via torsion angle dynamics using DYANA program with experimental NOE and hydrogen bond distance constraints and phi and chi1 dihedral angle constraints. The final structures were energy minimised with ECEPP/2 potential in FANTOM program. Binase consists of three alpha-helices in N-terminal part (residues 6-16, 26-32 and 41-44), five-stranded antiparallel beta-sheet in C-terminal part (residues 51-55, 70-75, 86-90, 94-99 and 104-108) and two-stranded parallel beta-sheet (residues 22-24 and 49-51). Three loops (residues 36-39, 56-67 and 76-83), which play significant role in biological functioning of binase, are flexible in solution. The differences between binase and barnase spatial structures in solution explain the differences in thermostability of binase, barnase and their hybrids.

Atomic solvation parameters for proteins in a membrane environment. Application to transmembrane alpha-helices.

Several sets of atomic solvation parameters imitating: (i) nonpolar environment of hydrocarbon core of a membrane, (ii) aqueous solution, and (iii) weakly-polar solvents have been developed. The parameters have been incorporated into the ECEPP/2 and CHARMM force fields and employed in non-restrained Monte Carlo and molecular dynamics simulations of membrane-spanning alpha-helical peptides (segment A of bacteriorhodopsin, melittin). Through these simulations, the structure and energetics of the helices have been examined as a function of the solvation term in the potential energy function. For the peptides under study, the set (i) of atomic solvation parameters reveals good retention of the alpha-helical conformation. By contrast, the simulations in vacuum or with the parameters imitating a polar solvent (sets (ii) or (iii)) show fast helix destabilization and tight packing of the structure accompanied by significant decreasing of the surface area accessible to solvent. Increased helical propensity for amino acid residues, population of side-chain rotamers as well as hydrogen-bonding pattern in nonpolar membrane-like environment agree well with available experimental and computational data. The problems related to further applications of the membrane-mimicking sets of atomic solvation parameters to simulations of membrane proteins and peptides are addressed.

Three-dimensional structure of toxin OSK1 from Orthochirus scrobiculosus scorpion venom.

A 600 MHz 1H NMR study of toxin OSK1, blocker of small-conductance Ca2+-activated K+ channels, is presented. The unambiguous sequential assignment of all the protons of the toxin was obtained using TOCSY, DQF-COSY, and NOESY experiments at pH 3.0 (10, 30, and 45 degrees C) in aqueous solution. 3J(N alpha), 3J(alphabeta) vicinal spin coupling constants were determined in high-resolution spectra. The cross-peak volumes in NOESY spectra and the coupling constants were used to define the local structure of the protein by the program HABAS and to generate torsion angle and interproton distance constraints for the program DIANA. Hydrogen-deuterium exchange rates of amide protons showed possible locations of hydrogen bonds. The hydrogen bond acceptors and disulfide bridges between residues 8-28, 14-33, and 18-35 were determined when analyzing distance distribution in preliminary DIANA structures. All constraints were used to obtain a set of 30 structures by DIANA. The resulting rms deviations over 30 structures are 1.30 A for the heavy atoms and 0.42 A for the backbone heavy atoms. The structures were refined by constrained energy minimization using the SYBYL program. Their analysis indicated the existence of the alpha-helix (residues 10-21) slightly distorted at the Cys14 residue, two main strands of the antiparallel beta-sheet (24-29, 32-38), and the extended fragment (2-6). The motif is stabilized by the disulfide bridges in the way, common to all known scorpion toxins. Using the fine spatial toxin structure, alignment of the homologues, mutagenesis analysis, and comparison of scorpion toxin family functions, we delineate some differences significant for the toxin specificity.

Three-dimensional structure of ectatomin from Ectatomma tuberculatum ant venom.

Two-dimensional 1H NMR techniques were used to determine the spatial structure of ectatomin, a toxin from the venom of the ant Ectatomma tuberculatum. Nearly complete proton resonance assignments for two chains of ectatomin (37 and 34 amino acid residues, respectively) were obtained using 2D TOCSY, DQF-COSY and NOESY experiments. The cross-peak volumes in NOESY spectra were used to define the local structure of the protein and generate accurate proton-proton distance constraints employing the MARDIGRAS program. Disulfide bonds were located by analyzing the global fold of ectatomin, calculated with the distance geometry program DIANA. These data, combined with data on the rate of exchange of amide protons with deuterium, were used to obtain a final set of 20 structures by DIANA. These structures were refined by unrestrained energy minimization using the CHARMm program. The resulting rms deviations over 20 structures (excluding the mobile N- and c-termini of each chain) are 0.75 A for backbone heavy atoms, and 1.25 A for all heavy atoms. The conformations of the two chains are similar. Each chain consists of two alpha-helices and a hinge region of four residues; this forms a hairpin structure which is stabilized by disulfide bridges. The hinge regions of the two chains are connected together by a third disulfide bridge. Thus, ectatomin forms a four-alpha-helical bundle structure.

Toxic principle of selva ant venom is a pore-forming protein transformer.

Ectatomin (Ea) is a newly isolated main toxic component of Ectatomma tuberculatum ant venom. Structural and electrophysiological studies were performed with purified Ea. The protein consists of two homologous polypeptide chains (37 and 34 residues) and forms a four alpha-helix bundle in aqueous solution. On insertion into artificial bilayer membranes, two Ea molecules form an ion pore. Our results suggest that the 'inside-out' mechanism of pore formation requires a significant movement of Ea helical parts. The pore formation in the cell membrane might well explain the toxic activity of Ea, not excluding at the same time its intracellular activities.

Three-dimensional structure of proteolytic fragment 163-231 of bacterioopsin determined from nuclear magnetic resonance data in solution.

546 NOESY cross-peak volumes were measured in the two-dimensional NOESY spectrum of proteolytic fragment 163-231 of bacterioopsin in organic solution. These data and 42 detected hydrogen bonds were applied for determining the peptide spatial structure. The fold of the polypeptide chain was determined by local structure analysis, a distance geometry approach and systematic search for energetically allowed side-chain rotamers which are consistent with experimental NOESY cross-peak volumes. The effective rotational correlation time of 6 ns for the molecule was evaluated from optimization of the local structure to meet NOE data and from the dependence on mixing time of the NiH/Ci alpha H cross-peak volumes of the residues in alpha-helical conformation. The resulting structure has two well defined alpha-helical regions, 168-191 and 198-227, with root-mean-square deviation 44 pm and 69 pm, respectively, between the backbone atoms in 14 final energy refined conformations. The alpha-helices correspond to transmembrane segments F and G of bacteriorhodopsin. The segment F contains proline 186, which introduces a kink of about 25 degrees with a disruption of the hydrogen bond with the NH group of the following residue. The segments are connected by a flexible loop region 192-197. Torsion angles chi 1 are unequivocally defined for 62% of side chains in the alpha-helices but half of them differ from electron cryo-microscopy (ECM) model of bacteriorhodopsin, apparently because of the low resolution of ECM. Nevertheless, the F and G segments can be packed as in the ECM model and with side-chain conformations consistent with all NMR data in solution.