Interactions of the pleckstrin homology domain with phosphatidylinositol phosphate and membranes: characterization via molecular dynamics simulations

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Polyphosphoinositide-Binding Domains: Insights from Peripheral Membrane and Lipid-Transfer Proteins

Within eukaryotic cells, biochemical reactions need to be organized on the surface of membrane compartments that use distinct lipid constituents to dynamically modulate the functions of integral proteins or influence the selective recruitment of peripheral membrane effectors. As a result of these complex interactions, a variety of human pathologies can be traced back to improper communication between proteins and membrane surfaces; either due to mutations that directly alter protein structure or as a result of changes in membrane lipid composition. Among the known structural lipids found in cellular membranes, phosphatidylinositol (PtdIns) is unique in that it also serves as the membrane-anchored precursor of low-abundance regulatory lipids, the polyphosphoinositides (PPIn), which have restricted distributions within specific subcellular compartments. The ability of PPIn lipids to function as signaling platforms relies on both non-specific electrostatic interactions and the selective stereospecific recognition of PPIn headgroups by specialized protein folds. In this chapter, we will attempt to summarize the structural diversity of modular PPIn-interacting domains that facilitate the reversible recruitment and conformational regulation of peripheral membrane proteins. Outside of protein folds capable of capturing PPIn headgroups at the membrane interface, recent studies detailing the selective binding and bilayer extraction of PPIn species by unique functional domains within specific families of lipid-transfer proteins will also be highlighted. Overall, this overview will help to outline the fundamental physiochemical mechanisms that facilitate localized interactions between PPIn lipids and the wide-variety of PPIn-binding proteins that are essential for the coordinate regulation of cellular metabolism and membrane dynamics.

Insights into binding of cholera toxin to GM1 containing membrane

Interactions of cholera toxin (CT) with membrane are associated with the massive secretory diarrhea seen in Asiatic cholera. Ganglioside GM1 has been shown to be responsible for the binding of the B subunit of cholera toxin (CT-B), which then helps CT to pass through the membrane, but the exact mechanism remains to be explored. In this work, we have carried out atomistic scale molecular dynamics simulation to investigate the structural changes of CT upon membrane binding and alteration in membrane structure and dynamics. Starting from the initial structure where the five units of B subunit bind with five GM1, only three of five units remain bound and the whole CT is tilted such that the three binding units are deeper in the membrane. The lipids that are in contact with those units of the CT-B behave differently from the rest of the lipids. Altogether, our results demonstrate the atomistic interaction of CT with GM1 containing lipid membrane and provide a probable mechanism of the early stage alteration of lipid structure and dynamics, which can make a passage for penetration of CT on membrane surface.

The effects of a protein osmolyte on the stability of the integral membrane protein glycerol facilitator

Osmolytes are naturally occurring molecules used by a wide variety of organisms to stabilize proteins under extreme conditions of temperature, salinity, hydrostatic pressure, denaturant concentration, and desiccation. The effects of the osmolyte trimethylamine N-oxide (TMAO) as well as the influence of detergent head group and acyl chain length on the stability of the Escherichia coli integral membrane protein glycerol facilitator (GF) tetramer to thermal and chemical denaturation by sodium dodecyl sulphate (SDS) are reported. TMAO promotes the association of the normally tetrameric α-helical protein into higher order oligomers in dodecyl-maltoside (DDM), but not in tetradecyl-maltoside (TDM), lyso-lauroylphosphatidyl choline (LLPC), or lyso-myristoylphosphatidyl choline (LMPC), as determined by dynamic light scattering (DLS); an octameric complex is particularly stable as indicated by SDS polyacrylamide gel electrophoresis. TMAO increases the heat stability of the GF tetramer an average of 10 °C in the 4 detergents and also protects the protein from denaturation by SDS. However, it did not promote re-association to the tetramer when added to SDS-dissociated protein. TMAO also promotes the formation of rod-like detergent micelles, and DLS was found to be useful for monitoring the structure of the protein and the redistribution of detergent during thermal dissociation of the protein. The protein is more thermally stable in detergents with the phosphatidylcholine head group (LLPC and LMPC) than in the maltoside detergents. The implications of the results for osmolyte mechanism, membrane protein stability, and protein-protein interactions are discussed.

Interactions of peripheral proteins with model membranes as viewed by molecular dynamics simulations

Many cellular signalling and related events are triggered by the association of peripheral proteins with anionic lipids in the cell membrane (e.g. phosphatidylinositol phosphates or PIPs). This association frequently occurs via lipid-binding modules, e.g. pleckstrin homology (PH), C2 and four-point-one, ezrin, radixin, moesin (FERM) domains, present in peripheral and cytosolic proteins. Multiscale simulation approaches that combine coarse-grained and atomistic MD simulations may now be applied with confidence to investigate the molecular mechanisms of the association of peripheral proteins with model bilayers. Comparisons with experimental data indicate that such simulations can predict specific peripheral protein–lipid interactions. We discuss the application of multiscale MD simulation and related approaches to investigate the association of peripheral proteins which contain PH, C2 or FERM-binding modules with lipid bilayers of differing phospholipid composition, including bilayers containing multiple PIP molecules.

Counterion-mediated pattern formation in membranes containing anionic lipids

Most lipid components of cell membranes are either neutral, like cholesterol, or zwitterionic, like phosphatidylcholine and sphingomyelin. Very few lipids, such as sphingosine, are cationic at physiological pH. These generally interact only transiently with the lipid bilayer, and their synthetic analogs are often designed to destabilize the membrane for drug or DNA delivery. However, anionic lipids are common in both eukaryotic and prokaryotic cell membranes. The net charge per anionic phospholipid ranges from -1 for the most abundant anionic lipids such as phosphatidylserine, to near -7 for phosphatidylinositol 3,4,5 trisphosphate, although the effective charge depends on many environmental factors. Anionic phospholipids and other negatively charged lipids such as lipopolysaccharides are not randomly distributed in the lipid bilayer, but are highly restricted to specific leaflets of the bilayer and to regions near transmembrane proteins or other organized structures within the plane of the membrane. This review highlights some recent evidence that counterions, in the form of monovalent or divalent metal ions, polyamines, or cationic protein domains, have a large influence on the lateral distribution of anionic lipids within the membrane, and that lateral demixing of anionic lipids has effects on membrane curvature and protein function that are important for biological control.

Preferred orientations of phosphoinositides in bilayers and their implications in protein recognition mechanisms

Phosphoinositides (PIPs), phosphorylated derivatives of phosphatidylinositol (PI), are essential regulatory lipids involved in various cellular processes, including signal transduction, membrane trafficking, and cytoskeletal remodeling. To gain insight into the protein-PIPs recognition process, it is necessary to study the inositol ring orientation (with respect to the membrane) of PIPs with different phosphorylation states. In this study, 8 PIPs (3 PIP, 2 PIP2, and 3 PIP3) with different phosphorylation and protonation sites have been separately simulated in two mixed bilayers (one with 20% phosphatidylserine (PS) lipids and another with PS lipids switched to phosphatidylcholine (PC) lipids), which roughly correspond to yeast membranes. Uniformity of the bilayer properties including hydrophobic thickness, acyl chain order parameters, and heavy atom density profiles is observed in both PS-contained and PC-enriched membranes due to the same hydrophobic core composition. The relationship between the inositol ring orientation (tilt and rotation angles) and its solvent-accessible surface area indicates that the orientation is mainly determined by its solvation energy. Different PIPs exhibit a clear preference in the inositol ring rotation angle. Surprisingly, a larger proportion of PIPs inositol rings stay closer to the surface of PS-contained membranes compared to PC-enriched ones. Such a difference is rationalized with the formation of more hydrogen bonds between the PS/PI headgroups and the PIPs inositol rings in PS-contained membranes. This hydrogen bond network could be functionally important; thus, the present results can potentially add important and detailed features into the existing protein-PIPs recognition mechanism.

Interactions of the auxilin-1 PTEN-like domain with model membranes result in nanoclustering of phosphatidyl inositol phosphates

Auxilin-1 is a neuron-specific membrane-binding protein involved in a late stage of clathrin-mediated endocytosis. It recruits Hsc70, thus initiating uncoating of the clathrin-coated vesicles. Interactions of auxilin-1 with the vesicle membrane are crucial for this function and are mediated via an N-terminal PTEN-like domain. We have used multiscale molecular dynamics simulations to probe the interactions of the auxilin-1 PTEN-like domain with lipid bilayers containing differing phospholipid composition, including bilayers containing phosphatidyl inositol phosphates. Our results suggest a novel, to our knowledge, model for the auxilin/membrane encounter and subsequent interactions. Negatively charged lipids (especially PIP2) enhance binding of auxilin to lipid bilayers and facilitate its correct orientation relative to the membrane. Mutations in three basic residues (R301E/R307E/K311E) of the C2 subdomain of the PTEN-like domain perturbed its interaction with the bilayer, changing its orientation. The interaction of membrane-bound auxilin-1 PTEN-like domain with negatively charged lipid headgroups results in nanoclustering of PIP2 molecules in the adjacent bilayer leaflet.

When fat is not bad: the regulation of actin dynamics by phospholipid signaling molecules

The actin cytoskeleton plays a key role in the plant morphogenesis and is involved in polar cell growth, movement of subcellular organelles, cell division, and plant defense. Organization of actin cytoskeleton undergoes dynamic remodeling in response to internal developmental cues and diverse environmental signals. This dynamic behavior is regulated by numerous actin-binding proteins (ABPs) that integrate various signaling pathways. Production of the signaling lipids phosphatidylinositol 4,5-bisphosphate and phosphatidic acid affects the activity and subcellular distribution of several ABPs, and typically correlates with increased actin polymerization. Here we review current knowledge of the inter-regulatory dynamics between signaling phospholipids and the actin cytoskeleton in plant cells.

Counterion-mediated cluster formation by polyphosphoinositides

Polyphosphoinositides (PPI) and in particular PI(4,5)P2, are among the most highly charged molecules in cell membranes, are important in many cellular signaling pathways, and are frequently targeted by peripheral polybasic proteins for anchoring through electrostatic interactions. Such interactions between PIP2 and proteins containing polybasic stretches depend on the physical state and the lateral distribution of PIP2 within the inner leaflet of the cell's lipid bilayer. The physical and chemical properties of PIP2 such as pH-dependent changes in headgroup ionization and area per molecule as determined by experiments together with molecular simulations that predict headgroup conformations at various ionization states have revealed the electrostatic properties and phase behavior of PIP2-containing membranes. This review focuses on recent experimental and computational developments in defining the physical chemistry of PIP2 and its interactions with counterions. Ca(2+)-induced changes in PIP2 charge, conformation, and lateral structure within the membrane are documented by numerous experimental and computational studies. A simplified electrostatic model successfully predicts the Ca(2+)-driven formation of PIP2 clusters but cannot account for the different effects of Ca(2+) and Mg(2+) on PIP2-containing membranes. A more recent computational study is able to see the difference between Ca(2+) and Mg(2+) binding to PIP2 in the absence of a membrane and without cluster formation. Spectroscopic studies suggest that divalent cation- and multivalent polyamine-induced changes in the PIP2 lateral distribution in model membrane are also different, and not simply related to the net charge of the counterion. Among these differences is the capacity of Ca(2+) but not other polycations to induce nm scale clusters of PIP2 in fluid membranes. Recent super resolution optical studies show that PIP2 forms nanoclusters in the inner leaflet of a plasma membrane with a similar size distribution as those induced by Ca(2+) in model membranes. The mechanisms by which PIP2 forms nanoclusters and other structures inside a cell remain to be determined, but the unique electrostatic properties of PIP2 and its interactions with multivalent counterions might have particular physiological relevance.

Defining the membrane-associated state of the PTEN tumor suppressor protein

Phosphatase and tensin-homolog deleted on chromosome 10 (PTEN) is a tumor-suppressor protein that regulates phosphatidylinositol 3-kinase (PI3-K) signaling by binding to the plasma membrane and hydrolyzing the 3' phosphate from phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P3) to form phosphatidylinositol (4,5)-bisphosphate (PI(4,5)P2). Several loss-of-function mutations in PTEN that impair lipid phosphatase activity and membrane binding are oncogenic, leading to the development of a variety of cancers, but information about the membrane-associated state of PTEN remains sparse. We have modeled a membrane-associated state of the truncated PTEN structure bound to PI(3,4,5)P3 via multiscale molecular dynamics simulations. We show that the location of the membrane-binding surface agrees with experimental observations and is robust to changes in lipid composition. The level of membrane interaction is substantially reduced in the phosphatase domain for the triple mutant R161E/K163E/K164E, in line with experimental results. We observe clustering of anionic lipids around the C2 domain in preference to the phosphatase domain, suggesting that the C2 domain is involved in nonspecific interactions with negatively charged lipid headgroups. Finally, our simulations suggest that the oncogenicity of the R335L mutation may be due to a reduction in the interaction of the mutant PTEN with anionic lipids.

Cancer associated E17K mutation causes rapid conformational drift in AKT1 pleckstrin homology (PH) domain

BACKGROUND

AKT1 (v-akt murine thymoma viral oncogene homologue 1) kinase is one of the most frequently activated proliferated and survival pathway of cancer. Recently it has been shown that E17K mutation in the Pleckstrin Homology (PH) domain of AKT1 protein leads to cancer by amplifying the phosphorylation and membrane localization of protein. The mutant has shown resistance to AKT1/2 inhibitor VIII drug molecule. In this study we have demonstrated the detailed structural and molecular consequences associated with the activity regulation of mutant protein.

METHODS

The docking score exhibited significant loss in the interaction affinity to AKT1/2 inhibitor VIII drug molecule. Furthermore, the molecular dynamics simulation studies presented an evidence of rapid conformational drift observed in mutant structure.

RESULTS

There was no stability loss in mutant as compared to native structure and the major cation-π interactions were also shown to be retained. Moreover, the active residues involved in membrane localization of protein exhibited significant rise in NHbonds formation in mutant. The rise in NHbond formation in active residues accounts for the 4-fold increase in the membrane localization potential of protein.

CONCLUSION

The overall result suggested that, although the mutation did not induce any stability loss in structure, the associated pathological consequences might have occurred due to the rapid conformational drifts observed in the mutant AKT1 PH domain.

GENERAL SIGNIFICANCE

The methodology implemented and the results obtained in this work will facilitate in determining the core molecular mechanisms of cancer-associated mutations and in designing their potential drug inhibitors.

Emerging methodologies to investigate lipid-protein interactions

Cellular membranes are composed of hundreds of different lipids, ion channels, receptors and scaffolding complexes that act as signalling and trafficking platforms for processes fundamental to life. Cellular signalling and membrane trafficking are often regulated by peripheral proteins, which reversibly interact with lipid molecules in highly regulated spatial and temporal fashions. In most cases, one or more modular lipid-binding domain(s) mediate recruitment of peripheral proteins to specific cellular membranes. These domains, of which more than 10 have been identified since 1989, harbour structurally selective lipid-binding sites. Traditional in vitro and in vivo studies have elucidated how these domains coordinate their cognate lipids and thus how the parent proteins associate with membranes. Cellular activities of peripheral proteins and subsequent physiological processes depend upon lipid binding affinities and selectivity. Thus, the development of novel sensitive and quantitative tools is essential in furthering our understanding of the function and regulation of these proteins. As this field expands into new areas such as computational biology, cellular lipid mapping, single molecule imaging, and lipidomics, there is an urgent need to integrate technologies to detail the molecular architecture and mechanisms of lipid signalling. This review surveys emerging cellular and in vitro approaches for studying protein-lipid interactions and provides perspective on how integration of methodologies directs the future development of the field.

Biophysical and computational studies of membrane penetration by the GRP1 pleckstrin homology domain

The pleckstrin homology (PH) domain of the general receptor for phosphoinositides 1 (GRP1) exhibits specific, high-affinity, reversible binding to phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P(3)) at the plasma membrane, but the nature and extent of the interaction between this bound complex and the surrounding membrane environment remains unclear. Combining equilibrium and nonequilibrium molecular dynamics (MD) simulations, NMR spectroscopy, and monolayer penetration experiments, we characterize the membrane-associated state of GRP1-PH. MD simulations show loops flanking the binding site supplement the interaction with PI(3,4,5)P(3) through multiple contacts with the lipid bilayer. NMR data show large perturbations in chemical shift for these loop regions on binding to PI(3,4,5)P(3)-containing DPC micelles. Monolayer penetration experiments and further MD simulations demonstrate that mutating hydrophobic residues to polar residues in the flanking loops reduces membrane penetration. This supports a "dual-recognition" model of binding, with specific GRP1-PH-PI(3,4,5)P(3) interactions supplemented by interactions of loop regions with the lipid bilayer.

Multiscale simulations suggest a mechanism for integrin inside-out activation

Integrins are large cell-surface adhesion receptors that can be activated to a high affinity state by the formation of an intracellular complex between the integrin β-subunit tail, the membrane, and talin. The F2 and F3 subdomains of the talin head play a key role in formation of this complex. Here, activation of the integrin αIIb/β3 dimer by the talin head domain was probed using multiscale molecular dynamics simulations. A number of novel insights emerge from these studies, including (i) the importance of the integrin αIIb subunit F992 and F993 residues in stabilizing the "off" state of the αIIb/β3 dimer, (ii) a crucial role for negatively charged groups in the F2-F3/membrane interaction, (iii) binding of the talin F2-F3 domain to negatively charged lipid headgroups in the membrane induces a reorientation of the β transmembrane (TM) domain, (iv) an increase in the tilt angle of the β TM domain relative to the bilayer normal helps to destabilize the α/β TM interaction and promote a scissor-like movement of the integrin TM helices. These results, combined with various published experimental observations, suggest a model for the mechanism of inside-out activation of integrins by talin.

Capturing spontaneous partitioning of peripheral proteins using a biphasic membrane-mimetic model

Membrane binding of peripheral proteins, mediated by specialized anchoring domains, is a crucial step for their biological function. Computational studies of membrane insertion, however, have proven challenging and largely inaccessible, due to the time scales required for the complete description of the process, mainly caused by the slow diffusion of the lipid molecules composing the membrane. Furthermore, in many cases, the nature of the membrane "anchor", i.e., the part of the protein that inserts into the membrane, is also unknown. Here, we address some of these issues by developing and employing a simplified representation of the membrane by a biphasic solvent model which we demonstrate can be used efficiently to capture and describe the process of hydrophobic insertion of membrane anchoring domains in all-atom molecular dynamics simulations. Applying the model, we have studied the insertion of the anchoring domain of a coagulation protein (the GLA domain of human protein C), starting from multiple initial configurations varying with regard to the initial orientation and height of the protein with respect to the membrane. In addition to efficiently and consistently identifying the "keel" region as the hydrophobic membrane anchor, within a few nanoseconds each configuration simulated showed a convergent height (2.20 ± 1.04 Å) and angle with respect to the interface normal (23.37 ± 12.48°). We demonstrate that the model can produce the same results as those obtained from a full representation of a membrane, in terms of both the depth of penetration and the orientation of the protein in the final membrane-bound form with an order of magnitude decrease in the required computational time compared to previous models, allowing for a more exhaustive search for the correct membrane-bound configuration.

Study of the Alzheimer's Aβ40 peptide in SDS micelles using molecular dynamics simulations

The interaction of the Alzheimer's amyloid beta peptide, Aβ40, with sodium dodecyl sulfate (SDS) micelles, together with the self-assembly of SDS molecules around the peptide from an initial random distribution were studied using atomistic and coarse-grained (CG) molecular dynamics simulations. In atomistic simulations, the peptide structure in the micelle was characterized by two helical regions connected through a short hinge. The initial structure of the system was shown to affect the simulation results. The atomistic self-assembly of SDS molecules resulted in a 38-molecule micelle around the peptide, along with some globules and individual molecules. Coarse-grained simulation results, however, did not show such a difference, and at the end of all CG simulations, a complete 60-molecule micelle was obtained, with the peptide located at the interface of the micelle with water. The obtained CG radial density profiles and SDS micelle size and shape properties were identical for all CG simulations.

A molecular dynamics investigation of lipid bilayer perturbation by PIP2

Phosphoinositides like phosphatidylinositol 4,5-bisphosphate (PIP(2)) are negatively charged lipids that play a pivotal role in membrane trafficking, signal transduction, and protein anchoring. We have designed a force field for the PIP(2) headgroup using quantum mechanical methods and characterized its properties inside a lipid bilayer using molecular dynamics simulations. Macroscopic properties such as area/headgroup, density profiles, and lipid order parameters calculated from these simulations agree well with the experimental values. However, microscopically, the PIP(2) introduces a local perturbation of the lipid bilayer. The average PIP(2) headgroup orientation of 45 degrees relative to the bilayer normal induces a unique, distance-dependent organization of the lipids that surround PIP(2). The headgroups of these lipids preferentially orient closer to the bilayer normal. This perturbation creates a PIP(2) lipid microdomain with the neighboring lipids. We propose that the PIP(2) lipid microdomain enables the PIP(2) to function as a membrane-bound anchoring molecule.

Cause and effect of melittin-induced pore formation: a computational approach

Melittin embedded in a palmitoyl oleyl phosphatidylcholine bilayer at a high peptide/lipid ratio (1:30) was simulated in the presence of explicit water and ions. The simulation results indicate the incipience of an ion-permeable water pore through collective membrane perturbation by bound peptides. The positively charged residues of melittin not only act as "anchors" but also disrupt the membrane, leading to cell lysis. A detailed analysis of the lipid tail order parameter profile depicts localized membrane perturbation. The lipids in the vicinity of the aqueous cavity adopt a tilted conformation, which allows local bilayer thinning. The prepore thus formed can be considered as the melittin-induced structural defects in the bilayer membrane. Because of the strong cationic nature, the melittin-induced prepore exhibits selectivity toward anions over cations. As Cl(-) ions entered into the prepore, they are electrostatically entrapped by positively charged residues located at its wall. The confined motion of the Cl(-) ions in the membrane interior is obvious from calculated diffusion coefficients. Moreover, reorientation of the local lipids occurs in such a way that few lipid heads along with peptide helices can line the surface of the penetrating aqueous phase. The flipping of lipids argued in favor of melittin-induced toroidal pore over a barrel-stave mechanism. Thus, our result provides atomistic level details of the mechanism of membrane disruption by antimicrobial peptide melittin.

PX- and FYVE-mediated interactions with membranes: simulation studies

Molecular dynamics simulations have been used to explore the interactions of two PI(3)P-binding domains with their PI ligands and with a phospholipid bilayer. Three simulations each of the EEA1-FYVE domain and the p40(phox)-PX domain have been compared: with the protein in an apo state, with a bound Ins(1,3)P(2) molecule, and bound to a PI(3)P molecule embedded in a lipid bilayer. Two main questions were addressed in analysis of the simulations: (i) the location of these domains relative to the lipid bilayer and (ii) their interactions with the lipids, both specific interactions via bound PI(3)P and nonspecific interactions with bilayer phospholipids. Both domains underwent a decrease in dynamic flexibility on binding to the ligand and to the membrane, this being more pronounced for the FYVE domain. Compared to their starting locations [docked to a membrane-inserted PI(3)P molecule], each of the domains penetrated more deeply into the lipid bilayer. For FYVE, nonspecific protein-lipid interactions were formed mainly by the N-terminal hydrophobic region of the protein. For PX, both the alpha1-alpha2 and the beta1-beta2 regions penetrated the bilayer. There appeared to be more marked dynamic fluctuations in hydrogen bonds between basic side chains and PI(3)P for FYVE than for PX, but for both domains, such interactions were maintained throughout the simulations. The simulations agree well with available biophysical data, suggesting this computational method may be used to predict protein-bilayer interactions for other PI-binding proteins.