The Search for Low Energy Conformational Families of Small Peptides: Searching for Active Conformations of Small Peptides in the Absence of a Known Receptor

Breast cancer is the most common cancer among women. Tamoxifen is the preferred drug for estrogen receptor-positive breast cancer treatment, yet many of these cancers are intrinsically resistant to tamoxifen or acquire resistance during treatment. Therefore, scientists are searching for breast cancer drugs that have different molecular targets. Previous work revealed that 8-mer and cyclic 9-mer peptides inhibit breast cancer in mouse and rat model systems, interacting with an unknown receptor, while peptides smaller than eight amino acids did not inhibit breast cancer. We have shown that the use of replica exchange molecular dynamics predicts structure and dynamics of active peptides, leading to the discovery of smaller peptides with full biological activity. These simulations identified smaller peptide analogs with a conserved turn, a β-turn formed in the larger peptides. These analogs inhibit estrogen-dependent cell growth in a mouse uterine growth assay, a test showing reliable correlation with human breast cancer inhibition. We outline the computational methods that were tried and used with the experimental information that led to the successful completion of this research.

Accurate Predictions of Water Cluster Formation, (H2O)n=2-10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Doubly-charged ethane ions: Solution to the dilemma of stability predicted by theory and instability observed in experiment

Potential energy curves have been computed for [C2H6]2+ ions and the results used to interpret the conspicuous absence of these ions in 2E mass spectra and in charge-stripping experiments. The energies and structures of geometry-optimized ground-state singlet and excited-state triplet [C2H6]2+ ions have been determined along with energies for different decomposition barriers and dissociation asymptotes. Although singlet and triplet [C2H6]2+ ions can exist as stable entities, they possess low energy barriers to decomposition. Vertical Franck-Condon transitions, involving electron impact ionization of ethane as well as charge-stripping collisions of [C2H6]+ ions, produce [C2H6]2+ ions which promptly dissociate since they are formed with energies in excess of various decomposition barriers. Appearance energies computed for doubly-charged ethane fragment ions are in accordance with experimental values.

Hydration of OCS with One to Four Water Molecules in Atmospheric and Laboratory Conditions

Carbonyl sulfide is the most abundant sulfur gas in the atmosphere. We have used MP2 and CCSD(T) theory to study the structures and thermochemistries of carbonyl sulfide interacting with one to four water molecules. We have completed an extensive search for clusters of OCS(H2O)n, where n = 1−4. We located three dimers, two trimers, five tetramers, and four pentamers with the MP2/aug-cc-pVDZ method. In each of the complexes with two or more waters, OCS preferentially interacts with low-energy water clusters. Our results match current theoretical and experimental literature, showing correlation with available geometries and frequencies for the OCS(H2O) species. The CCSD(T)/aug-cc-pVTZ thermochemical values combined with the average amount of OCS and the saturated concentration of H2O in the troposphere, lead to the prediction of 106 OCS(H2O) clusters·cm−3 and 102 OCS(H2O)2 clusters·cm−3 at 298 K. We predict the structures of OCS(H2O)n, n = 1−4 that should predominate in a low-temperature molecular beam and identify specific infrared vibrations that can be used to identify these different clusters.

Structures of Cage, Prism, and Book Isomers of Water Hexamer from Broadband Rotational Spectroscopy

Theory predicts the water hexamer to be the smallest water cluster with a three-dimensional hydrogen-bonding network as its minimum energy structure. There are several possible low-energy isomers, and calculations with different methods and basis sets assign them different relative stabilities. Previous experimental work has provided evidence for the cage, book, and cyclic isomers, but no experiment has identified multiple coexisting structures. Here, we report that broadband rotational spectroscopy in a pulsed supersonic expansion unambiguously identifies all three isomers; we determined their oxygen framework structures by means of oxygen-18–substituted water (H218O). Relative isomer populations at different expansion conditions establish that the cage isomer is the minimum energy structure. Rotational spectra consistent with predicted heptamer and nonamer structures have also been identified.

Hydration of the Bisulfate Ion: Atmospheric Implications

Using molecular dynamics configurational sampling combined with ab initio energy calculations, we determined the low energy isomers of the bisulfate hydrates. We calculated the CCSD(T) complete basis set (CBS) binding electronic and Gibbs free energies for 53 low energy isomers of HSO4–(H2O)n=1–6 and derived the thermodynamics of adding waters sequentially to the bisulfate ion and its hydrates. Comparing the HSO4–/H2O system to the neutral H2SO4/H2O cluster, water binds more strongly to the anion than it does to the neutral molecules. The difference in the binding thermodynamics of HSO4–/H2O and H2SO4/H2O systems decreases with increasing number of waters. The thermodynamics for the formation of HSO4–(H2O)n=1–5 is favorable at 298.15 K, and that of HSO4–(H2O)n=1–6 is favorable for T < 273.15 K. The HSO4– ion is almost always hydrated at temperatures and relative humidity values encountered in the troposphere. Because the bisulfate ion binds more strongly to sulfuric acid than it does to water, it is expected to play a role in ion-induced nucleation by forming a strong complex with sulfuric acid and water, thus facilitating the formation of a critical nucleus.

Accurate Predictions of Water Cluster Formation, (H2O)n=2−10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Broadband Fourier transform rotational spectroscopy for structure determination: The water heptamer

Over the recent years chirped-pulse, Fourier-transform microwave (CP-FTMW) spectrometers have chan- ged the scope of rotational spectroscopy. The broad frequency and large dynamic range make possible structural determinations in molecular systems of increasingly larger size from measurements of heavy atom (13C, 15N, 18O) isotopes recorded in natural abundance in the same spectrum as that of the parent isotopic species. The design of a broadband spectrometer operating in the 2–8 GHz frequency range with further improvements in sensitivity is presented. The current CP-FTMW spectrometer performance is benchmarked in the analyses of the rotational spectrum of the water heptamer, (H2O)7, in both 2– 8 GHz and 6–18 GHz frequency ranges. Two isomers of the water heptamer have been observed in a pulsed supersonic molecular expansion. High level ab initio structural searches were performed to pro- vide plausible low-energy candidates which were directly compared with accurate structures provided from broadband rotational spectra. The full substitution structure of the most stable species has been obtained through the analysis of all possible singly-substituted isotopologues (H218O and HDO), and a least-squares rm(1) geometry of the oxygen framework determined from 16 different isotopic species compares with the calculated O–O equilibrium distances at the 0.01 Å level.

Broadband Fourier transform rotational spectroscopy for structure determination: The water heptamer

Over the recent years chirped-pulse, Fourier-transform microwave (CP-FTMW) spectrometers have changed the scope of rotational spectroscopy. The broad frequency and large dynamic range make possible structural determinations in molecular systems of increasingly larger size from measurements of heavy atom (C-13, N-15, O-18) isotopes recorded in natural abundance in the same spectrum as that of the parent isotopic species. The design of a broadband spectrometer operating in the 2-8 GHz frequency range with further improvements in sensitivity is presented. The current CP-FTMW spectrometer performance is benchmarked in the analyses of the rotational spectrum of the water heptamer, (H2O)(7), in both 2-8 GHz and 6-18 GHz frequency ranges. Two isomers of the water heptamer have been observed in a pulsed supersonic molecular expansion. High level ab initio structural searches were performed to provide plausible low-energy candidates which were directly compared with accurate structures provided from broadband rotational spectra. The full substitution structure of the most stable species has been obtained through the analysis of all possible singly-substituted isotopologues ((H2O)-O-18 and HDO), and a least-squares r(m)((1)) geometry of the oxygen framework determined from 16 different isotopic species compares with the calculated O-O equilibrium distances at the 0.01 angstrom level. (C) 2013 Elsevier B.V. All rights reserved.