Pseudorotation motion in tetrahydrofuran: an ab initio study

Disentangling the Conformational Space and Structural Preferences of Tetrahydrofurfuryl Alcohol Using Rotational Spectroscopy and Computational Chemistry

The influence of the hydroxymethyl (CH2OH) group on the tetrahydrofuran (THF) ring structure was investigated by disentangling the gas phase conformational landscape of the sugar analogue tetrahydrofurfuryl alcohol (THFA). By combining rotational spectroscopy (6-20 GHz) and quantum chemical calculations, transitions corresponding to two stable conformers of THFA and their 13C isotopologues were observed and assigned in the rotational spectrum. The positions of the C atoms were precisely determined to unambiguously distinguish between nearly isoenergetic pairs of conformers that differ in their ring configurations: envelope (E) versus twist (T). The rotational spectrum confirms that the E ring geometry is favoured when the CH2OH fragment lies gauche (-) to the THF backbone (OCCO ~-60°) whereas the T form is more stable for the gauche (+) alignment of the substituent (OCCO ~+60°). The observed spectral intensities suggest that conformational relaxation of the THF geometry (E↔T) to the more stable form readily occurs within the pairs of g- and g+ conformers which is consistent with the low barriers (1.5-1.7 kJ mol-1) for conversion determined via transition state calculations. Insights into the intramolecular hydrogen bonding and other weak interactions stabilizing the lowest energy structures of THFA were derived and rationalized using non-covalent interaction analyses.

Deciphering the conformational preference and ionization dynamics of tetrahydrofuran: Insights from advanced spectroscopic techniques

Tetrahydrofuran (THF) has garnered significant attention due to its pivotal role in biological and chemical processes. The diverse array of conformations exhibited by THF profoundly impacts its reactivity and interactions with other molecules. Understanding these conformational preferences is crucial for comprehending its molecular behavior. In this study, we utilize infrared (IR) resonant vacuum ultraviolet photoionization/mass-analyzed threshold ionization (VUV-PI/MATI) mass spectroscopies to capture distinctive vibrational spectra of individual conformers, namely, "twisted" and "bent," within THF. Our conformer-specific vibrational spectra provide valuable insights into the relative populations of these two conformers. The analysis reveals that the twisted (C2) conformer is more stable than the bent (CS) conformer by 17 ± 15 cm-1. By precisely tuning the VUV photon energy to coincide with vibrational excitation via IR absorption, we selectively ionize specific conformers, yielding two-photon IR + VUV-PI/MATI spectra corresponding to the twisted and bent conformers. This investigation conclusively affirms that both the twisted and bent conformers coexist in the neutral state, while only the twisted conformer exists in the cationic state. These findings not only bridge gaps in existing knowledge but also provide profound insights into the behavior of this pivotal molecule in the realms of biology and medicine.

Photodissociation dynamics of tetrahydrofuran at 193 nm

Tetrahydrofuran (THF), a cyclic ether with the chemical formula C4H8O, can be considered the simplest analog of the deoxyribose backbone component of deoxyribonucleic acid (DNA). As such, it provides a useful model for probing the photochemistry of such biomolecular motifs. We present a velocity-map imaging study into the ultraviolet dissociation of THF at a wavelength of 193 nm. Excitation to the S1 state occurs via a 3s ← n transition involving a lone-pair electron on the oxygen atom, and has been shown by other authors to result in rapid ring opening via cleavage of one of the C-O bonds to form a ring-opened C4H8O diradical, followed by C-C bond cleavage over a longer timescale to form either OCH2 + C3H6 products (Channel 1a), HOCH2 + C2H5 products (Channel 1b), or OCH2CH2 + C2H4 products (Channel 2). The C2H4O products formed via Channel 2 are unstable on the timescale of our experiment and dissociate further to form CH3 and CHO. We also observe a number of minor products resulting from H or H2 loss from the primary photofragments. The speed distributions observed for all photofragments are broad, indicating excitation of a range of rotational and vibrational states of the products. The angular distributions of the photofragments show an interesting speed dependence: the slowest products have almost isotropic angular distributions, but the magnitude of the recoil anisotropy increases monotonically with photofragment speed. The fastest products exhibit highly anisotropic angular distributions, with the recoil anisotropy parameter β approaching its limiting value of -1 (-0.75 for Channel 1 and -0.5 for Channel 2). This behaviour is attributed to the range of timescales over which the diradical intermediate dissociates into the observed photofragments. Rapid dissociation leads to fast photofragments which retain the correlation between the transition dipole moment for the S1 ← S0 excitation (which lies perpendicular to the ring) and the photofragment velocities (which lie predominantly in the plane of the ring). Slow dissociation results in a high degree of energy redistribution into internal modes, slower photofragments, and loss of correlation between the photofragment velocities and the transition dipole. The higher barrier associated with dissociation via Channel 2 suggests somewhat longer lifetimes for the diradical intermediate and is consistent with a corresponding reduction in the maximum observed value for β.

Potential Energy Surfaces Sampled in Cremer-Pople Coordinates and Represented by Common Force Field Functionals for Small Cyclic Molecules

The complex conformations of the cyclic moieties impact the physical and chemical properties of molecules. In this work, we chose 22 molecules of four-, five-, and six-membered rings and performed a thorough conformational sampling using Cremer-Pople coordinates. With consideration of symmetries, we obtained a total of 1504 conformational structures for four-membered, 5576 for five-membered, and 13509 for six-membered rings. All well-known and many less well-known conformers for each molecule were identified. We represented the potential energy surfaces (PESs) by fitting the data to common analytical force field (FF) functional forms. We found that the general features of PESs can be described by the essential FF functional forms; however, the accuracy of representation can be improved remarkably by including the torsion-bond and torsion-angle coupling terms. The best fit yields R-squared (R²) values close to 1.0 and mean absolute errors in energy less than 0.3 kcal/mol.

Low-temperature heat capacity and pseudorotation in 2 methyltetrahydrofuran

The heat capacity and phase transitions of 2-methyltetrahydrofuran in the temperature range from 7 to 350 K were measured by adiabatic calorimetry. The smoothed molar thermodynamic functions in the condensed...

Bis-Spiro-Oxetane and Bis-Spiro-Tetrahydrofuran Pyrroline Nitroxide Radicals: Synthesis and Electron Spin Relaxation Studies

Synthesis of bis-spiro-oxetane and bis-spiro-tetrahydrofuran pyrroline nitroxide radicals relies on the Mitsunobu reaction-mediated double cyclizations of N-Boc protected pyrroline tetraols. Structures of the nitroxide radicals are supported by X-ray crystallography. In a trehalose/sucrose matrix at room temperature, the bis-spiro-oxetane nitroxide radical possesses electron spin coherence time, Tm ≈ 0.7 μs. The observed enhanced Tm is most likely associated with strong hydrogen bonding of oxetane moieties to the trehalose/sucrose matrix.

Study on the Gas-Phase Reactivity of Charged Pyridynes

The reactivities of three isomeric, charged ortho-pyridynes, the 1,2-, 2,3-, and 3,4-didehydropyridinium cations, were examined in the gas phase using Fourier-transform ion cyclotron resonance (FT-ICR) mass spectrometry. The structures of selected product ions were probed using collision-activated dissociation (CAD) experiments in a linear quadrupole ion trap (LQIT) mass spectrometer. Mechanisms based on quantum chemical calculations are proposed for the formation of all major products. The products of the reactions of the charged ortho-pyridynes in the gas phase were found to closely resemble those formed upon reactions of neutral ortho-arynes in solution, but the mechanisms of these reactions exhibit striking differences. Additionally, no radical reactions were observed for any of the charged ortho-pyridynes examined, in contrast to previous proposals that ortho-benzyne can occasionally react via radical mechanisms. Finally, the relative reactivities of those charged gaseous ortho-pyridynes that yielded similar product distributions were found to be affected mainly by the (calculated) vertical electron affinities of the dehydrocarbon sites, which suggests that the reactivity of these species is controlled by polar effects.

Conformational diversity of the THF molecule in N2 matrix by means of FTIR matrix isolation experiment and Car-Parrinello molecular dynamics simulations

Tetrahydrofuran (THF) is a widely used chemical compound, in particular as a solvent in organic and inorganic synthesis. The THF molecule has also an interesting property, namely, undergoes pseudorotation, similar to the case of the cyclopentane. Low energy difference between the envelope (C_s symmetry) and twisted (C₂ symmetry) conformations of the THF molecule leads to the interconversion between the two conformers. We study the influence of the molecular environment (N₂) on the C_s-C₂ equilibrium of tetrahydrofuran in the THF@N₂ system utilizing nitrogen matrix isolation infrared spectroscopy. We observe a different ratio between envelope (C_s) and twisted (C₂) conformations with respect to a change of the temperature. FTIR experimental studies are supported by the results of the static density functional theory calculations and Car-Parrinello molecular dynamics simulations. We focus on the dynamics of the pseudorotation process, in particular, the lifetime of the THF conformations and their mutual rearrangements. On the basis of the THF@N₂ matrix model, with explicit nitrogen molecules, the anharmonic infrared spectra are generated from the Fourier transformation of the dipole moment autocorrelation function.

VUV excited-state dynamics of cyclic ethers as a function of ring size

The ultrafast vacuum ultraviolet (VUV) Rydberg state dynamics of cyclic ethers reflects ring strain.

Determination of the cationic conformational structure of tetrahydrothiophene by one-photon MATI spectroscopy and Franck–Condon fitting

This study affords more accurate information regarding the pseudorotational PES describing the conformational interconversion of THT.

Single-walled carbon nanotubes in tetrahydrofuran solution: microsolvation from first-principles calculations

The molecular interactions between the commonly used solvent tetrahydrofuran (THF) and single-walled carbon nanotubes (SWCNT) are studied using density functional theory calculations and Car-Parrinello molecular dynamics simulations. The competitive interplay between THF-THF and THF-SWCNT interactions via C-H⋯O and C-H⋯π hydrogen bonds is analyzed in detail. The binding energies for different global and local energy minima configurations of THF monomers, dimers, trimers, and tetramers on SWCNT(10,0) were determined. The adsorbed species are analyzed in terms of their coordination to the surface via weak hydrogen bonds of the C-H⋯π type and in terms of their ability to form intermolecular C-H⋯O hydrogen bonds, which are responsible for the self-aggregation of THF molecules and a possible dimerization or tetramerization process. A special focus is put on the pseudorotation of the THF molecules at finite temperatures and on the formation of blue-shifting hydrogen bonds.

Potential-Energy Surfaces for Ring-Puckering Motions of Flexible Cyclic Molecules through Cremer-Pople Coordinates: Computation, Analysis, and Fitting

Ring-puckering motion in 12 flexible cyclic molecules is investigated by calculation and analysis of two-dimensional potential-energy surfaces (PESs) using the so-called ring-puckering coordinates proposed by Cremer and Pople. The PESs are calculated by means of density-functional theory using a B2PLYP-D3BJ exchange-correlation functional with a maug-cc-pVTZ basis set, and results are compared to the available experimental and theoretical data. Special care is devoted to the aspect of symmetry in such two-dimensional PESs, which are here reported for the first time also for molecules whose planar form has symmetry lower than D_{5 h} or C_{2 v}. The issue of PES fitting and that of solving the nuclear dynamics using ring-puckering coordinates are also addressed. Analytical formulations of the computed PESs using suitable functional forms with a limited set of parameters are provided.

Relation between the NMR data and the pseudorotational free-energy profile for oxolane

The conformation of five-membered furanose rings is a crucial issue for the structural analysis of many biologically-relevant molecules, including DNA and RNA. Oxolane can be treated as a prototypical furanose, composed only of saturated unsubstituted ring. In spite of its structural simplicity, providing the accurate quantitative description of the oxolane conformational features remains a great challenge for both the experimental and theoretical techniques. Here we show the method of recovering the free-energy profiles describing the conformational equilibrium in the oxolane ring (i.e. pseudorotation) based on the experimentally-inferred NMR data ([Formula: see text] coupling constants). The results remain in agreement with the quantum-mechanical-based molecular dynamics simulations and emphasize the large contributions of all ring conformers, even those located at the free-energy barriers. This includes the significant populations of limiting 3T2/2T3 and OE/EO shapes. Our findings provide another example of a poor applicability of the two-state model, which is routinely applied to analyze the NMR data in terms of population of different ring conformers.

The systematic influence of solvent on the conformational features of furanosides

The endo- and exo-anomeric effects are the two most recognizable stereoelectronic effects exhibited by carbohydrates. Their presence relies on the interactions between ring substituent(s) and ring oxygen atoms. Here, we report the finding of a new effect that partially controls the conformational properties of furanose rings and can be ascribed to the influence of the solvent on the electronic structure of a molecule. In contrast to anomeric effects, it is not dependent on either presence or absence of ring substituents. Its origins lie in a solvent-induced flux of atomic charges that involves atoms forming the furanose ring. This systematically changes the energy of the whole molecular system and alters the ring-distortion free energies by ∼2.5-6.5 kJ mol-1, favoring the geometries close to the twist 3T2/2T3 conformers and disfavoring the envelope OE/EO-like shapes. This intriguing effect has never been reported before, although it is expected to exist in all furanose rings. Along with more recognized stereoelectronic effects, this phenomenon contributes to a wide applicability of the two-state (north vs. south) model of pseudorotation in furanosides and, in the case of extremely flexible furanose rings, may change the preferred conformation type in comparison with the gas-phase-oriented predictions.

Exhaustive exploration of the conformational landscape of mono- and disubstituted five-membered rings by DFT and MP2 calculations

The conformational landscape of 22 different non, mono-, and disubstituted compounds with a five-membered ring was thoroughly explored by ab initio (MP2) and DFT (B3LYP and M06-2X) methods with the 6-311+G** basis set. Our results showed that the conformational preference of these compounds was governed mainly by the specific characteristics of the substituents, with a minor influence of the level of theory employed. After a detailed analysis of the computational data, we found an interesting preference of the electronegative substituents to take pseudo-axial positions, whereas alkyl groups preferred adopting the pseudo-equatorial locations. Such preferences were pronounced with MP2 and M06-2X and underestimated by B3LYP. Despite each level of theory affording different landscapes in many cases, as a general trend, we noticed that M06-2X afforded much higher correlation with the MP2 results than B3LYP.

The conformational landscape of 22 different non, mono-, and disubstituted compounds with a five-membered ring was thoroughly explored by ab initio (MP2) and DFT (B3LYP and M06-2X) methods with the 6-311+G** basis set.

Direct Observation of Tetrahydrofuranyl and Tetrahydropyranyl Peroxy Radicals via Cavity Ring-Down Spectroscopy

Room-temperature cavity ring-down (CRD) spectra of the à ← X̃ electronic transition of tetrahydrofuranyl peroxy (THFOO^•) and tetrahydropyranyl peroxy (THPOO^•) radicals were recorded. The peroxy radicals were produced by Cl-initiated oxidation of tetrahydrofuran and tetrahydropyran. Quantum chemical calculations of the lowest-energy conformers of all regioisomers of these two peroxy radicals have been carried out to aid the spectral simulation. Conformational identification and vibrational assignment were achieved by comparing the experimentally obtained spectra to the simulated ones. The absence of α-THPOO^• absorption peaks in the CRD spectrum is attributed to ring opening due to its weak C_α'O bond.

Hydrogen bonding and dominant conformations of hydrated sugar analogue complexes using tetrahydrofurfuryl alcohol as the model sugar molecule

Water molecules, which serve as both hydrogen bond donors and acceptors, have been found to influence the conformational landscape of gas-phase phenyl-β-d-glucopyranoside.

Vibrational spectroscopy of the mass-selected tetrahydrofurfuryl alcohol monomers and its dimers in gas phase using IR depletion and VUV single photon ionization

Tetrahydrofurfuryl alcohol (THFA, C₅H₁₀O₂) is a close chemical analog of the sugar rings present in the phosphate-deoxyribose backbone structure of the nucleic acids. In present report, the infrared (IR) spectra of the size-selected THFA monomer and its dimer have been investigated in a pulsed supersonic jet using infrared-vacuum ultraviolet (VUV) ionization. Herein, the laser light at 118nm wavelength served as the source of "soft" ionization in a time-of-flight mass spectrometer. The IR features for the monomers located at 3622cm^-1 can be assigned to the intramolecular hydrogen bonding stretch vibrations mainly referring to A and C conformers. Compared with the monomer, however, characteristic peaks for the dimer centered at 3415 and 3453cm^-1, red shifted 207 and 169cm^-1, respectively, were associated with the intermolecular hydrogen bonding stretch vibrations. Combined with the quantum-chemical calculations, the dimer in the gas phase preferred cyclic AC conformer stabled by forming two strong intermolecular hydrogen bonds, which shown the high hydrogen bond selectivity in the cluster. The conclusions drawn from the role played in the conformational flexibility by the hydroxyl and ether groups may be extended to other biomolecules.

Modelling charge transfer processes in C2+-tetrahydrofuran collision for ion-induced radiation damage in DNA building blocks

Investigations of collision-induced processes involving carbon ions and molecules of biological interest, in particular DNA building blocks, are crucial to model the effect of radiation on cells in order to improve medical treatments for cancer therapy. Using carbon ions appears to be one of the most efficient ways to increase biological effectiveness to damage cancerous cells by irradiating deep-seated tumors. Therefore, interest in accurate calculations to understand fundamental processes occurring in ion-molecule collision systems has been growing recently. In this context, the charge transfer process in the collision of C²⁺(1s²2s²) ions with the heterocyclic sugar moiety building block tetrahydrofuran (THF) was studied in order to interpret the mechanisms occurring at the molecular level. The molecular structure properties of THF were obtained by means of ab initio quantum chemistry methods. The role of the conformational structure and the orientation of the THF molecule in collision with C²⁺ ions are particularly discussed. Anisotropic effects of the process dynamics in the collision energy ranging from eV to keV by means of semiclassical treatment are also presented and compared to previous experimental and theoretical investigations. A detailed analysis of the obtained cross sections points out an increase in these values by three orders of magnitude by a change of the THF symmetry from C_2v to C_s in collision with C²⁺, which determines a more efficient charge transfer in this case.

Conformational structures of the tetrahydrofuran cation determined using one-photon mass-analyzed threshold ionization spectroscopy

One-photon vacuum ultraviolet mass-analyzed threshold ionization (VUV-MATI) spectroscopy was used to characterize the essential conformations of tetrahydrofuran (THF) and thus determine the stereochemistry of the furanose ring constituting the backbones of DNA and RNA.

Additive and Classical Drude Polarizable Force Fields for Linear and Cyclic Ethers

Empirical force field parameters consistent with the CHARMM additive and classical Drude based polarizable force fields are presented for linear and cyclic ethers. Initiation of the optimization process involved validation of the aliphatic parameters based on linear alkanes and cyclic alkanes. Results showed the transfer to cyclohexane to yield satisfactory agreement with target data; however, in the case of cyclopentane direct transfer of the Lennard-Jones parameters was not sufficient due to ring strain, requiring additional optimization of these parameters for this molecule. Parameters for the ethers were then developed starting with the available aliphatic parameters, with the nonbond parameters for the oxygens optimized to reproduce both gas- and condensed-phase properties. Nonbond parameters for the polarizable model include the use of an anisotropic electrostatic model on the oxygens. Parameter optimization emphasized the development of transferable parameters between the ethers of a given class. The ether models are shown to be in satisfactory agreement with both pure solvent and aqueous solvation properties, and the resulting parameters are transferable to test molecules. The presented force field will allow for simulation studies of ethers in condensed phase and provides a basis for ongoing developments in both additive and polarizable force fields for biological molecules.

Dynamics of tetrahydrofuran as minority component in a mixture with poly(2-(dimethylamino)ethyl methacrylate): A neutron scattering and dielectric spectroscopy investigation

We have investigated a mixture of poly(2-(dimethylamino)ethyl methacrylate) (PDMAEMA) and tetrahydrofuran (THF) (70 wt. % PDMAEMA/30 wt. % THF) by combining dielectric spectroscopy and quasielastic neutron scattering (QENS) on a labelled sample, focusing on the dynamics of the THF molecules. Two independent processes have been identified. The "fast" one has been qualified as due to an internal motion of the THF ring leading to hydrogen displacements of about 3 Å with rather broadly distributed activation energies. The "slow" process is characterized by an Arrhenius-like temperature dependence of the characteristic time which persists over more than 9 orders of magnitude in time. The QENS results evidence the confined nature of this process, determining a size of about 8 Å for the volume within which THF hydrogens' motions are restricted. In a complementary way, we have also investigated the structural features of the sample. This study suggests that THF molecules are well dispersed among side-groups nano-domains in the polymer matrix, ruling out a significant presence of clusters of solvent. Such a good dispersion, together with a rich mobility of the local environment, would prevent cooperativity effects to develop for the structural relaxation of solvent molecules, frustrating thereby the emergence of Vogel-Fulcher-like behavior, at least in the whole temperature interval investigated.

Robustness in the fitting of molecular mechanics parameters

Automated methods for force field parametrization have attracted renewed interest of the community, but the robustness issues associated with the often ill-conditioned nature of parameter optimization have been vastly underappreciated in the recent literature. For this reason, this article offers a detailed description of the origin and nature of these issues. This includes a discussion of the restrained electrostatic potential fit (RESP) charge model, which does contain explicit robustness-enhancing measures albeit not in the context of bonded parameters, and which forms an inspiration for the present work. It is also discussed how all the bonded parameters in a Class I force field can be simultaneously fit using the linear least squares (LLS) procedure, and a novel restraining strategy is presented that overcomes robustness issues in the LLS fitting of bonded parameters while minimally impacting the fitted values of well-behaved parameters. Two variants of this methodology are then validated through a number of case studies, including the fitting of bond-charge increments, which illustrates the method's potential for robustly solving general LLS problems beyond force field parametrization.

Reaction of chlorine radical with tetrahydrofuran: a theoretical investigation on mechanism and reactivity in gas phase

Reaction of chlorine (Cl) radical with heterocyclic saturated ether, tetrahydrofuran has been studied. The detailed reactivity and mechanism of this reaction is analyzed using hybrid density functional theory (DFT), B3LYP and BB1K methods, and aug-cc-pVTZ basis set. To explore the mechanism of the reaction of tetrahydrofuran with Cl radical, four possible sites of hydrogen atom (H) abstraction pathways in tetrahydrofuran were analyzed. The barrier height and rate constants are calculated for the four H-abstraction channels. The BB1K calculated rate constant for α-axial H-abstraction is comparable with the experimentally determined rate constant. It reflects that α-axial H-abstraction is the main degradation pathway of tetrahydrofuran with Cl radical. DFT-based reactivity descriptors are also calculated and these values describe α-axial H-abstraction as the main reaction channel.

Photoelectron and electron momentum spectroscopy of tetrahydrofuran from a molecular dynamical perspective

The results of experimental studies of the valence electronic structure of tetrahydrofuran employing He I photoelectron spectroscopy as well as Electron Momentum Spectroscopy (EMS) have been reinterpreted on the basis of Molecular Dynamical simulations employing the classical MM3 force field and large-scale quantum mechanical simulations employing Born-Oppenheimer Molecular Dynamics in conjunction with the dispersion corrected ωB97XD exchange-correlation functional. Analysis of the produced atomic trajectories demonstrates the importance of thermal deviations from the lowest energy path for pseudorotation, in the form of considerable variations of the ring-puckering amplitude. These deviations are found to have a significant influence on several outer-valence electron momentum distributions, as well as on the He I photoelectron spectrum.

ELECTRON DENSITY ANALYSIS OF WEAK VAN DER WAALS COMPLEXES

The nature of weak van der Waals interactions in different complexes of some atmospheric molecules such as CO 2, N 2 O , and N 2 was examined. Ab initio calculation was carried out at MP2 level of theory using Dunning's aug-cc-pVTZ basis set. Bader's theory of atoms in molecules (AIM) was employed to analyze electron density and to characterize the nature and properties of van der Waals interactions. A set of criteria, having been proposed in the context of AIM theory, was examined for these complexes. In spite of the parameter kinetics energy, per electron density is expected to be greater than unity for closed-shell interactions; we obtained values less than unity for many of these polyatomic systems. A set of limitations has also been outlined for the values of two AIM quantities: total energy density, H(r), and Laplacian of electron density, ∇2ρ, which correspond to different bond natures.

CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields

The widely used CHARMM additive all-atom force field includes parameters for proteins, nucleic acids, lipids, and carbohydrates. In the present article, an extension of the CHARMM force field to drug-like molecules is presented. The resulting CHARMM General Force Field (CGenFF) covers a wide range of chemical groups present in biomolecules and drug-like molecules, including a large number of heterocyclic scaffolds. The parametrization philosophy behind the force field focuses on quality at the expense of transferability, with the implementation concentrating on an extensible force field. Statistics related to the quality of the parametrization with a focus on experimental validation are presented. Additionally, the parametrization procedure, described fully in the present article in the context of the model systems, pyrrolidine, and 3-phenoxymethylpyrrolidine will allow users to readily extend the force field to chemical groups that are not explicitly covered in the force field as well as add functional groups to and link together molecules already available in the force field. CGenFF thus makes it possible to perform "all-CHARMM" simulations on drug-target interactions thereby extending the utility of CHARMM force fields to medicinally relevant systems.