Supplementary data are available at Bioinformatics online.

Supplementary data are available at Bioinformatics online.

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.

Bacterial cell wall is targeted by many antibiotics. Among them are lantibiotics, which realize their function via interaction with plasma membrane lipid-II molecule — a chemically conserved part of the cell wall synthesis pathway. To investigate structural and dynamic properties of this molecule, we have performed a series of nearly microsecond-long molecular dynamics simulations of lipid-II and some of its analogs in zwitterionic single component and charged mixed simulated phospholipid bilayers (the reference and the mimic of the bacterial plasma membrane, respectively). Extensive analysis revealed that lipid-II forms a unique “amphiphilic pattern” exclusively on the surface of the simulated bacterial membrane (and not in the reference one). We hypothesize that many lantibiotics exploit the conserved features of lipid-II along with characteristic modulation of the bacterial membrane as the “landing site”. This putative recognition mechanism opens new opportunities for studies on lantibiotics action and design of novel armament against resistant bacterial strains.

The new program DASHA is ah efficient implementation of common data processing steps for the protein internal dynamic analysis. The "model-free" parameters and their uncertainties (Lipari G., Szabo A.: J. Am. Chem. Soc. 104, 4546-4559 (1982) can be calculated from an arbitrary combination of experimental data sets (i.e. heteronuclear ~H-~~N or ZH-13C relaxation times and NOE vatues at different spectrometer frequencies). Anisotropy of the molecular rotational diffusion could be also taken into account without introduction of the new adjustable parameters into the spectral density function J(~), provided the structure of the molecule is known. Parameters of chemical (conformational) exchange can be estimated from the CPMG spin-lock frequency dependences (Bloom et al.:

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.

Incorporation of beta-sheet proteins into membrane is studied theoretically for the first time, and the results are validated by the direct experimental data. Using Monte Carlo simulations with implicit membrane, we explore spatial structure, energetics, polarity, and mode of insertion of two cardiotoxins with different membrane-destabilizing activity. Both proteins, classified as P- and S-type cardiotoxins, are found to retain the overall "three-finger" fold interacting with membrane core and lipid/water interface by the tips of the "fingers" (loops). The insertion critically depends upon the structure, hydrophobicity, and electrostatics of certain regions. The simulations reveal apparently distinct binding modes for S- and P-type cardiotoxins via the first loop or through all three loops, respectively. This rationalizes an earlier empirical classification of cardiotoxins into S- and P-type, and provides a basis for the analysis of experimental data on their membrane affinities. Accomplished with our previous simulations of membrane alpha-helices, the computational method may be used to study partitioning of proteins with diverse folds into lipid bilayers.

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.