Conformational features of linear and cyclic enkephalins. A computational study

Exploring Conformational Preferences of Leu-enkephalin Using the Conformational Search and Double-Hybrid DFT Energy Calculations

The conformational preferences of Leu-enkephalin (Leu-Enk) were explored by the conformational search and density functional theory (DFT) calculations. By a combination of low-energy conformers of each residue, the initial structures of the neutral Leu-Enk were generated and optimized using the ECEPP3 force field in the gas phase. These structures were reoptimized at the HF/3-21G(d) and M06-2X levels of theory with 6-31G(d) and 6-31+G(d) basis functions. We finally located the 139 structures with the relative energy <10 kcal mol-1 in the gas phase, from which the structures of the corresponding zwitterionic Leu-Enk were generated and reoptimized at the M06-2X/6-31+G(d) level of theory using the implicit solvation model based on density (SMD) in water. The conformational preferences of Leu-Enk were analyzed using Gibbs free energies corrected by single-point energies calculated at the double-hybrid DSD-PBEP86-D3BJ/def2-TZVP level of theory in the gas phase and in water. The neutral Leu-Enk dominantly adopted a folded structure in the gas phase stabilized by three H-bonds with a βII'-bend-like motif at the Gly3-Phe4 sequence and a close contact between the side chains of Phe4 and Leu5. The zwitterionic Leu-Enk exhibited a folded structure in water stabilized by three H-bonds with double β-bends such as a βII' bend at the Gly2-Gly3 sequence and a βI bend at the Gly3-Phe4 sequence. The calculated ensemble-averaged distance between CGly2 α and CLeu5 α of the zwitterionic Leu-Enk in water is consistent with the value estimated from the simulated annealing using the distance constraints derived from nuclear Overhauser effect spectroscopy (NOESY) spectra in water. Interestingly, the preferred conformations of the neutral and zwitterionic Leu-Enk are new folded structures not predicted by earlier computational studies. According to the refined model of the zwitterionic Leu-Enk bound to δ-opioid receptor (δOR), there were favorable interactions of the terminal charged groups of Leu-Enk with the side chains of charged residues of δOR as well as a favorable CAryl···H interaction of the Phe4 residue of Leu-Enk with Trp284 of δOR. Hence, these favorable interactions would induce the folded structure of the zwitterionic Leu-Enk with double β-bends isolated in water into the "bioactive conformation" like an extended structure when binding to δOR.

MOLS sampling and its applications in structural biophysics

This review describes the MOLS method and its applications. This computational method has been developed in our laboratory primarily to explore the conformational space of small peptides and identify features of interest, particularly the minima, i.e., the low energy conformations. A systematic "brute-force" search through the vast conformational space for such features faces the insurmountable problem of combinatorial explosion, whilst other techniques, e.g., Monte Carlo searches, are somewhat limited in their region of exploration and may be considered inexhaustive. The MOLS method, on the other hand, uses a sampling technique commonly employed in experimental design theory to identify a small sample of the conformational space that nevertheless retains information about the entire space. The information is extracted using a technique that is a variant of the self-consistent mean field technique, which has been used to identify, for example, the optimal set of side-chain conformations in a protein. Applications of the MOLS method to understand peptide structure, predict the structures of loops in proteins, predict three-dimensional structures of small proteins, and arrive at the best conformation, orientation, and positions of a small molecule ligand in a protein receptor site have all yielded satisfactory results.

Reaction of naphthalene-2,3-dicarboxaldehyde with enkephalins for LC-fluorescence and LC-MS analysis: conformational studies by molecular modeling and H/D exchange mass spectrometry

A new labeling method compatible with both laser-induced fluorescence (LIF) and MS detection for enkephalins, which uses naphthalene-2,3-dicarboxaldehyde (NDA) and a new nucleophilic agent (N,N-dimethylaminoethanethiol) is described. When the derivative is separated via reverse phase HPLC and detected via MS, two different peaks with similar exact mass but different fluorescence and fragmentation properties are obtained. To interpret these results, molecular modeling and H/D exchange mass spectrometry studies were investigated to test the hypothesis that the peak obtained by LC/LIF/MS analysis depends on the site of protonation of the labeled enkephalins. The peptides labeled with NDA and N,N-dimethylaminoethanethiol were separated on a reverse phase C18 column with a gradient of aqueous 0.1% formic acid and acetonitrile. In mass spectrometry, two peaks are observed with the same exact mass for each molecule while only one peak is detected using fluorescence. Tandem mass spectrometry experiments of ion m/z 809.5 were performed on each chromatographic peak; the first peak (which is not observed by LIF detection) gives a fragment corresponding to the loss of the aminothiol side chain while no fragmentation is observed on the second peak, which was detected by fluorescence. The hypothesis is that each peak represents the labeled enkephalin with different sites of protonation. According to this hypothesis, three fundamental conformations that were closed to the unlabeled leucine-enkephalin were obtained by molecular modeling: a beta-turn like conformation with two hydrogen bonds, a 3(10)-helix with an H bond, and finally, the extended form without any intramolecular interactions. H/D exchange mass spectrometry experiments with D(2)O and d(2-)formic acid as eluent was used to determine which conformation is involved in each peak.

Exploring the conformational space of protein loops using a mean field technique with MOLS sampling

We have recently developed a computational technique that uses mutually orthogonal Latin square sampling to explore the conformational space of oligopeptides in an exhaustive manner. In this article, we report its use to analyze the conformational spaces of 120 protein loop sequences in proteins, culled from the PDB, having the length ranging from 5 to 10 residues. The force field used did not have any information regarding the sequences or structures that flanked the loop. The results of the analyses show that the native structure of the loop, as found in the PDB falls at one of the low energy points in the conformational landscape of the sequences. Thus, a large portion of the structural determinants of the loop may be considered intrinsic to the sequence, regardless of either adjacent sequences or structures, or the interactions that the atoms of the loop make with other residues in the protein or in neighboring proteins.

Conformational study of Met-enkephalin based on the ECEPP force fields

We report a computational study of the small peptide Met-enkephalin based on the ECEPP/2 and ECEPP/3 force fields using the basin paving method. We have located a new global minimum when using the ECEPP/3 force field with peptide angles omega fixed at 180 degrees. With this new result, we can conclude that the lowest energy configurations of Met-enkephalin predicted based on all four versions of ECEPP have a classic gamma-turn centered at residue Gly3 and a beta-turn at residues Gly3-Phe4. However, minor differences between the structures also exist.

Protein structure prediction using mutually orthogonal Latin squares and a genetic algorithm

We combine a new, extremely fast technique to generate a library of low energy structures of an oligopeptide (by using mutually orthogonal Latin squares to sample its conformational space) with a genetic algorithm to predict protein structures. The protein sequence is divided into oligopeptides, and a structure library is generated for each. These libraries are used in a newly defined mutation operator that, together with variation, crossover, and diversity operators, is used in a modified genetic algorithm to make the prediction. Application to five small proteins has yielded near native structures.

Conformational studies on enkephalins using the MOLS technique

Conformational studies of two linear enkephalin molecules, Met-enkephalin and Leu-enkephalin, have been carried using the mutually orthogonal Latin squares (MOLS) technique with the ECEPP/3 force field. This technique was developed recently in our laboratory to perform an unbiased search of the conformational space of peptides and to locate low energy conformations. The present study identified all the folds predicted by other studies, and in addition picked up other energetically favorable structures. The results suggest that the peptide backbone exists as a mixture of folded and unfolded forms (approximately 50% each). The study also provides information on the distribution of the low energy conformations that we have classified on the basis of structural motifs, backbone hydrogen-bonding patterns, and root mean square deviations in atomic positions.

Enantioselective chromatography and molecular modeling of novel aryloxyaminopropan-2-ols with the alkyl carbamate function

A series of different racemic aryloxyaminopropan-2-ol derivatives 1a-d-3a-d with potential beta-adrenergic blocking effects related to propanolol 4 and atenolol 5 was resolved by HPLC using Chiralcel OD-H and Chiralpak AD as chiral stationary phases. Mobile phases consisted of a hexane/alcohol (propan-2-ol or ethanol) mixture doped with a modifier (DEA or TFA). The retention behavior of the compounds depended on the position of the carbamate attached to the aryloxy moiety and on the length of the alkyl residue in the carbamate. Enantiomers of the title compounds were baseline separated with the separation factors alpha and resolutions R(s) varying in the range of 1.34-4.55 and 1.50-10.65, respectively. The chromatographic systems developed can be used for the determination of the enantiomeric purity of the title compounds. Molecular modelling using empirical molecular mechanics and ab initio quantum chemistry methods provided low-energy structures in which sites of potential interactions responsible for retention behavior and chiral recognition could be identified.

Investigation of the complexation of (+)-catechin by beta-cyclodextrin by a combination of NMR, microcalorimetry and molecular modeling techniques

(+)-Catechin is a polyphenolic compound of natural origin that presents anti-oxidant properties of interest for therapeutics or cosmetics uses. Preliminary studies on inclusion into cyclodextrin cavities yielded contradictory results both for the quantitative (affinity constant) and qualitative description of the interaction. By a combination of several experimental and theoretical methods, the present study resolved the previous ambiguities about the interaction between (+)-catechin and beta-cyclodextrin. Thermodynamic data measured by isothermal titration calorimetry demonstrate that the binding is enthalpy driven. Excellent agreement has been obtained for the measurement of the association constant by NMR and microcalorimetry. The several docking modes obtained by systematic docking studies have been compared to intermolecular contacts measured by NMR and the overall geometry of the complex can be proposed.

Enhanced sampling of the molecular potential energy surface using mutually orthogonal latin squares: application to peptide structures

The computational identification of the optimal three-dimensional fold of even a small peptide chain from its sequence, without reference to other known structures, is a complex problem. There have been several attempts at solving this by sampling the potential energy surface of the molecule in a systematic manner. Here we present a new method to carry out the sampling, and to identify low energy conformers of the molecule. The method uses mutually orthogonal Latin squares to select (of the order of) n(2) points from the multidimensional conformation space of size m(n), where n is the number of dimensions (i.e., the number of conformational variables), and m specifies the fineness of the search grid. The sampling is accomplished by first calculating the value of the potential energy function at each one of the selected points. This is followed by analysis of these values of the potential energy to obtain the optimal value for each of the n-variables separately. We show that the set of the n-optimal values obtained in this manner specifies a low energy conformation of the molecule. Repeated application of the method identifies other low energy structures. The computational complexity of this algorithm scales as the fourth power of the size of the molecule. We applied this method to several small peptides, such as the neuropeptide enkephalin, and could identify a set of low energy conformations for each. Many of the structures identified by this method have also been previously identified and characterized by experiment and theory. We also compared the best structures obtained for the tripeptide (Ala)(3) by the present method, with those obtained by an exhaustive grid search, and showed that the algorithm is successful in identifying all the low energy conformers of this molecule.