Salvatore Profeta

Bio: Salvatore Profeta is an academic researcher from Florida State University. The author has contributed to research in topics: Ab initio & Ab initio quantum chemistry methods. The author has an hindex of 15, co-authored 26 publications receiving 4798 citations. Previous affiliations of Salvatore Profeta include University of California, San Francisco & Monsanto.

A new force field for molecular mechanical simulation of nucleic acids and proteins

Abstract: We present the development of a force field for simulation of nucleic acids and proteins. Our approach began by obtaining equilibrium bond lengths and angles from microwave, neutron diffraction, and prior molecular mechanical calculations, torsional constants from microwave, NMR, and molecular mechanical studies, nonbonded parameters from crystal packing calculations, and atomic charges from the fit of a partial charge model to electrostatic potentials calculated by ab initio quantum mechanical theory. The parameters were then refined with molecular mechanical studies on the structures and energies of model compounds. For nucleic acids, we focused on methyl ethyl ether, tetrahydrofuran, deoxyadenosine, dimethyl phosphate, 9-methylguanine-l-methylcytosine hydrogen-bonded complex, 9-methyladenine-l-methylthymine hydrogen-bonded complex, and 1,3-dimethyluracil base-stacked dimer. Bond, angle, torsional, nonbonded, and hydrogen-bond parameters were varied to optimize the agreement between calculated and experimental values for sugar pucker energies and structures, vibrational frequencies of dimethyl phosphate and tetrahydrofuran, and energies for base pairing and base stacking. For proteins, we focused on 4>,'lt maps of glycyl and alanyl dipeptides, hydrogen-bonding interactions involving the various protein polar groups, and energy refinement calculations on insulin. Unlike the models for hydrogen bonding involving nitrogen and oxygen electron donors, an adequate description of sulfur hydrogen bonding required explicit inclusion of lone pairs. There are two fundamental problems in simulating the struc­ tural and energetic properties of molecules: the first is how to choose an analytical been placed E(R) which correctly describes the energy of the system in terms of its 3N degrees of freedom. The second is how the simulation can search or span conforma­ tional space (R) in order to answer questions posed by the scientist interested in the properties of the system. For complex systems, solution to the first problem are an es­ sential first step in attacking the second problem, and thus, considerable effort has been placed in developing analytical functions that are simple enough to allow one to simulate the properties of complex molecules yet accurate enough to obtain meaningful estimates for structures and energies. In the case of the structures and thermodynamic stabilities of saturated hydrocarbons in inert solvents or the gas phase, the first problem has been essentially solved by molecular mechanics ap­ proaches of Allinger, I Ermer and Lifson,2 and their co-workers. However, for polar and ionic molecules in condensed phases, unsolved questions remain as to the best form of the analytical function E(R). In the area of proteins and peptides, seminal work has come from the Scheraga 3 and Lifson 4 schools. The Scheraga group has used both crystal packing (intermolecular) and con­ formational properties of peptides to arrive at force fields ECEPP, UNECEPP, and EPEN for modeling structural and thermodynamic properties of peptides and proteins. Levitt, using the energy refinement software developed in the Lifson group, has proposed a force field for proteins based on calculations on lysozyme,S and Gelin and Karplus have adapted this software along with many parameters from the Scheraga studies to do molecular dynamics

Origins of Selectivity for the [2+2] Cycloaddition of α,β-unsaturated Ketones within a Porous Self-assembled Organic Framework

Abstract: This article studies the origins of selectivity for the [2+2] cycloadditions of α,β-unsaturated ketones within a porous crystalline host. The host, formed by the self-assembly of a bis-urea macrocycle, contains accessible channels of ∼6 A diameter and forms stable inclusion complexes with a variety of cyclic and acyclic α,β-unsaturated ketone derivatives. Host 1 crystals provide a robust confined reaction environment for the highly selective [2+2] cycloaddition of 3-methyl-2-cyclopentenone, 2-cyclohexenone, and 2-methyl-2-cyclopentenone, forming their respective exo head-to-tail dimers in high conversion. The products are readily extracted from the self-assembled host and the crystalline host can be efficiently recovered and reused. Molecular modeling studies indicate that the origin of the observed selectivity is due to the excellent match between the size and shape of these guests to dimensions of the host channel and to the preorganization of neighboring enones into favorable reaction geometries. Small...

Mechanism of Initiation in the Phillips Ethylene Polymerization Catalyst: Redox Processes Leading to the Active Site

Abstract: The detailed mechanism by which ethylene polymerization is initiated by the inorganic Phillips catalyst (Cr/SiO2) without recourse to an alkylating cocatalyst remains one of the great unsolved mysteries of heterogeneous catalysis. Generation of the active catalyst starts with reduction of CrVI ions dispersed on silica. A lower oxidation state, generally accepted to be CrII, is required to activate ethylene to form an organoCr active site. In this work, a mesoporous, optically transparent monolith of CrVI/SiO2 was prepared using sol–gel chemistry in order to monitor the reduction process spectroscopically. Using in situ UV–vis spectroscopy, we observed a very clean, stepwise reduction by CO of CrVI first to CrIV, then to CrII. Both the intermediate and final states show XANES consistent with these oxidation state assignments, and aspects of their coordination environments were deduced from Raman and UV–vis spectroscopies. The intermediate CrIV sites are inactive toward ethylene at 80 °C. The CrII sites, wh...

Synthesis, Topoisomerase I Inhibitory Activity, and in Vivo Evaluation of 11-Azacamptothecin Analogs

Abstract: A series of analogs based on a novel template, 11-aza-(20S)-camptothecin, were obtained from total synthesis and tested as potential anticancer drugs in the topoisomerase I enzyme cleavable complex assay. The parent compound 11-aza-(20S)-camptothecin (8) was derived from a Friedlander condensation between the known aminopyridine derivative 3-(3-amino-4-picolylidene)-p-toluidine and optically active tricyclic ketone 7. Compound 8 had activity approximately twice that of (20S)-camptothecin in the calf thymus topoisomerase I cleavable complex assay. Compounds were prepared wherein the 11-aza nitrogen atom was quaternized as either the corresponding N-oxide or methyl iodide. Compounds with quaternized N-11 showed improved water solubility and were equipotent to the clinically investigated camptothecin analog topotecan in the cleavable complex assay. These compounds were evaluated in vivo in nude mice bearing HT-29 human colon carcinoma xenografts. The analog 11-aza-(20S)-camptothecin 11-N-oxide was found to significantly retard tumor growth when compared to untreated controls. Finally, 7,10-disubstituted 11-azacamptothecin analogs were synthesized using Pd(0) coupling reactions of 10-bromo-7-alkyl-11-aza-(20S)-camptothecins 19 and 20, which in turn were available from a Friedlander condensation of the novel bromopyridine derivatives 17a and 17b with 7. Among the 10-substituted series, a number of analogs displayed extremely high in vitro potency against topoisomerase I and improved aqueous solubility. A significant number of the compounds were found to be active in whole cell cytotoxicity assays and several were evaluated in nude mice bearing the HT-29 tumor xenografts. The most effective of these proved to be (S)-11-aza-7-ethyl-10-(aminohydroximinomethyl)camptothecin trifluoracetic acid salt (27), a potent topoisomerase I inhibitor which demonstrated excellent efficacy in both short term and in extended in vivo assays. A comparison between in vitro enzyme data and in vivo data from nude mouse studies in other compounds in this series revealed a poor overall correlation between topoisomerase inhibition in vitro and antitumor efficacy in vivo.

Scalable molecular dynamics with NAMD

Abstract: NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This article, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

Development and testing of a general amber force field.

Abstract: We describe here a general Amber force field (GAFF) for organic molecules. GAFF is designed to be compatible with existing Amber force fields for proteins and nucleic acids, and has parameters for most organic and pharmaceutical molecules that are composed of H, C, N, O, S, P, and halogens. It uses a simple functional form and a limited number of atom types, but incorporates both empirical and heuristic models to estimate force constants and partial atomic charges. The performance of GAFF in test cases is encouraging. In test I, 74 crystallographic structures were compared to GAFF minimized structures, with a root-mean-square displacement of 0.26 A, which is comparable to that of the Tripos 5.2 force field (0.25 A) and better than those of MMFF 94 and CHARMm (0.47 and 0.44 A, respectively). In test II, gas phase minimizations were performed on 22 nucleic acid base pairs, and the minimized structures and intermolecular energies were compared to MP2/6-31G* results. The RMS of displacements and relative energies were 0.25 A and 1.2 kcal/mol, respectively. These data are comparable to results from Parm99/RESP (0.16 A and 1.18 kcal/mol, respectively), which were parameterized to these base pairs. Test III looked at the relative energies of 71 conformational pairs that were used in development of the Parm99 force field. The RMS error in relative energies (compared to experiment) is about 0.5 kcal/mol. GAFF can be applied to wide range of molecules in an automatic fashion, making it suitable for rational drug design and database searching.

All-atom empirical potential for molecular modeling and dynamics studies of proteins.

Abstract: New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent−solvent, solvent−solute, and solute−solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volume...

Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids

Abstract: The parametrization and testing of the OPLS all-atom force field for organic molecules and peptides are described. Parameters for both torsional and nonbonded energetics have been derived, while the bond stretching and angle bending parameters have been adopted mostly from the AMBER all-atom force field. The torsional parameters were determined by fitting to rotational energy profiles obtained from ab initio molecular orbital calculations at the RHF/6-31G*//RHF/6-31G* level for more than 50 organic molecules and ions. The quality of the fits was high with average errors for conformational energies of less than 0.2 kcal/mol. The force-field results for molecular structures are also demonstrated to closely match the ab initio predictions. The nonbonded parameters were developed in conjunction with Monte Carlo statistical mechanics simulations by computing thermodynamic and structural properties for 34 pure organic liquids including alkanes, alkenes, alcohols, ethers, acetals, thiols, sulfides, disulfides, a...