Ground and Excited State Intramolecular Proton Transfer in Salicylic Acid: an Ab Initio Electronic Structure Investigation

Energetics of the ground and excited state intramolecular proton transfer in salicylic acid have been studied by ab initio molecular orbital calculations using the 6-31G** basis set at the restricted Hartree-Fock (RHF) and configuration interaction-single excitation (CIS) levels and also using the semiempirical method AM1 at the RHF level as well as with single and pair doubles excitation configuration interaction spanning eight frontier orbitals (PECI = 8). The ab initio potential energy profile for intramolecular proton transfer in the ground state reveals a single minimum corresponding to the primary form, in the first excited singlet state, however, there are two minima corresponding to the primary and tautomeric forms, separated by a barrier of similar to 6 kcal/mol, thus accounting for dual emission in salicylic acid. Electron density changes with electronic excitation and tautomerism indicate no zwitterion formation. Changes in spectral characteristics with change in pH, due to protonation and deprotonation of salicylic acid, are also accounted for, qualitatively. Although the AM1 calculations suggest a substantial barrier for proton transfer in the ground as well as the first excited state of SA, it predicts the transition wavelength in near quantitative accord with the experimental results for salicylic acid and its protonated and deprotonated forms.

A comparative study of the energetics, structures, and mechanisms of the HCN H HNC and LiCN <-> LiNC isomerizations

The potential energy surfaces of the HCN<->HNC and LiCN<->LiNC isomerization processes were determined by ab initio theory using fully optimized triple-zeta double polarization types of basis sets. Both the MP2 corrections and the QCISD level of calculations were performed to correct for the electron correlation. Results show that electron correlation has a considerable influence on the energetics and structures. Analysis of the intramolecular bond rearrangement processes reveals that, in both cases, H (or Li+) migrates in an almost elliptic path in the plane of the molecule. In HCN<->HNC, the migrating hydrogen interacts with the in-plane pi,pi* orbitals of CN, leading to a decrease in the C-N bond order. In LiCN<->LiNC, Li+ does not interact with the corresponding pi,pi* orbitals of CN.

A Doubly curved Quadrilateral finite-element for the analysis of laminated anisotropic thin shells of revoution

The details of development of the stiffness matrix for a doubly curved quadrilateral element suited for static and dynamic analysis of laminated anisotropic thin shells of revolution are reported. Expressing the assumed displacement state over the middle surface of the shell as products of one-dimensional first order Hermite polynomials, it is possible to ensure that the displacement state for the assembled set of such elements, is geometrically admissible. Monotonic convergence of total potential energy is therefore possible as the modelling is successively refined. Systematic evaluation of performance of the element is conducted, considering various examples for which analytical or other solutions are available.

Mean-field results of a lattice-gas model of multilayer adsorption

A lattice-gas model of multilayer adsorption has been solved in the mean-field approximation by a different numerical method. Earlier workers obtained a single solution for all values of temperature and pressure. In the present work, multiple solutions have been obtained in certain regions of temperature and pressure which give rise to bysteresis in the adsorption isotherm. In addition, we have obtained a parameter which behaves like an order parameter for the transition. The potential-energy function shows a double minimum in the region of bysteresis and a single maximum elsewhere.

A conformational approach to the study of the dynamics of enzyme inhibition: studies on thermolysin

The preferred conformations of β-phenylpropionyl-Image -phenylalanine (β-PPP) and N-carbobenzoxy-L-phenylalanine (Cbz-Phe), two inhibitors of thermolysin, have been determined by computing potential energy using empirial potential energy functions. Of the 15 to 20 conformations that are favoured for each of these inhibitors only a few have the right conformation to reach the active site of the enzyme. The conformer of β-PPP that initiates binding with the enzyme is different from the bound one, while for Cbz-Phe the bound and initiating conformers are quite similar. Thus, β-PPP favours the ‘induced fit’ model while Cbz-Phe follows the ‘lock and key’ model of binding. The inhibitors differ in their alignment at the active site.

Theoretical studies on the conformations of higher gangliosides

The possible conformations of higher gangliosides (GD3, GT1a. GT1b, GQ1b) have been determined by computing their potential energy using semi-empirical potential functions. The favoured conformation of the disialic acid fragment in these gangliosides is independent of its position (internal or terminal). The favoured conformations of these gangliosides have also been correlated to their biological activity. The results suggest that tetanus toxin and sendai virus may have a large binding site which can accommodate at least four sugar residues.

Theoretical studies on the conformations of higher gangliosides

The possible conformations of higher gangliosides (GD3, GT1a. GT1b, GQ1b) have been determined by computing their potential energy using semi-empirical potential functions. The favoured conformation of the disialic acid fragment in these gangliosides is independent of its position (internal or terminal). The favoured conformations of these gangliosides have also been correlated to their biological activity. The results suggest that tetanus toxin and sendai virus may have a large binding site which can accommodate at least four sugar residues.

X-Pro peptides. A theoretical study of the hydrogen bonded conformations of (alpha-aminoisobutyryl-L -prolyl)n sequences

Intramolecularly hydrogen bonded conformations of (Aib-Pro)n sequences have been analysed theoretically. Both 4-1 (C10 and 3-1 (C7 hydrogen bonded regular structures are shown to be stereochemically feasible. Conformational energies for the helical structures have been estimated using classical potential energy methods. Both C10 and C7 conformations have very similar energies. Pyrrolidine ring puckering has a pronounced effect on the energies, and only Cv-endo puckered Pro residues can be accommodated. The theoretical calculations using spectroscopic data suggest that the recently proposed novel 310 helical conformation for benzyloxycarbonyl(Aib-Pro)4-methyl ester is in solution, is indeed energetically and stereochemically favourable.

Complex carbohydrates: 1. Conformational studies on some oligosaccharides related to N-glycosyl proteins which interact with concanavalin A

Empirical potential energy calculations have been carried out to determine the preferred conformations of some oligosaccharides having the trimannosidic core structure (Man3GlcNAc2) and which interact with concanavalin A. In the minimum energy conformations for the trimannosidic core the mannose residue on the Man α(1–6) arm comes close to one of the N-acetylglucosamine residues of the core. The addition of N-acetylglucosamine residues to the terminal mannose residues does not alter the preferred conformation of the trimannosidic core although it alters the relative preference of some of the higher energy conformations. The minimum energy conformation broadly agrees with available X-ray data. The presence of a bisecting N-acetylglucosamine residue on the middle mannose does not push the trimannosidic core to any new conformation but it does alter the relative preference for a particular conformation.

Theoretical studies on penicillins, penicillin sulphones and their 1-oxa-1-dethia and 1-carba-1-dethia analogues: binding specificities of transeptidases and penicillinases

Empirical potential energy calculations have been carried out to determine the preferred conformations of penicillins and penicillin sulphones and their 1-oxa-1-dethia and 1-carba-1-dethia analogues. With the exception of 1-oxa-1-dethia penicillins, all the other compounds favour C2 and the C3 puckered conformations of their five-membered rings. Replacement of C2 methyl groups by hydrogen atoms as in bisnorpenicillin V or oxidation of sulphur in position 1 as in sulphones, makes the C3 puckered form much less favourable. Addition of an amino-acyl group at the C6 atom, however, makes the C3 puckered form more favoured in penicillin G or V and in 1-carba-1-dethia penicillins. Through the replacement of the sulphur atom at position 1 by an oxygen atom or by a -CH2 group increases the non-planarity of the lactam peptide bond, it significantly affects the relative disposition of the C3 carboxyl group with respect to the β-lactam ring. These conformational differences have been correlated with the biological activities of these compounds. The present study suggests that the conformation of the bicyclic ring system may be more important for initial binding with the crosslinking enzyme(s) involved in the biosynthesis of bacterial cell-wall peptidoglycan and that the mode of binding is influenced by the nature of the side-group at the C6 atom. These studies predict, in agreement with experimental results, that the 1-oxa-1-dethia penicillin nulceus is an inhibitor of penicillianses. The study also suggests that the stereospecificities of the crosslinking enzyme(s) and penicillinases are very similar with regard to the nature of the side-group at the 6 atom and the confirmation of the bicyclic ring system. However, the confirmational requirement for the bicyclic ring system appears to be more specific in the former enzyme than in the latter.

X-Pro peptides. A theoretical study of the hydrogen bonded conformations of (a-aminoisobutyryl-L -prolyl)n sequences

Intramolecularly hydrogen bonded conformations of (Aib-Pro)n sequences have been analysed theoretically. Both 4�1 (C10 and 3�1 (C7 hydrogen bonded regular structures are shown to be stereochemically feasible. Conformational energies for the helical structures have been estimated using classical potential energy methods. Both C10 and C7 conformations have very similar energies. Pyrrolidine ring puckering has a pronounced effect on the energies, and only Cγ-endo puckered Pro residues can be accommodated. The theoretical calculations using spectroscopic data suggest that the recently proposed novel 310 helical conformation for benzyloxycarbonyl(Aib-Pro)4-methyl ester is in solution, is indeed energetically and stereochemically favourable.

A study of interaction between two He+ ions in terms of the electronic forces

The potential energy curve of the He2+2 system dissociating into two He+ ions is examined in terms of the electronic force exerted on each nucleus as a function of the internuclear separation. The results are compared with the process of bond-formation in H2 from the separated atoms.

A Case Study of the Conformation of Poly(alpha-aminoisobutyric acid): alpha- or 310-Helix

The relative stabilities of a- and Blo-helical structures for polymers of a-aminoisobutyric acid (Aib) have been worked out, using the classical potential energy functions. To make a comparative study, we have used Buckingham "6-exp" and Kitaigorodsky's potential functions. Conformational analysis of the dipeptide segment with Aib residue indicates the necessity for nonplanar distortion of the peptide unit, which is a common feature in the observed crystal structures with Aib residues. In the range of Aw -10 to +loo studied, a-helical conformations are preferred in the region -3" < Aw < +loo, and Blo-helical conformations are preferred in the region -3" > Aw > -10'. Minimum energy conformations for right-handed structures are found in the +ue region of Aw and correspondingly for left-handed structures in the -ue region of Aw. For Aw - 6", a-helical structures have four- or near fourfold symmetry with h - 1.5 A. Such a helix with n = 4 and h = 1.5 A is termed an a'-helix. This structure is found to be consistent with the electron diffraction data of Malcolm3 and energetically more favorable than the standard 310-helix.

Intersoliton forces in weak-coupling quantum field theories

We offer a procedure for evaluating the forces exerted by solitons of weak-coupling field theories on one another. We illustrate the procedure for the kink and the antikink of the two-dimensional φ4 theory. To do this, we construct analytically a static solution of the theory which can be interpreted as a kink and an antikink held a distance R apart. This leads to a definition of the potential energy U(R) for the pair, which is seen to have all the expected features. A corresponding evaluation is also done for U(R) between a soliton and an antisoliton of the sine-Gordon theory. When this U(R) is inserted into a nonrelativistic two-body problem for the pair, it yields a set of bound states and phase shifts. These are found to agree with exact results known for the sine-Gordon field theory in those regions where U(R) is expected to be significant, i.e., when R is large compared to the soliton size. We take this agreement as support that our procedure for defining U(R) yields the correct description of the dynamics of well-separated soliton pairs. An important feature of U(R) is that it seems to give strong intersoliton forces when the coupling constant is small, as distinct from the forces between the ordinary quanta of the theory. We suggest that this is a general feature of a class of theories, and emphasize the possible relevance of this feature to real strongly interacting hadrons.

Ab initio molecular dynamics simulations of excited hydrogen halides and methyl halides

Using excited-state ab initio molecular dynamics simulations employing the complete-active-space self-consistent-field approach, we study the mechanism of photodissociation in terms of time evolution of structure, kinetic energy, charges and potential energy for the first excited state of hydrogen halides and methyl halides. Although the hydrogen halides and methyl halides are similar in the photodissociation mechanism, their dynamics are slightly different. The presence of the methyl group causes delay in photodissociation as compared to hydrogen halides.