Chiroptical properties of S-proline conformational isomers

The chiroptical properties of S-proline conformational isomers are examined on a theoretical model in which electronic wave functions are obtained from semiempirical molecular orbital calculations. The CNDO/S molecular orbital model is used to perform SCF-MO calculations on ground state electronic structure and excited states are constructed in the virtual orbital-configuration interaction approximation. Electronic rotatory strengths and dipole strengths are calculated directly from the complete (but approximate) molecular electronic wave functions. Zwitterionic, cationic, and anionic S-proline structures are studied twotypes of conformational variables are represented in the calculations: (1) pyrrolidine ring conformation; and (2) rotation about the C^α-COO^? bond. Rotatory strengths are found to be somewhat sensitive to rotational isomerism about the C^α-COO^? bond, but are found to be rather insensitive to conformational changes within the pyrrolidine ring. The CD spectrum of zwitterionic S-proline down to ～160 nm appears to be well accounted for by the theoretically calculated results if conformational preferences with respect to rotation about the C^α-COO^? bond can be assumed to exist in solution media. Furthermore, spectra-structure correlations are offered for the anionic and cationic forms of S-proline in solution.     

Short-term learning in conformational analysis

A method for learning short-term rules of conformational analysis is introduced. The technique works by discovering problems during the building of a conformation in Cartesian coordinate space, and builds an abstract critic suitable for reasoning in abstract symbolic space. The methods not only afford speed increases ranging from 1.0- to 2.3-fold in WIZARD (analysis completed in 100% to 43% of original run time), but can be modified to provide similar increases in other programs that use “template joining” and distance geometry. These methods also provide the basis for a longterm learning project.     

Modeling protein conformational transitions by a combination of coarse-grained normal mode analysis and robotics-inspired methods

Background

Obtaining atomic-scale information about large-amplitude conformational transitions in proteins is a challenging problem for both experimental and computational methods. Such information is, however, important for understanding the mechanisms of interaction of many proteins.

Methods

This paper presents a computationally efficient approach, combining methods originating from robotics and computational biophysics, to model protein conformational transitions. The ability of normal mode analysis to predict directions of collective, large-amplitude motions is applied to bias the conformational exploration performed by a motion planning algorithm. To reduce the dimension of the problem, normal modes are computed for a coarse-grained elastic network model built on short fragments of three residues. Nevertheless, the validity of intermediate conformations is checked using the all-atom model, which is accurately reconstructed from the coarse-grained one using closed-form inverse kinematics.

Results

Tests on a set of ten proteins demonstrate the ability of the method to model conformational transitions of proteins within a few hours of computing time on a single processor. These results also show that the computing time scales linearly with the protein size, independently of the protein topology. Further experiments on adenylate kinase show that main features of the transition between the open and closed conformations of this protein are well captured in the computed path.

Conclusions

The proposed method enables the simulation of large-amplitude conformational transitions in proteins using very few computational resources. The resulting paths are a first approximation that can directly provide important information on the molecular mechanisms involved in the conformational transition. This approximation can be subsequently refined and analyzed using state-of-the-art energy models and molecular modeling methods.

     

Chiroptical properties of phospholipids

The chiroptical properties of phospholipids were investigated by optical rotatory dispersion and circular dichroism. The spectra of phosphatidylcholine varied with the solvent used. The sign of the circular dichroism effects differed in some cases from previously reported. The results show that it is possible to distinguish different types of phospholipids from each other. The results also show that optical rotatory dispersion and circular dichroism might be a very useful tool for the determination of the configuration of various phospholipids.     

NMR-based conformational analysis of sphingomyelin in bicelles

Sphingomyelin (SM) is a common sphingolipid in mammalian membranes and is known to be substantially involved in cellular events such as the formation of lipid rafts. Despite its biological significance, conformation of SM in a membrane environment remains unclear because the noncrystalline property and anisotropic environment of lipid bilayers hampers the application of X-ray crystallography and NMR measurements. In this study, to elucidate the conformation of SM in membranes, we utilized bicelles as a substitute for a lipid bilayer membrane. First, we demonstrated through (31)P NMR, (2)H NMR, and dynamic light scattering experiments that SM forms both oriented and isotropic bicelles by changing the ratio of SM/dihexanoyl phosphatidylcholine. Then, we determined the conformation of SM in isotropic bicelles on the basis of coupling constants and NOE correlations in (1)H NMR and found that the C2-C6 and amide groups of SM take a relatively rigid conformation in bicelles.     

Computational conformational analysis of cyclohexaglycyl

Conformational hyperspace searching procedures have been used in conjunction with molecular mechanics calculations to characterize a series of minimum energy cyclohexaglycyl conformations (MECs). The key MECs are identical to those found by Scheraga m an earlier study using different computational techniques. The gas phase global MEC is stabilized by numerous intramolecular hydrogen bonds between (solvent) inaccessible amide groups; in contrast to the solution and crystalline phase global MEC where all of the amide groups are accessible for solvent interactions and/or intermolecular hydrogen bonding. The conformational hyperspace searching procedures are also extremely useful in protein model building studies.     

Computation techniques in the conformational analysis of carbohydrates

A growing number of modern studies of carbohydrates is devoted to spatial mechanisms of their participation in the cell recognition processes and directed design of inhibitors of these processes. Any progress in this field is impossible without the development of theoretical conformational analysis of carbohydrates. In this review, we generalize literature data on the potentialities of using of different molecular-mechanic force fields, the methods of quantum mechanics, and molecular dynamics to study the conformation of glycoside bond. A possibility of analyzing the reactivity of carbohydrates with the computation techniques is also discussed in brief.     

Theoretical conformational analysis of phospholipids bilayers

We present a computational approach describing the conformation of lipid molecules (1-2-dipalmitoyl-sn-glycero-3 phosphocholine (DPPC)) organized in bilayers. The classical semi-empirical method used in peptide conformational analysis has been extended successfully to lipids. The excellent agreement between our theoretical predictions and recent experimental data on the molecular organization of lipid bilayers suggests that the method could be a valuable tool in the lipid conformational analysis but also in the prediction of orientation and mode of insertion of amphiphilic molecules into the lipid bilayer.     

Computation techniques in the conformational analysis of carbohydrates

A growing number of modern studies of carbohydrates is devoted to spatial mechanisms of their participation in the cell recognition processes and directed design of inhibitors of these processes. No progress in this field is possible without the development of theoretical conformational analysis of carbohydrates. In this review, we generalize literature data on the potentialities of using various molecular-mechanic force fields, the methods of quantum mechanics, and molecular dynamics to study the conformation of glycoside linkage. A possibility of analyzing the reactivity of carbohydrates with the computation techniques is also discussed in brief.     

Chiroptical properties of diamino carboxylic acids

Diamino carboxylic acids have recently come to the attention of scientists working in the field of early life and its development. These are the monomers of a hypothetic early form of genetic material, the so-called Peptide Nucleic Acid (PNA) (Nielson et al., Proc Natl Acad Sci USA 2000;97:3868-3871). Since all biopolymers rely on a specific handedness of their building blocks, the question of symmetry breaking occurs in diamino acids and PNA in the same way as in amino acids and proteins. One possible mechanism for triggering this, is asymmetric photochemistry in interstellar/circumstellar matter by means of circularly polarized light (Bailey et al., Science 2005;281:672-674; Bailey, Orig Life Evol Biosphere 2001;21:167-183; Buscherm?hle, Astrophys J 2005;624:821-826; Meierhenrich, Angew Chem Int Ed Engl 2005;44:5630-5634). Here we have measured the CD-spectra of four chiral diamino carboxylic acids, three of which were found in the Murchison meteorite (Meierhenrich, Proc Natl Acad Sci USA 2004;101:9182-9186). The spectra show a uniform peak at 200 nm. These results and additional quantum mechanical calculations of the involved molecular orbitals support the assumption that the process of symmetry breaking in diamino acids does not depend significantly on the length of the side chain. This means that one process alone could suffice to lead to symmetry breaking in all four measured diamino carboxylic acids and might even to some extent be transferable to monoamino acids, the monomers of proteins.     

Characterization of the conformational domains of bradykinin by computational methods

The AMBER 4.0 force field was used to perform a characterization of the conformational profile of the nonapeptide bradykinin. A thorough conformational search was carried out using molecular dynamics as sampling technique, by computing cycles of high (900 K) and low (300 K) temperature trajectories. A total of 2400 minima were generated and subsequently clustered using the root-mean-square of the backbone dihedral angles as criterium. After the use of a tolerance value of 20deg;, the conformations were clustered in 233 unique conformations with energies up to 40 kcal/mol above the lowest minimum. The analysis of the low-energy conformations indicate that the peptide exhibits a high tendency to adopt a β-turn at the C-terminus and a propensity to adopt a bent structure at the N-terminus. These results are in agreement with the experimental evidence reported in the literature and provide detailed information necessary to understand the conformational preferences of the peptide.     

Solution conformational study of Scyliorhinin I analogues with conformational constraints by two-dimensional NMR and theoretical conformational analysis.

Two analogues of Scyliorhinin I (Scyl), a tachykinin with N-MeLeu in position 8 and a 1,5-disubstituted tetrazole ring between positions 7 and 8, introduced in order to generate local conformational constraints, were synthesized using the solid-phase method. Conformational studies in water and DMSO-d6 were performed on these peptides using a combination of the two-dimensional NMR technique and theoretical conformational analysis. The algorithm of conformational search consisted of the following three stages: (i) extensive global conformational analysis in order to find all low-energy conformations; (ii) calculation of the NOE effects and vicinal coupling constants for each of the low energy conformations; (iii) determining the statistical weights of these conformations by means of a nonlinear least-squares procedure, in order to obtain the best fit of the averaged simulated spectrum to the experimental one. In both solvents the three-dimensional structure of the analogues studied can be interpreted only in terms of an ensemble of multiple conformations. For [MeLeu8]Scyl, the C-terminal 6-10 fragment adopts more rigid structure than the N-terminal one. In the case of the analogue with the tetrazole ring in DMSO-d6 the three-dimensional structure is characterized by two dominant conformers with similar geometry of their backbones. They superimpose especially well (RMSD = 0.28 A) in the 6-9 fragments. All conformers calculated in both solvents superimpose in their C-terminal fragments much better than those of the first analogue. The results obtained indicate that the introduction of the tetrazole ring into the Scyl molecule rigidifies its structure significantly more than that of MeLeu.