Hydrophobic interactions play a key role in the folding and maintenance of the 3-dimensional structure of proteins, as well as in the binding of ligands (e.g. drugs) to protein targets. Therefore, quantitative assessment of spatial hydrophobic (lipophilic) properties of these molecules is indispensable for the development of efficient computational methods in drug design. One possible solution to the problem lies in application of a concept of the 3-dimensional molecular hydrophobicity potential (MHP). The formalism of MHP utilizes a set of atomic physicochemical parameters evaluated from octanol-water partition coefficients (log P) of numerous chemical compounds. It permits detailed assessment of the hydrophobic and/or hydrophilic properties of various parts of molecules and may be useful in analysis of protein-protein and protein-ligand interactions. This review surveys recent applications of MHP-based techniques to a number of biologically relevant tasks. Among them are: (i) Detailed assessment of hydrophobic/hydrophilic organization of proteins; (ii) Application of this data to the modeling of structure, dynamics, and function of globular and membrane proteins, membrane-active peptides, etc. (iii) Employment of the MHP-based criteria in docking simulations for ligands binding to receptors. It is demonstrated that the application of the MHP-based techniques in combination with other molecular modeling tools (e.g. Monte Carlo and molecular dynamics simulations, docking, etc.) permits significant improvement to the standard computational approaches, provides additional important insights into the intimate molecular mechanisms driving protein assembling in water and in biological membranes, and helps in the computer-aided drug discovery process.

The major representatives of Elapidae snake venom, cytotoxins (CTs), share similar three-fingered fold and exert diverse range of biological activities against various cell types. CT-induced cell death starts from the membrane recognition process, whose molecular details remain unclear. It is known, however, that the presence of anionic lipids in cell membranes is one of the important factors determining CT-membrane binding. In this work, we therefore investigated specific interactions between one of the most abundant of such lipids, phosphatidylserine (PS), and CT 4 of Naja kaouthia using a combined, experimental and modeling, approach. It was shown that incorporation of PS into zwitterionic liposomes greatly increased the membrane-damaging activity of CT 4 measured by the release of the liposome-entrapped calcein fluorescent dye. The CT-induced leakage rate depends on the PS concentration with a maximum at approximately 20% PS. Interestingly, the effects observed for PS were much more pronounced than those measured for another anionic lipid, sulfatide. To delineate the potential PS binding sites on CT 4 and estimate their relative affinities, a series of computer simulations was performed for the systems containing the head group of PS and different spatial models of CT 4 in aqueous solution and in an implicit membrane. This was done using an original hybrid computational protocol implementing docking, Monte Carlo and molecular dynamics simulations. As a result, at least three putative PS-binding sites with different affinities to PS molecule were delineated. Being located in different parts of the CT molecule, these anion-binding sites can potentially facilitate and modulate the multi-step process of the toxin insertion into lipid bilayers. This feature together with the diverse binding affinities of the sites to a wide variety of anionic targets on the membrane surface appears to be functionally meaningful and may adjust CT action against different types of cells.

Abstract: 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.