Structure and dynamics of β-lactoglobulin in complex with dodecyl sulfate and laurate: A molecular dynamics study

β-Lactoglobulin and Glycodelin: Two Sides of the Same Coin?

The two lipocalins, β-lactoglobulin (βLg) and glycodelin (Gd), are possibly the most closely related members of the large and widely distributed lipocalin family, yet their functions appear to be substantially different. Indeed, the function of β-lactoglobulin, a major component of ruminant milk, is still unclear although neonatal nutrition is clearly important. On the other hand, glycodelin has several specific functions in reproduction conferred through distinct, tissue specific glycosylation of the polypeptide backbone. It is also associated with some cancer outcomes. The glycodelin gene, PAEP, reflecting one of its names, progestagen-associated endometrial protein, is expressed in many though not all primates, but the name has now also been adopted for the β-lactoglobulin gene (HGNC, www.genenames.org). After a general overview of the two proteins in the context of the lipocalin family, this review considers the properties of each in the light of their physiological functional significance, supplementing earlier reviews to include studies from the past decade. While the biological function of glycodelin is reasonably well defined, that of β-lactoglobulin remains elusive.

Energetic and structural effects of the Tanford transition on ligand recognition of bovine β-lactoglobulin

Bovine β-lactoglobulin, an abundant protein in whey, is a promising nanocarrier for peroral administration of drug-like hydrophobic molecules, a process that involves transit through the different acidic conditions of the human digestive tract. Among the several pH-induced conformational rearrangements that this lipocalin undergoes, the Tanford transition is particularly relevant. This transition, which occurs with a midpoint around neutral pH, involves a conformational change of the E-F loop that regulates accessibility to the primary binding site. The effect of this transition on the ligand binding properties of this protein has scarcely been explored. In this study, we carried out an energetic and structural characterization of β-lactoglobulin molecular recognition at pH values above and below the zone in which the Tanford transition occurs. The combined analysis of crystallographic, calorimetric, and molecular dynamics data sheds new light on the interplay between self-association, ligand binding, and the Tanford pre- and post-transition conformational states, revealing novel aspects underlying the molecular recognition mechanism of this enigmatic lipocalin.

Structural mechanism of the Tanford transition of bovine β-lactoglobulin through microsecond molecular dynamics simulations

X-ray analysis has provided structural data about a pH-driven conformational change in β-lactoglobulin (BLG) known as the Tanford transition, which occurs at around pH 7 and involves the EF loop, which acts as a lid closing the internal cavity of the protein below pH 7 and opening it above pH 7. NMR studies using wild-type BLG have encountered problems trying to explain the Tanford transition, however, they have provided important insight using a dimeric BLG mutant, revealing that the opening and closure of the EF loop consists of two types of motions in the microsecond and milliseconds timescales. This provides valuable information indicating that the dimeric state is a good model to study the Tanford transition, although the understanding of this structural change is still lacking at the atomic level. We performed microsecond molecular dynamics (MD) simulations starting from different conformations of BLG in the monomeric and dimeric state, with protonated and deprotonated E89, in order to explore the Tanford transition. Our results provide structural information for the transition from the closed to the open conformation in BLG and show it occurs in the dimeric state in the microsecond timescale, in line with the fast motion observed through NMR experiments. In addition, MD simulations coupled to MMGBSA approach indicated that the most populated conformer of BLG in the open state is able to bind ligands with similar affinity to that of BLG at neutral pH obtained through crystallographic experiments.Communicated by Ramaswamy H. Sarma.

Effects of Polar Head Nature and Tail Length of Single-Chain Lipids on the Conformational Stability of β-Lactoglobulin

Interaction between β-lactoglobulin and single-chain lipids, differing for either the length of the aliphatic chain or the molecular properties of the headgroup, was investigated at neutral and acidic pH to determine the impact on the thermal stability of the protein. Differential scanning calorimetry results with different fatty acids (from C10:0 to C18:0) show a correlation of both melting temperature and unfolding enthalpy of the protein with the ligand binding affinity, and the maximum effect was found for palmitic acid (PLM). The influence of the lipid polar head was investigated by comparing PLM with lyso-palmitoylphosphatidylcholine (LPC), which possesses the same aliphatic chain. At neutral pH, the stabilizing effect of LPC is less favorable compared to PLM. However, fluorescence results revealed that LPC can bind into the protein calyx even at acidic pH, at variance with fatty acids. Molecular dynamics simulations indicated that this difference is due to the ability of the polar head of LPC to interact with the protein loop that regulates the shift (Tanford transition) between open and closed state of the binding site of β-lactoglobulin. The results provide a rationale for how a ligand has the ability to access the protein active site at acidic conditions by overcoming the Tanford transition, and they demonstrate that β-lactoglobulin can deliver ligands with tailored properties of the polar head in a wide pH range.

In Silico Characterization of the Binding Modes of Surfactants with Bovine Serum Albumin

The binding interactions of the surfactants: anionic sodium dodecyl sulphate (SDS), cationic cetyltrimethylammonium bromide (CTAB), non-ionic octyl glucoside (OG), and zwitterionic 3-[Hexadecyl(dimethyl)ammonio]-1-propanesulfonate (HPS), with bovine serum albumin (BSA) were investigated by computer simulation. The results disclosed that the surfactants bound stably between hydrophobic subdomain IIA and IIIA where tryptophan-213 residue, an important intrinsic fluorophore in BSA is housed. The interactions of the surfactants with the BSA were electrostatic and hydrophobic interactions. The head-groups of SDS, HPS and OG formed hydrogen bonds with the BSA, while that of CTAB was shielded from intermolecular hydrogen-bonding due to intervening methyl groups. Subsequently, molecular dynamics (MD) simulation of the protein-surfactant complexes revealed that hydrogen bonds formed by OG were stronger than those of SDS and HPS. However, the decomposed force-field energies showed that OG had the least interaction energy with the BSA. In addition to MD simulation, it was found by density functional theory (DFT) that the differences in the coulomb interaction energies can be attributed to charge distribution in the surfactants. Overall, free energies calculated by linear interaction energy (LIE) proved that the binding of each surfactant was dominated by differences between van der Waals interactions in bound and free states.

β-Lactoglobulin as a potential carrier for bioactive molecules

In this study, the interaction and binding behavior of anesthetic tetracaine (TET) with bovine β-lactoglobulin (LGB) isoform A and a mixture of isoforms A and B were investigated under varying environmental conditions (pH, ionic strength, concentration, LGB-TET complex molar ratio). A wide range of analytical techniques (dynamic light scattering (DLS), electrophoretic mobility, UV-Vis spectroscopy, circular dichroism (CD), quartz crystal microbalance (QCM-D) were used to analyze the physicochemical properties of the complexes in bulk solution and on the surface of gold. The experiments revealed that TET, which is amphiphilic, could bind with LGB not only in the β-barrel but also onto the surface. The zeta potential of the LGB becomes more positively charged upon interaction with TET due to electrostatic interaction of the amino group present in the TET structure. Changes in the zeta potential values are mostly visible above pH 6 for all tested systems. CD spectra show that interaction with the ligand does not change the secondary structure of the protein. The physicochemical properties of LGB-TET complex were examined in an adsorbed state on a gold surface using the QCM-D method. Additionally, molecular docking was used to evaluate the most likely binding site for TET with LGB.

Molecular Dynamics Simulation for Mechanism Elucidation of Food Processing and Safety: State of the Art

Molecular dynamics (MD) simulation is a useful technique to study the interaction between molecules and how they are affected by various processes and processing conditions. This review summarizes the application of MD simulations in food processing and safety, with an emphasis on the effects that emerging nonthermal technologies (for example, high hydrostatic pressure, pulsed electric field) have on the molecular and structural characteristics of foods and biomaterials. The advances and potential projection of MD simulations in the science and engineering aspects of food materials are discussed and focused on research work conducted to study the effects of emerging technologies on food components. It is expected by showing key case studies that it will stir novel developments as a valuable tool to study the effects of emerging food technologies on biomaterials. This review is useful to food researchers and the food industry, as well as researchers and practitioners working on flavor and nutraceutical encapsulations, dietary carbohydrate product developments, modified starches, protein engineering, and other novel food applications.

β-Lactoglobulin: An efficient nanocarrier for advanced delivery systems

Thanks to the progress of nanotechnology there are several agent-delivery systems that can be selected to achieve rapid and specific delivery of a wide variety of biologically active agents. Consequently, the manipulation and engineering of biopolymers has become one of the most exciting subjects for those who study delivery systems on the nanoscale. In this regard, both nanoparticle formation and a carrier role have been observed in the case of the globular milk whey protein, β-lactoglobulin (β-LG), setting it apart from many other proteins. To date, many efforts adopting different approaches have created β-LG nanoparticles useful in forming delivery systems for various agents with specific targets. In this review, the potential of β-LG to play the role of an efficient and diverse carrier protein, as well as its ability to form a well-targeted nano-scale delivery system is discussed.

Exploring the conformational and binding properties of unphosphorylated/phosphorylated monomeric and trimeric Bcl-2 through docking and molecular dynamics simulations

B-cell lymphoma (Bcl-2) is commonly associated with the progression and preservation of cancer and certain lymphomas; therefore, it is considered as a biological target against cancer. Nevertheless, evidence of all its structural binding sites has been hidden because of the lack of a complete Bcl-2 model, given the presence of a flexible loop domain (FLD), which is responsible for its complex behavior. FLD region has been implicated in phosphorylation, homotrimerization, and heterodimerization associated with Bcl-2 antiapoptotic function. In this contribution, homology modeling, molecular dynamics (MD) simulations in the microsecond (µs) time-scale and docking calculations were combined to explore the conformational complexity of unphosphorylated/phosphorylated monomeric and trimeric Bcl-2 systems. Conformational ensembles generated through MD simulations allowed for identifying the most populated unphosphorylated/phosphorylated monomeric conformations, which were used as starting models to obtain trimeric complexes through protein-protein docking calculations, also submitted to µs MD simulations. Principal component analysis showed that FLD represents the main contributor to total Bcl-2 mobility, and is affected by phosphorylation and oligomerization. Subsequently, based on the most representative unphosphorylated/phosphorylated monomeric and trimeric Bcl-2 conformations, docking studies were initiated to identify the ligand binding site of several known Bcl-2 inhibitors to explain their influence in homo-complex formation and phosphorylation. Docking studies showed that the different conformational states experienced by FLD, such as phosphorylation and oligomerization, play an essential role in the ability to make homo and hetero-complexes. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 393-413, 2016.

Searching the conformational complexity and binding properties of HDAC6 through docking and molecular dynamic simulations

Histone deacetylases (HDACs) are a family of proteins involved in the deacetylation of histones and other non-histones substrates. HDAC6 belongs to class II and shares similar biological functions with others of its class. Nevertheless, its three-dimensional structure that involves the catalytic site remains unknown for exploring the ligand recognition properties. Therefore, in this contribution, homology modeling, 100-ns-long Molecular Dynamics (MD) simulation and docking calculations were combined to explore the conformational complexity and binding properties of the catalytic domain 2 from HDAC6 (DD2-HDAC6), for which activity and affinity toward five different ligands have been reported. Clustering analysis allowed identifying the most populated conformers present during the MD simulation, which were used as starting models to perform docking calculations with five DD2-HDAC6 inhibitors: Cay10603 (CAY), Rocilinostat (RCT), Tubastatin A (TBA), Tubacin (TBC), and Nexturastat (NXT), and then were also submitted to 100-ns-long MD simulations. Docking calculations revealed that the five inhibitors bind at the DD2-HDAC6 binding site with the lowest binding free energy, the same binding mode is maintained along the 100-ns-long MD simulations. Overall, our results provide structural information about the molecular flexibility of apo and holo DD2-HDAC6 states as well as insight of the map of interactions between DD2-HDAC6 and five well-known DD2-HDAC6 inhibitors allowing structural details to guide the drug design. Finally, we highlight the importance of combining different theoretical approaches to provide suitable structural models for structure-based drug design.

Energetic and conformational features linked to the monomeric and dimeric states of bovine BLG

Bovine β-lactoglobulin (BLG) belong to the lipocalin family. This is a group of proteins involved in the binding and transporting of hydrophobic molecules. Experimental and theoretical reports have stated its complex structural behavior in solution, with coupled effects between homodimerization and ligand recognition. Nonetheless, structural evidence at the atomic level about the cause of this coupled effect has not been reported to date. To address this issue microsecond molecular dynamics (MD) simulations were combined with the molecular mechanics generalized Born surface area (MM/GBSA) approach, clustering analysis and principal component analysis (PCA), to explore the conformational complexity of BLG protein-protein self-association and palmitic acid (PLM) or dodecyl sulfate (SDS) ligand recognition in the monomeric and dimeric state. MD simulations, coupled to the MM/GBSA method, revealed that dimerization exerts contrasting effects on the ligand-binding capacity of BLG. Protein dimerization decreases PLM affinity, promoting dimer association. For SDS the dimeric state increases affinity, enhancing dimer dissociation. MD simulations based on PCA revealed that while few differences in the conformational subspace are observed between the free and bound monomer and dimer coupling for PLM, substantial changes are observed between the free and bound monomer and dimer coupling for SDS.

Understanding the molecular basis of agonist/antagonist mechanism of GPER1/GPR30 through structural and energetic analyses

The G-protein coupled receptors (GPCRs) represent the largest superfamily of membrane proteins in charge to pass the cell signaling after binding with their cognate ligands to the cell interior. In breast cancer, a GPCR named GPER1 plays a key role in the process of growth and the proliferation of cancer cells. In a previous study, theoretical methods were applied to construct a model of GPER1, which later was submitted to molecular dynamics (MD) simulations to perform a docking calculation. Based on this preceding work, it is known that GPER1 is sensitive to structural differences in its binding site. However, due to the nature of that past study, conformational changes linked to the ligand binding were not observed. Therefore, in this study, in order to explore the conformational changes coupled to the agonist/antagonist binding, MD simulations of about 0.25μs were performed for the free and bound states, summarizing 0.75μs of MD simulation in total. For the bound states, one agonist (G-1) and antagonist (G-15) were chosen since is widely known that these two molecules cause an impact on GPER1 mobility. Based on the conformational ensemble generated through MD simulations, we found that despite G-1 and G-15 being stabilized by similar map of residues, the structural differences between both ligands impact the hydrogen bond pattern not only at the GPER1 binding site but also along the seven-helix bundle, causing significant differences in the conformational mobility along the extracellular and cytoplasmic domain, and to a lesser degree in the curvatures of helix 2, helix 3 and helix 7 between the free and bound states, which is in agreement with reported literature, and might be linked to microscopic characteristics of the activated-inactivated transition. Furthermore, binding free energy calculations using the MM/GBSA method for the bound states, followed by an alanine scanning analysis allowed us to identify some important residues for the complex stabilization.

Energetic and flexibility properties captured by long molecular dynamics simulations of a membrane-embedded pMHCII-TCR complex

Although crystallographic data have provided important molecular insight into the interactions in the pMHC-TCR complex, the inherent features of this structural approach cause it to only provide a static picture of the interactions. While unbiased molecular dynamics simulations (UMDSs) have provided important information about the dynamic structural behavior of the pMHC-TCR complex, most of them have modeled the pMHC-TCR complex as soluble, when in physiological conditions, this complex is membrane bound; therefore, following this latter UMDS protocol might hamper important dynamic results. In this contribution, we performed three independent 300 ns-long UMDSs of the pMHCII-TCR complex anchored in two opposing membranes to explore the structural and energetic properties of the recognition of pMHCII by the TCR. The conformational ensemble generated through UMDSs was subjected to clustering and Cartesian principal component analyses (cPCA) to explore the dynamical behavior of the pMHCII-TCR association. Furthermore, based on the conformational population sampled through UMDSs, the effective binding free energy, per-residue free energy decomposition, and alanine scanning mutations were explored for the native pMHCII-TCR complex, as well as for 12 mutations (p1-p12MHCII-TCR) introduced in the native peptide. Clustering analyses and cPCA provide insight into the rocking motion of the TCR onto pMHCII, together with the presence of new electrostatic interactions not observed through crystallographic methods. Energetic results provide evidence of the main contributors to the pMHC-TCR complex formation as well as the key residues involved in this molecular recognition process.

Factors affecting the interactions between beta-lactoglobulin and fatty acids as revealed in molecular dynamics simulations

Beta-lactoglobulin (BLG), a bovine dairy protein, is a promiscuously interacting protein that can bind multiple hydrophobic ligands. Fatty acids (FAs), common hydrophobic molecules bound to BLG, are important sources of fuel for life because they yield large quantities of ATP when metabolized. The binding affinity increases with the length of the ligands, indicating the importance of the van der Waals (vdW) interactions between the hydrocarbon tail and the hydrophobic calyx of BLG. An exception to this rule is caprylic acid (OCA) which is two-carbon shorter but has a stronger binding affinity than capric acid. Theoretical calculations in the current literature are not accurate enough to shed light on the underlying physics of this exception. The computed affinity values are greater for longer fatty acids without respect for the caprylic exception and those values are generally several orders of magnitude away from the experimental data. In this work, we used hybrid steered molecular dynamics to accurately compute the binding free energies between BLG and the five saturated FAs of 8 to 16 carbon atoms. The computed binding free energies agree well with experimental data not only in rank but also in absolute values. We gained insights into the exceptional behavior of caprylic acid in the computed values of entropy and electrostatic interactions. We found that the electrostatic interaction between the carboxyl group of caprylic acid and the two amino groups of K60/69 in BLG is much stronger than the vdW force between the OCA's hydrophobic tail and the BLG calyx. This pulls OCA to the top of the beta barrel where it is easier to fluctuate, giving rise to greater entropy of OCA at the binding site.

The differences in binding 12-carbon aliphatic ligands by bovine β-lactoglobulin isoform A and B studied by isothermal titration calorimetry and X-ray crystallography

Isoforms A (LGB-A) and B (LGB-B) of bovine lactoglobulin, the milk protein, differ in positions 64 (D↔G) and 118 (V↔A). Interactions of LGB-A and LGB-B with sodium dodecyl sulfate (SDS), dodecyltrimethylammonium chloride (DTAC) and lauric acid (LA), 12-carbon ligands possessing differently charged polar groups, were investigated using isothermal titration calorimetry and X-ray crystallography, to study the proton linkage phenomenon and to distinguish between effects related to different isoforms and different ligand properties. The determined values of ΔS and ΔH revealed that for all ligands, binding is entropically driven. The contribution from enthalpy change is lower and shows strong dependence on type of buffer that indicates proton release from the protein varying with protein isoform and ligand type and involvement of LA and Asp64 (in isoform A) in this process. The ligand affinities for both isoforms were arranged in the same order, DTAC < LA < SDS, and were systematically lower for variant B. The entropy change of the complexation process was always higher for isoform A, but these values were compensated by changes in enthalpy, resulting in almost identical ΔG for complexes of both isoforms. The determined crystal structures showed that substitution in positions 64 and 118 did not influence the overall structure of LGB complexes. The chemical character of the ligand polar group did not affect the position of its aliphatic chain in protein β-barrel, indicating a major role of hydrophobic interactions in ligand binding that prevailed even with the repulsion between positively charged DTAC and lysine residues located at binding site entrance.

Spectroscopic and theoretical investigation of oxali-palladium interactions with β-lactoglobulin

The possibility of using a small cheap dairy protein, β-lactoglobulin (β-LG), as a carrier for oxali-palladium for drug delivery was studied. Their binding in an aqueous solution at two temperatures of 25 and 37°C was investigated using spectroscopic techniques in combination with a molecular docking study. Fluorescence intensity changes showed combined static and dynamic quenching during β-LG oxali-palladium binding, with the static mode being predominant in the quenching mechanism. The binding and thermodynamic parameters were determined by analyzing the results of quenching and those of the van't Hoff equation. According to obtained results the binding constants at two temperatures of 25 and 37°C are 3.3×10(9) M(-1) and 18.4×10(6) M(-1) respectively. Fluorescence resonance energy transfer (FRET) showed that the experimental results and the molecular docking results were coherent. An absence change of β-LG secondary structure was confirmed by the CD results. Molecular docking results agreed fully with the experimental results since the fluorescence studies also revealed the presence of two binding sites with a negative value for the Gibbs free energy of binding of oxali-palladium to β-LG. Furthermore, molecular docking and experimental results suggest that the hydrophobic effect plays a critical role in the formation of the oxali-palladium complex with β-LG. This agreement between molecular docking and experimental results implies that docking studies may be a suitable method for predicting and confirming experimental results, as shown in this study. Hence, the combination of molecular docking and spectroscopy methods is an effective innovative approach for binding studies, particularly for pharmacophores.

Sequence-based Gaussian network model for protein dynamics

Motivation: Gaussian network model (GNM) is widely adopted to analyze and understand protein dynamics, function and conformational changes. The existing GNM-based approaches require atomic coordinates of the corresponding protein and cannot be used when only the sequence is known.

Results: We report, first of its kind, GNM model that allows modeling using the sequence. Our linear regression-based, parameter-free, sequence-derived GNM (L-pfSeqGNM) uses contact maps predicted from the sequence and models local, in the sequence, contact neighborhoods with the linear regression. Empirical benchmarking shows relatively high correlations between the native and the predicted with L-pfSeqGNM B-factors and between the cross-correlations of residue fluctuations derived from the structure- and the sequence-based GNM models. Our results demonstrate that L-pfSeqGNM is an attractive platform to explore protein dynamics. In contrast to the highly used GNMs that require protein structures that number in thousands, our model can be used to study motions for the millions of the readily available sequences, which finds applications in modeling conformational changes, protein–protein interactions and protein functions.

Contact:zerozhua@126.com

Supplementary information:Supplementary data are available at Bioinformatics online.

Ligand binding and self-association cooperativity of β-lactoglobulin

Unlike most small globular proteins, lipocalins lack a compact hydrophobic core. Instead, they present a large central cavity that functions as the primary binding site for hydrophobic molecules. Not surprisingly, these proteins typically exhibit complex structural dynamics in solution, which is intricately modified by intermolecular recognition events. Although many lipocalins are monomeric, an increasing number of them have been proven to form oligomers. The coupling effects between self-association and ligand binding in these proteins are largely unknown. To address this issue, we have calorimetrically characterized the recognition of dodecyl sulfate by bovine β-lactoglobulin, which forms weak homodimers at neutral pH. A thermodynamic analysis based on coupled-equilibria revealed that dimerization exerts disparate effects on the ligand-binding capacity of β-lactoglobulin. Protein dimerization decreases ligand affinity (or, reciprocally, ligand binding promotes dimer dissociation). The two subunits in the dimer exhibit a positive, entropically driven cooperativity. To investigate the structural determinants of the interaction, the crystal structure of β-lactoglobulin bound to dodecyl sulfate was solved at 1.64 Å resolution.

β-lactoglobulin's conformational requirements for ligand binding at the calyx and the dimer interphase: a flexible docking study

β-lactoglobulin (BLG) is an abundant milk protein relevant for industry and biotechnology, due significantly to its ability to bind a wide range of polar and apolar ligands. While hydrophobic ligand sites are known, sites for hydrophilic ligands such as the prevalent milk sugar, lactose, remain undetermined. Through the use of molecular docking we first, analyzed the known fatty acid binding sites in order to dissect their atomistic determinants and second, predicted the interaction sites for lactose with monomeric and dimeric BLG. We validated our approach against BLG structures co-crystallized with ligands and report a computational setup with a reduced number of flexible residues that is able to reproduce experimental results with high precision. Blind dockings with and without flexible side chains on BLG showed that: i) 13 experimentally-determined ligands fit the calyx requiring minimal movement of up to 7 residues out of the 23 that constitute this binding site. ii) Lactose does not bind the calyx despite conformational flexibility, but binds the dimer interface and an alternate Site C. iii) Results point to a probable lactolation site in the BLG dimer interface, at K141, consistent with previous biochemical findings. In contrast, no accessible lysines are found near Site C. iv) lactose forms hydrogen bonds with residues from both monomers stabilizing the dimer through a claw-like structure. Overall, these results improve our understanding of BLG's binding sites, importantly narrowing down the calyx residues that control ligand binding. Moreover, our results emphasize the importance of the dimer interface as an insufficiently explored, biologically relevant binding site of particular importance for hydrophilic ligands. Furthermore our analyses suggest that BLG is a robust scaffold for multiple ligand-binding, suitable for protein design, and advance our molecular understanding of its ligand sites to a point that allows manipulation to control binding.

Molecular dynamics of a thermostable multicopper oxidase from Thermus thermophilus HB27: structural differences between the apo and holo forms

Molecular dynamic (MD) simulations have been performed on Tth-MCO, a hyperthermophilic multicopper oxidase from thermus thermophilus HB27, in the apo as well as the holo form, with the aim of exploring the structural dynamic properties common to the two conformational states. According to structural comparison between this enzyme and other MCOs, the substrate in process to electron transfer in an outer-sphere event seems to transiently occupy a shallow and overall hydrophobic cavity near the Cu type 1 (T1Cu). The linker connecting the β-strands 21 and 24 of the second domain (loop (β21–β24)_D2) has the same conformation in both states, forming a flexible lid at the entrance of the electron-transfer cavity. Loop (β21–β24)_D2 has been tentatively assigned a role occluding the access to the electron-transfer site. The dynamic of the loop (β21–β24)_D2 has been investigated by MD simulation, and results show that the structures of both species have the same secondary and tertiary structure during almost all the MD simulations. In the simulation, loop (β21–β24)_D2 of the holo form undergoes a higher mobility than in the apo form. In fact, loop (β21–β24)_D2 of the holo form experiences a conformational change which enables exposure to the electron-transfer site (open conformation), while in the apo form the opposite effect takes place (closed conformation). To confirm the hypothesis that the open conformation might facilitate the transient electron-donor molecule occupation of the site, the simulation was extended another 40 ns with the electron-donor molecule docked into the protein cavity. Upon electron-donor molecule stabilization, loops near the cavity reduce their mobility. These findings show that coordination between the copper and the protein might play an important role in the general mobility of the enzyme, and that the open conformation seems to be required for the electron transfer process to T1Cu.