Geoffrey P. F. Wood

A restricted-open-shell complete-basis-set model chemistry.

A restricted-open-shell model chemistry based on the complete basis set-quadratic Becke3 (CBS-QB3) model is formulated and denoted ROCBS-QB3. As the name implies, this method uses spin-restricted wave functions, both for the direct calculations of the various components of the electronic energy and for extrapolating the correlation energy to the complete-basis-set limit. These modifications eliminate the need for empirical corrections that are incorporated in standard CBS-QB3 to compensate for spin contamination when spin-unrestricted wave functions are used. We employ an initial test set of 19 severely spin-contaminated species including doublet radicals and both singlet and triplet biradicals. The mean absolute deviation (MAD) from experiment for the new ROCBS-QB3 model (3.6+/-1.5 kJ mol(-1)) is slightly smaller than that of the standard unrestricted CBS-QB3 version (4.8+/-1.5 kJ mol(-1)) and substantially smaller than the MAD for the unrestricted CBS-QB3 before inclusion of the spin correction (16.1+/-1.5 kJ mol(-1)). However, when applied to calculate the heats of formation at 298 K for the moderately spin-contaminated radicals in the G2/97 test set, ROCBS-QB3 does not perform quite as well as the standard unrestricted CBS-QB3, with a MAD from experiment of 3.8+/-1.6 kJ mol(-1) (compared with 2.9+/-1.6 kJ mol(-1) for standard CBS-QB3). ROCBS-QB3 performs marginally better than standard CBS-QB3 for the G2/97 set of ionization energies with a MAD of 4.1+/-0.1 kJ mol(-1) (compared with 4.4+/-0.1 kJ mol(-1)) and electron affinities with a MAD of 3.9+/-0.2 kJ mol(-1) (compared with 4.3+/-0.2 kJ mol(-1)), but the differences in MAD values are comparable to the experimental uncertainties. Our overall conclusion is that ROCBS-QB3 eliminates the spin correction in standard CBS-QB3 with no loss in accuracy.

Nature of Glycine and Its α-Carbon Radical in Aqueous Solution: A Theoretical Investigation

Quantum chemistry calculations and classical molecular dynamics simulations have been used to examine the equilibria in solution between the neutral and zwitterionic forms of glycine and also of the glycyl radical. The established preference (by 30 kJ mol(-1)) for the zwitterion of glycine was confirmed by both the quantum chemical calculations and the classical molecular dynamics simulations. The best agreement with experiment was derived from thermodynamic integration calculations of explicitly solvated systems, which gives a free energy difference of 36.6 ± 0.6 kJ mol(-1). In contrast, for the glycyl radical in solution, the neutral form is preferred, with a calculated free energy difference of 54.8 ± 0.6 kJ mol(-1). A detailed analysis of the microsolvation environments of each species was carried out by evaluating radial distribution functions and hydrogen bonding patterns. This analysis provides evidence that the change in preference between glycine and glycyl radical is due to the inherent gas-phase stability of the neutral α-carbon radical rather than to any significant difference in the solvation behavior of the constituent species.

Secondary Structure Assignment of Amyloid-β Peptide Using Chemical Shifts

The distinct conformational dependence of chemical shifts caused by α-helices and β-sheets renders NMR chemical shift analysis a powerful tool for the structural determination of proteins. However, the time scale of NMR experiments can make a secondary structure assignment of highly flexible peptides or proteins, which may be converting between conformational substates, problematic. For instance the amyloid-β monomer, according to NMR chemical shifts, adopts a predominately random coil structure in aqueous solution (with <3% α-helical content). Molecular dynamics simulations, on the other hand, suggest that α-helical content can be significant (10-25%). In this paper, we explore the possible reasons for this discrepancy and show that the different results from experiments and theory are not necessarily mutually exclusive but may reflect a general problem of secondary structure assignment of conformationally flexible biomolecules.

Reactions of benzene and 3-methylpyrrole with the •OH and •OOH radicals: an assessment of contemporary density functional theory methods.

A high-level quantum chemistry investigation has been carried out for the addition and abstraction reactions by the radicals (•)OH and (•)OOH to and from the model alkenes 3-methylpyrrole and benzene. These models were chosen to reflect the functionalities contained in the side chain of the amino acid tryptophan. The W1BD procedure was used to calculate benchmark barriers and reaction energies for the smaller model system of (•)OOH addition to ethylene. It was found that the CBS-QB3 methodology compares best with the W1BD benchmark, demonstrating a mean absolute deviation (MAD) from W1BD of 3.9 kJ mol(-1). For the reactions involving the (•)OH radical and benzene or 3-methylpyrrole, addition is favored over abstraction in all cases. In particular the CBS-QB3 calculations suggest a barrierless addition reaction of the (•)OH radical to position two of 3-methylpyrrole. For the analogous addition and abstraction reactions involving the (•)OOH radical, the same order of reactivity was found, albeit with higher barriers. A number of other processes involving the addition of the (•)OOH radical were also investigated. The main findings of these studies determined that the initial (•)OOH barrier of stepwise addition to 3-methylpyrrole (+18.8 kJ mol(-1)) is significantly smaller than the concerted addition barrier (+71.5 kJ mol(-1)). This conclusion contrasts starkly with the situation for ethylene in which it is well established that the concerted process has the smaller barrier. A considerable variety of contemporary density functional theory procedures have been tested to examine their accuracy in predicting the CBS-QB3 results. It was found that the best overall performing method was UBMK with an MAD of 7.3 kJ mol(-1). A number of other functionals additionally performed well. They included UM06, RM06, UXYG3 and RXYG3, all of which have MADs of less than 8 kJ mol(-1).

Modeling β-Scission Reactions of Peptide Backbone Alkoxy Radicals: Backbone C-C Bond Fission.

To model the C-C β-scission reactions of backbone peptide alkoxy radicals, enthalpies and barriers for the fragmentation of four substituted alkoxy radicals have been calculated with a variety of ab initio molecular orbital theory and density functional theory procedures. The high-level methods examined include CBS-QB3, variants of the G3 family, and W1. Simpler methods include HF, MP2, QCISD, B3-LYP, BMK, and MPW1K with a range of basis sets. We find that good accuracy can be achieved with the G3(MP2)//B3-LYP and G3X(MP2)-RAD methods. Lower-cost methods producing reasonable results are single-point energy calculations with UB3-LYP/6-311+G(3df,2p), RB3-LYP/6-311+G(3df,2p), UBMK/6-311+G(3df,2p), and RBMK/6-311+G(3df,2p) on geometries optimized with UB3-LYP/6-31G(d) or UBMK/6-31G(d). Heats of formation at 0 K for the alkoxy radicals and their fragmentation products were also calculated. We predict ΔfH0 values for the alkoxy radicals of -71.4 ((•)OCH2CH [Formula: see text] O), -102.5 ((•)OCH(CH3)CH [Formula: see text] O), -176.6 ((•)OCH(CH3)C(NH2) [Formula: see text] O), and -264.6 ((•)OC(CH3)(NHCH [Formula: see text] O)CH [Formula: see text] O) kJ mol(-)(1). For the fragmentation products NH2C((•)) [Formula: see text] O and CH( [Formula: see text] O)NHC(CH3) [Formula: see text] O, we predict ΔfH0 values of -5.9 kJ mol(-)(1) and -352.8 kJ mol(-)(1).

Molecular Investigation of the Mechanism of Non-Enzymatic Hydrolysis of Proteins and the Predictive Algorithm for Susceptibility

A number of potential degradation routes can limit the shelf life of a biotherapeutic. While these are experimentally measurable, the tests to do so require a substantial investment in both time and material, resources rarely available early in the drug development process. To address the potential degradation route of non-enzymatic hydrolysis, we performed a molecular modeling analysis, together with an experimental study, to gain detailed insight into the reaction. On the basis of the mechanism, an algorithm for predicting the likely cleavage sites of a protein has been created. This algorithm measures four key properties during a molecular dynamics simulation, which relate to the key steps of the hydrolysis mechanism, in particular the rate-determining step (which can vary depending on the local environment). The first two properties include the secondary structure and the surface exposure of the amide bond, both of which help detect if the addition of the proton to the amide bond is possible. The second two properties relate to whether the side chain can cyclize and form a furane ring. These two properties are the orientation of the side chain relative to the amide bond and the number of hydrogen bonds between the side chain and the surrounding protein. Overall, the algorithm performs well at identifying reactive versus nonreactive bonds. The algorithm correctly classifies nearly 90% of all amide bonds following an aspartic or glutamic acid residue as reactive or nonreactive.

Mechanistic Insights into Radical-Mediated Oxidation of Tryptophan from ab Initio Quantum Chemistry Calculations and QM/MM Molecular Dynamics Simulations.

An assessment of the mechanisms of (•)OH and (•)OOH radical-mediated oxidation of tryptophan was performed using density functional theory calculations and ab initio plane-wave Quantum Mechanics/Molecular Mechanics (QM/MM) molecular dynamics simulations. For the (•)OH reactions, addition to the pyrrole ring at position 2 is the most favored site with a barrierless reaction in the gas phase. The subsequent degradation of this adduct through a H atom transfer to water was intermittently observed in aqueous-phase molecular dynamics simulations. For the (•)OOH reactions, addition to the pyrrole ring at position 2 is the most favored pathway, in contrast to the situation in the model system ethylene, where concerted addition to the double bond is preferred. From the (•)OOH position 2 adduct QM/MM simulations show that formation of oxy-3-indolanaline occurs readily in an aqueous environment. The observed transformation starts from an initial rupture of the O-O bond followed by a H atom transfer with the accompanying loss of an (•)OH radical to solution. Finally, classical molecular dynamics simulations were performed to equate observed differential oxidation rates of various tryptophan residues in monoclonal antibody fragments. It was found that simple parameters derived from simulation correlate well with the experimental data.

Solvation of the Glycyl Radical

The effect of adding explicit water molecules to the neutral (N) and zwitterionic (Z) forms of the glycyl radical has been examined. The results show that a minimum of three water molecules is required to stabilize the Z radical as a local minimum, with an energy gap of 123 kJ mol-1 between the N and Z forms at this point, in favor of the N form. Increasing the number of water molecules to ∼20 leads to a converged Z-N energy difference of ∼50 kJ mol-1 still in favor of the N form, even though the radical is not considered fully solvated from a structural point of view. Thus, energetic convergence is determined mainly by solvation of the polar functional groups, and a complete coverage of the entire molecule is not necessary. Because aqueous closed-shell glycine exists as a zwitterion while aqueous glycyl radical prefers the neutral form, the conversion between the two necessitates a change along the hydrogen-abstraction reaction pathway. In this regard, the transition structure for α-hydrogen abstraction by the ·OH radical has greater resemblance to glycine than to the glycyl radical. Overall, the barrier for hydrogen abstraction from Z glycine is larger than that from the N isomer, and this might act to provide some protection against radical damage to the free amino acid in the (aqueous) biological environment.

Dynamics and thermodynamics of Ibuprofen conformational isomerism at the crystal/solution interface.

Conformational flexibility of molecules involved in crystal growth and dissolution is rarely investigated in detail and usually considered to be negligible in the formulation of mesoscopic models of crystal growth. In this work, we set out to investigate the conformational isomerism of ibuprofen as it approaches and is incorporated in the morphologically dominant {100} crystal face, in a range of different solvents: water, 1-butanol, toluene, cyclohexanone, cyclohexane, acetonitrile, and trichloromethane. To this end, we combine extensive molecular dynamics and well-tempered metadynamics simulations to estimate the equilibrium distribution of conformers, compute conformer-conformer transition rates, and extract the characteristic relaxation time of the conformer population in solution, adsorbed at the solid/liquid interface, incorporated in the crystal in contact with the mother solution, and in the crystal bulk. We find that, while the conformational equilibrium distribution is weakly dependent on the solvent, relaxation times are instead significantly affected by it. Furthermore, differences in the relaxation dynamics are enhanced on the crystal surface, where conformational transitions become slower and specific conformational transition pathways are hindered. This leads us to observe that the dominant mechanisms of conformational transition can also change significantly moving from the bulk solution to the crystal interface, even for a small molecule with limited conformational flexibility such as ibuprofen. Our findings suggest that understanding conformational flexibility is key to provide an accurate description of the solid/liquid interface during crystal dissolution and growth, and therefore, its relevance should be systematically assessed in the formulation of mesoscopic growth models.