NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Relativistic quark spin coupling effects in the nucleon electromagnetic form factors

We investigate the effect of different forms of relativistic spin coupling of constituent quarks in the nucleon electromagnetic form factors. The four-dimensional integrations in the two-loop Feynman diagram are reduced to the null-plane, such that the light-front wave function is introduced in the computation of the form factors. The neutron charge form factor is very sensitive to different choices of spin coupling schemes, once its magnetic moment is fitted to the experimental value. The scalar coupling between two quarks is preferred by the neutron data, when a reasonable fit of the proton magnetic momentum is found. (C) 2000 Elsevier Science B.V.

Isotropic and anisotropic spin-spin interactions and a quantum phase transition in a dinuclear Cu(II) compound

We report electron-paramagnetic resonance (EPR) studies at similar to 9.5 GHz (X band) and similar to 34 GHz (Q band) of powder and single-crystal samples of the compound Cu(2)[TzTs](4) [N-thiazol-2-yl-toluenesulfonamidatecopper(II)], C(40)H(36)Cu(2)N(8)O(8)S(8), having copper(II) ions in dinuclear units. Our data allow determining an antiferromagnetic interaction J(0)=(-113 +/- 1) cm(-1) (H(ex)=-J(0)S(1)center dot S(2)) between Cu(II) ions in the dinuclear unit and the anisotropic contributions to the spin-spin coupling matrix D (H(ani)=S(1)center dot D center dot S(2)), a traceless symmetric matrix with principal values D/4=(0.198 +/- 0.003) cm(-1) and E/4=(0.001 +/- 0.003) cm(-1) arising from magnetic dipole-dipole and anisotropic exchange couplings within the units. In addition, the single-crystal EPR measurements allow detecting and estimating very weak exchange couplings between neighbor dinuclear units, with an estimated magnitude parallel to J(')parallel to=(0.060 +/- 0.015) cm(-1). The interactions between a dinuclear unit and the ""environment"" of similar units in the structure of the compound produce a spin dynamics that averages out the intradinuclear dipolar interactions. This coupling with the environment leads to decoherence, a quantum phase transition that collapses the dipolar interaction when the isotropic exchange coupling with neighbor dinuclear units equals the magnitude of the intradinuclear dipolar coupling. Our EPR experiments provide a new procedure to follow the classical exchange-narrowing process as a shift and collapse of the line structure (not only as a change of the resonance width), which is described with general (but otherwise simple) theories of magnetic resonance. Using complementary procedures, our EPR measurements in powder and single-crystal samples allow measuring simultaneously three types of interactions differing by more than three orders of magnitude (between 113 cm(-1) and 0.060 cm(-1)).

Long-range spin-qubit interaction mediated by microcavity polaritons

We study the optically induced coupling between spins mediated by polaritons in a planar microcavity. In the strong-coupling regime, the vacuum Rabi splitting introduces anisotropies in the spin coupling. Moreover, due to their photonlike mass, polaritons provide an extremely long spin coupling range. This suggests the realization of two-qubit all-optical quantum operations within tens of picoseconds with spins localized as far as hundreds of nanometers apart.

Computation of (3)J(HH) coupling constants with a combination of density functional theory and semiempirical calculations. Application to complex molecules

This work evaluates the efficiency of economic levels of theory for the prediction of (3)J(HH) spin-spin coupling constants, to be used when robust electronic structure methods are prohibitive. To that purpose, DFT methods like mPW1PW91. B3LYP and PBEPBE were used to obtain coupling constants for a test set whose coupling constants are well known. Satisfactory results were obtained in most of cases, with the mPW1PW91/6-31G(d,p)//B3LYP/6-31G(d,p) leading the set. In a second step. B3LYP was replaced by the semiempirical methods PM6 and RM1 in the geometry optimizations. Coupling constants calculated with these latter structures were at least as good as the ones obtained by pure DFT methods. This is a promising result, because some of the main objectives of computational chemistry - low computational cost and time, allied to high performance and precision - were attained together. (C) 2012 Elsevier B.V. All rights reserved.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP/6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP/6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7 +/- 0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm. (C) 2008 American Institute of Physics.

Solution structure of robustoxin, the lethal neurotoxin from the funnel-web spider Atrax robustus

The solution structure of robustoxin, the lethal neurotoxin from the Sydney funnel-web spider Atrax robustus, has been determined from 2D H-1 NMR data, Robustoxin is a polypeptide of 42 residues cross-linked by four disulphide bonds, the connectivities of which were determined from NMR data and trial structure calculations to be 1-15, 8-20, 14-31 and 16-42 (a 1-4/2-6/3-7/5-8 pattern), The structure consists of a small three-stranded, anti-parallel beta-sheet and a series of interlocking gamma-turns at the C-terminus. It also contains a cystine knot, thus placing it in the inhibitor cystine knot motif family of structures, which includes the omega-conotoxins and a number of plant and animal toxins and protease inhibitors. Robustoxin contains three distinct charged patches on its surface, and an extended loop that includes several aromatic and non-polar residues, Both of these structural features may play a role in its binding to the voltage-gated sodium channel. (C) 1997 Federation of European Biochemical Societies.

Structure determination of the three disulfide bond isomers of alpha-conotoxin GI: A model for the role of disulfide bonds in structural stability

The three possible disulfide bonded isomers of alpha-conotoxin GI have been selectively synthesised and their structures determined by H-1 NMR spectroscopy. alpha-Conotoxin GI derives from the venom of Conus geographus and is a useful neuropharmacological tool as it selectively binds to the nicotinic acetylcholine receptor (nAChR), a ligand-gated ion channel involved in nerve signal transmission. The peptide has the sequence ECCNPACGRHYSC-NH2, and the three disulfide bonded isomers are referred to as GI(2-7;3-13), GI(2-13;3-7) and GI(2-3;7-13). The NMR structure for the native isomer GI(2-7;3-13) is of excellent quality, with a backbone pairwise RMSD of 0.16 Angstrom for a family of 35 structures, and comprises primarily a distorted 3(10),, helix between residues 5 to 11. The two non-native isomers exhibit multiple conformers in solution, with the major populated forms being different in structure both from each other and from the native form. Structure-activity relationships for the native GI(2-7;3-13) as well as the role of the disulfide bonds on folding and stability of the three isomers are examined. It is concluded that the disulfide bonds in alpha-conotoxin GI play a crucial part in determining both the structure and stability of the peptide. A trend for increased conformational heterogeneity was observed in the order of GI(2-7;3-13) < GI(2-13;3-7) < GI(2-3;7-13). It was found that the peptide bond joining Cys2 to Cys3 in GI(2-3;7-13) is predominantly trans, rather than cis as theoretically predicted. These structural data are used to interpret the varying nAChR binding of the non-native forms. A model for the binding of native GI(2-7;3-13) to the mammalian nAChR is proposed, with an alpha-subunit binding face made up of Cys2, Asn4, Pro5, Ala6 and Cys7 and a selectivity face, comprised of Arg9 and His10. These two faces orient the molecule between the alpha and delta subunits of the receptor. The structure of the CCNPAC sequence of the native GI(2-7;3-13) is compared to the structure of the identical sequence from the toxic domain of heat-stable enterotoxins, which forms part of the receptor binding region of the enterotoxins, but which has a different disulfide connectivity. (C) 1998 Academic Press Limited.

Three-dimensional solution structure of alpha-conotoxin MII by NMR spectroscopy: Effects of solution environment on helicity

alpha-Conotoxin MII, a 16-residue polypeptide from the venom of the piscivorous cone snail Conus magus, is a potent and highly specific blocker of mammalian neuronal nicotinic acetylcholine receptors composed of alpha 3 beta 2 subunits. The role of this receptor type in the modulation of neurotransmitter release and its relevance to the problems of addiction and psychosis emphasize the importance of a structural understanding of the mode of interaction of MII with the alpha 3 beta 2 interface. Here we describe the three-dimensional solution structure of MIT determined using 2D H-1 NMR spectroscopy. Structural restraints consisting of 376 interproton distances inferred from NOEs and 12 dihedral restraints derived from spin-spin coupling constants were used as input for simulated annealing calculations and energy minimization in the program X-PLOR. The final set of 20 structures is exceptionally well-defined with mean pairwise rms differences over the whole molecule of 0.07 Angstrom for the backbone atoms and 0.34 Angstrom for all heavy atoms. MII adopts a compact structure incorporating a central segment of alpha-helix and beta-turns at the N- and C-termini. The molecule is stabilized by two disulfide bonds, which provide cross-links between the N-terminus and both the middle and C-terminus of the structure. The susceptibility of the structure to conformational change was examined using several different solvent conditions. While the global fold of MII remains the same, the structure is stabilized in a more hydrophobic environment provided by the addition of acetonitrile or trifluoroethanol to the aqueous solution. The distribution of amino acid side chains in MII creates distinct hydrophobic and polar patches on its surface that may be important for the specific interaction with the alpha 3 beta 2 neuronal nAChR. A comparison of the structure of MII with other neuronal-specific alpha-conotoxins provides insights into their mode of interaction with these receptors.

Role of the 6-20 disulfide bridge in the structure and activity of epidermal growth factor

Two synthetic analogues of murine epidermal. growth factor, [Abu6, 20] mEGF4-48 (where Abu denotes amino-butyric acid) and [G1, M3, K21, H40] mEGF1-48, have been investigated by NMR spectroscopy. [Abu6, 20] mEGF4-48 was designed to determine the contribution of the 6-20 disulfide bridge to the structure and function of mEGF The overall structure of this analogue was similar to that of native mEGF, indicating that the loss of the 6-20 disulfide bridge did not affect the global fold of the molecule. Significant structural differences were observed near the N-terminus, however, with the direction of the polypeptide chain between residues four and nine being altered such that these residues were now located on the opposite face of the main beta-sheet from their position in native mEGF Thermal denaturation experiments also showed that the structure of [Abu6, 20] mEGF4-48 was less stable than that of mEGF. Removal of this disulfide bridge resulted in a significant loss of both mitogenic activity in Balb/c 3T3 cells and receptor binding on A431 cells compared with native mEGF and mEGF4-48, implying that the structural changes in [Abu6, 20] mEGF4-48, although limited to the N-terminus, were sufficient to interfere with receptor binding. The loss of binding affinity probably arose mainly from steric interactions of the dislocated N-terminal region with part of the receptor binding surface of EGF [G1, M3, K21, H40] mEGF1-48 was also synthesized in order to compare the synthetic polypeptide with the corresponding product of recombinant expression. Its mitogenic activity in Balb/c 3T3 cells was similar to that of native mEGF and analysis of its H-1 chemical shifts suggested that its structure was also very similar to native.

Conotoxin TVIIA, a novel peptide from the venom of Conus tulipa - 2. Three-dimensional solution structure

The three-dimensional solution structure of conotoxin TVIIA, a 30-residue polypeptide from the venom of the piscivorous cone snail Conus tulipa, has been determined using 2D H-1 NMR spectroscopy. TVIIA contains six cysteine residues which form a 'four-loop' structural framework common to many peptides from Conus venoms including the omega-, delta-, kappa-, and mu O-conotoxins. However, TVIIA does not belong to these well-characterized pharmacological classes of conotoxins, but displays high sequence identity with conotoxin GS, a muscle sodium channel blocker from Conus geographus. Structure calculations were based on 562 interproton distance restraints inferred from NOE data, together with 18 backbone and nine side-chain torsion angle restraints derived from spin-spin coupling constants. The final family of 20 structures had mean pairwise rms differences over residues 2-27 of 0.18 +/- 0.05 Angstrom for the backbone atoms and 1.39 +/- 0.33 Angstrom for all heavy atoms. The structure consists of a triple-stranded, antiparallel beta sheet with +2x, -1 topology (residues 7-9, 16-20 and 23-27) and several beta turns. The core of the molecule is formed by three disulfide bonds which form a cystine knot motif common to many toxic and inhibitory polypeptides. The global fold, molecular shape and distribution of amino-acid sidechains in TVIIA is similar to that previously reported for conotoxin GS, and comparison with other four-loop conotoxin structures provides further indication that TVIIA and GS represent a new and distinct subgroup of this structural family. The structure of TVIIA determined in this study provides the basis for determining a structure-activity relationship for these molecules and their interaction with target receptors.

Stereochemistry of cyclopentane derivatives from (2,3)J(CH) dependence on dihedral angle (theta H-C-C-X)

The (2,3)J(CH) dependence on dihedral angle (theta H-C-C-X) for cyclopentane derivatives was investigated. We observed that the combined use of experimentally obtained (2,3)J(CH) values and the theoretically determined dihedral angles between the corresponding nuclei can be used to infer the relative stereochemistry of the ring substituents in cyclopentane derivatives. There is a good correlation between the magnitude of (3)J(CH) and the dihedral angle between the hydrogen and the coupled carbon (R-2 = 0.88). Copyright (C) 2008 John Wiley & Sons, Ltd.

Linearization of a naturally occurring circular protein maintains structure but eliminates hemolytic activity

Cyclotides are a recently discovered family of disulfide rich proteins from plants that contain a circular protein backbone. They are exceptionally stable, as exemplified by their use in native medicine of the prototypic cyclotide kalata B1. The peptide retains uterotonic activity after the plant from which it is derived is boiled to make a medicinal tea. The circular backbone is thought to be in part responsible for the stability of the cyclotides, and to investigate its role in determining structure and biological activity, an acyclic derivative, des-(24-28)-kalata B1, was chemically synthesized and purified. This derivative has five residues removed from the 29-amino acid circular backbone of kalata B1 in a loop region corresponding to a processing site in the biosynthetic precursor protein. Two-dimensional NMR spectra of the peptide were recorded, assigned, and used to identify a series of distance, angle, and hydrogen bonding restraints. These were in turn used to determine a representative family of solution structures. Of particular interest was a determination of the structural similarities and differences between des-(2428)-kalata B1 and native kalata B1. Although the overall three-dimensional fold remains very similar to that of the native circular protein, removal of residues 24-28 of kalata B1 causes disruption of some structural features that are important to the overall stability. Furthermore, loss of hemolytic activity is associated with backbone truncation and linearization.

Desulfovibrio gigas ferredoxin II: redox structural modulation of the [3Fe–4S] cluster

J Biol Inorg Chem (2006) 11: 307–315 DOI 10.1007/s00775-005-0077-2

Fases e criticalidade no modelo ashkin - teller de três cores

The usual Ashkin-Teller (AT) model is obtained as a superposition of two Ising models coupled through a four-spin interaction term. In two dimension the AT model displays a line of fixed points along which the exponents vary continuously. On this line the model becomes soluble via a mapping onto the Baxter model. Such richness of multicritical behavior led Grest and Widom to introduce the N-color Ashkin-Teller model (N-AT). Those authors made an extensive analysis of the model thus introduced both in the isotropic as well as in the anisotropic cases by several analytical and computational methods. In the present work we define a more general version of the 3-color Ashkin-Teller model by introducing a 6-spin interaction term. We investigate the corresponding symmetry structure presented by our model in conjunction with an analysis of possible phase diagrams obtained by real space renormalization group techniques. The phase diagram are obtained at finite temperature in the region where the ferromagnetic behavior is predominant. Through the use of the transmissivities concepts we obtain the recursion relations in some periodical as well as aperiodic hierarchical lattices. In a first analysis we initially consider the two-color Ashkin-Teller model in order to obtain some results with could be used as a guide to our main purpose. In the anisotropic case the model was previously studied on the Wheatstone bridge by Claudionor Bezerra in his Master Degree dissertation. By using more appropriated computational resources we obtained isomorphic critical surfaces described in Bezerra's work but not properly identified. Besides, we also analyzed the isotropic version in an aperiodic hierarchical lattice, and we showed how the geometric fluctuations are affected by such aperiodicity and its consequences in the corresponding critical behavior. Those analysis were carried out by the use of appropriated definitions of transmissivities. Finally, we considered the modified 3-AT model with a 6-spin couplings. With the inclusion of such term the model becomes more attractive from the symmetry point of view. For some hierarchical lattices we derived general recursion relations in the anisotropic version of the model (3-AAT), from which case we can obtain the corresponding equations for the isotropic version (3-IAT). The 3-IAT was studied extensively in the whole region where the ferromagnetic couplings are dominant. The fixed points and the respective critical exponents were determined. By analyzing the attraction basins of such fixed points we were able to find the three-parameter phase diagram (temperature £ 4-spin coupling £ 6-spin coupling). We could identify fixed points corresponding to the universality class of Ising and 4- and 8-state Potts model. We also obtained a fixed point which seems to be a sort of reminiscence of a 6-state Potts fixed point as well as a possible indication of the existence of a Baxter line. Some unstable fixed points which do not belong to any aforementioned q-state Potts universality class was also found