Search

Distributed modeling approach may have much better performance and accuracy compared with lumped-parameter hydrologic models. The main goals of this research are: investigating the possibility of combining distributed hydrological models with an one-dimensional hydraulic model and simulating waterways in large watersheds with limited hydrological and hydraulic data. Then performing sensitivity analysis on different parameters in order to identify the parameters containing the major influences on results. In the current research, an innovative approach in Gridded Surface Subsurface Hydrologic Analysis (GSSHA) model, the cross-sections of all 414 waterways in the 3450 km2 Karvandar watershed, used for flow routing calculations, are uniquely extracted. Then, the effect of three essential factors are evaluated. These factors are accuracy of the digital topographic model, cell size of grid network, and density of streams, on the results of GSSHA model simulations. This watershed is located in southeastern Iran, has a dry climate with limited available hydrological data. Results showed that peak discharges obtained from the GSSHA model, developed based on a DEM with a spatial resolution of 12.5 m, are slightly (< 4%) lower than the corresponding values ​​in the GSSHA model with a 30 m DEM resolution. This fact confirms that the use of the topographic model with a lower spatial resolution has no substantial effects on the accuracy of simulation. Also, the peak discharges increased significantly (44% to 57%) by increasing the density of waterways in the GSSHA model. Furthermore, results showed that peak discharge obtained from three models with grid cell sizes of 100, 150, and 200 m (base model), are close together. Comparing with two models of coarser grids (250 and 300 m), significant differences observed, which indicated that the grids larger than 200 m could induce substantial errors in results. [ABSTRACT FROM AUTHOR]

In this work, the role of Lewis acid–base (LAB) interaction on the stereoselectivity of the Diels–Alder (DA) reaction has been studied by DFT in gas and solution (dichloromethane) phases. The calculations were performed at the B3LYP/6-311G++ (d, p) level. Two different series of DA reactions were investigated: (1)—three mono-substituted cyclopentadienes + dimethyl(vinyl)borane; (2)—five α,β-unsaturated carbonyl compounds + cyclopenta-2,4-dien-1-yldimethylborane. The reacting diene and dienophile pairs were chosen to restrict LAB interaction to the exo reaction pathway. It was found that in some of the examined cases, the favorable LAB interaction is so strong that it can lead to a completely exo-selective DA reaction. Furthermore, a novel multistep synthetic method was hypothesized for preparing exo cycloadduct with near 100% stereoselectivity. Our results can open up new avenues toward the rational design of exo-selective DA reactions for synthesizing novel bioorganic compounds. [ABSTRACT FROM AUTHOR]

Designing drug delivery systems for therapeutic compounds whose receptors are located in the cytosol of cells is challenging as a bilayer cell membrane is negatively charged. The newly designed drug delivery systems should assist the mentioned drugs in passing the membrane barriers and achieving their targets. This study concentrated on developing novel ionic liquids (ILs) that interact effectively with cell membranes. These ILs are based on glucose-containing choline and are expected to be non-toxic. The binding energies of the known pharmaceutically active ionic liquids were calculated at the B3LYP/6-311++G(d,p) level in the gas phase and compared with those of our newly designed carbohydrate-based ionic liquids. Subsequently, we employed MD simulations to obtain information about the interactions of these known and designed ILs with the cell membrane. In our approach, we adopted QM and MD studies and illustrated that there could be a link between the QM and MD results. [ABSTRACT FROM AUTHOR]

DNA polymerase η (polη) is of significant value for designing new families of anticancer drugs. This protein takes a role in many stages of the cell cycle, including DNA replication, translesion DNA synthesis, and the repairing process of DNA. According to many studies, a high level of expression of polη in most cases has been associated with low rates of patients' survival, regardless of considering the stage of tumor cells. Thus, the design of new drugs with fewer side effects to inhibit polη in cancerous cells has attracted attention in recent years. This project aims to design and explore the alternative inhibitors for polη, which are based on carbohydrates and amino acids. In terms of physicochemical properties, they are similar to the traditional anticancer drugs such as Cytarabine (cytosine arabinose). These alternative inhibitors are supposed to disrupt the DNA replication process in cancerous cells and prevent the tumor cells from mitosis. These newly designed structures, which are based on natural products, are expected to be non-toxic and to have the same chemotherapeutic impact as the traditional agents. The combinatorial use of quantum mechanics studies and molecular dynamic simulation has enabled us to precisely predict the inhibition mechanism of the newly designed structure, which is based on carbohydrates and amino acids, and compare it with that of the traditional chemotherapeutic drugs such as Cytarabine. Our results suggest that the inhibitors containing the natural building blocks of amino acid and carbohydrate could be considered alternative drugs for Cytarabine to block polη. [ABSTRACT FROM AUTHOR]

We have investigated the pathways of stability for NHCs derived from azole, di-azole, n-tetrazole, and ab-tetrazole (1a, 2a, 3a, and 4a, respectively), at the M06/6-311++G** level of theory. Optimization and vibrational frequency calculations of ground states (GS) and transition states (TS) are performed to identify Gibbs free energies and nature of stationary points, respectively. Two possible pathways of stability for 1a-4a are compared and contrasted which entail dimerization through hydrogen bonding (HB) and covalent bonding (CB). The CB pathway comprises head to head (HH) and head to tail (HT) dimerizations. Plausible reaction profiles are illustrated for 1a-4a along with the mechanism of each dimerization. Structures 1a-3a show one possibility for HB while 4a represents two possibilities. Structures 1a and 4a display two HH dimers while 2a and 3a show one. Structures 1a-4a undergo HT dimerizations to yield three possible dimers which include trans, cis, and [2+3] isomers. Interestingly, for all 1a-4a, HB dimerization turns out as the most favorable stability pathway for showing no barrier of reaction. Structures 4b and 4c indicate the highest stability with respect to their initial 4a compared to remaining HB dimers 1b-3b. In addition, the 1,2-H shift appears as a possible rearrangement for 1a-4a to yield their corresponding tautomers (1i, 2h, 3h, and 4k, respectively). The reaction profile of this rearrangement indicates that 1a-4a favor HB dimerization pathway more than 1,2-H shift, in terms of kinetic and thermodynamic. [ABSTRACT FROM AUTHOR]

In recent years, there has been increasing interest in damaged DNA and RNA nucleobases. These damaged nucleobases can cause DNA mutation, resulting in various diseases such as cancer. Alkylating agents are mutagenic and carcinogenic in a variety of prokaryotic and eukaryotic organisms. The present study employs density functional theory (DFT/B3LYP) with the 6-311++G(d,p) basis set to investigate the effect of chemical damage in O-alkyl pyrimidines such as O-methylthymine, O-methylcytosine and O-methylthymine. We compared the intrinsic properties, such as proton affinities, gas phase acidities, equilibrium tautomerization and nucleobase pair's hydrogen bonding properties, of these molecules with those in the normal nucleobases thymine and cytosine. The results are of interest for chemical reasons and also possibly for biological purposes since biological media can be quite non-polar. Furthermore, we found that N1-H of O-methylthymine is less acidic than N1-H of thymine, suggesting that alkyl DNA glycosylase enzyme cannot discriminate this damaged nucleobase from a normal thymine nucleobase. This result indicates that the conjugated base anion of O-methylthymine would be a worse leaving group and O-methylthymine is repaired in genome by demethylation rather than enzyme-catalyzed excision at N1. [ABSTRACT FROM AUTHOR]

In the present density functional theory study, we have compared intrinsic properties of non-natural nucleobases ( acA, acG, acC, and acT nucleobases) such as proton affinities, gas phase acidities, tautomerization, and hydrogen-bonding properties with those in normal Watson-Crick nucleobases (A, G, C, T nucleobases). The hydrogen-bonding interactions in non-natural and Watson-Crick base pairs were studied at B3LYP/6-311++G (d,p) level regarding their geometries, energies, and topological features of the electron density. The quantum theory of atoms-in-molecule (QTAIM) and natural bond orbital (NBO) analyses were employed to elucidate the interaction characteristics in base pairs. The electron density ρ( r) as well as its Laplacian $$ \nabla^{2} $$ ρ( r) at the hydrogen bond critical point predicted by QTAIM is strongly correlated with hydrogen bond structural parameter and the second-order perturbation energies in NBO scheme. Our results show that most of hydrogen bonds in normal Watson-Crick and non-natural base pairs must be considered as electrostatic interactions. Results of calculations revealed that energetic values of hydrogen bonds in T- A base pair are more than those in ac T- ac A base pair, while values of hydrogen bonds in C- G base pair and ac C- ac G base pair are almost the same. These results confirmed stability order of stabilization energies of these base pairs. [ABSTRACT FROM AUTHOR]

Glutathione tripeptide (γ-glutamyl-cysteinyl-glycine) is a flexible molecule and its conformational energy landscape is strongly influenced by forming intramolecular hydrogen bond, its charge and the environment. This study employs DFT-B3LYP method with the 6-31+G (d,p) basis set to carry out conformational analysis of neutral, zwitterionic, cationic, and anionic forms of glutathione. In analyzing the structural characteristics of these structures, intramolecular hydrogen bonds were identified and characterized in details by topological parameters such as electron density ρ( r) and Laplacian of electron density $$ \nabla^{2} $$ρ( r) from Bader's atom in molecules theory. Charge transfer energies based on natural bond orbital analysis are also considered to interpret these intramolecular hydrogen bonds. Our results show that these hydrogen bonds are partially electrostatic and partially covalent in nature, in which the covalent contribution increases as the stabilization energy of hydrogen bond increases. Furthermore, the back bone and side chain (Ramachandran map) orientations of various ionic forms of glutathione have been studied and conformation of each constitution of glutathione tripeptide (i.e., Glu, Cys, and Gly moieties) was determined. In most species side chain conformation were found to be hindered gauche-gauche orientation by intramolecular hydrogen bonds. [ABSTRACT FROM AUTHOR]

Thiotepa ( N,N′, N″-triethylenethiophosphoramide) and its major metabolite (Tepa) as trifunctional alkylating agents has recently been used in cancer therapy. In vivo and vitro studies show the possible pathways of alkylation of DNA by Thiotepa and Tepa. Two pathways are suggested, but the main pathway of mechanism remains unclear. In pathway 1, forming cross-links with DNA molecules can be carried out via two different mechanisms. In first mechanism, these agents undergo the ring opening reaction which is initiated by protonating aziridine, which then becomes the main target of nucleophilic attack by the N7-Guanine of DNA. The second probable mechanism is ring opening of aziridyl group by nucleophilic attack of N7-Guanine without initial protonation. Thiotepa and Tepa in pathway 2, act as a cell penetrating carrier for aziridine, which is released via hydrolysis. The released aziridine can form a cross-link with N7-Guanine. In this study, we calculated the activation free energy and kinetic rate constant for alkylating the Guanine via the first pathway to determine the most precise mechanism by applying density functional theory using B3LYP method. We carried out geometrical optimizations with the conductor-like polarizable continuum model to account for the solvent effect, and the results were compared with those in the gas phase. Hyperconjugation stabilization factors that affect on stability of generated transition state were investigated by natural bond order analysis. Furthermore, quantum theory of atoms in molecules analysis was performed to extract the bond critical points properties because the electron densities can be considered as a good description of the strength of different types of interactions. [ABSTRACT FROM AUTHOR]

Etheno-damaged DNA adducts such as 3, N-ethenocytosine, N,3-ethenoguanine, and 1, N-ethenoguanine are associated with carcinogenesis and cell death. These inevitable damages are counteracted by glycosylase enzymes, which cleave damaged nucleobases from DNA. Escherichia coli alkyl purine DNA glycosylase is the enzyme responsible for excising damaged etheno adducts from DNA in humans. In an effort to understand the intrinsic properties of these molecules, we examined gas-phase acidity values and proton affinities (PA) of multiple sites of these molecules as well as equilibrium tautomerization and base pairing properties by quantum mechanical calculations. We also used calculations to compare the acidic and basic properties of these etheno adduct with those of the normal bases-cytosine and guanine nucleobases. We hypothesize that alkyl DNA glycosylase may cleave certain damaged nucleobases as anions and that the active site may take advantage of a nonpolar environment to favor deprotonated cytosine or guanine as a leaving group versus damaged nucleobases. [ABSTRACT FROM AUTHOR]

The relative stabilities and noncovalent interactions of the six low-lying energy tautomers of cytosine nucleobase with some biological anions (such as F, Cl, and CN) have been investigated in gas phase by density functional theory (DFT) method in conjunction with 6-311++G (d,p) atomic basis set. Furthermore, to systematically investigate all possible tautomerisms from cytosine induced by proton transfer, we describe a study of structural tautomer interconversion in the gas phase and in a continuum solvent using DFT calculation. We carried out geometrical optimizations with the integral equation formalism of polarizable continuum (IEF-PCM) model to account for the solvent effect, and the results were compared with those in the gas phase. The result of calculation revealed that anions bind mostly in a bidentate manner via hydrogen bond, and stabilization energies of these complexes are larger than those in the case when anions bind in a monodentate manner. The quantum theory of atoms in molecules (QTAIM), natural bonding orbital (NBO) and energy decomposition analysis (EDA) have also been applied to understand the nature of hydrogen bond interactions in these complexes. NBO analysis reveals that the interaction patterns between the anions and the tautomers are σ-type interaction between lone pairs and $$ \sigma_{{_{{{\text{N}}--{\text{H}}}} }}^{*} $$, $$ \sigma_{{_{{{\text{O}}--{\text{H}}}} }}^{*} $$ and $$ \sigma_{{_{{{\text{H}}--{\text{F}}}} }}^{*} $$ antibonding orbitals. Also, according to these theories, the interactions are found to be partially electrostatic and partially covalent. EDA results identify that these bonds have less than 35% covalent character and more than 65% electrostatic, and the covalent character increases in different anions in the order F > CN > Cl. On the other hand, orbital interaction energies of complexes of F anion are more than those of Cl and CN complexes. The lower orbital interaction energies in complexes of Cl and CNanions imply less charge transfer and stronger ionic bond character. Furthermore, relationship between the orbital interaction energy ( E) with hydrogen bonding energy ( E) and the electron density (ρ( r)) with hydrogen bonding energy of F, Cl and CN complexes have also been investigated. [ABSTRACT FROM AUTHOR]

N,N′,N″-triethylenethiophosphoramide (Thiotepa) and its oxo analogue (Tepa) as the major metabolite are trifunctional alkylating agents with a broad spectrum of antitumor activity. In vivo and vitro studies show alkylation of DNA by Thiotepa and Tepa can follow two pathways, but it remains unclear which pathway represents the precise mechanism of action. In pathway 1, these agents are capable of forming cross-links with DNA molecules via two different mechanisms. In the first mechanism, the ring opening reaction is initiated by protonating the aziridine, which then becomes the primary target of nucleophilic attack by the N7-Guanine. The second one is a direct nucleophilic ring opening of aziridyl group. Thiotepa and Tepa in pathway 2, act as a cell penetrating carrier for aziridine, which is released via hydrolysis. The released aziridine can form a cross-link with N7-Guanine. In this study, we calculated the activation free energy and kinetic rate constant for hydrolysis of these agents and explored interaction of aziridine with Guanine to predict the most probable mechanism by applying density functional theory (DFT) using B3LYP method. In addition, solvent effect was introduced using the conductor-like polarizable continuum model (CPCM) in water, THF and diethylether. Hyperconjugation stabilization factors that have an effect on stability of generated transition state were investigated by natural bond order (NBO) analysis. Furthermore, quantum theory of atoms in molecules (QTAIM) analysis was performed to extract the bond critical points (BCP) properties, because the electron densities can be considered as a good description of the strength of different types of interactions. [ABSTRACT FROM AUTHOR]

The cation radicals of DNA constituents generated by the ionizing radiation initiate the alteration of the bases, which is one main type of cytotoxic DNA lesions. These cation radical spices are known for their role in producing nucleic acid strand break, and it is important to identify the cation radical formation at particular atomic site in these molecules so that the major pathway for the nucleic acid damage may be trapped. In the present study, we explored theoretically energetic, structural, and electronic properties of the possible radicals formed via proton atom abstraction at various sites of sugar part of deoxycytidine cation radical by employing density functional theory at B3LYP/6-311++G (d,p) level. The computation revealed 0.0-22.6 kcal/mol energy disparity in these radicals. Radical-centered carbon increases the extent of bonding with its adjacent atoms. This tendency should be important in predicting the reactivity of sugar-based radicals. Based on DFT calculations, sugar radicals of deoxycytidine have following stability order: raH1′ > raH2′ > raH4′ > raH3′ > raH5′ > raO5′H > raO3′H. Furthermore, influence of cation radical formation on acidities of multiple sites in deoxycytidine nucleosides was investigated. For instance upon cation radical formation, ΔH of O3′H and O5′H sites of deoxycytidine varies from 348.6 and 351.5 to 228.8 and 227.5 kcal/mol, respectively. [ABSTRACT FROM AUTHOR]

Influence of the addition of water molecules ( n = 1-6) on the interaction energy between Li, Na, K cations and indole molecule as tryptophan amino acid residue is considered at MP2(FULL)/6-311++G(d,p)//B3LYP/6-311++G(d,p) levels of theory. The calculations suggest that the size of cation and the number of water molecules are two important factors that affect the interaction energy between the hydrated metal cation and indole molecule. The strength of cation-π interactions get substantially reduced when the metal ion is solvated or the size of metal cation increases. Quantum theory of atoms in molecules analysis of cation-π interaction indicates that there is a correlation between the electron density (ρ( r)) in the cage critical points generated upon complexation and the distance between metal cation and centroid of phenyl ring in indole molecule. [ABSTRACT FROM AUTHOR]

The B3LYP/6-311++G (d,p) density functional approach was used to study the gas-phase metal affinities of Guanosine (ribonucleoside) for the Li, Na, K, Mg, Ca, Zn, and Cu cations. In this study we determine coordination geometries, binding strength, absolute metal ion affinities, and free energies for the most stable products. We have also compared the results for Guanosine, with our previously reported results for 2′-Deoxyguanosine. Based on the results, it is obvious that MIA is strongly dependent on the charge-to-size ratio of the cation. Guanosine interacts more strongly with Zn than do with Mg, Ca, and Cu and therefore stronger interactions lead to higher MIA. In both free molecules and their complexes, the Syn orientation of the base is stabilized by an intramolecular O5′-H···N3 hydrogen bond and the anti orientation of the base is stabilized by an intramolecular C-H···O hydrogen bond formed between the (C8-H8) and the O5′ atom of the sugar moiety. It is also interesting to mention that linear correlation between calculated MIA values and the atomic numbers ( Z) of the metal ions of Li, Na, and K were found. Furthermore, the influences of metal cationization on the strength of the N-glycosidic bond, torsion angles, angle of pseudorotation ( P), and intramolecular C-H···O and O-H···O hydrogen bonds have been studied. Natural bond orbital (NBO) analysis was performed to calculate the charge transfer and natural population analysis of the complexes. Quantum theory of atoms in molecules (QTAIM) was also applied to determine the nature of interactions. [ABSTRACT FROM AUTHOR]

Hydrogen bonding interactions between thymine nucleobase and 2′-deoxythymidine nucleoside (dT) with some biological anions such as F (fluoride), Cl (chloride), OH (hydroxide), and NO (nitrate) have been explored theoretically. In this study, complexes have been studied by density functional theory (B3LYP method and 6-311++G (d,p) basis set). The relevant geometries, energies, and characteristics of hydrogen bonds (H-bonds) have been systematically investigated. There is a correlation between interaction energy and proton affinity for complexes of thymine nucleobase. The nature of all the interactions has been analyzed by means of the natural bonding orbital (NBO) and quantum theory atoms in molecules (QTAIM) approaches. Donors, acceptors, and orbital interaction energies were also calculated for the hydrogen bonds. Excellent correlations between structural parameter (δ R) and electron density topological parameter ( ρ) as well as between E(2) and ρ have been found. It is interesting that hydrogen bonds with anions can affect the geometry of thymine and 2′-deoxythymidine molecules. For example, these interactions can change the bond lengths in thymine nucleobase, the orientation of base unit with respect to sugar ring, the furanose ring puckering, and the C1′-N1 glycosidic linkage in dT nucleoside. Thus, it is necessary to obtain a fundamental understanding of chemical behavior of nucleobases and nucleosides in presence of anions. [ABSTRACT FROM AUTHOR]