On the stability of single- and double-stranded DNA helices: The application of the PPT-MCF method on large fragments of DNA

The pseudo-polarization tensor mutually consistent field (PPT -MCF ) method recently introduced [1] has been applied to study the stacking interactions between the nucleotide bases in large periodic B-DNA fragments. The effects on the global and local binding properties caused by replacing one base in the periodic sequence by another base are investigated. The increase in the stability for comparable fragments owing to this base substitution is further enforced in the case of periodic alternating helices. The most important results are that the stacking interaction between two bases is slowly converging with the interbase distance and that the average contribution per base to the binding energy is repulsive. Furthermore, the energetical properties of double helix models in B- and Z-DNA configurations, respectively, consisting of up to five base pairs have been compared. It turns out that the G C G C sequence in Z-DNA is significantly more stable than either in periodic or periodic alternating B-DNA. In these cases the average energy contribution of a single Watson–Crick-type base pair is predicted also to be positive. From the calculations it follows that the double helix is not stabilized owing to the hydrogen bonding between the bases belonging to both strands, in contradiction to most other investigations.     

Stability and Unfolding Mechanism of the N‐terminal β‐Hairpin from [2Fe‐2S] Ferredoxin I by Molecular Dynamics Simulations

The stability and unfolding mechanism of the N‐terminal β‐hairpin of the [2Fe‐2S] ferredoxin I from the blue‐green alga Aphanothece sacrum in pure methanol, 40% (v/v) methanol‐water, and pure water systems were investigated by 10 ns molecular dynamics simulations under periodic boundary conditions. The β‐hairpin was mostly in its native‐like state in pure methanol, whereas it unfolds dramatically following the ‘zip‐up’ mechanism when it was placed in pure water. Both interstrand and inside‐turn hydrogen bonds account for the stability of the β‐hairpin in its native‐like conformation, whereas hydrophobic interactions among nonpolar side chains are responsible for maintaining its stable loop‐like intermediate structures in 40% (v/v) methanol‐water. Reducing solvent polarity seems to increase the stability of the β‐hairpin in its native‐like structure. Methanol is likely to mimic the partially hydrophobic environment around the N‐terminal β‐hairpin by the subsequent α‐helix.     

Relative contributions of desolvation, inter- and intramolecular interactions to binding affinity in protein kinase systems

In several previous studies, we performed sensitivity analysis to gauge the relative importance of different atomic partial charges in determining protein-ligand binding. In this work, we gain further insights by decomposing these results into three contributions: desolvation, intramolecular interactions, and intermolecular interactions, again based on a Poisson continuum electrostatics model. Three protein kinase-inhibitor systems have been analyzed: CDK2-deschloroflavopiridol, PKA-PKI, and LCK-PP2. Although our results point out the importance of specific intermolecular interactions to the binding affinity, they also reveal the remarkable contributions from the solvent-mediated intramolecular interactions in some cases. Thus, it is necessary to look beyond analyzing protein-ligand interactions to understand protein-ligand recognition or to gain insights into designing ligands and proteins. In analyzing the contributions of the three components to the overall binding free energy, the PKA-PKI system with a much larger ligand was found to behave differently from the other two systems with smaller ligands. In the former case, the intermolecular interactions are very favorable, and together with the favorable solvent-mediated intramolecular interactions, they overcome the large desolvation penalties to give a favorable electrostatics contribution to the overall binding affinity. On the other hand, the other two systems with smaller ligands only present modest intermolecular interactions and they are not or are only barely sufficient to overcome the desolvation penalty even with the aid of the favorable intramolecular contributions. As a result, the binding affinity of these two systems do not or only barely benefit from electrostatics contributions.     

Interstrand side chain--side chain interactions in a designed beta-hairpin: significance of both lateral and diagonal pairings

The contributions of interstrand side chain-side chain contacts to beta-sheet stability have been examined with an autonomously folding beta-hairpin model system. RYVEV(D)PGOKILQ-NH2 ((D)P = D-proline, O = ornithine) has previously been shown to adopt a beta-hairpin conformation in aqueous solution, with a two-residue loop at D-Pro-Gly. In the present study, side chains that display interstrand NOEs (Tyr-2, Lys-9, and Leu-11) are mutated to alanine or serine, and the conformational impact of the mutations is assessed. In the beta-hairpin conformation Tyr-2 and Leu-11 are directly across from one another (non-hydrogen bonded pair). This "lateral" juxtaposition of two hydrophobic side chains appears to contribute to beta-hairpin conformational stability, which is consistent with results from other beta-sheet model studies and with statistical analyses of interstrand residue contacts in protein crystal structures. Interaction between the side chains of Tyr-2 and Lys-9 also stabilizes the beta-hairpin conformation. Tyr-2/Lys-9 is a "diagonal" interstrand juxtaposition because these residues are not directly across from one another in terms of the hydrogen bonding registry between the strands. This diagonal interaction arises from the right-handed twist that is commonly observed among beta-sheets. Evidence of diagonal side chain-side chain contacts has been observed in other autonomously folding beta-sheet model systems, but we are not aware of other efforts to determine whether a diagonal interaction contributes to beta-sheet stability.     

Polymers in an external periodic field: the effect of the excluded volume interactions

By developing and making use of the "transfer operator" formalism, we calculate the number density and average Flory end-to-end distance of the polymers placed in an external periodic field. The considered mathematical problem is of immediate relevance for such realistic physical systems as the homopolymers immersed in the host structure of alternating layers that have different affinities for homopolymers (e.g., lamellar microphases of copolymers, ripple morphology of the mixed brush, and lipidwater systems). In contrast to the conventional ground state dominance approximation, the developed method makes it possible to calculate the characteristic size (Flory radius R(F)) of the polymers in the direction of applied external periodic field, with the effect of the excluded volume taken into account. The excluded volume interactions are shown to qualitatively change the behavior of R(F) as a function of the reduced field strength theta relative to the case of ideal Gaussian polymers. In particular, in the limit of strong fields theta>1 the average Flory radius R(F) is found to saturate to its minimal value, which is calculated as a function of the excluded volume parameter u. This finding is in distinct contrast to the result for the Flory radius R(F) in the case of ideal polymers where R(F) approaches zero as the interaction parameter theta increases.     

Folding of a stable DNA motif involves a highly cooperative network of interactions

Hairpins are structural elements that play important roles in the folding and function of RNA and DNA. The extent of cooperativity in folding is an important aspect of the RNA folding problem. We reasoned that an investigation into the origin of cooperativity might be best carried out on a stable nucleic acid system with a limited number of interactions, such as a stable DNA hairpin loop. The stable d(cGNAg) hairpin loop motif (closing base pair in lower case; loop in upper case; N = A, C, G, or T) is stabilized through only three interactions: two loop-loop hydrogen bonds in a sheared GA base pair and a loop-closing base pair interaction. Herein, we investigate this network of interactions and test whether the loop-loop and loop-closing base pair interactions communicate. Thermodynamic measurements of nucleotide analogue substituted oligonucleotides were used to probe the additivity of the interactions. On the basis of double mutant cycles, all interactions were found to be nonadditive and interdependent, suggesting that loop-loop and loop-closing base pair interactions form in a highly cooperative manner. When double mutant cycles were repeated in the absence of the other interaction, nonadditivity was significantly reduced suggesting that coupling is indirect and requires all three interactions in order to be optimal. A cooperative network of interactions helps explain the structural and energetic bases of stability in certain DNA hairpins and paves the way for similar studies in more complex nucleic acid systems.     

Dynamics of noncovalent interactions in all-α and all-β class proteins: implications for the stability of amyloid aggregates

A fully folded functional protein is stabilized by several noncovalent interactions. When a protein undergoes conformational motions, the existing noncovalent interactions may be maintained. They may also break or new interactions may be formed. Knowledge of the dynamical nature of the different types of noncovalent interactions is extremely important to understand the structural stability, function, and folding of a protein. There are experimental limitations to investigate the dynamics of different noncovalent interactions simultaneously in a biomolecule. We have carried out molecular dynamics simulations on four different proteins, two belonging to all-α class proteins and the other two are representatives of all-β class proteins. The dynamical nature of eight different noncovalent interactions was studied by monitoring the maximum residence time (MRT) and lifetime (LT). The conventional hydrogen bonds are the dominant interactions in all four proteins, and the majority of those formed between the main-chain atoms were maintained during most of the simulation time with MRT greater than 10 ns. Such interactions with more than 1 ns lifetime provide stability to the secondary structures, and hence they are responsible for the overall stability of the protein. The weak C-H···O hydrogen bond is the next major type of interactions. However, a large number of such interactions are observed between the main-chain atoms only in all-β proteins as interstrand interactions, and, surprisingly, they are observed during most part of the simulation although their average lifetime is only about 20 to 30 ps. The strong cation···π and salt-bridge interactions are present few in number. However, in many cases they are almost uninterrupted indicating the higher strength of these interactions. Four other interactions involving the π-electron cloud of aromatic rings are very small in number, and, in many cases, their presence is not maintained throughout the simulation. Our results clearly indicate that the weak C-H···O interactions between the main-chain atoms are the distinguishing factor between the all-α and all-β class of proteins, and these interstrand interactions can provide additional stability to all-β protein structures. Based on these results, we hypothesize that such weak C-H···O interstrand interactions could play a major role in providing stability to amyloid type of aggregates that are responsible for the pathological state of many proteins.     

ICMMM2D: an accurate method to include planar dielectric interfaces via image charge summation

An extension of the MMM2D method, called ICMMM2D, is presented to deal with the electrostatic interactions in partially periodic systems that are confined along one direction by two planar dielectric interfaces. The method handles all electrostatic interactions that are induced by the presence of dielectric interfaces by an image charge method. Our method accurately treats repeated image charges under the planar dielectric interfaces as well as the periodic images that are due to the periodic boundary conditions along the other two directions. The scaling of the method with the number of charges, N, is still N5/3, and the overhead involved approximately doubles the CPU time compared to the original MMM2D method. The errors are fully under control and the error bounds can be preset up to computer accuracy.     

Structure and conformational behavior of biopolymers by density functional calculations employing periodic boundary conditions. I. The case of polyglycine, polyalanine, and poly-alpha-aminoisobutyric acid in vacuo.

Fully quantum mechanical calculations exploiting periodic boundary conditions (PBC) have been applied to the study of four different regular structures (alpha- and 3(10)-helix, fully extended and repeated gamma-turns) of the infinite polypeptides of glycine, alanine, and alpha-aminoisobutyric acid (Aib) in vacuo. alpha-Helix is predicted to be the most stable conformer for polyalanine and polyglycine, being stabilized over the 3(10)-helix mainly by more favorable dipole-dipole interactions. Contrary to previous suggestions, steric effects and hydrogen-bond strengths are comparable for both helix structures. 3(10)-Helix is preferred for poly-Aib, since in this case alpha-helix is strongly distorted due to unfavorable intrachain repulsions. Extended structures and repeated gamma-turns are much less stable than helix structures for all of the polypeptides examined, mainly due to the absence of favorable long-range interactions. The optimized geometries are in good agreement with the available experimental data and reveal a remarkable dependence on the nature of the residue forming the polypeptides; at the same time the electronic and structural parameters of each residue strongly depend on the secondary structure of the polypeptides.     

How do size-expanded DNA nucleobases enhance duplex stability? Computational analysis of the hydrogen-bonding and stacking ability of xDNA bases

Computational chemistry (B3LYP, MP2) is used to study the properties of size-expanded DNA nucleobases generated by inserting a benzene spacer into the natural nucleobases. Although the addition of the spacer does not significantly affect the hydrogen-bonding properties of natural nucleobases, the orientation of the base about the glycosidic bond necessary for Watson-Crick binding is destabilized, which could have implications for the selectivity of expanded bases, as well as the stability of expanded duplexes. Consideration of the (stacked) binding energies in the preferred relative orientation of natural and expanded nucleobases aligned according to their centers of mass reveals that the stacking within natural dimers can be increased by up to 50% upon expansion of one nucleobase and up to 90% upon expansion of two nucleobases. The implications of these findings to the stability of expanded duplexes were revealed by considering simplified models of natural and mixed duplexes composed of four nucleobases. Although intra- and interstrand interactions within double helices are typically less than those predicted when nucleobases are stacked according to their centers of mass, some nucleobases utilize their full stacking potential within double helices, where both intra- and interstrand interactions can be significant. Most importantly, increasing the size of nucleobases within the duplex significantly increases both intra- and interstrand stacking interactions. Specifically, some interactions are double the magnitude of the corresponding intrastrand interactions in natural helices, and even greater increases in interstrand interactions are sometimes found. Thus, our work suggests that mixed duplexes composed of natural bases hydrogen bound to expanded bases may exploit the increase in the inherent stacking ability of the expanded bases in more than one way and thereby afford duplexes with greater stability than natural DNA.     

Development of the cyclic cluster model formalism for Kohn-Sham auxiliary density functional theory methods

The development of the cyclic cluster model (CCM) formalism for Kohn-Sham auxiliary density functional theory (KS-ADFT) methods is presented. The CCM is a direct space approach for the calculation of perfect and defective systems under periodic boundary conditions. Translational symmetry is introduced in the CCM by integral weighting. A consistent weighting scheme for all two-center and three-center interactions appearing in the KS-ADFT method is presented. For the first time, an approach for the numerical integration of the exchange-correlation potential within the cyclic cluster formalism is derived. The presented KS-ADFT CCM implementation was applied to covalent periodic systems. The results of cyclic and molecular cluster model (MCM) calculations for trans-polyacetylene, graphene, and diamond are discussed as examples for systems periodic in one, two, and three dimensions, respectively. All structures were optimized. It is shown that the CCM results represent the results of MCM calculations in the limit of infinite molecular clusters. By analyzing the electronic structure, we demonstrate that the symmetry of the corresponding periodic systems is retained in CCM calculations. The obtained geometric and electronic structures are compared with available data from the literature.     

Interstrand cross-link of DNA by covalently linking a pair of abasic sites

A pair of apurinic/apyrimidinic sites formed in DNA has been covalently connected with bis(aminooxy) derivatives. The efficacy of the interstrand cross-link is associated with the structural tethering of two aminooxy groups. The interstrand cross-link constructed stable DNA scaffolds for enzyme alignment.     

New aspects of the repair and genotoxicity of psoralen photoinduced lesions in DNA.

Several approaches are described aiming at a better understanding of the genotoxicity of psoralen photoinduced lesions in DNA. Psoralens can photoinduce different types of photolesions including 3,4- and 4',5'-monoadducts and interstrand cross-links, oxidative damage (in the case of 3-carbethoxypsoralen (3-CPs)) and even pyrimidine dimers (in the case of 7-methylpyrido(3,4-c)psoralen (MePyPs)). The characterization and detection of different types of lesions has been essential for the analysis of their possible contributions to genotoxicity. For example, oxidative damage photoinduced by 3-CPs can be detected by the formamidopyrimidine glycosylase (FPG) protein. Furthermore, it is shown how the presence of MePyPs induced monoadducts may interfere with the photoreactivation of concomitantly induced pyrimidine dimers, how the ratio of monoadducts and interstrand cross-links (CL) affects the occurrence of double-strand breaks during the repair of photolesions and genotoxicity. In vitro treatment of yeast plasmids, followed by transformation, also indicates that the repair of photoadducts on exogenous DNA differs for 8-methoxy-psoralen (8-MOP) induced mono- and diadducts and for monoadducts alone. The recombinational rad52 dependent pathway is not needed for the repair of 8-MOP induced monoadducts. The results obtained suggest that the genotoxic effects of psoralens are conditioned by the nature, number, ratio and sequence distribution of the photolesions induced in DNA.     

The stability of columns comprising alternating triphenylene and hexaphenyltriphenylene molecules: variations in the structure of the hexaphenyltriphenylene component

Previous investigations of complementary polytopic interaction (CPI) columnar mesophases, in which the columns are built up of alternating hexaalkoxytriphenylene (HAT) and hexaphenyltriphenylene (PTP) molecules, concentrated mainly on the effect of variations in the structure of the HAT component. This investigation is concerned with the effect of variations in the structure of the PTP component and, in particular, variations in the position of an alkoxy side chain in the phenyl ring. Stable columnar mesophases are obtained when a hexyloxy substituent is placed in the meta- or para-position but not in the ortho-position. In the case of the meta- and para-substituted systems, the two-component CPI columnar phases are stable over a considerably larger temperature range than the one-component HAT systems. The evidence suggests that unfavourable PTP/PTP stacking is as much a driving force for the formation of these mixed stacks as is favourable PTP/HAT stacking, but both need to be explained in terms of the sum of atomically dispersed van der Waals and coulombic interactions. On cooling from the isotropic into the Col_h phase, the columnar phase based on a 1:1 mixture of hexakis(hexyloxy)tripenylene and the meta-hexyloxy-substituted PTP gives an unusual texture consisting of 'viking-axe'-shaped structures.     

Arrangement and nature of intermolecular hydrogen bonding in complex biomolecular systems: modeling the vitamin C---L-alanine interaction

Vitamin C is known as an essential dietary supplement and implicated in diverse biological processes. We present here a theoretical study on the nature of hydrogen bonding of vitamin C in biological systems. For this reason, the complexes of vitamin C (VC) with neutral and zwitterionic L-alanine (as the simplest chiral amino acid) were studied at the MP2/6-311++G(d,p) level of theory. In the gas phase, neutral L-alanine leads to more stable complexes than the zwitterionic forms while the reverse is true in the aqueous phase. The complexes are formed via two hydrogen bond interactions, which result in a ring-like hydrogen-bonded networks. The nature of H-bonds was characterized in terms of natural bond orbital and quantum theory of atoms in molecule analyses (QTAIM). The H-bonds in the studied complexes were electrostatic in nature; however, in the case of shorter and directional H-bonds and ionic interactions, contributions of covalent character were also non-negligible. Natural energy decomposition analysis of hydrogen-bonded complexes reveals that the charge transfer and electrical components are the largest contributors for the interaction energies of complexes. Natural resonance theory analysis suggests higher resonance weight for charge-assisted interactions of vitamin C---alanine (zwitterionic) complexes, where the total interaction energy is considerably higher than that of neutral alanine.     

Conformational analysis—CXIV : Molecular mechanics studies of sulfoxides

The conformational characteristics of dimethyl sulfoxide and a number of 6-membered ring sulfoxides have been studied. The thiane, 1,3- and 1,4-dithiane and 1,3,5-trithiane ring systems with various oxide substitutions have been examined. It is found that the chair forms are more stable than the twist or boat in all cases. The energy profiles of the twist-boat manifold are in many cases highly unusual, and quite different from anything so far known experimentally. The axial-equatorial preference of the oxygen is highly variable, depending on the steric and electrostatic interactions found in the particular case.     

Dielectric response of modified Hubbard models with neutral-ionic and Peierls transitions

The dipole P(F) of systems with periodic boundary conditions in a static electric field F is applied to one-dimensional Peierls-Hubbard models for organic charge-transfer (CT) salts. Exact results for P(F) are obtained for finite systems of N=14 and 16 sites that are almost converged to infinite chains in deformable lattices subject to a Peierls transition. The electronic polarizability per site, alpha(el)=(partial differential P/partial differential F)0, of rigid stacks with alternating transfer integrals t(1+/-delta) diverges at the neutral-ionic transition for delta=0 but remains finite for delta>0 in dimerized chains. The Peierls or dimerization mode couples to charge fluctuations along the stack and results in large vibrational contributions alpha(vib) that are related to partial differential P/ partial differential delta and that peak sharply at the Peierls transition. The extension of P(F) to correlated electronic states yields the dielectric response kappa of models with neutral-ionic or Peierls transitions, where kappa peaks >100 are found with parameters used previously for variable ionicity rho and vibrational spectra of CT salts. The calculated kappa accounts for the dielectric response of CT salts based on substituted TTF's (tetrathiafulvalene) and substituted CA's (chloranil). The role of lattice stiffness appears clearly in models: soft systems have a Peierls instability at small rho and continuous crossover to large rho, while stiff stacks such as TTF-CA have a first-order transition with discontinuous rho that is both a neutral-ionic and Peierls transition. The transitions are associated with tuning the electronic ground state of insulators via temperature or pressure in experiments, or via model parameters in calculations.     

Duplex structure of a minimal nucleic acid

The crystal structure of an 8-mer (S)-GNA duplex is presented. As a tool for phasing, the anomalous diffraction of two copper(II) ions within two artificial metallo-base pairs was employed. The duplex structure confirms a canonical Watson-Crick base pairing scheme of GNA with antiparallel strands. The duplex secondary structure is distinct from canonical A- and B-form nucleic acids and can be described as a right-handed helical ribbon wrapped around the helix axis, resulting in a large hollow core. Most intriguingly, neighboring base pairs slide strongly against each other, resulting in extensive interstrand base-base hydrophobic interactions along with unusual hydrophobic intrastrand interactions of nucleobases with their backbone. These results reveal how a minimal nucleic acid backbone can support highly stable Watson-Crick-like duplex formation.     

Isotropic periodic sum: a method for the calculation of long-range interactions

This work presents an accurate and efficient approach to the calculation of long-range interactions for molecular modeling and simulation. This method defines a local region for each particle and describes the remaining region as images of the local region statistically distributed in an isotropic and periodic way, which we call isotropic periodic images. Different from lattice sum methods that sum over discrete lattice images generated by periodic boundary conditions, this method sums over the isotropic periodic images to calculate long-range interactions, and is referred to as the isotropic periodic sum (IPS) method. The IPS method is not a lattice sum method and eliminates the need for a reciprocal space sum. Several analytic solutions of IPS for commonly used potentials are presented. It is demonstrated that the IPS method produces results very similar to that of Ewald summation, but with three major advantages, (1) it eliminates unwanted symmetry artifacts raised from periodic boundary conditions, (2) it can be applied to potentials of any functional form and to fully and partially homogenous systems as well as finite systems, and (3) it is more computationally efficient and can be easily parallelized for multiprocessor computers. Therefore, this method provides a general approach to an efficient calculation of long-range interactions for various kinds of molecular systems.