Roman G. Efremov

Spatial Structure of the Dimeric Transmembrane Domain of the Growth Factor Receptor ErbB2 Presumably Corresponding to the Receptor Active State

Proper lateral dimerization of the transmembrane domains of receptor tyrosine kinases is required for biochemical signal transduction across the plasma membrane. The spatial structure of the dimeric transmembrane domain of the growth factor receptor ErbB2 embedded into lipid bicelles was obtained by solution NMR, followed by molecular dynamics relaxation in an explicit lipid bilayer. ErbB2 transmembrane segments associate in a right-handed α-helical bundle through the N-terminal tandem GG4-like motif Thr652-X3-Ser656-X3-Gly660, providing an explanation for the pathogenic power of some oncogenic mutations.

Unique dimeric structure of BNip3 transmembrane domain suggests membrane permeabilization as a cell death trigger.

BNip3 is a prominent representative of apoptotic Bcl-2 proteins with rather unique properties initiating an atypical programmed cell death pathway resembling both necrosis and apoptosis. Many Bcl-2 family proteins modulate the permeability state of the outer mitochondrial membrane by forming homo- and hetero-oligomers. The structure and dynamics of the homodimeric transmembrane domain of BNip3 were investigated with the aid of solution NMR in lipid bicelles and molecular dynamics energy relaxation in an explicit lipid bilayer. The right-handed parallel helix-helix structure of the domain with a hydrogen bond-rich His-Ser node in the middle of the membrane, accessibility of the node for water, and continuous hydrophilic track across the membrane suggest that the domain can provide an ion-conducting pathway through the membrane. Incorporation of the BNip3 transmembrane domain into an artificial lipid bilayer resulted in pH-dependent conductivity increase. A possible biological implication of the findings in relation to triggering necrosis-like cell death by BNip3 is discussed.

PLATINUM: a web tool for analysis of hydrophobic/hydrophilic organization of biomolecular complexes

UNLABELLED The PLATINUM (Protein-Ligand ATtractions Investigation NUMerically) web service is designed for analysis and visualization of hydrophobic/hydrophilic properties of biomolecules supplied as 3D-structures. Furthermore, PLATINUM provides a number of tools for quantitative characterization of the hydrophobic/hydrophilic match in biomolecular complexes e.g. in docking poses. These complement standard scoring functions. The calculations are based on the concept of empirical Molecular Hydrophobicity Potential (MHP). AVAILABILITY The PLATINUM web tool as well as detailed documentation and tutorial are available free of charge for academic users at http://model.nmr.ru/platinum/. PLATINUM requires Java 5 or higher and Adobe Flash Player 9. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Spatial Structure and pH-dependent Conformational Diversity of Dimeric Transmembrane Domain of the Receptor Tyrosine Kinase EphA1

Eph receptors are found in a wide variety of cells in developing and mature tissues and represent the largest family of receptor tyrosine kinases, regulating cell shape, movements, and attachment. The receptor tyrosine kinases conduct biochemical signals across plasma membrane via lateral dimerization in which their transmembrane domains play an important role. Structural-dynamic properties of the homodimeric transmembrane domain of the EphA1 receptor were investigated with the aid of solution NMR in lipid bicelles and molecular dynamics in explicit lipid bilayer. EphA1 transmembrane segments associate in a right-handed parallel α-helical bundle, region (544-569)2, through the N-terminal glycine zipper motif A550X3G554X3G558. Under acidic conditions, the N terminus of the transmembrane helix is stabilized by an N-capping box formed by the uncharged carboxyl group of Glu547, whereas its deprotonation results in a rearrangement of hydrogen bonds, fractional unfolding of the helix, and a realignment of the helix-helix packing with appearance of additional minor dimer conformation utilizing seemingly the C-terminal GG4-like dimerization motif A560X3G564. This can be interpreted as the ability of the EphA1 receptor to adjust its response to ligand binding according to extracellular pH. The dependence of the pKa value of Glu547 and the dimer conformational equilibrium on the lipid head charge suggests that both local environment and membrane surface potential can modulate dimerization and activation of the receptor. This makes the EphA1 receptor unique among the Eph family, implying its possible physiological role as an “extracellular pH sensor,” and can have relevant physiological implications.

Left-Handed Dimer of EphA2 Transmembrane Domain: Helix Packing Diversity among Receptor Tyrosine Kinases

The Eph receptor tyrosine kinases and their membrane-bound ephrin ligands control a diverse array of cell-cell interactions in the developing and adult organisms. During signal transduction across plasma membrane, Eph receptors, like other receptor tyrosine kinases, are involved in lateral dimerization and subsequent oligomerization presumably with proper assembly of their single-span transmembrane domains. Spatial structure of dimeric transmembrane domain of EphA2 receptor embedded into lipid bicelle was obtained by solution NMR, showing a left-handed parallel packing of the transmembrane helices (535-559)(2). The helices interact through the extended heptad repeat motif L(535)X(3)G(539)X(2)A(542)X(3)V(546)X(2)L(549) assisted by intermolecular stacking interactions of aromatic rings of (FF(557))(2), whereas the characteristic tandem GG4-like motif A(536)X(3)G(540)X(3)G(544) is not used, enabling another mode of helix-helix association. Importantly, a similar motif AX(3)GX(3)G as was found is responsible for right-handed dimerization of transmembrane domain of the EphA1 receptor. These findings serve as an instructive example of the diversity of transmembrane domain formation within the same family of protein kinases and seem to favor the assumption that the so-called rotation-coupled activation mechanism may take place during the Eph receptor signaling. A possible role of membrane lipid rafts in relation to Eph transmembrane domain oligomerization and Eph signal transduction was also discussed.

A Solvent Model for Simulations of Peptides in Bilayers. I. Membrane-Promoting α-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.

Human copper-transporting ATPase ATP7B (the Wilson's disease protein): biochemical properties and regulation.

Wilsons disease protein (WNDP) is a product of a gene ATP7B that is mutated in patients with Wilsons disease, a severe genetic disorder with hepatic and neurological manifestations caused by accumulation of copper in the liver and brain. In a cell, WNDP transports copper across various cell membranes using energy of ATP-hydrolysis. Copper regulates WNDP at several levels, modulating its catalytic activity, posttranslational modification, and intracellular localization. This review summarizes recent studies on enzymatic function and copper-dependent regulation of WNDP. Specifically, we describe the molecular architecture and major biochemical properties of WNDP, discuss advantages of the recently developed functional expression of WNDP in insect cells, and summarize the results of the ligand-binding studies and molecular modeling experiments for the ATP-binding domain of WNDP. In addition, we speculate on how copper binding may regulate the activity and intracellular distribution of WNDP, and what role the human copper chaperone Atox1 may play in these processes.

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.

Spatial structure of zervamicin IIB bound to DPC micelles: implications for voltage-gating.

Zervamicin IIB is a 16-amino acid peptaibol that forms voltage-dependent ion channels with multilevel conductance states in planar lipid bilayers and vesicular systems. The spatial structure of zervamicin IIB bound to dodecylphosphocholine micelles was studied by nuclear magnetic resonance spectroscopy. The set of 20 structures obtained has a bent helical conformation with a mean backbone root mean square deviation value of approximately 0.2 A and resembles the structure in isotropic solvents (Balashova et al., 2000. NMR structure of the channel-former zervamicin IIB in isotropic solvents. FEBS Lett 466:333-336). The N-terminus represents an alpha-helix, whereas the C-terminal part has a mixed 3(10)/alpha(R) hydrogen-bond pattern. In the anisotropic micelle environment, the bending angle on Hyp10 (23 degrees) is smaller than that (47 degrees) in isotropic solvents. In the NOESY (Nuclear Overhauser Effect Spectroscopy) spectra, the characteristic attenuation of the peptide signals by 5- and 16-doxylstearate relaxation probes indicates a peripheral mode of the peptaibol binding to the micelle with the N-terminus immersed slightly deeper into micelle interior. Analysis of the surface hydrophobicity reveals that the zervamicin IIB helix is amphiphilic and well suited to formation of a tetrameric transmembrane bundle, according to the barrel-stave mechanism. The results are discussed in a context of voltage-driven peptaibol insertion into membrane.

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, α‐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.