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.

A study of interaction in the lowest singlet and triplet states of H2

The potential energy curves of the ground state and the first excited state of H2 are examined in terms of the electronic force acting on each nucleus. The results reveal the detailed course of events that occur when two hydrogen atoms with parallel and antiparallel electron spins approach one another from a large internuclear separation.

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.

Multi-channel curve crossing problems: analytically solvable model

We have proposed a general method for finding the exact analytical solution for the multi-channel curve crossing problem in the presence of delta function couplings. We have analysed the case where aa potential energy curve couples to a continuum (in energy) of the potential energy curves.

Buckled nano rod - a two state system: quantum effects on its dynamics

We consider a suspended elastic rod under longitudinal compression. The compression can be used to adjust potential energy for transverse displacements from the harmonic to the double well regime. The two minima in potential energy curve describe two possible buckled states. Using transition state theory (TST) we have calculated the rate of conversion from one state to other. If the strain epsilon = 4 epsilon c the simple TST rate diverges. We suggest a method to correct this divergence for quantum calculations. We also find that zero point energy contributions can be quite large so that single mode calculations can lead to large errors in the rate.

Dynamics of Barrierless and Activated Chemical Reactions in a Dispersive Medium within the Fractional Diffusion Equation Approach

Barrierless chemical reactions have often been modeled as a Brownian motion on a one-dimensional harmonic potential energy surface with a position-dependent reaction sink or window located near the minimum of the surface. This simple (but highly successful) description leads to a nonexponential survival probability only at small to intermediate times but exponential decay in the long-time limit. However, in several reactive events involving proteins and glasses, the reactions are found to exhibit a strongly nonexponential (power law) decay kinetics even in the long time. In order to address such reactions, here, we introduce a model of barrierless chemical reaction where the motion along the reaction coordinate sustains dispersive diffusion. A complete analytical solution of the model can be obtained only in the frequency domain, but an asymptotic solution is obtained in the limit of long time. In this case, the asymptotic long-time decay of the survival probability is a power law of the Mittag−Leffler functional form. When the barrier height is increased, the decay of the survival probability still remains nonexponential, in contrast to the ordinary Brownian motion case where the rate is given by the Smoluchowski limit of the well-known Kramers' expression. Interestingly, the reaction under dispersive diffusion is shown to exhibit strong dependence on the initial state of the system, thus predicting a strong dependence on the excitation wavelength for photoisomerization reactions in a dispersive medium. The theory also predicts a fractional viscosity dependence of the rate, which is often observed in the reactions occurring in complex environments.