Simulating structural and thermodynamic properties of carcinogen-damaged DNA

A pair of stereoisomeric covalent adducts to guanine in double-stranded DNA, derived from the reaction of mutagenic and tumorigenic metabolites of benzo[a]pyrene, have been well characterized structurally and thermodynamically. Both high-resolution NMR solution structures and an array of thermodynamic data are available for these 10S (+)- and 10R (-)-trans-anti -[BP]-N(2)-dG adducts in double-stranded deoxyoligonucleotides. The availability of experimentally well-characterized duplexes containing these two stereoisomeric guanine adducts provides an opportunity for evaluating the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method for computing thermodynamic properties from molecular dynamics ensembles. We have carried out 3-ns molecular dynamics simulations, using NMR solution structures as the starting models for the 10S (+)- and 10R (-)-trans-anti-dG adducts in a DNA duplex 11-mer using AMBER 6.0. We employed the MM-PBSA method to compute the free energies, enthalpies, and entropies of the two adducts. Our complete thermodynamic analysis agrees quite well with the full experimental thermodynamic characterization of these adducts, showing essentially equal stabilities of the two adducts. We also calculated the nuclear Overhauser effect (NOE) distances from the molecular dynamics trajectories, and compared them against the experimental NMR-derived NOE distances. Our results showed that the simulated structures are in good agreement with the NMR experimental NOE data. Furthermore, the molecular dynamics simulations provided new structural and biological insights. Specifically, the puzzling observation that the BP aromatic ring system in the 10S (+)-trans-anti-dG adduct is more exposed to the aqueous solvent than the 10R (-)-trans-anti-dG adduct, is rationalized in terms of the adduct structures. The structural and thermodynamic features of these stereoisomeric adducts are also discussed in relation to their reported low susceptibilities to nucleotide excision repair.     

Structural requirements of the fructan-lipid interaction

Fructans are a group of fructose-based oligo- and polysaccharides. They are proposed to be involved in membrane protection of plants during dehydration. In accordance with this hypothesis, they show an interaction with hydrated lipid model systems. However, the structural requirements for this interaction are not known both with respect to the fructans as to the lipids. To get insight into this matter, the interaction of several inulins and levan with lipids was investigated using a monomolecular lipid system or the MC 540 probe in a bilayer system. MD was used to get conformational information concerning the polysaccharides. It was found that levan-type fructan interacted comparably with model membranes composed of glyco- or phospholipids but showed a preference for lipids with a small headgroup. Furthermore, it was found that there was an inulin chain-length-dependent interaction with lipids. The results also suggested that inulin-type fructan had a more profound interaction with the membrane than levan-type fructan. MD simulations indicated that the favorable conformation for levan is a helix, whereas inulin tends to form random coil structures. This suggests that flexibility is an important determinant for the fructan-lipid interaction.     

Sequence context dependence of tandem guanine:adenine mismatch conformations in RNA: a continuum solvent analysis

Guanine:adenine (G:A) mismatches and in particular tandem G:A (tG:A) mismatches are frequently observed in biological RNA molecules and can serve as sites for tertiary interaction, metal binding and protein recognition. Depending on the surrounding sequence tG:A mismatches can adopt different basepairing topologies. In the sequence context (5'-) GGAC (tandem G:A in bold) a face-to-face (imino or Watson-Crick-like) pairing is preferred whereas in the CGAG context, G and A adopt a sheared arrangement. Systematic conformational searches with a generalized Born continuum model and molecular dynamics simulations including explicit water molecules and ions have been used to generate face-to-face and sheared tG:A mismatches in both CGAG and GGAC sequence contexts. Conformations from both approaches were evaluated using the same force field and a Poisson-Boltzmann continuum solvent model. Although the substate analysis predicted the sheared arrangement to be energetically preferred in both sequence contexts, a significantly greater preference of the sheared form was found for the CGAG context. In agreement with the experimental observation, the analysis of molecular dynamics trajectories indicated a preference of the sheared form in the case of the CGAG-context and a favorization of the face-to-face form in the case of the GGAC context. The computational studies allowed to identify energetic contributions that stabilize or destabilize the face-to-face and sheared tandem mismatch topologies. The calculated nonpolar solvation and Lennard-Jones packing interaction were found to stabilize the sheared topology independent of the sequence context. Electrostatic contributions are predicted to make the most significant contribution to the sequence context dependence on the structural preference of tG:A mismatches.     

Structure-function relationship studies of bovine parathyroid hormone [bPTH(1-34)] analogues containing alpha-amino-iso-butyric acid (Aib) residues

The N-terminal 1-34 fragments of the parathyroid hormone (PTH) and parathyroid hormone-related protein (PTHrP) elicit the full spectrum of bone-related biological activities of the intact native sequences. It has been suggested that the structural elements essential for bioactivity are two helical segments located at the N-terminal and C-terminal sequences, connected by hinges or flexible points around positions 12 and 19. In order to assess the relevance of the local conformation around Gly(12) upon biological function, we synthesized and characterized the following PTH(1-34) analogues containing Aib residues: (I) A-V-S-E-I-Q-F-nL-H-N-Aib-G-K-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-Y-NH(2) ([Nle(8,18), Aib(11), Nal(23),Tyr(34)]bPTH(1-34)-NH(2)); (II) A-V-S-E-I-Q-F-nL-H-N-L-Aib-K-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-Y-NH(2) ([Nle(8,18), Aib(12),Nal(23),Tyr(34)]bPTH(1-34)-NH(2)); (III) A-V-S-E-I-Q-F-nL-H-N-L-G-Aib-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-Y-NH(2) ([Nle(8,18), Aib(13), Nal(23),Tyr(34)]bPTH(1-34)-NH(2)); (IV) A-V-S-E-I-Q-F-nL-H-N-Aib-Aib-K-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-YNH(2) ([Nle(8,18), Aib(11,12), Nal(23),Tyr(34)]bPTH(1-34)-NH(2)); (V) A-V-S-E-I-Q-F-nL-H-N-L-Aib-Aib-H-L-S-S-nL-E-R-V-E-Nal-L-R-K-K-L-Q-D-V-H-N-Y-NH(2) ([Nle(8,18), Aib(12,13),Nal(23),Tyr(34)]bPTH(1-34)-NH(2)). (nL= Nle; Nal= L-(2-naphthyl)-alanine; Aib= alpha-amino-isobutyric acid.) The introduction of Aib residues at position 11 in analogue I or at positions 11 and 12 in analogue IV resulted in a 5-20-fold lower efficacy and a substantial loss of binding affinity compared to the parent compound [Nle(8,18), Nal(23),Tyr(34)]bPTH(1-34)-NH(2). Both binding affinity and adenylyl cyclase stimulation activity are largely restored when the Aib residues are introduced at position 12 in analogue II, 13 in analogue III, and 12-13 in analogue V. The conformational properties of the analogues in aqueous solution containing dodecylphosphocholine micelles were studied by CD, two-dimensional (2D) NMR and computer simulations. The results indicated the presence of two helical segments in all analogues, located at the N-terminal and C-terminal sequences. Insertion of Aib residues at positions 12 and 13, or of Aib dyads at positions 11-12 and 12-13, enhances the stability of the N-terminal helix of all analogues. In all analogues the Aib residues are included in the helical segments. These results confirmed the importance of the helical structure in the N-terminal activation domain, as well as of the presence of the Leu(11) hydrophobic side chain in the native sequence, for PTH-like bioactivity.     

Molecular dynamics simulations of A. T-rich oligomers: sequence-specific binding of Na+ in the minor groove of B-DNA

Molecular dynamics simulations have been employed to probe the sequence-specific binding of sodium ions to the minor groove of B-DNA of three A. T-rich oligomers having identical compositions but different orders of the base pairs: C(AT)(4)G, CA(4)T(4)G, and CT(4)A(4)G. Recent experimental investigations, either in crystals or in solution, have shown that monovalent cations bind to DNA in a sequence-specific mode, preferentially in the narrow minor groove regions of uninterrupted sequences of four or more adenines (A-tracts), replacing a water molecule of the ordered hydration structure, the hydration spine. Following this evidence, it has been hypothesized that in A-tracts these events may be responsible for structural peculiarities such as a narrow minor groove and a curvature of the helix axis. The present simulations confirm a sequence specificity of the binding of sodium ions: Na(+) intrusions in the first layer of hydration of the minor groove, with long residence times, up to approximately 3 ns, are observed only in the minor groove of A-tracts but not in the alternating sequence. The effects of these intrusions on the structure of DNA depend on the ion coordination: when the ion replaces a water molecule of the spine, the minor groove becomes narrower. Ion intrusions may also disrupt the hydration spine modifying the oligomer structure to a large extent. However, in no case intrusions were observed to locally bend the axis toward the minor groove. The simulations also show that ions may reside for long time periods in the second layer of hydration, particularly in the wider regions of the groove, often leading to an opening of the groove.     

Conformational analysis of endomorphin-2 by molecular dynamics methods

Endomorphin-2 (EM2, H-Tyr-Pro-Phe-Phe-NH(2)) is a highly potent and selective mu-opioid receptor agonist. A conformational analysis of EM2 was carried out by simulated annealing (SA) and molecular dynamics (MD) methods. Molecular modeling was conducted on both neutral (N-terminal NH(2)) and charged (N-terminal NH(3) (+)) molecules. Based on the results of NMR investigations showing an equilibrium mixture of cis and trans Tyr(1)-Pro(2) peptide bonds for EM2 in solution, simulations were performed with restrained cis-Pro and trans-Pro peptide bonds, too. A separate SA study with unrestrained Pro peptide bonds was also conducted. Preferred conformational states are presented in Ramachandran plots. The g(+), g(-), and trans populations of the aromatic amino acid residue side chains were determined in chi(1) space. The distances between the N-terminal N atom and the other backbone N and O atoms, and the distances between the centers of the aromatic rings and the Pro(2) ring, were determined. The energy distribution of the structures obtained by SA was calculated. The preferred secondary structural elements were different kinds of beta-turns, an inverse gamma-turn located in the N-terminal region, and regular and inverse gamma-turns located in the C-terminal region. These turns were stabilized by intramolecular H-bonds and bifurcated H-bonds.     

Molecular dynamics simulation studies of cytochrome b5 from outer mitochondrial and microsomal membrane

Two forms of cytochrome b(5) have been identified, associated with the outer membrane of liver mitochondria (OM cyt b(5)) and with the membrane of the endoplasmic reticulum (microsomal, Mc cyt b(5)). These proteins have very similar structures, but differ significantly in physical properties, with the OM cyt b(5) exhibiting a more negative reduction potential, higher stability, and stronger interactions with the heme. We perform molecular dynamics simulations to probe the structures and fluctuations of the two proteins in solution, to help explain the observed physical differences. We find that the structures of the two proteins, highly similar in the crystal, differ in position of a surface loop involving residues 49-51 in solution. Hydrophobic residues Ala-18, Ile-32, Leu-36, and Leu-47 tend to cluster together on the surface of rat OM cyt b(5), blocking water access to the protein interior. In bovine Mc cyt b(5), two of these positions, Ser-18 and Arg-47, are occupied by hydrophilic residues. This leads to breaking the hydrophobic cluster and allowing the protein to occupy a more open conformation. A measure of this structural transition is the opening of a cleft on the protein surface, which is 5 A wider in the OM cyt b(5) simulation compared to the Mc form. The OM protein also appears to have a more compact hydrophobic core in its beta-sheet region. These effects may be used to explain observed stability differences between the two proteins.     

Brownian dynamics simulation of helix-capping motifs

Helix-capping motifs are believed to play an important role in stabilizing alpha-helices and defining helix start and stop signals. We performed microsecond scale Brownian dynamics simulations to study ten XAAD sequences, with X = (A,E,I,L,N,Q,S,T,V,Y), to examine their propensity to form helix capping motifs and correlate these results with those obtained from analyzing a structural database of proteins. For the widely studied capping box motif S**D, where the asterisk can be any amino acid residue, the simulations suggested that one of the two hydrogen bonds proposed earlier as a stabilizing factor might not be as important. On the other hand, side-chain interactions between the capping residue and the third residue downstream on the polypeptide chain might also play a role in stabilizing this motif. These results are consistent with explicit-solvent molecular dynamics simulations of two capping box motifs found in the proteins BPTI and alpha-dendrotoxin. Principal component analysis of the SAAD trajectory showed that the first three principal components, after those corresponding to translational-rotational motion were removed, accounted for more than half of the conformational fluctuations. The first component separated the conformational space into two parts with the all-helical conformation and the capping box motif lying largely in one part. The second component, on the other hand, could be used to describe conformational transitions between the all-helical form and the capping box motif.     

Interaction of pollen-specific actin-depolymerizing factor with actin

We have examined the interaction of recombinant lily pollen ADF, LlADF1, with actin and found that whilst it bound both G- and F-actin, it had a much smaller effect on the polymerization and depolymerization rate constants than the maize vegetative ADF, ZmADF3. An antiserum specific to pollen ADF, antipADF, was raised and used to localize pollen ADF in daffodil--a plant in which massive reorganizations of the actin cytoskeleton have been seen to occur as pollen enters and exits dormancy. We show, for the first time, an ADF decorating F-actin in cells that did not result from artificial increase in ADF concentration. In dehydrated pollen this ADF : actin array is replaced by actin : ADF rodlets and aggregates of actin, which presumably act as a storage form of actin during dormancy. In germinated pollen ADF has no specific localization, except when an adhesion is made at the tip where actin and ADF now co-localize. These activities of pollen ADF are discussed with reference to the activities of ZmADF3 and other members of the ADF/cofilin group of proteins.