De novo design of a pentameric coiled-coil: decoding the motif for tetramer versus pentamer formation in water-soluble phospholamban.

Water-soluble phospholamban (WSPLB) is a designed, water-soluble analogue of the pentameric membrane protein phospholamban (PLB), which contains the same core and interhelical residues as PLB, with only the solvent-exposed positions mutated. WSPLB contains the same secondary and quaternary structure as PLB. The hydrophobic cores of PLB and WSPLB contain Leu and Ile at the a- and d-positions of a heptad repeat (abcdefg) from residues 31-52, while residues 21-30 are rich in polar amino acids at these positions. While the full-length WSPLB forms pentamers in solution, truncated peptides lacking residues 21-30 are largely tetrameric. Thus, truncation of residues 1-20 promotes a switch from pentamer to tetramer formation. Here, the motifs for WSPLB pentamerization were elucidated by characterizing a series of peptides, which were progressively truncated in this polar 'switch' region. When fully present, the 'switch' region promotes pentamer formation in WSPLB, by destabilizing a more stable tetrameric species which exists in its absence. We find that the burial of hydrogen bonding residues from 21 to 30 drives WSPLB from a tetramer to a pentamer, with direct implications for coiled-coil design.     

X-ray structure of a water-soluble analog of the membrane protein phospholamban: sequence determinants defining the topology of tetrameric and pentameric coiled coils

Phospholamban (PLB) is a pentameric transmembrane protein that regulates the Ca(2+)-dependent ATPase SERCA2a in sarcoplasmic reticulum membranes. We previously described the computational design of a water-soluble variant of phospholamban, WSPLB, which reproduced many of the structural and functional properties of the native membrane-soluble protein. While the full-length WSPLB forms a pentamer in solution, a truncated variant forms very stable tetramers. To obtain insight into the tetramer-pentamer cytoplasmic switch, we solved the crystal structure of the truncated construct, WSPLB 21-52. This peptide has a heptad sequence repeat with Leu residues at a- and Ile at d-positions from residues 31-52. The crystal structure revealed that WSPLB 21-52 adopted an antiparallel tetrameric coiled coil. This topology contrasts with the parallel topology of an analogue of the coiled-coil of GCN4 with the same Leu(a) Ile(d) repeat. Analysis of these structures revealed how the nature of the partially exposed residues at e- and g-positions influence the topology formed by the bundle. We also constructed a model for the pentameric form of PLB using the coiled-coil parameters derived from a single monomer in the tetrameric structure. This model suggests that both buried and interfacial hydrogen bonds are important for stabilizing the parallel pentamer.     

(alphaMe)Nva: stereoselective syntheses and preferred conformations of selected model peptides.

Using different stereoselective chemical and chemoenzymatic approaches we synthesized the chiral, Calpha-methylated alpha-amino acid L-(alphaMe)Nva with a short, linear side-chain. A set of terminally protected model peptides to the pentamer level containing either (alphaMe)Nva or Nva in combination with Ala and/or Aib was prepared using solution methods and characterized fully. Two (alphaMe)Nva peptides were also synthesized using side-chain hydrogenation of the corresponding Calpha-methyl, Calpha-allylglycine (Mag) peptides. A detailed solution and crystal-state conformational analysis based on FT-IR absorption, 1H NMR and X-ray diffraction techniques allowed us to define that: (i) (alphaMe)Nva is an effective beta-turn and 3(10)-helix former; and (ii) the relationship between (alphaMe)Nva chirality and the screw sense of the turn/helix formed is that typical of protein amino acids, i.e. L-(alphaMe)Nva induces the preferential formation of right-handed folded structures. In more general terms, this study reinforced previous conclusions that peptides based on alpha-amino acids with a Calpha-methyl substituent and a Calpha-linear alkyl substituent are characterized by a strong tendency to fold into turn and helical structures.     

The protonation state of the catalytic aspartates in plasmepsin II

Assigning the correct protonation state to the catalytic residues is essential for a realistic modelling of an enzyme's active site. Plasmepsins are pharmaceutically relevant aspartic proteases involved in haemoglobin degradation by Plasmodium spp. In aspartic proteases, one of the two catalytic aspartates is protonated, while the other is negatively charged. Here, multiple explicit-water molecular dynamics simulations of plasmepsin II, uncomplexed and with a hydroxypropylamine peptidomimetic inhibitor, indicate that protonation of Asp214 favours a stable active site structure. Moreover, the protonation state of the catalytic aspartate has a strong influence on a linear chain of hydrogen bonds with the adjacent side chains.     

Solution structures by 1H NMR of the novel cyclic trypsin inhibitor SFTI-1 from sunflower seeds and an acyclic permutant

SFTI-1 is a recently discovered cyclic peptide trypsin inhibitor from sunflower seeds comprising 14 amino acid residues. It is the most potent known Bowman-Birk inhibitor and the only naturally occurring cyclic one. The solution structure of SFTI-1 has been determined by 1H-NMR spectroscopy and compared with a synthetic acyclic permutant. The solution structures of both are remarkably similar. The lowest energy structures from each family of 20 structures of cyclic and acyclic SFTI-1 have an rmsd over the backbone and heavy atoms of 0.29 A and 0.66 A, respectively. The structures consist of two short antiparallel beta-strands joined by an extended loop containing the active site at one end. Cyclic SFTI-1 also has a hairpin turn completing the cycle. Both molecules contain particularly stable arrangements of cross-linking hydrogen bonds between the beta-strands and a single disulfide bridge, making them rigid and well defined in solution. These stable arrangements allow both the cyclic and acyclic variants of SFTI-1 to inhibit trypsin with very high potencies (0.5 nM and 12.1 nM, respectively). The cyclic nature of SFTI-1 appears to have evolved to provide higher trypsin inhibition as well as higher stability. The solution structures are similar to the crystal structure of the cyclic inhibitor in complex with trypsin. The lack of a major conformational change upon binding suggests that the structure of SFTI-1 is rigid and already pre-organized for maximal binding due to minimization of entropic losses compared to a more flexible ligand. These properties make SFTI-1 an ideal platform for the design of small peptidic pharmaceuticals or pesticides.     

Unfolding of an alpha-helix in water.

We describe a 1 ns molecular dynamics simulation of an 18-residue peptide (corresponding to a portion of the H helix of myoglobin) in water. The initial helical conformation progressively frays to a more disordered structure, with the loss of internal secondary structure generally proceeding from the C-terminus toward the N-terminus. Although a variety of mechanisms are involved in the breaking of helical hydrogen bonds, the formation of transient turn structures, with i----i + 3 hydrogen bonds, and bifurcated hydrogen-bond structures intermediate between alpha and turn or 3(10) structures is a common motif. In some cases a single water molecule is inserted into an internal hydrogen bond, but it is also common to have several water molecules involved in transient intermediates. Overall, the results provide new information about the detailed mechanisms by which helices are made and broken in aqueous solution.     

Structure of the 1-36 amino-terminal fragment of human phospholamban by nuclear magnetic resonance and modeling of the phospholamban pentamer

The structure of a 36-amino-acid-long amino-terminal fragment of phospholamban (phospholamban[1-36]) in aqueous solution containing 30% trifluoroethanol was determined by nuclear magnetic resonance. The peptide, which comprises the cytoplasmic domain and six residues of the transmembrane domain of phospholamban, assumes a conformation characterized by two alpha-helices connected by a turn. The residues of the turn are Ile18, Glu19, Met20, and Pro21, which are adjacent to the two phosphorylation sites Ser16 and Thr17. The proline is in a trans conformation. The helix comprising amino acids 22-36 is well determined (the root mean square deviation for the backbone atoms, calculated for a family of 18 nuclear magnetic resonance structures is 0.57 A). Recently, two molecular models of the transmembrane domain of phospholamban were proposed in which a symmetric homopentamer is composed of a left-handed coiled coil of alpha-helices. The two models differ by the relative orientation of the helices. The model proposed by,Simmerman et al. (H.K. Simmerman, Y.M. Kobayashi, J.M. Autry, and L.R. Jones, 1996, J. Biol. Chem. 271:5941-5946), in which the coiled coil is stabilized by a leucine-isoleucine zipper, is similar to the transmembrane pentamer structure of the cartilage oligomeric membrane protein determined recently by x-ray (V. Malashkevich, R. Kammerer, V Efimov, T. Schulthess, and J. Engel, 1996, Science 274:761-765). In the model proposed by Adams et al. (P.D. Adams, I.T. Arkin, D.M. Engelman, and A.T. Brunger, 1995, Nature Struct. Biol. 2:154-162), the helices in the coiled coil have a different relative orientation, i.e., are rotated clockwise by approximately 50 degrees. It was possible to overlap and connect the structure of phospholamban[1-36] derived in the present study to the two transmembrane pentamer models proposed. In this way two models of the whole phospholamban in its pentameric form were generated. When our structure was connected to the leucine-isoleucine zipper model, the inner side of the cytoplasmic domain of the pentamer (where the helices face one another) was lined by polar residues (Gln23, Gln26, and Asn30), whereas the five Arg25 side chains were on the outer side. On the contrary, when our structure was connected to the other transmembrane model, in the inner side of the cytoplasmic domain of the pentamer, the five Arg25 residues formed a highly charged cluster.     

Characterization of crystal water molecules in a high-affinity inhibitor and hematopoietic prostaglandin D synthase complex by interaction energy studies

Hematopoietic prostaglandin D synthase (H-PGDS) is one of the two enzymes that catalyze prostaglandin D₂ synthesis and a potential therapeutic target of allergic and inflammatory responses. To reveal key molecular interactions between a high-affinity ligand and H-PGDS, we designed and synthesized a potent new inhibitor (K_D: 0.14?nM), determined the crystal structure in complex with human H-PGDS, and quantitatively analyzed the ligand–protein interactions by the fragment molecular orbital calculation method. In the cavity, 10 water molecules were identified, and the interaction energy calculation indicated their stable binding to the surface amino acids in the cavity. Among them, 6 water molecules locating from the deep inner cavity to the peripheral part of the cavity contributed directly to the ligand binding by forming hydrogen bonding interactions. Arg12, Gly13, Gln36, Asp96, Trp104, Lys112 and an essential co-factor glutathione also had strong interactions with the ligand. A strong repulsive interaction between Leu199 and the ligand was canceled out by forming a hydrogen bonding network with the adjacent conserved water molecule. Our quantitative studies including crystal water molecules explained that compounds with an elongated backbone structure to fit from the deep inner cavity to the peripheral part of the cavity would have strong affinity to human H-PGDS.     

Rational design of carbonitrile-carboxaldehyde cation receptor models: probing the nature of the heteroatom–metal interaction

Amino acids are constituents of proteins and enzymes which take part almost in all metabolic reactions. Glutamic acid, with an ability to form a negatively charged side chain, plays a major role in intra and intermolecular interactions of proteins, peptides, and enzymes. An exhaustive conformational analysis has been performed for all eight possible forms at B3LYP/cc-pVTZ level. All possible neutral, zwitterionic, protonated, and deprotonated forms of glutamic acid structures have been investigated in solution by using polarizable continuum model mimicking water as the solvent. Nine families based on the dihedral angles have been classified for eight glutamic acid forms. The electrostatic effects included in the solvent model usually stabilize the charged forms more. However, the stability of the zwitterionic form has been underestimated due to the lack of hydrogen bonding between the solute and solvent; therefore, it is observed that compact neutral glutamic acid structures are more stable in solution than they are in vacuum. Our calculations have shown that among all eight possible forms, some are not stable in solution and are immediately converted to other more stable forms. Comparison of isoelectronic glutamic acid forms indicated that one of the structures among possible zwitterionic and anionic forms may dominate over the other possible forms. Additional investigations using explicit solvent models are necessary to determine the stability of charged forms of glutamic acid in solution as our results clearly indicate that hydrogen bonding and its type have a major role in the structure and energy of conformers.     

High-resolution structures of scytalone dehydratase-inhibitor complexes crystallized at physiological pH

Scytalone dehydratase is a molecular target of inhibitor design efforts aimed at preventing the fungal disease caused by Magnaporthe grisea. A method for cocrystallization of enzyme with inhibitors at neutral pH has produced several crystal structures of enzyme-inhibitor complexes at resolutions ranging from 1.5 to 2.2 A. Four high resolution structures of different enzyme-inhibitor complexes are described. In contrast to the original X-ray structure of the enzyme, the four new structures have well-defined electron density for the loop region comprising residues 115-119 and a different conformation between residues 154 and 160. The structure of the enzyme complex with an aminoquinazoline inhibitor showed that the inhibitor is in a position to form a hydrogen bond with the amide of the Asn131 side chain and with two water molecules in a fashion similar to the salicylamide inhibitor in the original structure, thus confirming design principles. The aminoquinazoline structure also allows for a more confident assignment of donors and acceptors in the hydrogen bonding network. The structures of the enzyme complexes with two dichlorocyclopropane carboxamide inhibitors showed the two chlorine atoms nearly in plane with the amide side chain of Asn131. The positions of Phe53 and Phe158 are significantly altered in the new structures in comparison to the two structures obtained from crystals grown at acidic pH. The multiple structures help define the mobility of active site amino acids critical for catalysis and inhibitor binding.     

The preferred conformation of the tripeptide Ala-Phe-Ala in water is an inverse gamma-turn: implications for protein folding and drug design

Recent studies have provided evidence that peptides as short as tripeptides do adopt preferred conformations. Here we report that the tripeptide Ala-Phe-Ala (AFA) in aqueous solution preferentially forms an inverse gamma-turn. Circular dichroism (CD) indicated the presence of a predominant turn structure, and Fourier transform infrared (FTIR) bands suggested the presence of a gamma-turn forming a bifurcated H-bond with the solvent molecules. The high-resolution structure was obtained by a combined use of NMR spectroscopy and calculations. On the basis of 30 unambiguous ROESY-derived distance restraints (including the Halpha-NH NOE between Ala(1) and Ala(3) and a hydrogen bond between the CO group of Ala(1) and the NH group of Ala(3)), calculations clearly demonstrated the presence of an inverse gamma-turn centered on Phe(2). From NOE data, we estimated a mole fraction for the gamma-turn of 0.65. Since for AFA an extended beta-strand was also reported [Eker, F., Griebenow, K., Cao, X., Nafie, L. A., and Schweitzer-Stenner, R. (2004) Proc. Natl. Acad. Sci. U.S.A. 101, 10054-10059], we investigated the possibility that gamma-turn and beta-strand may represent two major conformations. By using a best-fit procedure that calculated experimental NOEs as weighted averages of the effects originating from both structures, we were able to calculate with good accuracy the backbone NOEs at 280 K in terms of the two limiting conformers, yielding a mole fraction for the gamma-turn and beta-strand conformations of 0.60 and 0.40, respectively, in good agreement with those found by NOE data. The implication of the existence of a preferred conformation by a small structural element is discussed in the context of the nucleation of protein folding events and the design of small peptide and peptidomimetic drugs.     

Molecular packing and intermolecular interactions in two structural polymorphs of N-palmitoylethanolamine, a type 2 cannabinoid receptor agonist

The molecular structure, packing properties, and intermolecular interactions of two structural polymorphs of N-palmitoylethanolamine (NPEA) have been determined by single-crystal X-ray diffraction. Polymorphs alpha and beta crystallized in monoclinic space group P2(1)/c and orthorhombic space group Pbca, respectively. In both polymorphs, NPEA molecules are organized in a tail-to-tail manner, resembling a bilayer membrane. Although the molecular packing in polymorph alpha is similar to that in N-myristoylethanolamine and N-stearoylethanolamine, polymorph beta is a new form. The acyl chains in both polymorphs are tilted by approximately 35 degrees with respect to the bilayer normal, with their hydrocarbon moieties packed in an orthorhombic subcell. In both structures, the hydroxy group of NPEA forms two hydrogen bonds with the hydroxy groups of molecules in the opposite leaflet, resulting in extended, zig-zag type H-bonded networks along the b-axis in polymorph alpha and along the a-axis in polymorph beta. Additionally, the amide N-H and carbonyl groups of adjacent molecules are involved in N-H...O hydrogen bonds that connect adjacent molecules along the b-axis and a-axis, respectively, in alpha and beta. Whereas in polymorph alpha the L-shaped NPEA molecules in opposite layers are arranged to yield a Z-like organization, in polymorph beta one of the two NPEA molecules is rotated 180 degrees , leading to a W-like arrangement. Lattice energy calculations indicate that polymorph alpha is more stable than polymorph beta by approximately 2.65 kcal/mol.     

Thermodynamic stability of beta-peptide helices and the role of cyclic residues

Beta-peptides are emerging as an attractive class of peptidomimetic molecules. In contrast to naturally occurring alpha-peptides, short oligomers of beta-amino acids (comprising just 4-6 monomers) exhibit stable secondary structures that make them amenable for quantitative, concerted experimental and theoretical studies of the effects of particular chemical interactions on structure. In this work, molecular simulations are used to study the thermodynamic stability of helical conformations formed by beta-peptides containing varying proportions of acyclic (beta(3)) and cyclic (ACH) residues. More specifically, several beta-peptides differing only in their content of cyclic residues are considered in this work. Previous computational studies of beta-peptides have relied mostly on energy minimization of molecular dynamics simulations. In contrast, our study relies on density-of-states based Monte Carlo simulations to calculate the free energy and examine the stability of various folded structures of these molecules along a well-defined order parameter. By resorting to an expanded-ensemble formalism, we are able to determine the free energy required to unfold specific molecules, a quantity that could be measured directly through single-molecule force spectroscopy. Simulations in both implicit and explicit solvents have permitted a systematic study of the role of cyclic residues and electrostatics on the stability of secondary structures. The molecules considered in this work are shown to exhibit stable H-14 helical conformations and, in some cases, relatively stable H-12 conformations, thereby suggesting that solvent quality may be used to manipulate the hydrogen-bonding patterns and structure of these peptides.     

New reagents, reactions, and peptidomimetics for drug design.

It has been a major focus in our laboratories to prepare novel reagents and peptidomimetic structures for drug design. We have designed and prepared novel guanidinylation reagents that can be employed in solution or as solid phase reagents. We and others have utilized the reagent 3-(diethoxyphosphoryloxy)-1,2,3-benzotriazin-4(3H)-one (DEPBT) for amide bond formation to couple sterically hindered structures. These couplings proceed with remarkably strong resistance to racemization. In the area of peptidomimetics, we have incorporated novel building blocks to create biologically active compounds. These building blocks include thioether and alkylamine bridges, beta-methylated, and beta,beta-dimethylated amino acid residues. These mimetic structures have been incorporated into specific target molecules such as opioids to obtain cyclic peptidomimetics with potent and selective biological activity.     

Crystal structure of an archaeal pentameric riboflavin synthase in complex with a substrate analog inhibitor: stereochemical implications

Whereas eubacterial and eukaryotic riboflavin synthases form homotrimers, archaeal riboflavin synthases from Methanocaldococcus jannaschii and Methanothermobacter thermoautrophicus are homopentamers with sequence similarity to the 6,7-dimethyl-8-ribityllumazine synthase catalyzing the penultimate step in riboflavin biosynthesis. Recently it could be shown that the complex dismutation reaction catalyzed by the pentameric M. jannaschii riboflavin synthase generates riboflavin with the same regiochemistry as observed for trimeric riboflavin synthases. Here we present crystal structures of the pentameric riboflavin synthase from M. jannaschii and its complex with the substrate analog inhibitor, 6,7-dioxo-8-ribityllumazine. The complex structure shows five active sites located between adjacent monomers of the pentamer. Each active site can accommodate two substrate analog molecules in anti-parallel orientation. The topology of the two bound ligands at the active site is well in line with the known stereochemistry of a pentacyclic adduct of 6,7-dimethyl-8-ribityllumazine that has been shown to serve as a kinetically competent intermediate. The pentacyclic intermediates of trimeric and pentameric riboflavin synthases are diastereomers.     

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.     

Triple helical structure and stabilization of collagen-like molecules with 4(R)-hydroxyproline in the Xaa position

In this study, we examine the relationships between the structure and stability of five related collagen-like molecules that have hydroxyproline residues occupying positions not observed in vertebrate collagen. Two of the molecules contain valine or threonine and form stable triple helices in water. Three of the molecules contain allo-threonine (an enantiomer of threonine), serine, or alanine, and are not stable. Using molecular dynamics simulation methods, we examine possible explanations for the stability difference, including considering the possibility that differences in solvent shielding of the essential interchain hydrogen bonds may result in differences in stability. By comparing the structures of threonine- and allo-threonine-containing molecules in six polar and nonpolar solvation conditions, we find that solvent shielding is not an adequate explanation for the stability difference. A closer examination of the peptides shows that the structures of the unstable molecules are looser, having weaker intermolecular hydrogen bonds. The weakened hydrogen bonds result from extended Yaa residue Psi-angles that prevent optimal geometry. The Phi-Psi-maps of the relevant residues suggest that each residue's most favorable Psi-angle determines the corresponding collagen-like molecule's stability. Additionally, we propose that these molecules illustrate a more general feature of triple-helical structures: interchain hydrogen bonds are always longer and weaker than ideal, so they are sensitive to relatively small changes in molecular structure. This sensitivity to small changes may explain why large stability differences often result from seemingly small changes in residue sequence.     

Helix stability of oligoglycine,oligoalanine, and oligo‐β‐alanine dodecamers reflected by hydrogen‐bond persistence

Helices are important structural/recognition elements in proteins and peptides. Stability and conformational differences between helices composed of α‐ and β‐amino acids as scaffolds for mimicry of helix recognition has become a theme in medicinal chemistry. Furthermore, helices formed by β‐amino acids are experimentally more stable than those formed by α‐amino acids. This is paradoxical because the larger sizes of the hydrogen‐bonding rings required by the extra methylene groups should lead to entropic destabilization. In this study, molecular dynamics simulations using the second‐generation force field, AMOEBA (Ponder, J.W., et al., Current status of the AMOEBA polarizable force field. J Phys Chem B, 2010. 114 (8): p. 2549–64.) explored the stability and hydrogen‐bonding patterns of capped oligo‐β‐alanine, oligoalanine, and oligoglycine dodecamers in water. The MD simulations showed that oligo‐β‐alanine has strong acceptor+2 hydrogen bonds, but surprisingly did not contain a large content of 3₁₂‐helical structures, possibly due to the sparse distribution of the 3₁₂‐helical structure and other structures with acceptor+2 hydrogen bonds. On the other hand, despite its backbone flexibility, the β‐alanine dodecamer had more stable and persistent <3.0 Å hydrogen bonds. Its structure was dominated more by multicentered hydrogen bonds than either oligoglycine or oligoalanine helices. The 3₁ (PII) helical structure, prevalent in oligoglycine and oligoalanine, does not appear to be stable in oligo‐β‐alanine indicating its competition with other structures (stacking structure as indicated by MD analyses). These differences are among the factors that shape helical structural preferences and the relative stabilities of these three oligopeptides. Proteins 2014; 82:3043–3061. © 2014 Wiley Periodicals, Inc.     

Effect of Coordinated Ions on Structure and Flexibility of Parallel G-quadruplexes: A Molecular Dynamics Study

Abstract Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na(+) or K(+) ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 ? from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na(+) counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na(+) ions, within the quadruplex cavity, are more mobile than coordinated K(+) ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.