PAPQMD/AM1 Parametrization of the bonded term of aromatic biomolecules

The suitability of the PAPQMD/AM1 (Program for Approximate Parametrization from Quantum Mechanical Data/based on AM1 calculations) strategy to provide force-field parameters for large heteroaromatic compounds was studied. For this purpose, PAPQMD/AMI stretching and bending parameters for adenine, cytosine, thymine, guanine, and uracil were computed and compared with experimentally derived force-field parameters. Furthermore, equilibrium geometries and vibrational spectra for the five nucleic acid bases obtained from molecular mechanical calculations using both PAPQMD/AM1 and AMBER (Assisted Model Building with Energy Refinement) All Atoms force fields were compared with the available experimental data. © 1994 John Wiley & Sons, Inc.     

Molecular Parameterisation and the Theoretical Calculation of Electrode Potentials

The difference in reduction potentials between ortho and para-benzoquinones has been calculated. The employs gas phase ab initio and semi-empirical computations in combination with free energy perturbation theory applied to gas and solution phase Monte Carlo simulations. The effects on calculated results of altering solute electrostatic parameterisation in solution phase simulations is examined. Atom centred charges derived from the molecular electrostatic potentials, MEPs, from optimised ab initio wavefunctions and charges generated by consideration of hydrogen bonded complexes are considered. Parameterisation of hydroxyl torsions in hydroquinone molecules is treated in a physically realistic manner. The coupled torsional system of the ortho-hydrobenzoquinone molecule is described by a potential energy surface calculated using gas phase AM1 semi-empirical computations rather than the simple torsional energy functions frequently employed in such calculations. Calculated differences in electrode potentials show that the electrostatic interactions of quinone and hydroquinone molecules in aqueous solution are not well described by atom centred charges derived from ab initio calculated MEPs. Moreover, results in good agreement with the experimental reduction potential difference can be obtained by employing high level ab initio calculations and solution phase electrostatic parameters developed by consideration of hydrogen bonded complexes.     

Electrostatic interactions play an essential role in DNA repair and cold-adaptation of Uracil DNA glycosylase

Life has adapted to most environments on earth, including low and high temperature niches. The increased catalytic efficiency and thermoliability observed for enzymes from organisms living in constantly cold regions when compared to their mesophilic and thermophilic cousins are poorly understood at the molecular level. Uracil DNA glycosylase (UNG) from cod (cUNG) catalyzes removal of uracil from DNA with an increased k_cat and reduced K_m relative to its warm-active human (hUNG) counterpart. Specific issues related to DNA repair and substrate binding/recognition (K_m) are here investigated by continuum electrostatics calculations, MD simulations and free energy calculations. Continuum electrostatic calculations reveal that cUNG has surface potentials that are more complementary to the DNA potential at and around the catalytic site when compared to hUNG, indicating improved substrate binding. Comparative MD simulations combined with free energy calculations using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method show that large opposing energies are involved when forming the enzyme-substrate complexes. Furthermore, the binding free energies obtained reveal that the Michaelis-Menten complex is more stable for cUNG, primarily due to enhanced electrostatic properties, suggesting that energetic fine-tuning of electrostatics can be utilized for enzymatic temperature adaptation. Energy decomposition pinpoints the residual determinants responsible for this adaptation. Figure Electrostatic isosurfaces of cod uracil DNA glycosylase in complex with double stranded DNA     

A computational study of expanded heterocyclic nucleosides in DNA

The first molecular dynamics study of a series of heterospacer-expanded tricyclic bases in DNA using modified force field parameters in AMBER is detailed. The expanded purine nucleoside monomers have been designed to probe the effects of a heteroaromatic spacer ring on the structure, function, and dynamics of the DNA helix. The heterobase scaffold has been expanded with a furan, pyrrole, or thiophene spacer ring. This structural modification increases the polarizability of the bases and provides an additional hydrogen bond donor with the amine hydrogen of the pyrrole ring or hydrogen bond acceptor with the furan or thiophene ring free electron pairs. The polarizability of the expanded bases were determined by AM1 calculations and the results of the MD simulations of 20-mers predict that the modified curvature of the expanded base leads to a much larger major groove, while the effect on the minor groove is negligible. Overall, the structure resembles A-DNA. MD simulations of 10-mers suggest that the balance between base pairing vs. base stacking and intercalation can be shifted towards the latter due to the increased surface area and polarizability of the expanded bases.     

Modeling of Supramolecular Properties of Molecular Tweezers, Clips, and Bowls

The electrostatic potential surface (EPS) is calculated for molecular tweezers, clips, and bowls at different levels of theory (semiempirical AM1, ab initio HF/6-31G*, and density functional theory pBP/DN**). According to these calculations, the molecular electrostatic potential (MEP) on the concave side of the molecular tweezers and clips is suprisingly negative for hydrocarbons. This finding seems to be a general phenomenon in nonconjugated ?-electron systems with concave-convex topology and it explains the receptor properties of the molecular tweezers and clips. Analogous calculations performed for the conjugated aromatic molecular bowls show different results. The DFT calculations predict that in these systems the more negative MEP lies on the concave side similar to the findings for the nonconjugated molecular tweezer- and clip-systems, whereas the AM1 calculation leads to the opposite result that the MEP is more negative on convex side of the bowl-systems.     

The electrostatic potential of B-DNA

Electrostatic potentials around DNA are obtained by solving the nonlinear Poisson-Boltzmann (PB) equation. The detailed charge distribution of the DNA and the different polarizabilities of the macromolecule and solvent are included explicitly in the calculations. The PB equation is solved using extensions of a finite difference approach applied previously to proteins. Electrical potentials and ion concentrations are compared to those obtained with simpler models. It is found that the shape of the dielectric boundary between the macromolecule and solvent has significant effects on the calculated potentials near the surface, particularly in the grooves. Sequence-specific patterns are found, the most surprising result being the existence of positive regions of potential near the bases in both the major and minor grooves. The effect of solvent and ionic atmosphere screening of phosphate-phosphate repulsions is studied, and an effective dielectric function, appropriate for molecular mechanics simulations, is derived.     

The electrostatic molecular potential of tRNAPhe. IV. The potentials and steric accessibilities of sites associated with the bases.

The sites of the 76 nucleic acid bases of tRNAPhe potentially reactive towards electrophiles are studied by calculations on the associated molecular electrostatic potentials and the static steric accessibilities. Each of these sites is treated in its environment within the macromolecule. The influence of various schemes of screening by countercations of the backbone phosphates on the electrostatic potentials is investigated. The possible significance of the potentials and accessibilities in connection with observed chemical reactivities is discussed.     

Binding free energy calculations of adenosine deaminase inhibitors

The interactions between four inhibitors and adenosine deaminase (ADA) were examined by calculating their binding free energies after molecular dynamics simulations. A bonded model was used to represent the electrostatic potentials of the zinc coordination site. The charge distribution of the model was derived by using a two-stage electrostatic potential fitting calculations. The calculated binding free energies are in good agreement with the experimental data and the ranking of binding affinities is well reproduced. Notably, our findings suggest that non-polar contributions play an important role for ADA-inhibitor interactions.     

Modeling Lanthanide Complexes: Towards the Theoretical Design of Light Conversion Molecular Devices

Theoretical techniques have been developed and/or improved to predict the molecular structure of lanthanide complexes which were used to calculate their electronic properties, in particular, their electronic spectra and energy levels necessary to calculate the rates of energy transfer from the ligands to the metal ion. The molecular structure has been obtained by the SMLC/AM1 (Sparkle Model for the Calculation of Lanthanide Complexes – Austin Model 1) model where the lanthanide ion is simulated by a sparkle implemented into the AM1 Hamiltonian used to perform a HF-SCF (Hartree-Fock Self-Consistent Field) calculation. The previous implementation of the SMLC/AM1 model (sparkle/1) involving only two parameters has been generalized to be consistent with the AM1 Hamiltonian and the new model (sparkle/2) significantly improved the prediction of molecular structures of Eu(III) complexes. For the electronic spectra and energy level calculations of the lanthanide complexes the model replaces the metal ion by a point charge with the ligands held in their positions as determined by the SMLC/AM1 model, and uses a INDO/S-CI (intermediate neglect of differential overlap/spectroscopic-configuration interaction) model. A preliminary study of the solvent effects on the absorption spectra of the free ligand is also presented. For the ligand-lanthanide ion energy transfer Fermi's golden rule is used with the multipolar and exchange mechanisms being implemented and tested for several complexes. These theoretical techniques have been applied to several complexes yielding very good results when compared to experimental data as well as predictions for the molecular and electronic structures and the relative contributions of the mechanisms for the energy transfer rates. This revised version was published online in June 2006 with corrections to the Cover Date.     

A Mesoscale Model of DNA and Its Renaturation

A mesoscale model of DNA is presented (3SPN.1), extending the scheme previously developed by our group. Each nucleotide is mapped onto three interaction sites. Solvent is accounted for implicitly through a medium-effective dielectric constant and electrostatic interactions are treated at the level of Debye-Hückel theory. The force field includes a weak, solvent-induced attraction, which helps mediate the renaturation of DNA. Model parameterization is accomplished through replica exchange molecular dynamics simulations of short oligonucleotide sequences over a range of composition and chain length. The model describes the melting temperature of DNA as a function of composition as well as ionic strength, and is consistent with heat capacity profiles from experiments. The dependence of persistence length on ionic strength is also captured by the force field. The proposed model is used to examine the renaturation of DNA. It is found that a typical renaturation event occurs through a nucleation step, whereby an interplay between repulsive electrostatic interactions and colloidal-like attractions allows the system to undergo a series of rearrangements before complete molecular reassociation occurs.     

Electrostatic Potential of B-DNA: Effect of Interionic Correlations

Modified Poisson-Boltzmann (MPB) equations have been numerically solved to study ionic distributions and mean electrostatic potentials around a macromolecule of arbitrarily complex shape and charge distribution. Results for DNA are compared with those obtained by classical Poisson-Boltzmann (PB) calculations. The comparisons were made for 1:1 and 2:1 electrolytes at ionic strengths up to 1 M. It is found that ion-image charge interactions and interionic correlations, which are neglected by the PB equation, have relatively weak effects on the electrostatic potential at charged groups of the DNA. The PB equation predicts errors in the long-range electrostatic part of the free energy that are only ∼1.5 kJ/mol per nucleotide even in the case of an asymmetrical electrolyte. In contrast, the spatial correlations between ions drastically affect the electrostatic potential at significant separations from the macromolecule leading to a clearly predicted effect of charge overneutralization.     

Explicit-solvent molecular dynamics simulations of a DNA tetradecanucleotide duplex: lattice-sum versus reaction-field electrostatics

Five long-timescale (10 ns) explicit-solvent molecular dynamics simulations of a DNA tetradecanucleotide dimer are performed using the GROMOS 45A4 force field and the simple-point-charge water model, in order to investigate the effect of the treatment of long-range electrostatic interactions as well as of the box shape and size on the structure and dynamics of the molecule (starting from an idealised B-DNA conformation). Long-range electrostatic interactions are handled using either a lattice-sum (LS) method (particle–particle–particle–mesh; one simulation performed within a cubic box) or a cutoff-based reaction-field (RF) method (four simulations, with long-range cutoff distances of 1.4 or 2.0 nm and performed within cubic or truncated octahedral periodic boxes). The overall double-helical structure, including Watson–Crick (WC) base-pairing, is well conserved in the simulation employing the LS scheme. In contrast, the WC base-pairing is nearly completely disrupted in the four simulations employing the RF scheme. These four simulations result in highly distorted compact (cutoff distance of 1.4 nm) or extended (cutoff distance of 2 nm) structures, irrespective of the shape and size of the computational box. These differences observed between the two schemes seem correlated with large differences in the radial distribution function between charged entities (backbone phosphate groups and sodium counterions) within the system.     

Electrostatics of Polymorphic DNA

The molecular electrostatic potential (MEP) and the molecular electrostatic field (MEF) are associated with significantly different patterns of distribution in the nucleic acids and their constituents. In particular, a) while the values of the minimal potentials at the reactive sites of the bases or at the phosphates increase manyfold when going from the subunits to the double helix, the values of the field undergo only very small changes under the same circumstances and b) while the deepest potentials are located in the grooves of the double helix, the greatest fields are concentrated on the phosphates of the backbone. They are also influenced differently by such environmental factors as counterion screening: while the absolute values of the potentials are profoundly reduced, the fields are increased with respect to those of the unscreened acids. MEP and MEF also govern the electrostatics of interaction of DNA with different types of species. The MEP being of particular significance in this respect for interaction with cations and the MEF for the association with neutral dipolar molecules. A number of examples are given to illustrate the significance of this situation for different conformers of DNA.     

DNA minor groove electrostatic potential: influence of sequence-specific transitions of the torsion angle gamma and deoxyribose conformations

The structural adjustments of the sugar-phosphate DNA backbone (switching of the γ angle (O5′–C5′–C4′–C3′) from canonical to alternative conformations and/or C2′-endo → C3′-endo transition of deoxyribose) lead to the sequence-specific changes in accessible surface area of both polar and non-polar atoms of the grooves and the polar/hydrophobic profile of the latter ones. The distribution of the minor groove electrostatic potential is likely to be changing as a result of such conformational rearrangements in sugar-phosphate DNA backbone. Our analysis of the crystal structures of the short free DNA fragments and calculation of their electrostatic potentials allowed us to determine: (1) the number of classical and alternative γ angle conformations in the free B-DNA; (2) changes in the minor groove electrostatic potential, depending on the conformation of the sugar-phosphate DNA backbone; (3) the effect of the DNA sequence on the minor groove electrostatic potential. We have demonstrated that the structural adjustments of the DNA double helix (the conformations of the sugar-phosphate backbone and the minor groove dimensions) induce changes in the distribution of the minor groove electrostatic potential and are sequence-specific. Therefore, these features of the minor groove sizes and distribution of minor groove electrostatic potential can be used as a signal for recognition of the target DNA sequence by protein in the implementation of the indirect readout mechanism.     

Role of structure and dynamics of DNA with cisplatin and oxaliplatin adducts in various sequence contexts on binding of HMGB1a

Anti-cancer drugs, such as cisplatin and oxaliplatin, covalently bind to adjacent guanine bases in DNA to form intra-strand adducts. Differential recognition of drug–DNA adducts by the protein HMGB1a has been related to the differences in efficacy of these drugs in tumours. Additionally, the bases flanking the adduct (the sequence context) also have a marked effect on HMGB1a binding affinity. We perform atomistic molecular dynamics simulations of DNA with cisplatin and oxaliplatin adducts in four sequence contexts (AGGC, CGGA, TGGA and TGGT) in the absence and presence of HMGB1a. The structure of HMGB1a-bound drug–DNA molecules is independent of sequence and drug identity, confirming that differential recognition cannot be explained by the protein-bound structure. The differences in the static and conformational dynamics of the drug–DNA structures in the absence of the protein explain some but not all trends in differential binding affinity of HMGB1a. Since the minor groove width and helical bend of all drug–DNA molecules in the unbound state are lower than the protein-bound state, HMGB1a must actively deform the DNA during binding. The thermodynamic pathway between the unbound and protein-bound states could be an additional factor in the binding affinity of HMGB1a for drug–DNA adducts in various sequence contexts.     

Molecular similarity analysis between insect juvenile hormone and N, N-diethyl-m-toluamide (DEET) analogs may aid design of novel insect repellents

Molecular similarity analysis of stereoelectronic properties between natural insect juvenile hormone (JH), -a synthetic insect juvenile hormone mimic (JH-mimic, undecen-2-yl carbamate), and N, N-diethyl-m-toluamide (DEET) and its analogs reveals similarities that may aid the design of more efficacious insect repellents and give a better insight into the mechanism of repellent action. The study involves quantum chemical calculations using the AM1 semi-empirical computational method enabling a conformational search for the lowest and most abundant energy conformers of JH, JH-mimic, and 15 DEET compounds, followed by complete geometry optimization of the conformers. Similarity analyses of stereoelectronic properties such as structural parameters, atomic charges, dipole moments, molecular electrostatic potentials, and highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies were performed on JH, JH-mimic and the DEET compounds. The similarity of stereoelectronic attributes of the amide/ester moiety, the negative electrostatic potential regions beyond the van der Waals surface, and the large distribution of hydrophobic regions in the compounds appear to be the three important factors leading to a similar interaction with the JH receptor. The similarity of electrostatic profiles beyond the van der Waals surface is likely to play a crucial role in molecular recognition interaction with the JH receptor from a distance. This also suggests electrostatic bioisosterism of the amide group of the DEET compounds and JH-mimic and, thus, a model for molecular recognition at the JH receptor. The insect repellent property of the DEET analogs may thus be attributed to a conflict of complementarity for the JH receptor binding sites.     

Revisiting the planarity of nucleic acid bases: Pyramidilization at glycosidic nitrogen in purine bases is modulated by orientation of glycosidic torsion

We describe a novel, fundamental property of nucleobase structure, namely, pyramidilization at the N1/9 sites of purine and pyrimidine bases. Through a combined analyses of ultra-high-resolution X-ray structures of both oligonucleotides extracted from the Nucleic Acid Database and isolated nucleotides and nucleosides from the Cambridge Structural Database, together with a series of quantum chemical calculations, molecular dynamics (MD) simulations, and published solution nuclear magnetic resonance (NMR) data, we show that pyramidilization at the glycosidic nitrogen is an intrinsic property. This property is common to isolated nucleosides and nucleotides as well as oligonucleotides—it is also common to both RNA and DNA. Our analysis suggests that pyramidilization at N1/9 sites depends in a systematic way on the local structure of the nucleoside. Of note, the pyramidilization undergoes stereo-inversion upon reorientation of the glycosidic bond. The extent of the pyramidilization is further modulated by the conformation of the sugar ring. The observed pyramidilization is more pronounced for purine bases, while for pyrimidines it is negligible. We discuss how the assumption of nucleic acid base planarity can lead to systematic errors in determining the conformation of nucleotides from experimental data and from unconstrained MD simulations.     

Electrostatics of Polymorphic DNA

Abstract

The molecular electrostatic potential (MEP) and the molecular electrostatic field (MEF) are associated with significantly different patterns of distribution in the nucleic acids and their constituents. In particular, a) while the values of the minimal potentials at the reactive sites of the bases or at the phosphates increase manyfold when going from the subunits to the double helix, the values of the field undergo only very small changes under the same circumstances and b) while the deepest potentials are located in the grooves of the double helix, the greatest fields are concentrated on the phosphates of the backbone. They are also influenced differently by such environmental factors as counterion screening: while the absolute values of the potentials are profoundly reduced, the fields are increased with respect to those of the unscreened acids. MEP and MEF also govern the electrostatics of interaction of DNA with different types of species. The MEP being of particular significance in this respect for interaction with cations and the MEF for the association with neutral dipolar molecules. A number of examples are given to illustrate the significance of this situation for different conformers of DNA.