Extracting Ramachandran torsional angle distributions from 2D NMR data using Tikhonov regularization /

Solid state nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for studying structural and dynamical properties of disordered and partially ordered materials, such as glasses, polymers, liquid crystals, and biological materials. In particular, twodimensional( 2D) NMR methods such as ^^C-^^C correlation spectroscopy under the magicangle- spinning (MAS) conditions have been used to measure structural constraints on the secondary structure of proteins and polypeptides. Amyloid fibrils implicated in a broad class of diseases such as Alzheimer's are known to contain a particular repeating structural motif, called a /5-sheet. However, the details of such structures are poorly understood, primarily because the structural constraints extracted from the 2D NMR data in the form of the so-called Ramachandran (backbone torsion) angle distributions, g{^,'4)), are strongly model-dependent. Inverse theory methods are used to extract Ramachandran angle distributions from a set of 2D MAS and constant-time double-quantum-filtered dipolar recoupling (CTDQFD) data. This is a vastly underdetermined problem, and the stability of the inverse mapping is problematic. Tikhonov regularization is a well-known method of improving the stability of the inverse; in this work it is extended to use a new regularization functional based on the Laplacian rather than on the norm of the function itself. In this way, one makes use of the inherently two-dimensional nature of the underlying Ramachandran maps. In addition, a modification of the existing numerical procedure is performed, as appropriate for an underdetermined inverse problem. Stability of the algorithm with respect to the signal-to-noise (S/N) ratio is examined using a simulated data set. The results show excellent convergence to the true angle distribution function g{(j),ii) for the S/N ratio above 100.

Application of Reptation Quantum Monte Carlo and Related Methods to LiH

This work investigates mathematical details and computational aspects of Metropolis-Hastings reptation quantum Monte Carlo and its variants, in addition to the Bounce method and its variants. The issues that concern us include the sensitivity of these algorithms' target densities to the position of the trial electron density along the reptile, time-reversal symmetry of the propagators, and the length of the reptile. We calculate the ground-state energy and one-electron properties of LiH at its equilibrium geometry for all these algorithms. The importance sampling is performed with a single-determinant large Slater-type orbitals (STO) basis set. The computer codes were written to exploit the efficiencies engineered into modern, high-performance computing software. Using the Bounce method in the calculation of non-energy-related properties, those represented by operators that do not commute with the Hamiltonian, is a novel work. We found that the unmodified Bounce gives good ground state energy and very good one-electron properties. We attribute this to its favourable time-reversal symmetry in its target density's Green's functions. Breaking this symmetry gives poorer results. Use of a short reptile in the Bounce method does not alter the quality of the results. This suggests that in future applications one can use a shorter reptile to cut down the computational time dramatically.

FT-IR and MAS NMR analysis of montmorillonite K10 supported MF2 reagents and their activity as catalysis /

ZnF2, CdF2, and CUF2 have been adsorbed onto the surface of montmorillonite K10, and the infrared and 19F, 27 AI, and 29Si MAS NMR spectra of the reagents over a range of loadings have been obtained. CUF2 was observed to attack the Si02 layer and form the complex CuSiF6, Zn F2 tends to attack the aluminium oxide layer, in which Zn isomorphously replaces AI, and forms AIF3 and AIF4 - complexes. All the spectroscopic evidence ruled out the formation of any AI-F and/or Si-F free species as CdF2 is adsorbed on the surface of montmorillonite K10. The reactivity of MF2-K10 reagents towards Friedel-Crafts benzylation of benzene with benzyl chloride varied from one reagent to another. ZnF2-K10 was observed to be the most reactive and CUF2 was the least reactive.

The investigation of the reduction mechanism of the 5a-reductase enzyme of Penicillium decumbens /

The 5a-reductase of Penicillium decumbens ATCC 10436 was used as a model for the mammalian enzyme to investigate the mechanism of reduction of testosterone to 5adihydrotestosterone . The purpose of this study was to search for specific 5a-reductase inhibitors which antagonize prostate cancer . In a whole-cell biotransformation mode, this organism reduced testosterone (1) to 5a-dihydrosteroids (8) and 5aandrostane- 3, 17-dione (9) in yields of 28% and 37% respectively. Control experiments have shown that 5aandrostane- 3, 17-dione (9) can be produced from the corresponding alcohol (8) in a subsequent reaction separate from that catalysed by the 5a-reductase enzyme . Androst-4- ene-3, 17-dione (2) is reduced to give only (9) with a recovery of 80% The stereochemistry of the reduction was determined by 500 MHz ^H NMR analysis of the products resulting from the deuterium labelled substrates. The results were obtained by an analysis of the NOE difference spectra, double-quantum filtered phase sensitive COSY 2-D spectra, and ^^c-Ir 2-D shift correlation spectra of deuterium labelled products. According to the unambiguous assignment of the signals due to H-4a and H-4Ii in 5a-dihydro steroids, the NMR data show clearly that addition of hydrogen to the 4{5)K bond has occurred in a trans manner at positions 413 and 5a. To Study the reduction mechanism of this enzyme, several substrates were prepared as following; 3-methyleneandrost-4-en- 17fi-ol(3), androst-4-en-17i5-ol(5) , androst-4-en-3ii, 17fi-diol (6) and 4, 5ii-epoxyandrostane-3, 17-dione (7) . Results suggest that this enzyme system requires an oxygen atom at the 3-position of the steroid in order to bind the substrate. Furthermore, the mechanism of this 5a-reductase may proceed via direct addition of hydrogen at the 4,5 position without involvement of a carbonyl group as an intermediate.

Chelated fluoroboron cations : synthesis and multinuclear nmr studies

The preparation of chelated difluoroboron cations (DD)BF2+, where DD is a saturated polydentate tertiary-amine or polydentate aromatic ligand, has been systematically studied by using multinuclear solution and solid state nuclear magnetic resonance spectroscopy and fast atom bombardment mass spectrometry. Three new methods of synthesis of (DD)BF2+ cations are reported, and compared with the previous method of reacting a chelating donor with Et20.BF3. The methods most effective for aromatic donors such as 1,1O-phenanthroline are ineffective for saturated polydentate tertiary-amines like N,N,N' ,Nil ,Nil-pentamethyldiethylenetriamine. Polydentate tertiary-amine donors that form 5-membered rings upon bidentate chelation were found to chelate effectively when the BF2 source contained two leaving groups (a heavy halide and a Lewis base such as pyridine =pyr or isoxazole =ISOX), i.e., pyr.BF2X (X = CI or Br), ISOX.BF2X and (pyr)2BF2+. Those that would form 6membered rings upon chelation do not chelate by any of the four methods. Polydentate aromatic ligands chelate effectively when the BF2 source contained a weak Lewis base, e.g., ISOX.BF3, ISOX.BF2X and Et20.BF3. Bidentate chelation by polydentate tertiaryamine and aromatic donors leads to nmr parameters that are significantly different then their (D)2BF2+ relatives (D =monod~ntate t-amines or pyridines). The chelated haloboron cations (DD)BFCI+, and (DD)BFBr+ were generated from D.BFX2 adducts for all ligands that form BF2+ cations above. In addition, the (DD)BCI2+ and (DD)BBr2+ cations were formed from D.BX3 adducts by the chelating aromatic ligands, except for the aromatic ligand 1,8-bis(dimethylamino)naphthalene, which formed only the (DD)BF2+ cation, apparently due to its extreme steric hindrance. Chelation by a donor is a two-step reaction. For polydentate tertiary-amine ligands, the two rates appear to be very dependent on the two possible leaving groups on the central boron atom. The order of increasing ease of displacement for the donors was: pyr < Cl < Br < ISOX. The rate of chelation by polydentate aromatic ligands appears to be dependent on the displacement of the first ligand from the boron. The order of increasing ease of displacement for the donors was: pyr < CI < ISOX ~ Br < Et20.

NMR studies of mixed tetrahaloborates and some related boron trihalide complexes /|nGary John Schrobilgen. -- 260 St. Catharines, Ont. : [s. n.],

Nuclear magnetic resonance spectroscopy has been used to study donor-acceptor complexes of boron trifluoride with several ureas, tetramethylthiourea, tetramethylselenourea, and tetramethylquanidine as well as adducts of tetramethyl- -urea with BF2Cl, BFC1 2 , and BC1 3 - A large number of mixed tetrahaloborate ions, including some of the ternary ones such as BF2CIBr-,have been obtained by ligand exchange reactions and studied by NMR techniques. The bonding in these ions is of the same inherent interest as the bonding in the isoelectronic tetrahalomethanes which have been the subject of many detailed studies and have been involved in a controversy concerning the existence of and the nature of "fluorine hyperconjugation" or C-F P1T- Pn bonding_ Ligand exchange reactions also gave rise to the difluoroboron cation, (TMU)20BF2+o The difluoroboron cation has been observed in solutions of TMU-BF3 , and has been proposed as a possible intermediate for fluorine exchange reactions in BF3 adducts.

NMR studies of strong hydrogen bond to fluoride ion

Systems such as MF/diol (M = alkali metal) and }1F/carboxylic acid were subjected to IH, I9F and 13C nmr study to investigate the nature of the very strong H-bonding of fluoride ions with these systems. Evidence indicates a strong H-bond in diol-fluoride systems (~H ~ -(56) kJ mol-I) which is stronger than most 'typical' H-bonds (~H = -(12-40) kJ mol-I), but weaker than that reported for carboxylic acid-fluoride systems (~H ~ -(120) kJ mol-I). Approximate fluoride H-bonded shifts (o(OH)OHF) were evaluated for MF/diol systems from IH chemical shift measurements. No direct correlation was observed between I9F chemical shift and H-bond strength. Thermodynamic parameters were calculated from temperature dependent IH and 19F shifts. Preliminary studies of BUn 4NF-acetylacetone by I9F nmr were conducted at low temperatures and a possible Jmax (ca. 400 Hz) is reported for the fluoride ion H-bonded to acetylacetone. Highfield shift for non-protonated carbons and downfield shift for protonated carbons were observed in carboxylic acid/KF systems. Significant decreas$in I3C TI due to strong H-bonding to fluoride ions were also detected in both diol and carboxylic acid systems. Anomalous results were obtained, such as increasing NOE with increasing temperature in neat 1,2-ethanediol (values above the theoretical maximum of 1.988) and in 1,2-ethanediol/KF. The large 13C NOE's for carboxy carbons in neat carboxylic acids which are. further enhanced by the addition of KF are also unusual.

NMR studies of the exchange reactions of CH3CN.BX3 with excess CH3CN /|nJoseph Fogelman. -- 260 St. Catharines, Ont. : [s. n.],

Exchange reactions between molecular complexes and excess acid or base are well known and have been extensively surveyed in the literature(l). Since the exchange mechanism will, in some way involve the breaking of the labile donor-acceptor bond, it follows that a discussion of the factors relating to bonding in molecular complexes will be relevant. In general, a strong Lewis base and a strong Lewis acid form a stable adduct provided that certain stereochemical requirements are met. A strong Lewis base has the following characteristics (1),(2) (i) high electron density at the donor site. (ii) a non-bonded electron pair which has a low ionization potential (iii) electron donating substituents at the donor atom site. (iv) facile approach of the site of the Lewis base to the acceptor site as dictated by the steric hindrance of the substituents. Examples of typical Lewis bases are ethers, nitriles, ketones, alcohols, amines and phosphines. For a strong Lewis acid, the following properties are important:( i) low electron density at the acceptor site. (ii) electron withdrawing substituents. (iii) substituents which do not interfere with the close approach of the Lewis base. (iv) availability of a vacant orbital capable of accepting the lone electron pair of the donor atom. Examples of Lewis acids are the group III and IV halides such (M=B, AI, Ga, In) and MX4 - (M=Si, Ge, Sn, Pb). The relative bond strengths of molecular complexes have been investigated by:- (i) (ii) (iii) (iv) (v] (vi) dipole moment measurements (3). shifts of the carbonyl peaks in the IIIR. (4) ,(5), (6) .. NMR chemical shift data (4),(7),(8),(9). D.V. and visible spectrophotometric shifts (10),(11). equilibrium constant data (12), (13). heats of dissociation and heats of reactions (l~), (16), (17), (18), (19). Many experiments have bben carried out on boron trihalides in order to determine their relative acid strengths. Using pyridine, nitrobenzene, acetonitrile and trimethylamine as reference Lewis bases, it was found that the acid strength varied in order:RBx3 > BC1 3 >BF 3 • For the acetonitrile-boron trihalide and trimethylamine boron trihalide complexes in nitrobenzene, an-NMR study (7) showed that the shift to lower field was. greatest for the BB~3 adduct ~n~ smallest for the BF 3 which is in agreement with the acid strengths. If electronegativities of the substituents were the only important effect, and since c~ Br ,one would expect the electron density at the boron nucleus to vary as BF3

On the path integral formulation and the evaluation of quantum statistical averages

Four problems of physical interest have been solved in this thesis using the path integral formalism. Using the trigonometric expansion method of Burton and de Borde (1955), we found the kernel for two interacting one dimensional oscillators• The result is the same as one would obtain using a normal coordinate transformation, We next introduced the method of Papadopolous (1969), which is a systematic perturbation type method specifically geared to finding the partition function Z, or equivalently, the Helmholtz free energy F, of a system of interacting oscillators. We applied this method to the next three problems considered• First, by summing the perturbation expansion, we found F for a system of N interacting Einstein oscillators^ The result obtained is the same as the usual result obtained by Shukla and Muller (1972) • Next, we found F to 0(Xi)f where A is the usual Tan Hove ordering parameter* The results obtained are the same as those of Shukla and Oowley (1971), who have used a diagrammatic procedure, and did the necessary sums in Fourier space* We performed the work in temperature space• Finally, slightly modifying the method of Papadopolous, we found the finite temperature expressions for the Debyecaller factor in Bravais lattices, to 0(AZ) and u(/K/ j,where K is the scattering vector* The high temperature limit of the expressions obtained here, are in complete agreement with the classical results of Maradudin and Flinn (1963) .

Optimization of trial wave functions for use in quantum Monte Carlo with application to LiH

Methods for both partial and full optimization of wavefunction parameters are explored, and these are applied to the LiH molecule. A partial optimization can be easily performed with little difficulty. But to perform a full optimization we must avoid a wrong minimum, and deal with linear-dependency, time step-dependency and ensemble-dependency problems. Five basis sets are examined. The optimized wavefunction with a 3-function set gives a variational energy of -7.998 + 0.005 a.u., which is comparable to that (-7.990 + 0.003) 1 of Reynold's unoptimized \fin ( a double-~ set of eight functions). The optimized wavefunction with a double~ plus 3dz2 set gives ari energy of -8.052 + 0.003 a.u., which is comparable with the fixed-node energy (-8.059 + 0.004)1 of the \fin. The optimized double-~ function itself gives an energy of -8.049 + 0.002 a.u. Each number above was obtained on a Bourrghs 7900 mainframe computer with 14 -15 hrs CPU time.

Quantum Monte Carlo : some theoretical and numerical studies

In Part I, theoretical derivations for Variational Monte Carlo calculations are compared with results from a numerical calculation of He; both indicate that minimization of the ratio estimate of Evar , denoted EMC ' provides different optimal variational parameters than does minimization of the variance of E MC • Similar derivations for Diffusion Monte Carlo calculations provide a theoretical justification for empirical observations made by other workers. In Part II, Importance sampling in prolate spheroidal coordinates allows Monte Carlo calculations to be made of E for the vdW molecule var He2' using a simplifying partitioning of the Hamiltonian and both an HF-SCF and an explicitly correlated wavefunction. Improvements are suggested which would permit the extension of the computational precision to the point where an estimate of the interaction energy could be made~

Simultaneous extraction of order parameters and orientational distribution fuctions from 31[superscript]P NMR spectra of magnetically partially oriented phospholipid bilayers

Order parameter profiles extracted from the NMR spectra of model membranes are a valuable source of information about their structure and molecular motions. To al1alyze powder spectra the de-Pake-ing (numerical deconvolution) ~echnique can be used, but it assumes a random (spherical) dist.ribution of orientations in the sample. Multilamellar vesicles are known to deform and orient in the strong magnetic fields of NMR magnets, producing non-spherical orientation distributions. A recently developed technique for simultaneously extracting the anisotropies of the system as well as the orientation distributions is applied to the analysis of partially magnetically oriented 31p NMR spectra of phospholipids. A mixture of synthetic lipids, POPE and POPG, is analyzed to measure distortion of multilamellar vesicles in a magnetic field. In the analysis three models describing the shape of the distorted vesicles are examined. Ellipsoids of rotation with a semiaxis ratio of about 1.14 are found to provide a good approximation of the shape of the distorted vesicles. This is in reasonable agreement with published experimental work. All three models yield clearly non-spherical orientational distributions, as well as a precise measure of the anisotropy of the chemical shift. Noise in the experimental data prevented the analysis from concluding which of the three models is the best approximation. A discretization scheme for finding stability in the algorithm is outlined

Solid state NMR chemical shifts as an alternative to diffraction data in the determination of SIC polytypic structures

Silicon carbide, which has many polytypic modifications of a very simple and very symmetric structure, is an excellent model system for exploring, the relationship between chemical shift, long-range dipolar shielding, and crystal structure in network solids. A simple McConnell equation treatment of bond anisotropy effects in a poly type predicts chemical shifts for silicon and carbon sites which agree well with the experiment, provided that contributions from bonds up to 100 A are included in the calculation. The calculated chemical shifts depend on three factors: the layer stacking sequence, electrical centre of gravity, and the spacings between silicon and carbon layers. The assignment of peaks to lattice sites is proved possible for three polytypes (6H, 15R, and 3C). The fact that the calculated chemical shifts are very sensitive to layer spacings provides us a potential way to detennine and refine a crystal structure. In this work, the layer spacings of 6H SiC have been calculated and are within X-ray standard deviations. Under this premise, the layer spacings of 15R have been detennined. 29Si and 13C single crystal nmr studies of 6H SiC polytype indicate that all silicons and carbons are magnetically anisotropic. The relationship between a magnetic shielding tensor component and layer spacings has been derived. The comparisons between experimental and semi-empirical chemical shielding tensor components indicate that the paramagnetic shielding of silicon should be included in the single crystal chemical shift calculation.

Variational Monte Carlo estimation of the dissociation energy of CuH using correlated sampling

A new approach to treating large Z systems by quantum Monte Carlo has been developed. It naturally leads to notion of the 'valence energy'. Possibilities of the new approach has been explored by optimizing the wave function for CuH and Cu and computing dissociation energy and dipole moment of CuH using variational Monte Carlo. The dissociation energy obtained is about 40% smaller than the experimental value; the method is comparable with SCF and simple pseudopotential calculations. The dipole moment differs from the best theoretical estimate by about 50% what is again comparable with other methods (Complete Active Space SCF and pseudopotential methods).