Unusual conformational changes in 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase as revealed by X-ray crystallography and NMR

The crystal structure of Escherichia coli 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) in complex with MgADP has been determined at 1.5-A resolution with a crystallographic R factor of 0.191. The solution structure of HPPK in complex with Mg(2+) and beta,gamma-methyleneadenosine 5'-triphosphate (MgAMPPCP) has been determined using a simulated annealing protocol with 3,523 experimental NMR restraints. The root mean square deviation of the ensemble of 20 refined conformers that represent the solution structure from the mean coordinate set derived from them is 0.74 +/- 0.26 A for all backbone atoms and 0.49 +/- 0.22 A when residues Pro(14), Pro(44)-Gln(50), and Arg(84)-Pro(91) are excluded. Binding of MgADP causes significant changes in the conformation and dynamical property of three loops of HPPK that are involved in catalysis. A dramatic, unusual conformational change is that loop 3 moves away from the active center significantly with some residues moving by >17 A. The binding of MgADP also stabilizes loop 1 and loop 3 but makes loop 2 more mobile. Very similar conformational and dynamical changes are observed in the NMR solution structure of HPPK.MgAMPPCP. The conformational and dynamical changes may play important roles in both substrate binding and product release in the catalytic cycle.     

Mechanism of the conformational transitions in 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase as revealed by NMR spectroscopy

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the transfer of pyrophosphate from ATP to 6-hydroxymethyl-7,8-dihydropterin (HP), leading to the biosynthesis of folate cofactors. HPPK undergoes dramatic conformational changes during its catalytic cycle, and the conformational changes are essential for enzymatic catalysis. Thus, the enzyme is not only an attractive target for developing antimicrobial agents but also an excellent model system for studying the catalytic mechanism of enzymatic pyrophosphoryl transfer as well as the role of protein dynamics in enzymatic catalysis. In the present study, we report the NMR solution structures of the binary complex HPPK*MgAMPCPP and the ternary complex HPPK*MgAMPCPP*DMHP, where alpha,beta-methyleneadenosine triphosphate (AMPCPP) and 7, 7-dimethyl-6-hydroxypterin (DMHP) are the analogues of the substrates ATP and HP, respectively. The results suggest that the three catalytic loops of the binary complex of HPPK can assume multiple conformations in slow exchanges as evidenced by multiple sets of NMR signals for several residues in loops 2 and 3 and the very weak or missing NH cross-peaks for several residues in loops 1 and 3. However, the ternary complex shows only one set of NMR signals, and the cross-peak intensities are rather uniform, suggesting that the binding of the second substrate shifts the multiple conformations of the binary complex to an apparently single conformation of the ternary complex. The NMR behaviors and conformations of the binary complex HPPK*MgAMPCPP are significantly different from those of HPPK in complex with Mgbeta,gamma-methyleneadenosine triphosphate (MgAMPPCP). It is suggested that the conformational properties of the binary substrate complex HPPK*MgATP be represented by those of HPPK*MgAMPCPP, because MgAMPCPP is a better MgATP analogue for HPPK with respect to both binding affinity and bound conformation.     

Triggering loops and enzyme function: identification of loops that trigger and modulate movements

Enzyme function often involves a conformational change. There is a general agreement that loops play a vital role in correctly positioning the catalytically important residues. Nevertheless, predicting the functional loops and most importantly their role in enzyme function remains a difficult task. A major reason for this difficulty is that loops that undergo conformational change are frequently not well conserved in their primary sequence. beta1,4-Galactosyltransferase is one such enzyme. There, the amino acid sequence of a long loop that undergoes a large conformational change upon substrate binding is not well conserved. Our molecular dynamics simulations show that the large conformational change in the long loop is brought about by a second, interacting loop. Interestingly, while the structural change of the second loop is much smaller than that of the long loop, its sequence (particularly glycine residues) is highly conserved. We further examine the generality of the proposition that there are loops that trigger movements but nevertheless show little or no structural changes in crystals. We focus on two other enzymes, enolase and lipase. We chose these enzymes, since they too undergo conformational change upon ligand binding, however, they have different folds and different functions. Through multiple sets of simulations we show that the conformational change of the functional loop(s) is brought about through communication of flexibility by triggering loops that have several glycine residues. We further propose that similar to the conservation of common favorable fold types and structural motifs, evolution has also conserved common "skillful" mechanisms. Mechanisms may be conserved across different folds, sequences and functions, with adaptation to specific enzymatic roles.     

Is the critical role of loop 3 of Escherichia coli 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase in catalysis due to loop-3 residues arginine-84 and tryptophan-89? Site-directed mutagenesis, biochemical, and crystallographic studies

Deletion mutagenesis, biochemical, and X-ray crystallographic studies have shown that loop 3 of Escherichia coli 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is required for the assembly of the active center, plays an important role in the stabilization of the ternary complex of HPPK with MgATP and 6-hydroxymethyl-7,8-dihydropterin (HP), and is essential for catalysis. Whether the critical functional importance of loop 3 is due to the interactions between residues R84 and W89 and the two substrates has been addressed by site-directed mutagenesis, biochemical, and X-ray crystallographic studies. Substitution of R84 with alanine causes little changes in the dissociation constants and kinetic constants of the HPPK-catalyzed reaction, indicating that R84 is not important for either substrate binding or catalysis. Substitution of W89 with alanine increases the K(d) for the binding of MgATP by a factor of 3, whereas the K(d) for HP increases by a factor of 6, which is due to the increase in the dissociation rate constant. The W89A mutation decreases the rate constant for the chemical step of the forward reaction by a factor of 15 and the rate constant for the chemical step of the reverse reaction by a factor of 25. The biochemical results of the W89A mutation indicate that W89 contributes somewhat to the binding of HP and more significantly to the chemical step. The crystal structures of W89A show that W89A has different conformations in loops 2 and 3, but the critical catalytic residues are positioned for catalysis. When these results are taken together, they suggest that the critical functional importance of loop 3 is not due to the interactions of the R84 guanidinium group or the W89 indole ring with the substrates.     

Probing the binding mechanism of mercaptoguanine derivatives as inhibitors of HPPK by docking and molecular dynamics simulations

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a promising antimicrobial target involved in the folate biosynthesis pathway. Although, the results from crystallographic studies of HPPK have attracted a great interest in the design of novel HPPK inhibitors, the mechanism of action of HPPK due to inhibitor binding remains questionable. Recently, mercaptoguanine derivatives were reported to inhibit the pyrophosphoryl transfer mechanism of Staphylococcus aureus HPPK (SaHPPK). The present study is an attempt to understand the SaHPPK-inhibitors binding mechanism and to highlight the key residues that possibly involve in the complex formation. To decipher these questions, we used the state-of-the-art advanced insilico approach such as molecular docking, molecular dynamics (MD), molecular mechanics-generalized Born surface area approach. Domain cross correlation and principle component analysis were applied to the snapshots obtained from MD revealed that the compounds with high binding affinity stabilize the conformational dynamics of SaHPPK. The binding free energy estimation showed that the van der Waals and electrostatic interactions played a vital role for the binding mechanism. Additionally, the predicted binding free energy was in good agreement with the experimental values (R² = .78). Moreover, the free energy decomposition on per-residue confirms the key residues that significantly contribute to the complex formation. These results are expected to be useful for rational design of novel SaHPPK inhibitors.     

Disparate degrees of hypervariable loop flexibility control T-cell receptor cross-reactivity, specificity, and binding mechanism

αβ T-cell receptors (TCRs) recognize multiple antigenic peptides bound and presented by major histocompatibility complex molecules. TCR cross-reactivity has been attributed in part to the flexibility of TCR complementarity-determining region (CDR) loops, yet there have been limited direct studies of loop dynamics to determine the extent of its role. Here we studied the flexibility of the binding loops of the αβ TCR A6 using crystallographic, spectroscopic, and computational methods. A significant role for flexibility in binding and cross-reactivity was indicated only for the CDR3α and CDR3β hypervariable loops. Examination of the energy landscapes of these two loops indicated that CDR3β possesses a broad, smooth energy landscape, leading to rapid sampling in the free TCR of a range of conformations compatible with different ligands. The landscape for CDR3α is more rugged, resulting in more limited conformational sampling that leads to specificity for a reduced set of peptides as well as the major histocompatibility complex protein. In addition to informing on the mechanisms of cross-reactivity and specificity, the energy landscapes of the two loops indicate a complex mechanism for TCR binding, incorporating elements of both conformational selection and induced fit in a manner that blends features of popular models for TCR recognition.     

Extended molecular dynamics of a c-kit promoter quadruplex

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.     

Molecular motions and conformational changes of HPPK

6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) belongs to a class of catalytic enzymes involved in phosphoryl transfer and is a new target for the development of novel antimicrobial agents. In the present study, the fundamental consideration is to view the overall structure of HPPK as a network of interacting residues and to extract the most cooperative collective motions that define its global dynamics. A coarse-grained model, harmonically constrained according to HPPK's crystal structure is used. Four crystal structures of HPPK (one apo and three holo forms with different nucleotide and pterin analogs) are studied with the goal of providing insights about the function-dynamic correlation and ligand induced conformational changes. The dynamic differences are examined between HPPK's apo- and holo-forms, because they are involved in the catalytic reaction steps. Our results indicate that the palm-like structure of HPPK is nearly rigid, whereas the two flexible loops: L2 (residues 43-53) and L3 (residues 82-92) exhibit the most concerted motions for ligand recognition and presumably, catalysis. These two flexible loops are involved in the recognition of HPPKs nucleotide and pterin ligands, whereas the rigid palm region is associated with binding of these cognate ligands. Six domains of collective motions are identified, comprised of structurally close but not necessarily sequential residues. Two of these domains correspond to the flexible loops (L2 and L3), whereas the remaining domains correspond to the rigid part of the molecule.     

Interdependence of backbone flexibility,residue conservation,and enzyme function: a case study on beta1,4-galactosyltransferase-I

Beta1,4-galactosyltransferase-I (beta4Gal-T1) catalyzes the transfer of a galactose from UDP-galactose to N-acetylglucosamine. A recent crystal structure determination of the substrate-bound enzyme reveals a large conformational change, which creates binding sites for the oligosaccharide and alpha-lactalbumin, when compared to the ligand-free structure. The conformational changes take place in a 21-residue-long loop (I345-H365) and in a smaller loop containing a tryptophan residue (W314) flanked by glycines (Y311-G316; Trp loop). A series of molecular dynamics simulations carried out with an implicit solvent model and with explicit water successfully identify flexibility in the two loops and in another interacting loop. These observations are confirmed by limited proteolysis experiments that reveal an intrinsic flexibility of the long loop. The multiple simulation runs starting with the substrate-free structure show that the long loop moves toward its conformation in the ligand-bound structure; however, it gets stabilized in an intermediate position. The Trp loop moves in the opposite direction to that of the long loop, making contacts with residues in the long loop. Remarkably, when the Trp loop is restrained in its starting conformation, no large conformational change takes place in the long loop, indicating residue communication of flexibility. Sequence and structural analysis of the beta4Gal-T1 family with 37 known sequences reveals that in contrast to the unconserved long loop, which undergoes a much larger conformational change, the Trp loop including the glycines is highly conserved. These observations lead us to propose a new functional mechanism that may be conserved by evolution to perform a variety of functions.     

Molecular dynamics, crystallography and mutagenesis studies on the substrate gating mechanism of prolyl oligopeptidase

Altered prolyl oligopeptidase (PREP) activity is found in many common neurological and other genetic disorders, and in some cases PREP inhibition may be a promising treatment. The active site of PREP resides in an internal cavity; in addition to the direct interaction between active site and substrate or inhibitor, the pathway to reach the active site (the gating mechanism) must be understood for more rational inhibitor design and understanding PREP function. The gating mechanism of PREP has been investigated through molecular dynamics (MD) simulation combined with crystallographic and mutagenesis studies. The MD results indicate the inter-domain loop structure, comprised of 3 loops at residues, 189-209 (loop A), 577-608 (loop B), and 636-646 (loop C) (porcine PREP numbering), are important components of the gating mechanism. The results from enzyme kinetics of PREP variants also support this hypothesis: When loop A is (1) locked to loop B through a disulphide bridge, all enzyme activity is halted, (2) nicked, enzyme activity is increased, and (3) removed, enzyme activity is only reduced. Limited proteolysis study also supports the hypothesis of a loop A driven gating mechanism. The MD results show a stable network of H-bonds that hold the two protein domains together. Crystallographic study indicates that a set of known PREP inhibitors inhabit a common binding conformation, and this H-bond network is not significantly altered. Thus the domain separation, seen to occur in lower taxa, is not involved in the gating mechanism for mammalian PREP. In two of the MD simulations we observed a conformational change that involved the breaking of the H-bond network holding loops A and B together. We also found that this network was more stable when the active site was occupied, thus decreasing the likelihood of this transition.     

Flexibility and structural conservation in a c-KIT G-quadruplex

A quadruplex sequence from the promoter region of the c-KIT gene forms a stable quadruplex, as characterized by crystallographic and NMR methods. Two new crystal structures are reported here, together with molecular dynamics simulation studies on these quadruplex crystal structures and an NMR structure. The new crystal structures, each in a distinct space group and lattice packing arrangement, together with the existing structures, demonstrate that the c-KIT quadruplex fold does not change with differing environments, suggesting that quadruplex topological dynamism is not a general phenomenon. The single and dinucleotide loops in these structures show a high degree of conformational flexibility within the three crystal forms and the NMR ensemble, with no evidence of clustering to particular conformers. This is in accord with the findings of high loop flexibility from the molecular dynamics studies. It is suggested that intramolecular quadruplexes can be grouped into two broad classes (i) those with at least one single-nucleotide loop, often showing singular topologies even though loops are highly flexible, and (ii) with all loops comprising at least two nucleotides, leading to topological dynamism. The loops can have more stable and less dynamic base-stacked secondary structures.     

Essential roles of a dynamic loop in the catalysis of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the transfer of pyrophosphoryl group from ATP to 6-hydroxymethyl-7,8-dihydropterin (HP) following an ordered bi-bi mechanism with ATP as the first substrate. The rate-limiting step of the reaction is product release, and the complete active center is assembled and sealed only upon the binding of both ATP and HP. The assembly of the active center involves large conformational changes in three catalytic loops, among which loop 3 undergoes the most dramatic and unusual changes. To investigate the roles of loop 3 in catalysis, we have made a deletion mutant, which has been investigated by biochemical and X-ray crystallographic analysis. The biochemical data showed that the deletion mutation does not have significant effects on the dissociation constants or the rate constants for the binding of the first substrate MgATP or its analogues. The dissociation constant of HP for the mutant increases by a factor of approximately 100, which is due to a large increase in the dissociation rate constant. The deletion mutation causes a shift of the rate-limiting step in the reaction and a decrease in the rate constant for the chemical step by a factor of approximately 1.1 x 10(5). The crystal structures revealed that the deletion mutation does not affect protein folding, but the catalytic center of the mutant is not fully assembled even upon the formation of the ternary complex and is not properly sealed. The results together suggest that loop 3 is dispensable for the folding of the protein and the binding of the first substrate MgATP, but is required for the assembling and sealing of the active center. The loop plays an important role in the stabilization of the ternary complex and is critical for catalysis.     

Dynamics and stability of GCAA tetraloops with 2-aminopurine and purine substitutions

The contributions of various interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using the substitutions of functional groups. The guanosine (G) in the first tetraloop position or in the C-G closing base pair was replaced by 2-aminopurine (AP), and the individual tetraloop's adenosines (A) were replaced by purine (PUR). These substitutions eliminated particular hydrogen bonds thought to stabilize the GCAA tetraloop. For each substitution, molecular dynamics (MD) simulations were carried out in an aqueous solution with sodium counterions, using the CHARMM27 force field. The MD simulations showed that the substitutions in the first (G-->AP) and the third (A-->PUR) position of the GCAA tetraloop did not significantly influence the conformation of the hairpin. A long-lived bridging water molecule observed in the GCAA loop was present in both modified loops. The substitutions made in the last loop position (A-->PUR) or in the C-G base pair closing the tetraloop (G-->AP) to some extent influenced the loop structure and dynamics. These loops did not display the long-lived bridging water molecules. When the second A in the GCAA loop was replaced by PUR, the first A in the loop was observed in the anti or in the syn orientation about the glycosyl bond. The G to AP substitution in C-G base pair led to a change of their arrangement from the Watson-Crick to wobble. The MD simulations of the hairpin with C-AP wobble closing base pair showed increased conformational dynamics of the hairpin. The changes of hairpin formation free energy associated with the substitutions of individual bases were calculated by the free energy perturbation method. Our theoretical estimates suggest a larger destabilization for the G to AP substitutions in GCAA loop than for the substitutions of individual A's by PUR, which is in accordance with experimental tendency. The calculations predicted a similar free energy change for G to AP substitutions in the GCAA tetraloop and in the C-G closing base pair.     

Computer aided study of ligand binding with catalytic domain of Avian sarcoma virus integrase and its ligand binding loops.

We have carried out 1 nanosecond (ns) Molecular Dynamics (MD) simulations of the drug Y3 (4-acetylamino-5-hydroxynaphthalene- 2, 7-disulfonic acid) complexed with catalytic domain of Avian sarcoma virus Integrase (ASV-IN), both in vacuum and in the presence of explicit solvent. Starting models were obtained on the basis of PDB co- ordinates (1A5X) of ASV-IN-Y3 complex. Mn(2+) cation was present in the active site. To neutralize the positive charge in the presence of explicit solvent, eight Cl(-)anions were added. Energy Minimization (EM) and MD simulations, for both the systems, were carried out using Sander's module of AMBER5.0 with all atom force field. We also carried out 1 ns MD simulation on two flexible loops--L1 (Gly54-Gln62) and L2 (Trp138-Met155) playing crucial role in interaction of IN with the drug, under differing environmental conditions (vacuum, aqueous and organic solvent methanol). Comparison of the conformational changes in the loops, monomer and dimer is presented in the paper. Our results showed that the conformation of the loop region was closest to crystallographic data in case of monomer and constrained loops in aqueous environment. However, the dimer in vacuum was more stable than monomer. The beta sheet structure of the monomer in aqueous environment was unstable. Latter also took long time for equilibration. The box formed by loops L1 and L2 from two sub units (IINA and INB) of the dimer satisfies prerequisites for ligand recognition site and seems to be the functional biological unit.     

Conformational relaxation following hydride transfer plays a limiting role in dihydrofolate reductase catalysis

The catalytic cycle of an enzyme is frequently associated with conformational changes that may limit maximum catalytic throughput. In Escherichia coli dihydrofolate reductase, release of the tetrahydrofolate (THF) product is the rate-determining step under physiological conditions and is associated with an "occluded" to "closed" conformational change. In this study, we demonstrate that in dihydrofolate reductase the closed to occluded conformational change in the product ternary complex (E.THF.NADP (+)) also gates progression through the catalytic cycle. Using NMR relaxation dispersion, we have measured the temperature and pH dependence of microsecond to millisecond time scale backbone dynamics of the occluded E.THF.NADP (+) complex. Our studies indicate the presence of three independent dynamic regions, associated with the active-site loops, the cofactor binding cleft, and the C-terminus and an adjacent loop, which fluctuate into discrete conformational substates with different kinetic and thermodynamic parameters. The dynamics of the C-terminally associated region is pH-dependent (p K a < 6), but the dynamics of the active-site loops and cofactor binding cleft are pH-independent. The active-site loop dynamics access a closed conformation, and the accompanying closed to occluded rate constant is comparable to the maximum pH-independent hydride transfer rate constant. Together, these results strongly suggest that the closed to occluded conformational transition in the product ternary complex is a prerequisite for progression through the catalytic cycle and that the rate of this process places an effective limit on the maximum rate of the hydride transfer step.     

Conformational dynamics in loop swap mutants of homologous fibronectin type III domains

Fibronectin type III (FN-III) domains are autonomously folded modules found in a variety of multidomain proteins. The 10th FN-III domain from fibronectin (fnFN10) and the 3rd FN-III domain from tenascin-C (tnFN3) have 27% sequence identity and the same overall fold; however, the CC' loop has a different pattern of backbone hydrogen bonds and the FG loop is longer in fnFN10 compared to tnFN3. To examine the influence of length, sequence, and context in determining dynamical properties of loops, CC' and FG loops were swapped between fnFN10 and tnFN3 to generate four mutant proteins and backbone conformational dynamics on ps-ns and mus-ms timescales were characterized by solution (15)N-NMR spin relaxation spectroscopy. The grafted loops do not strongly perturb the properties of the protein scaffold; however, specific effects of the mutations are observed for amino acids that are proximal in space to the sites of mutation. The amino acid sequence primarily dictates conformational dynamics when the wild-type and grafted loop have the same length, but both sequence and context contribute to conformational dynamics when the loop lengths differ. The results suggest that changes in conformational dynamics of mutant proteins must be considered in both theoretical studies and protein design efforts.     

Cation-mediated interplay of loops in chaperonin-10

The ubiquitously occurring chaperonins consist of a large tetradecameric Chaperonin-60, forming a cylindrical assembly, and a smaller heptameric Chaperonin-10. For a functional protein folding cycle, Chaperonin-10 caps the cylindrical Chaperonin-60 from one end forming an asymmetric complex. The oligomeric assembly of Chaperonin-10 is known to be highly plastic in nature. In Mycobacterium tuberculosis, the plasticity has been shown to be modulated by reversible binding of divalent cations. Binding of cations confers rigidity to the metal binding loop, and also promotes stability of the oligomeric structure. We have probed the conformational effects of cation binding on the Chaperonin-10 structure through fluorescence studies and molecular dynamics simulations. Fluorescence studies show that cation binding induces reduced exposure and flexibility of the dome loop. The simulations corroborate these results and further indicate a complex landscape of correlated motions between different parts of the molecule. They also show a fascinating interplay between two distantly spaced loops, the metal binding "dome loop" and the GroEL-binding "mobile loop", suggesting an important cation-mediated role in the recognition of Chaperonin-60. In the presence of cations the mobile loop appears poised to dock onto the Chaperonin-60 structure. The divalent metal ions may thus act as key elements in the protein folding cycle, and trigger a conformational switch for molecular recognition.     

Molecular dynamics simulations of Guanine quadruplex loops: advances and force field limitations

A computational analysis of d(GGGGTTTTGGGG)(2) guanine quadruplexes containing either lateral or diagonal four-thymidine loops was carried out using molecular dynamics (MD) simulations in explicit solvent, locally enhanced sampling (LES) simulations, systematic conformational search, and free energy molecular-mechanics, Poisson Boltzmann, surface area (MM-PBSA) calculations with explicit inclusion of structural monovalent cations. The study provides, within the approximations of the applied all-atom additive force field, a qualitatively complete analysis of the available loop conformational space. The results are independent of the starting structures. Major conformational transitions not seen in conventional MD simulations are observed when LES is applied. The favored LES structures consistently provide lower free energies (as estimated by molecular-mechanics, Poisson Boltzmann, surface area) than other structures. Unfortunately, the predicted optimal structure for the diagonal loop arrangement differs substantially from the atomic resolution experiments. This result is attributed to force field deficiencies, such as the potential misbalance between solute-cation and solvent-cation terms. The MD simulations are unable to maintain the stable coordination of the monovalent cations inside the diagonal loops as reported in a recent x-ray study. The optimal diagonal and lateral loop arrangements appear to be close in energy although a proper inclusion of the loop monovalent cations could stabilize the diagonal architecture.