Probing the structure of the Ff bacteriophage major coat protein transmembrane helix dimer by solution NMR

The transmembrane (TM) segment of the major coat protein from Ff bacteriophage has been extensively studied as an example of dimerization in detergent and lipid bilayer systems. However, almost all the information regarding this interaction has been gained through mutagenesis studies, with little direct structural information being available. To this end solution NMR has the potential to provide new insights into structure of the dimer. In order to evaluate the utility of this approach we have studied a selectively ¹⁵N-labeled peptide containing the TM segment of MCP (MCP_TM) by solution NMR. This peptide was found to give rise to detergent concentration-dependent spectra that were assigned to monomeric and dimeric forms. The standard free energy of this interaction in SDS was estimated from these spectra and found to be consistent with weak but specific dimerization. In addition, similar spectra could be obtained in β-octyl glucoside with intermolecular paramagnetic relaxation experiments demonstrating a parallel arrangement of TM helices in the dimer. In both detergents backbone chemical shift differences between monomeric and dimeric forms of MCP_TM showed that the largest changes occur around its GXXXG motif. The resulting structural model is consistent with observations made for MCP mutants previously characterized in biological membranes, opening the door to detailed structural characterization of this form of MCP. These results also have general implications for the study of weakly interacting TM segments by solution NMR since the use of similar sample conditions should allow structural data to be accessed for oligomeric states from a wide range systems that undergo biologically relevant but weak associations in the membrane.     

Towards a structural understanding of the smallest known oncoprotein: investigation of the bovine papillomavirus E5 protein using solution-state NMR

The homodimeric E5 protein from bovine papillomavirus activates the platelet-derived growth factor β receptor through transmembrane (TM) helix-helix interactions leading to uncontrolled cell growth. Detailed structural information for the E5 dimer is essential if we are to uncover its unique mechanism of action. In vivo mutagenesis has been used to identify residues in the TM domain critical for dimerization, and we previously reported that a truncated synthetic E5 peptide forms dimers via TM domain interactions. Here we extend this work with the first application of high-resolution solution-state NMR to the study of the E5 TM domain in SDS micelles. Using selectively 15N-labelled peptides, we first probe sample homogeneity revealing two predominate species, which we interpret to be monomer and dimer. The equilibrium between the two states is shown to be dependent on detergent concentration, revealed by intensity shifts between two sets of peaks in 15N-(1)H HSQC experiments, highlighting the importance of sample preparation when working with these types of proteins. This information is used to estimate a free energy of association (ΔGx°=-3.05 kcal mol(-1)) for the dimerization of E5 in SDS micelles. In addition, chemical shift changes have been observed that indicate a more pronounced change in chemical environment for those residues expected to be at the dimer interface in vivo versus those that are not. Thus we are able to demonstrate our in vitro dimer is comparable to that defined in vivo, validating the biological significance of our synthetic peptide and providing a solid foundation upon which to base further structural studies. Using detergent concentration to modulate oligomeric state and map interfacial residues by NMR could prove useful in the study of other homo-oligomeric transmembrane proteins.     

Spatial structure and dimer--monomer equilibrium of the ErbB3 transmembrane domain in DPC micelles

In present work the interaction of two TM α-helices of the ErbB3 receptor tyrosine kinase from the ErbB or HER family (residues 639-670) was studied by means of NMR spectroscopy in a membrane-mimicking environment provided by the DPC micelles. The ErbB3 TM segment appeared to form a parallel symmetric dimer in a left-handed orientation. The interaction between TM spans is accomplished via the non-standard motif and is supported by apolar contacts of bulky side chains and by stacking of aromatic rings together with π-cation interactions of Phe and Arg side chains. The investigation of the dimer--monomer equilibrium revealed thermodynamic properties of the assembly and the presence of two distinct regimes of the dimerization at low and at high peptide/detergent ratio. It was found that the detergent in case of ErbB3 behaves not as an ideal solvent, thus affecting the dimer--monomer equilibrium. Such behavior may account for the problems occurring with the refolding and stability of multispan helical membrane proteins in detergent solutions. The example of ErbB3 allows us to conclude that the thermodynamic parameters of dimerization, measured in micelles for two different helical pairs, cannot be compared without the investigation of their dependence on detergent concentration.     

Without its N-finger, the main protease of severe acute respiratory syndrome coronavirus can form a novel dimer through its C-terminal domain

The main protease (M(pro)) of severe acute respiratory syndrome coronavirus (SARS-CoV) plays an essential role in the extensive proteolytic processing of the viral polyproteins (pp1a and pp1ab), and it is an important target for anti-SARS drug development. It was found that SARS-CoV M(pro) exists in solution as an equilibrium of both monomeric and dimeric forms, and the dimeric form is the enzymatically active form. However, the mechanism of SARS-CoV M(pro) dimerization, especially the roles of its N-terminal seven residues (N-finger) and its unique C-terminal domain in the dimerization, remain unclear. Here we report that the SARS-CoV M(pro) C-terminal domain alone (residues 187 to 306; M(pro)-C) is produced in Escherichia coli in both monomeric and dimeric forms, and no exchange could be observed between them at room temperature. The M(pro)-C dimer has a novel dimerization interface. Meanwhile, the N-finger deletion mutant of SARS-CoV M(pro) also exists as both a stable monomer and a stable dimer, and the dimer is formed through the same C-terminal-domain interaction as that in the M(pro)-C dimer. However, no C-terminal domain-mediated dimerization form can be detected for wild-type SARS-CoV M(pro). Our study results help to clarify previously published controversial claims about the role of the N-finger in SARS-CoV M(pro) dimerization. Apparently, without the N-finger, SARS-CoV M(pro) can no longer retain the active dimer structure; instead, it can form a new type of dimer which is inactive. Therefore, the N-finger of SARS-CoV M(pro) is not only critical for its dimerization but also essential for the enzyme to form the enzymatically active dimer.     

Rough energy landscapes in protein folding: dimeric E. coli Trp repressor folds through three parallel channels

The folding mechanism of the dimeric Escherichia coli Trp repressor (TR) is a kinetically complex process that involves three distinguishable stages of development. Following the formation of a partially folded, monomeric ensemble of species, within 5 ms, folding to the native dimer is controlled by three kinetic phases. The rate-limiting step in each phase is either a non-proline isomerization reaction or a dimerization reaction, depending on the final denaturant concentration. Two approaches have been employed to test the previously proposed folding mechanism of TR through three parallel channels: (1) unfolding double-jump experiments demonstrate that all three folding channels lead directly to native dimer; and (2) the differential stabilization of the transition state for the final step in folding and the native dimer, by the addition of salt, shows that all three channels involve isomerization of a dimeric species. A refined model for the folding of Trp repressor is presented, in which all three channels involve a rapid dimerization reaction between partially folded monomers followed by the isomerization of the dimeric intermediates to yield native dimer. The ensemble of partially folded monomers can be captured at equilibrium by low pH; one-dimensional proton NMR spectra at pH 2.5 demonstrate that monomers exist in two distinct, slowly interconverting conformations. These data provide a potential structural explanation for the three-channel folding mechanism of TR: random association of two different monomeric forms, which are distinguished by alternative packing modes of the core dimerization domain and the DNA-binding, helix-turn-helix, domain. One, perhaps both, of these packing modes contains non-native contacts.     

NMR analysis of the monomeric form of a mutant unliganded bovine neurophysin: comparison with the crystal structure of a neurophysin dimer

Determination of the structure of the unliganded monomeric state of neurophysin is central to an understanding of the allosteric relationship between neurophysin peptide-binding and dimerization. We examined this state by NMR, using the weakly dimerizing H80E mutant of bovine neurophysin-I. The derived structure, to which more than one conformer appeared to contribute, was compared with the crystal structure of the unliganded des 1-6 bovine neurophysin-II dimer. Significant conformational differences between the two proteins were evident in the orientation of the 3,10 helix, in the 50-58 loop, in beta-turns, and in specific intrachain contacts between amino- and carboxyl domains. However, both had similar secondary structures, in independent confirmation of earlier circular dichroism studies. Previously suggested interactions between the amino terminus and the 50-58 loop in the monomer were also confirmed. Comparison of the observed differences between the two proteins with demonstrated effects of dimerization on the NMR spectrum of bovine neurophysin-I, and preliminary investigation of the effects of dimerization on H80E spectra, allowed tentative distinction between the contributions of sequence and self-association differences to the difference in conformation. Regions altered by dimerization encompass most binding site residues, providing a potential explanation of differences in binding affinity between the unliganded monomeric and dimeric states. Differences between monomer and dimer states in turns, interdomain contacts, and within the interdomain segment of the 50-58 loop suggest that the effects of dimerization on intrasubunit conformation reflect the need to adjust the relative positions of the interface segments of the two domains for optimal interaction with the adjacent subunit and/or reflect the dual role of some residues as participants both at the interface and in interdomain contacts.     

A designed point mutant in fis1 disrupts dimerization and mitochondrial fission

Mitochondrial and peroxisomal fission are essential processes with defects resulting in cardiomyopathy and neonatal lethality. Central to organelle fission is Fis1, a monomeric tetratricopeptide repeat (TPR)-like protein whose role in assembly of the fission machinery remains obscure. Two nonfunctional, Saccharomyces cerevisiae Fis1 mutants (L80P or E78D/I85T/Y88H) were previously identified in genetic screens. Here, we find that these two variants in the cytosolic domain of Fis1 (Fis1ΔTM) are unexpectedly dimeric. A truncation variant of Fis1ΔTM that lacks an N-terminal regulatory domain is also found to be dimeric. The ability to dimerize is a property innate to the native Fis1ΔTM amino acid sequence as we find this domain is dimeric after transient exposure to elevated temperature or chemical denaturants and is kinetically trapped at room temperature. This is the first demonstration of a specific self-association in solution for the Fis1 cytoplasmic domain. We propose a three-dimensional domain-swapped model for dimerization that is validated by a designed mutation, A72P, which potently disrupts dimerization of wild-type Fis1. A72P also disrupts dimerization of nonfunctional variants, indicating a common structural basis for dimerization. The obligate monomer variant A72P, like the dimer-promoting variants, is nonfunctional in fission, consistent with a model in which Fis1 activity depends on its ability to interconvert between monomer and dimer species. These studies suggest a new functionally important manner in which TPR-containing proteins may reversibly self-associate.     

The monomer-dimer equilibrium of stromal cell-derived factor-1 (CXCL 12) is altered by pH, phosphate, sulfate, and heparin

Chemokines, like stromal cell-derived factor-1 (SDF1/CXCL12), are small secreted proteins that signal cells to migrate. Because SDF1 and its receptor CXCR4 play important roles in embryonic development, cancer metastasis, and HIV/AIDS, this chemokine signaling system is the subject of intense study. However, it is not known whether the monomeric or dimeric structure of SDF1 is responsible for signaling in vivo. Previous structural studies portrayed the SDF1 structure as either strictly monomeric in solution or dimeric when crystallized. Here, we report two-dimensional NMR, pulsed-field gradient diffusion and fluorescence polarization measurements at various SDF1 concentrations, solution conditions, and pH. These results demonstrate that SDF1 can form a dimeric structure in solution, but only at nonacidic pH when stabilizing counterions are present. Thus, while the previous NMR structural studies were performed under acidic conditions that strongly promote the monomeric state, crystallographic studies used nonacidic buffer conditions that included divalent anions shown here to promote dimerization. This pH-sensitive aggregation behavior is explained by a dense cluster of positively charged residues at the SDF1 dimer interface that includes a histidine side chain at its center. A heparin disaccharide shifts the SDF1 monomer-dimer equilibrium in the same manner as other stabilizing anions, suggesting that glycosaminoglycan binding may be coupled to SDF1 dimerization in vivo.     

The folding energy landscape of the dimerization domain of Escherichia coli Trp repressor: a joint experimental and theoretical investigation

Enhanced structural insights into the folding energy landscape of the N-terminal dimerization domain of Escherichia coli tryptophan repressor, [2-66]2 TR, were obtained from a combined experimental and theoretical analysis of its equilibrium folding reaction. Previous studies have shown that the three intertwined helices in [2-66]2 TR are sufficient to drive the formation of a stable dimer for the full-length protein, [2-107]2 TR. The monomeric and dimeric folding intermediates that appear during the folding reactions of [2-66]2 TR have counterparts in the folding mechanism of the full-length protein. The equilibrium unfolding energy surface on which the folding and dimerization reactions occur for [2-66]2 TR was examined with a combination of native-state hydrogen exchange analysis, pepsin digestion and matrix-assisted laser/desorption mass spectrometry performed at several concentrations of protein and denaturant. Peptides corresponding to all three helices in [2-66]2 TR show multi-layered protection patterns consistent with the relative stabilities of the dimeric and monomeric folding intermediates. The observation of protection exceeding that offered by the dimeric intermediate in segments from all three helices implies that a segment-swapping mechanism may be operative in the monomeric intermediate. Protection greater than that expected from the global stability for a single amide hydrogen in a peptide from the C-helix possibly and another from the A-helix may reflect non-random structure, possibly a precursor for segment swapping, in the urea-denatured state. Native topology-based model simulations that correspond to a funnel energy landscape capture both the monomeric and dimeric intermediates suggested by the HX MS data and provide a rationale for the progressive acquisition of secondary structure in their conformational ensembles.     

Polar residue tagging of transmembrane peptides

Studies that focus on packing interactions between transmembrane (TM) helices in membrane proteins would greatly benefit from the ability to investigate their association and packing interactions in multi-spanning TM domains. However, the production, purification, and characterization of such units have been impeded by their high intrinsic hydrophobicity. We describe the polar tagging approach to biophysical analysis of TM segment peptides, where incorporation of polar residues of suitable type and number at one or both peptide N- and C-termini can serve to counterbalance the apolar nature of a native TM segment, and render it aqueous-soluble. Using the native TM sequences of the human erythrocyte protein glycophorin A (GpA) and bacteriophage M13 major coat protein (MCP), properties of tags such as Lys, His, Asp, sarcosine, and Pro-Gly are evaluated, and general procedures for tagging a given TM segment are presented. Gel-shift assays on sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) establish that various tagged GpA TM segments spontaneously insert into micellar membranes, and exhibit native TM dimeric states. Sedimentation equilibrium analytical centrifugation is used to confirm that Lys-tagged GpA peptides retain the native dimer state. Two-dimensional nuclear magnetic resonance (NMR) spectroscopy studies on Lys-tagged TM MCP peptides selectively enriched with N-15 illustrate the usefulness of this system for evaluating monomer-dimer equilibria in micelle environments. The overall results suggest that polar-tagging of hydrophobic (TM) peptides approach constitutes a valuable tool for the study of protein-protein interactions in membranes.     

The domain-swapped dimer of cyanovirin-N is in a metastable folded state: reconciliation of X-ray and NMR structures

The structure of the potent HIV-inactivating protein cyanovirin-N was previously found by NMR to be a monomer in solution and a domain-swapped dimer by X-ray crystallography. Here we demonstrate that, in solution, CV-N can exist both in monomeric and in domain-swapped dimeric form. The dimer is a metastable, kinetically trapped structure at neutral pH and room temperature. Based on orientational NMR constraints, we show that the domain-swapped solution dimer is similar to structures in two different crystal forms, exhibiting solely a small reorientation around the hinge region. Mutation of the single proline residue in the hinge to glycine significantly stabilizes the protein in both its monomeric and dimeric forms. By contrast, mutation of the neighboring serine to proline results in an exclusively dimeric protein, caused by a drastic destabilization of the monomer.     

Structural role of PufX in the dimerization of the photosynthetic core complex of Rhodobacter sphaeroides

Monomeric and dimeric PufX-containing core complexes have been purified from membranes of wild-type Rhodobacter sphaeroides. Reconstitution of both samples by detergent removal in the presence of lipids leads to the formation of two-dimensional crystals constituted of dimeric core complexes. Two-dimensional crystals were further analyzed by cryoelectron microscopy and atomic force microscopy. A projection map at 26-A resolution reveals that core complexes assemble in an "S"-shaped dimeric complex. Each core complex is composed of one reaction center, 12 light-harvesting 1 alpha/beta-heterodimers, and one PufX protein. The light-harvesting 1 assemblies are open with a gap of density of approximately 30-A width and surround oriented reaction centers. A maximum density is found at the dimer junction. Based on the projection map, a model is proposed, in which the two PufX proteins are located at the dimer junction, consistent with the finding of dimerization of monomeric core complexes upon reconstitution. This localization of PufX in the core complex implies that PufX is the structural key for the dimer complex formation rather than a channel-forming protein for the exchange of ubiquinone/ubiquinol between the reaction center and the cytochrome bc1 complex.     

Structural insights into the intertwined dimer of fyn SH2

Src homology 2 domains are interaction modules dedicated to the recognition of phosphotyrosine sites incorporated in numerous proteins found in intracellular signaling pathways. Here we provide for the first time structural insight into the dimerization of Fyn SH2 both in solution and in crystalline conditions, providing novel crystal structures of both the dimer and peptide‐bound structures of Fyn SH2. Using nuclear magnetic resonance chemical shift analysis, we show how the peptide is able to eradicate the dimerization, leading to monomeric SH2 in its bound state. Furthermore, we show that Fyn SH2's dimer form differs from other SH2 dimers reported earlier. Interestingly, the Fyn dimer can be used to construct a completed dimer model of Fyn without any steric clashes. Together these results extend our understanding of SH2 dimerization, giving structural details, on one hand, and suggesting a possible physiological relevance of such behavior, on the other hand.     

Nonoxidative cadmium-dependent dimerization of Cd7-metallothionein from rabbit liver.

The effect of free Cd(II) ions on monomeric Cd7-metallothionein-2 (MT) from rabbit liver has been studied. Slow, concentration-dependent dimerization of this protein was observed by gel filtration chromatographic studies. The dimeric MT form, isolated by gel filtration, contains approximately two additional and more weakly bound Cd(II) ions per monomer. The incubation of MT dimers with complexing agents EDTA and 2-mercaptoethanol leads to the dissociation of dimers to monomers. The results of circular dichroism (CD) and electronic absorption studies indicate that the slow dimerization process is preceded by an initial rapid Cd-induced rearrangement of the monomeric Cd7-MT structure. The 113Cd NMR spectrum of the MT dimer revealed only four 113Cd resonances at chemical shift positions similar to those observed for the Cd4 cluster of the well-characterized monomeric 113Cd7-MT. This result suggests that on dimer formation major structural changes occur in the original three-metal cluster domain of Cd7-MT.     

Probing the initial stage of aggregation of the Abeta(10-35)-protein: assessing the propensity for peptide dimerization

Characterization of the early stages of peptide aggregation is of fundamental importance in elucidating the mechanism of the formation of deposits associated with amyloid disease. The initial step in the pathway of aggregation of the Abeta-protein, whose monomeric NMR structure is known, was studied through the simulation of the structure and stability of the peptide dimer in aqueous solution. A protocol based on shape complementarity was used to generate an assortment of possible dimer structures. The structures generated based on shape complementarity were evaluated using rapidly computed estimates of the desolvation and electrostatic interaction energies to identify a putative stable dimer structure. The potential of mean force associated with the dimerization of the peptides in aqueous solution was computed for both the hydrophobic and the electrostatic driven forces using umbrella sampling and classical molecular dynamics simulation at constant temperature and pressure with explicit solvent and periodic boundary conditions. The comparison of the two free energy profiles suggests that the structure of the peptide dimer is determined by the favorable desolvation of the hydrophobic residues at the interface. Molecular dynamics trajectories originating from two putative dimer structures indicate that the peptide dimer is stabilized primarily through hydrophobic interactions, while the conformations of the peptide monomers undergo substantial structural reorganization in the dimerization process. The finding that the phi-dimer may constitute the ensemble of stable Abeta(10-35) dimer has important implications for fibril formation. In particular, the expulsion of water molecules at the interface might be a key event, just as in the oligomerization of Abeta(16-22) fragments. We conjecture that events prior to the nucleation process themselves might involve crossing free energy barriers which depend on the peptide-peptide and peptide-water interactions. Consistent with existing experimental studies, the peptides within the ensemble of aggregated states show no signs of formation of secondary structure.     

Modulation of glycophorin A transmembrane helix interactions by lipid bilayers: molecular dynamics calculations

Starting from the glycophorin A dimer structure determined by NMR, we performed simulations of both dimer and monomer forms in explicit lipid bilayers with constant normal pressure, lateral area, and temperature using the CHARMM potential. Analysis of the trajectories in four different lipids reveals how lipid chain length and saturation modulate the structural and energetic properties of transmembrane helices. Helix tilt, helix-helix crossing angle, and helix accessible volume depend on lipid type in a manner consistent with hydrophobic matching concepts: the most relevant lipid property appears to be the bilayer thickness. Although the net helix-helix interaction enthalpy is strongly attractive, analysis of residue-residue interactions reveals significant unfavorable electrostatic repulsion between interfacial glycine residues previously shown to be critical for dimerization. Peptide volume is nearly conserved upon dimerization in all lipid types, indicating that the monomeric helices pack equally well with lipid as dimer helices do with one another. Enthalpy calculations indicate that the helix-environment interaction energy is lower in the dimer than in the monomer form, when solvated by unsaturated lipids. In all lipid environments there is a marked preference for lipids to interact with peptide predominantly through one rather than both acyl chains. Although our trajectories are not long enough to allow a full thermodynamic treatment, these results demonstrate that molecular dynamics simulations are a powerful method for investigating the protein-protein, protein-lipid, and lipid-lipid interactions that determine the structure, stability and dynamics of transmembrane alpha-helices in membranes.     

Tyrosine Sulfation of Chemokine Receptor CCR2 Enhances Interactions with Both Monomeric and Dimeric Forms of the Chemokine Monocyte Chemoattractant Protein-1 (MCP-1)

Chemokine receptors are commonly post-translationally sulfated on tyrosine residues in their N-terminal regions, the initial site of binding to chemokine ligands. We have investigated the effect of tyrosine sulfation of the chemokine receptor CCR2 on its interactions with the chemokine monocyte chemoattractant protein-1 (MCP-1/CCL2). Inhibition of CCR2 sulfation, by growth of expressing cells in the presence of sodium chlorate, significantly reduced the potency for MCP-1 activation of CCR2. MCP-1 exists in equilibrium between monomeric and dimeric forms. The obligate monomeric mutant MCP-1(P8A) was similar to wild type MCP-1 in its ability to induce leukocyte recruitment in vivo, whereas the obligate dimeric mutant MCP-1(T10C) was less effective at inducing leukocyte recruitment in vivo. In two-dimensional NMR experiments, sulfated peptides derived from the N-terminal region of CCR2 bound to both the monomeric and dimeric forms of wild type MCP-1 and shifted the equilibrium to favor the monomeric form. Similarly, MCP-1(P8A) bound more tightly than MCP-1(T10C) to the CCR2-derived sulfopeptides. NMR chemical shift mapping using the MCP-1 mutants showed that the sulfated N-terminal region of CCR2 binds to the same region (N-loop and β3-strand) of both monomeric and dimeric MCP-1 but that binding to the dimeric form also influences the environment of chemokine N-terminal residues, which are involved in dimer formation. We conclude that interaction with the sulfated N terminus of CCR2 destabilizes the dimerization interface of inactive dimeric MCP-1, thus inducing dissociation to the active monomeric state.     

Hepatitis C Virus RNA Replication and Virus Particle Assembly Require Specific Dimerization of the NS4A Protein Transmembrane Domain

Hepatitis C virus (HCV) NS4A is a single-pass transmembrane (TM) protein essential for viral replication and particle assembly. The sequence of the NS4A TM domain is highly conserved, suggesting that it may be important for protein-protein interactions. To test this hypothesis, we measured the potential dimerization of the NS4A TM domain in a well-characterized two-hybrid TM protein interaction system. The NS4A TM domain exhibited a strong homotypic interaction that was comparable in affinity to glycophorin A, a well-studied human blood group antigen that forms TM homodimers. Several mutations predicted to cluster on a common surface of the NS4A TM helix caused significant reductions in dimerization, suggesting that these residues form an interface for NS4A dimerization. Mutations in the NS4A TM domain were further examined in the JFH-1 genotype 2a replicon system; importantly, all mutations that destabilized NS4A dimers also caused defects in RNA replication and/or virus assembly. Computational modeling of NS4A TM interactions suggests a right-handed dimeric interaction of helices with an interface that is consistent with the mutational effects. Furthermore, defects in NS4A oligomerization and virus particle assembly of two mutants were rescued by NS4A A15S, a TM mutation recently identified through forward genetics as a cell culture-adaptive mutation. Together, these data provide the first example of a functionally important TM dimer interface within an HCV nonstructural protein and reveal a fundamental role of the NS4A TM domain in coordinating HCV RNA replication and virus particle assembly.     

Structural and energetic requirements for a second binding site at the dimeric β-lactoglobulin interface

Through experimental and theoretical approaches, it has been shown that bovine β-lactoglobulin (βlg) uses its hydrophobic cavity or calyx as the primary binding site for hydrophobic molecules, whereas the existence of a second ligand binding site at the dimeric interface has only been structurally identified for vitamin D3 (VD3). This binding exists even in the thermally denatured state, suggesting the prevalence of this secondary site. Although crystallographic experiments have suggested that VD3 can bind to both monomeric and dimeric states without significant structural differences, theoretical and experimental reports have proposed some structural requirements. Thus, in this study, based on known experimental data, the dynamic interaction of VD3 with the monomeric or dimeric forms of βlg was investigated through a protocol combining blind docking and 2 microsecond molecular dynamics simulations coupled with binding free energy and per-residue binding free energy decomposition analyses using the Molecular Mechanics Generalized Born Surface Area approach. Binding free energy calculations allowed us to estimate the energetic differences of coupling VD3 at the calyx and the dimeric interface for the monomeric or dimeric state, revealing that the dimeric structure is required to form a stable complex with VD3 at the dimeric interface. This also has an important impact on the dimerization process, whereas although the monomeric state also forms a stable complex with VD3 at the dimeric interface, the incorporation of the entropy component contributed to producing a marginally favorable binding free energy. Finally, the per-residue decomposition analysis provided energetic information about the most relevant residues in stabilizing the different systems.