Resonance energy transfer (RET) from anthrylvinyl-labeled phosphatidylcholine (AV-PC) or cardiolipin (AV-CL) to cytochrome c (cyt c) heme moiety was employed to assess the molecular-level details of protein interactions with lipid bilayers composed of PC with 2.5 (CL2.5), 5 (CL5), 10 (CL10), or 20 (CL20) mol % CL under conditions of varying ionic strength and lipid/protein molar ratio. Monte Carlo analysis of multiple data sets revealed a subtle interplay between 1), exchange of the neutral and acidic lipid in the protein-lipid interaction zone; 2), CL transition into the extended conformation; and 3), formation of the hexagonal phase. The switch between these states was found to be controlled by CL content and salt concentration. At ionic strengths ≥40 mM, lipid bilayers with CL fraction not exceeding 5 mol % exhibited the tendency to transform from lamellar to hexagonal phase upon cyt c adsorption, whereas at higher contents of CL, transition into the extended conformation seems to become thermodynamically favorable. At lower ionic strengths, deviations from homogeneous lipid distributions were observed only for model membranes containing 2.5 mol % CL, suggesting the existence of a certain surface potential critical for assembly of lipid lateral domains in protein-lipid systems that may subsequently undergo morphological transformations depending on ambient conditions. These characteristics of cyt c-CL interaction are of great interest, not only from the viewpoint of regulating cyt celectron transfer and apoptotic propensities, but also to elucidate the general mechanisms by which membrane functional activities can be modulated by protein-lipid interactions. © 2010 by the Biophysical Society.