Elongation cutoff technique: low-order scaling SCF method

The elongation cutoff technique at restricted Hartree-Fock (HF) level of theory in conventional type of calculations, i.e., with two electron integrals (TEI) stored on a disc, is presented for two model systems. It is demonstrated that the number of TEI in the elongation cutoff calculations increases linearly with the system size thus, allowing to extend the conventional type of calculations to bigger systems. The step CPU (central processing unit) time in the elongation cutoff calculations is much lower than in the HF reference calculations. Such behavior reduces significantly the prefactor in the quadratic scaling relation. The total CPU time in the elongation calculation is about 40% lower than in the conventional HF calculations or comparable to direct type of calculations with the quantum fast multipoles method employed. It is shown that by introducing the interaction radius one can obtain linear scaling in the SCF calculations. Figure: The structure of density matrix and total CPU timings for polyglycine clusters in the elongation cutoff calculations. The structure of density matrix and total CPU timings for polyglycine clusters in the elongation cutoff calculations     

Ab Initio and Molecular Mechanics Calculations of Various Substituted Ureas - Rotational Barriers and a New Parametrization for Ureas

High-level ab initio calculations have been performed on urea, methylurea, tetramethylurea and three isomers of dimethylurea to obtain accurate rotational barriers. Results of MP2(fc)/6-31 G(d) calculations are compared to those with lower basis sets and semiempirical calculations. The MM2(87) force field has been parameterized.     

Structure,energetics and properties of some molecules with potent anti-HIV activity: a theoretical study

Density functional theory (DFT; B3LYP) and Hartree–Fock (HF; 3-21G, 6-31G(d) and 6-311G(d,p)) calculations with complete geometry optimisations are carried out in the ground state on five 6-aminoquinolone derivatives, which have been proved to be highly effective in inhibiting HIV replication, to study their structures, energetics and HOMO–LUMO correlation with physiological action. The gas-phase calculations and single-point polarisable continuum model water-phase calculations show that the molecules are highly effective in inhibiting HIV replication, which is in excellent agreement with the experiment. Structural features, energies, charge densities and HOMO–LUMO correlation have been found to substantiate the experimental findings. Compound 4 (pyrazine) shows some special features in DFT calculations which are not found in HF calculations. In the present series, HF results are more reliable as expected.     

A new approach to possible substrate binding mechanisms for nitrile hydratase

We combined normal mode analysis (NMA) with cavity calculations as a method to get more insight into static crystal structures. We used nitrile hydratase (NHase) as a case study, and the crystal structure of a complex of Pseudonocardia thermophila NHase (1UGP) with n-butyric acid was chosen as a reference structure. The reference structure was compared with the other available NHase crystal structures. Cavity calculations of the static structures showed the entrances to the active site and also a possible function of the N-terminal in the substrate selection of the Co-type NHase. When NMA was combined with cavity calculations, a closing-opening passage was observed. Analysis of low frequency modes combined with cavity calculations led us to propose "breathing" and "flip-flop" mechanisms which might be a key part of the substrate binding mechanism.     

Divalent cations and the electrostatic potential around DNA: Monte Carlo and Poisson-Boltzmann calculations

The predictions of counterion condensation theory for divalent ions were tested by comparison with the results of Monte Carlo calculations on an all-atom model of DNA. Monovalent-divalent competition at the polyelectrolyte surface was investigated by varying the partial molar volume of divalent ions. To assess the viability of using Poisson-Boltzmann (PB) calculations for determining divalent ion concentrations at DNA surfaces, Monte Carlo (MC) calculations were compared with PB calculations using different models of the dielectric continuum. It was determined that, while standard PB calculations of divalent ion surface densities are about 25-30% below those predicted by MC techniques, and somewhat larger than errors previously determined for monovalent ions, errors due to the use of the mean-field approximation of PB theory are smaller than those arising from common assumptions regarding the dielectric continuum.     

Ab initio multiconfiguration reference perturbation theory calculations on the energetics of low-energy spin states of iron(III) porphyrins

Although a major goal of inorganic spectroscopy is to determine the energetics of the low-lying spin states of transition metal complexes, surprisingly little has been accomplished in this respect by means of accurate ab initio calculations. Against this context, we present ab initio multiconfiguration reference perturbation theory (CASPT2) calculations with large basis sets on the low-lying spin states of Fe(III)(P)Cl and [Fe(P)Cl](+) (P(2-)=porphinato). The CASPT2 results on the energetics of various low-lying spin states studied differ significantly, sometimes even dramatically, from those obtained from density functional theory calculations.     

An electron spin resonance investigation and molecular orbital calculation of the anion radical intermediate in the enzymatic cis-trans isomerization of furylfuramide, a nitrofuran derivative of ethylene

The enzymatic cis-trans isomerization of nitrofuran derivatives has been proposed to occur via the formation of a radical anion intermediate. ESR investigations, in conjunction with intermediate neglect of differential overlap (INDO) molecular orbital calculations, support this concept by demonstrating the enzymatic generation of cis and trans radical anions of 3-(5-nitro-2-furyl)-2-(2-furyl) acrylamide. The INDO calculations further indicate that the rotational barrier between the cis and trans anion radicals of this compound is only 5--10 kcal/mol, whereas a 70 kcal/mol barrier exists for the parent geometric isomers. Hyperfine splitting constants for the cis-trans conformers have been assigned on the basis of INDO calculations. Surprisingly, only the nitrogen hyperfine splitting of the nitro group is distinguishably different in the two conformers, a result which is not inconsistent with the INDO calculations.     

Conformational analysis of octa- and tetrabromo tetraphenylporphyrins and their Ni(II) and Tb(III) complexes

Molecular mechanics (MM) calculations were used to analyze the puckering of metalloporphyrins as a function of metal ion size and the position of substituents on the porphyrin periphery, on a three series of octa- and tetrabromo tetraphenylporphyrins: without metal, and with Ni(II), and Tb(III) as representative small and large metal ions, respectively. Molecular energy optimization calculations were carried out using the Consistent Force Field (CFF) program, with the parameters developed previously and new parameters for bromine atom. Normal-coordinate structural decomposition (NSD) analysis was performed on the equilibrium structures obtained by MM calculations. The conformers are also stereochemically characterized, compared with available X-ray structures and with the conformers obtained in our previous MM study using chloro instead of bromo beta-pyrrole substituents.     

An atomic-level mechanism for molybdenum nitrogenase. Part 1. Reduction of dinitrogen

The properties of the Fe and Mo sites of the iron-molybdenum cofactor of nitrogenase with respect to binding and activation of N(2) have been studied by molecular mechanics calculations on the local protein environment and by density functional theory (DFT) calculations on subsections of the cofactor. The DFT calculations indicate that the homocitrate ligand of the cofactor can become monodentate on reduction, allowing N(2) to bind at Mo. In addition, the neighboring Fe atom plays a crucial role in N(2) reduction by stabilizing the initial reduced N(2) species and by facilitating cleavage of the N-N bond. The various possible isomers for partially reduced N(2) intermediates have been compared by DFT, and a detailed model for the reduction of N(2) is developed based on these results, together with chemical precedents and the available biochemical data for nitrogenase.     

CD of nucleic acids: III. Calculated CD of RNAs from new A, U, G, and C transition-moment parameters

New adenine (A) and uracil (U) π → π* transition-moment parameters have been derived from a recently developed semiempirical procedure. Using conformational energy probabilities based on the Boltzmann equation, the new parameters were assigned by optimizing the calculated CD of cyclic nucleotides against measured CD. The derived A-and U-parameters (along with guanine and cytosine parameters derived previously by the same procedure) have been assessed in CD spectral calculations of some polyribonucleic acid sequences, in assumed A-class geometries. Comparisons have been made between CD spectra calculated from the newly derived parameters and those calculated from parameters obtained from a combination of crystal optical measurements and quantum-mechanical calculations. Although some spectral differences do occur, for the RNA sequences considered, no major disagreements were found in CD spectral signs and shapes, between measurements and calculations. Overall, the results indicate that the newly derived A-, U-, G-, and C-parameters show better agreement between theory and experiment than those used in previous nucleic acid CD calculations.     

Power calculations for a general class of family-based association tests: dichotomous traits

Using large-sample theory, we present a unified approach to power calculations for family-based association tests. Currently available methods for power calculations are restricted to special designs or require approximations or simulations. Our analytical approach to power calculations is broadly applicable in many settings. We discuss power calculations for two scenarios that have high practical relevance and in which power previously could only be assessed by simulation studies or by approximations: (1) studies using both affected and unaffected offspring and (2) studies with missing parental information. When the population prevalence is high, it can be worthwhile to genotype unaffected offspring. For many scenarios, high power can be achieved with reasonable sample sizes, even when no parental information is available.     

Comparison of three different ways of measuring distances between residences and high voltage power lines

The aim of this study was to evaluate whether distance data based on calculations by use of digitalized geographical information systems (GIS) and distance data based on measurements on 1:5000 maps agree sufficiently with on site distance measurements to be used as input to magnetic field calculations in epidemiological studies. The analysis were performed by use of weighted kappa (kappa(w)) statistical method described by Bland and Altman for comparison of measures of agreement. Map measurements showed better agreement with on site measurements than GIS calculations did. However, we consider both methods appropriate for use in larger epidemiological studies if the results are interpreted with caution. GIS calculations have the advantage of being both time and cost saving.     

An energetic correlation of ab initio and NMR studies of the 3'-gauche effect in 3'-substituted thymidines

A straightforward correlation of our experimental NMR findings on 3'-substituted thymidine derivatives with that of the ab initio calculations shows that (i) the delta Go298kNRM of N reversible S equilibrium in nucleoside can be predicted from the ab initio calculated delta ES-N obtained from 6-311++G** level of theory; (ii) the substituent-dependent steric and stereoelectronic effects on the bias of the two-state N reversible S equilibrium in nucleosides can also be predicted from the ab initio calculations with sufficiently large basis functions, and (iii) the necessity of mimicking the solvation behaviour of the experimental NMR measurement condition in the ab initio calculations of biomolecules is also emphasized.     

Electronic absorption spectroscopy and time-dependent density functional theory calculations on the nickel(II) complex of 1,4-bis(pyrrol-2-ylmethyleneamino)butane

The electronic structure of the 1,4-bis(pyrrol-2-ylmethyleneamino)butane nickel(II) complex has been studied using electronic absorption spectroscopy and density functional theory (DFT) calculations. The DFT optimised structure is in excellent agreement with the X-ray crystal structure of the complex and time-dependent DFT calculations have been used to probe the nature of the transitions observed in the electronic absorption spectrum.     

Theoretical study of the infrared spectrum of 5-phenyl-1,3,4-oxadiazole-2-thiol by using DFT calculations

We have carried out a structural and vibrational study for 5-phenyl-1,3,4-oxadiazole-2-thiol by using the infrared (IR) spectrum and theoretical calculations. For a complete assignment of the compound IR spectrum, density functional theory calculations were combined with Pulay's scaled quantum mechanical force field methodology in order to fit the theoretical wavenumber values to the experimental ones. An agreement between theoretical and available experimental results was found. The theoretical vibrational calculations allowed us to obtain a set of scaled force constants fitting the observed wavenumbers. The results were then used to predict the Raman spectra, for which there are no experimental data. The nature of the benzyl and oxadiazole rings was studied by means of natural bond order and atoms in molecules theory calculations. In addition, the frontier molecular (highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO)) orbitals were analysed and compared with those calculated for the oxadiazole molecule.     

Theoretical studies of the ATP hydrolysis mechanism of myosin

The ATP hydrolysis mechanism of myosin was studied using quantum chemical (QM) and molecular dynamics calculations. The initial model compound for QM calculations was constructed on the basis of the energy-minimized structure of the myosin(S1dc)-ATP complex, which was determined by molecular mechanics calculations. The result of QM calculations suggested that the ATP hydrolysis mechanism of myosin consists of a single elementary reaction in which a water molecule nucleophilically attacked gamma-phosphorus of ATP. In addition, we performed molecular dynamics simulations of the initial and final states of the ATP hydrolysis reaction, that is, the myosin-ATP and myosin-ADP.Pi complexes. These calculations revealed roles of several amino acid residues (Lys185, Thr186, Ser237, Arg238, and Glu459) in the ATPase pocket. Lys185 maintains the conformation of beta- and gamma-phosphate groups of ATP by forming the hydrogen bonds. Thr186 and Ser237 are coordinated to a Mg(2+) ion, which interacts with the phosphates of ATP and therefore contributes to the stabilization of the ATP structure. Arg238 and Glu459, which consisted of the gate of the ATPase pocket, retain the water molecule acting on the hydrolysis at the appropriate position for initiating the hydrolysis.     

Differences and similarities in binding of pyruvate and L-lactate in the active site of M4 and H4 isoforms of human lactate dehydrogenase

We present QM/MM calculations that show differences in geometries of active sites of M(4) and H(4) isoforms of human LDH ligated with oxamate, pyruvate or L-lactate. As the consequence of these differences, binding isotope effects of the methyl hydrogen atoms of pyruvate and l-lactate may be used to experimentally distinguish these isoforms. Based on the FEP calculations we argue that L-lactate is a better candidate for the experimental studies. Our calculations of energies of interactions of ligands with the active site residues provide explanation for the observed experimentally sensitivity to inhibition of the M(4) isoenzyme isoform and pinpoint the differences to interactions of the ligand with the histidine residue. We conclude that pyruvate interacts much stronger in the active site of H(4) than M(4) isoform and that the latter interactions are weaker than with water molecules in the aqueous solution.     

Calculation of transmission coefficients for some permeant molecules in human red cell and resting axolemma squid axon membranes

Quantum mechanical calculations of transmission coefficients for some permeant molecules across the human red cell and resting axolemma squid axon membranes are carried out. The calculations depend on (i) the molecular weight of the molecule and (ii) the depth and width of the potential well of the membrane. In most cases good agreement between calculated and experimental values is found.     

A comparison of the circular dichroism spectra of synthetic DNA sequences of the homopurine . homopyrimidine and mixed purine- pyrimidine types.

We have obtained the ultraviolet circular dichroism spectra of two repeating trinucleotide DNAs, poly [d(A-G-G).d(C-C-T)] and poly[d(A-A-G).d(C-T-T)], that have all purines on one strand and all pyrimidines on the other. These spectra, together with spectra of other synthetic polymers, can be combined to give 3 first-neighbor calculations of the spectrum of poly[d(A).d(T)] and 2 first-neighbor calculations of the spectrum of poly [d(G).d(C)]. The results show (1) that first-neighbor calculations utilizing only spectra of homopurine.homopyrimidine DNA sequences are no more accurate than are similar calculations that involve spectra of mixed purine-pyrimidine sequences, demonstrating that double-stranded homopurine.homopyrimidine sequences do not obviously belong to a special class of secondary conformations, and (2) that the wavelength region above 250 nm in the CD spectra of synthetic DNAs is least predictable from first-neighbor equations, probably because this region is especially sensitive to sequence-dependent conformational differences.