Torsional path integral Monte Carlo method for calculating the absolute quantum free energy of large molecules

A new technique for evaluating the absolute free energy of large molecules is presented. Quantum-mechanical contributions to the intramolecular torsions are included via the torsional path integral Monte Carlo (TPIMC) technique. Importance sampling schemes based on uncoupled free rotors and harmonic oscillators facilitate the use of the TPIMC technique for the direct evaluation of quantum partition functions. Absolute free energies are calculated for the molecules ethane, n-butane, n-octane, and enkephalin, and quantum contributions are found to be significant. Comparison of the TPIMC technique with the harmonic oscillator approximation and a variational technique is performed for the ethane molecule. For all molecules, the quantum contributions to free energy are found to be significant but slightly smaller than the quantum contributions to internal energy

Quantum free energies of the conformers of glycine on an ab initio potential energy surface

The torsional path integral Monte Carlo (TPIMC) technique is used to study the five lowest-energy conformers of glycine. The theoretical method provides an anharmonic and quantum-mechanical description of conformational free energy and is used for the first time with an ab initio potential energy surface. The 3-dimensional torsional potential energy surface of glycine was obtained at the MP2/6-311++G** level of theory and is optimized with respect to the non-torsional degrees of freedom. Calculated conformer populations compare well with those reported in recent matrix-isolation infrared spectroscopy experiments. An additional conformer, not yet observed, is predicted to be heavily populated in the thermal equilibria probed by experiment, and a new explanation for its elusiveness is provided. Quantum effects, such as zero point energy, are found to substantially alter conformer populations, and an algorithm for estimating the role of non-torsional vibrations in the conformational thermodynamics of a molecule is introduced

Torsional path integral Monte Carlo method for the quantum simulation of large molecules

A molecular application is introduced for calculating quantum statistical mechanical expectation values of large molecules at nonzero temperatures. The Torsional Path Integral Monte Carlo (TPIMC) technique applies an uncoupled winding number formalism to the torsional degrees of freedom in molecular systems. The internal energy of the molecules ethane, n-butane, n-octane, and enkephalin are calculated at standard temperature using the TPIMC technique and compared to the expectation values obtained using the harmonic oscillator approximation and a variational technique. All studied molecules exhibited significant quantum mechanical contributions to their internal energy expectation values according to the TPIMC technique. The harmonic oscillator approximation approach to calculating the internal energy performs well for the molecules presented in this study but is limited by its neglect of both anharmonicity effects and the potential coupling of intramolecular torsion

Torsional anharmonicity in the conformational thermodynamics of flexible molecules

We present an algorithm for calculating the conformational thermodynamics of large, flexible molecules that combines ab initio electronic structure theory calculations with a torsional path integral Monte Carlo (TPIMC) simulation. The new algorithm overcomes the previous limitations of the TPIMC method by including the thermodynamic contributions of non-torsional vibrational modes and by affordably incorporating the ab initio calculation of conformer electronic energies, and it improves the conventional ab initio treatment of conformational thermodynamics by accounting for the anharmonicity of the torsional modes. Using previously published ab initio results and new TPIMC calculations, we apply the algorithm to the conformers of the adrenaline molecule

Collision-induced conformational changes in glycine

We present quantum dynamical calculations on the conformational changes of glycine in collisions with the He, Ne, and Ar rare-gas atoms. For two conformer interconversion processes (III-->I and IV-->I), we find that the probability of interconversion is dependent on several factors, including the energy of the collision, the angle at which the colliding atom approaches the glycine molecule, and the strength of the glycine-atom interaction. Furthermore, we show that attractive interactions between the colliding atom and the glycine molecule catalyze conformer interconversion at low collision energies. In previous infrared spectroscopy studies of glycine trapped in rare-gas matrices and helium clusters, conformer III has been consistently observed, but conformer IV has yet to be conclusively detected. Because of the calculated thermodynamic stability of conformer IV, its elusiveness has been attributed to the IV-->I conformer interconversion process. However, our calculations present little indication that IV-->I interconversion occurs more readily than III-->I interconversion. Although we cannot determine whether conformer IV interconverts during experimental Ne- and Ar-matrix depositions, our evidence suggests that the conformer should be present in helium droplets. Anharmonic vibrational frequency calculations illustrate that previous efforts to detect conformer IV may have been hindered by the overlap of its IR-absorption bands with those of other conformers. We propose that the redshifted symmetric –CH2 stretch of conformer IV provides a means for its conclusive experimental detection

Reduced dimensionality spin-orbit dynamics of CH3 + HCl reversible arrow CH4 Cl on ab initio surfaces

A reduced dimensionality quantum scattering method is extended to the study of spin-orbit nonadiabatic transitions in the CH3 + HCl reversible arrow CH4 + Cl(P-2(J)) reaction. Three two-dimensional potential energy surfaces are developed by fitting a 29 parameter double-Morse function to CCSD(T)/IB//MP2/cc-pV(T+d)Z-dk ab initio data; interaction between surfaces is described by geometry-dependent spin-orbit coupling functions fit to MCSCF/cc-pV(T+d)Z-dk ab initio data. Spectator modes are treated adiabatically via inclusion of curvilinear projected frequencies. The total scattering wave function is expanded in a vibronic basis set and close-coupled equations are solved via R-matrix propagation. Ground state thermal rate constants for forward and reverse reactions agree well with experiment. Multi-surface reaction probabilities, integral cross sections, and initial-state selected branching ratios all highlight the importance of vibrational energy in mediating nonadiabatic transition. Electronically excited state dynamics are seen to play a small but significant role as consistent with experimental conclusions. (C) 2011 American Institute of Physics. [doi:10.1063/1.3592732

Positron emission tomography with f18-fluorodeoxyglucose in the staging and preoperative evaluation of malignant pleural mesothelioma

AbstractObjectives: The purpose of this study was to evaluate the utility of positron emission tomography with F18-fluorodeoxyglucose in the preoperative evaluation and staging of malignant mesothelioma in patients who were candidates for aggressive combined modality therapy. Methods: Eighteen consecutive patients with biopsy-proven malignant mesothelioma underwent positron emission tomographic scanning. The results of positron emission tomographic imaging were compared with results obtained by computed tomography, mediastinoscopy, thoracoscopy, and pathologic examination of surgical specimens. All patients fasted and received an average of 14.5 ± 2.7 mCi of F18-fluorodeoxyglucose for positron emission tomographic scanning. Attenuation-corrected whole-body and regional emission images of the chest and upper abdomen were acquired and formatted into transaxial, coronal, and sagittal images. Results: All primary malignant mesotheliomas accumulated F18-fluorodeoxyglucose, and the mean standardized uptake value was 7.6 (range, 3.33-14.85; n = 9). There were no false-negative results of positron emission tomography. Identification of occult extrathoracic metastases by positron emission tomography was the basis for excluding two patients from surgical therapy. There were two false-positive results of positron emission tomography: increased F18-fluorodeoxyglucose uptake in the contralateral chest that was negative by thoracoscopic biopsy (n = 1) and increased abdominal F18-fluorodeoxyglucose uptake after partial colectomy for diverticular disease (n = 1). Conclusions: Positron emission tomography can identify malignant pleural mesothelioma and appears to be a useful noninvasive staging modality for patients being considered for aggressive combined modality therapy. (J Thorac Cardiovasc Surg 2000;120:128-33

Cell-State-Specific Metabolic Dependency in Hematopoiesis and Leukemogenesis

The balance between oxidative and non-oxidative glucose metabolism is essential for a number of pathophysiological processes. By deleting enzymes that affect aerobic glycolysis with different potencies, we examine how modulating glucose metabolism specifically affects hematopoietic and leukemic cell populations. We find that deficiency in the M2 pyruvate kinase isoform (PKM2) reduces levels of metabolic intermediates important for biosynthesis and impairs progenitor function without perturbing hematopoietic stem cells (HSC), whereas lactate dehydrogenase-A (LDHA) deletion significantly inhibits the function of both HSC and progenitors during hematopoiesis. In contrast, leukemia initiation by transforming alleles putatively affecting either HSC or progenitors is inhibited in the absence of either PKM2 or LDHA, indicating that the cell state-specific responses to metabolic manipulation in hematopoiesis do not apply to the setting of leukemia. This finding suggests that fine-tuning the level of glycolysis may be therapeutically explored for treating leukemia while preserving HSC function.National Institutes of Health (U.S.) (Grants P30CA147882 and R01CA168653)Smith Family FoundationBurroughs Wellcome FundVirginia and D.K. Ludwig Fund for Cancer ResearchDamon Runyon Cancer Research Foundatio