Characterization of low-energy conformational domains for Met-enkephalin

Transition pathway and its free-energy profile: a protocol for protein folding simulations

We propose a protocol that provides a systematic definition of reaction coordinate and related free-energy profile as the function of temperature for the protein-folding simulation. First, using action-derived molecular dynamics (ADMD), we investigate the dynamic folding pathway model of a protein between a fixed extended conformation and a compact conformation. We choose the pathway model to be the reaction coordinate, and the folding and unfolding processes are characterized by the ADMD step index, in contrast to the common a priori reaction coordinate as used in conventional studies. Second, we calculate free-energy profile as the function of temperature, by employing the replica-exchange molecular dynamics (REMD) method. The current method provides efficient exploration of conformational space and proper characterization of protein folding/unfolding dynamics from/to an arbitrary extended conformation. We demonstrate that combination of the two simulation methods, ADMD and REMD, provides understanding on molecular conformational changes in proteins. The protocol is tested on a small protein, penta-peptide of met-enkephalin. For the neuropeptide met-enkephalin system, folded, extended, and intermediate sates are well-defined through the free-energy profile over the reaction coordinate. Results are consistent with those in the literature.

A simulation strategy for the atomistic modeling of flexible molecules covalently tethered to rigid surfaces: application to peptides

A computational strategy to model flexible molecules tethered to a rigid inert surface is presented. The strategy is able to provide uncorrelated relaxed microstructures at the atomistic level. It combines an algorithm to generate molecules tethered to the surface without atomic overlaps, a method to insert solvent molecules and ions in the simulation box, and a powerful relaxation procedure. The reliability of the strategy has been investigated by simulating two different systems: (i) mixed monolayers consisting of binary mixtures of long-chain alkyl thiols of different lengths adsorbed on a rigid inert surface and (ii) CREKA (Cys-Arg-Glu-Lys-Ala), a short linear pentapeptide that recognizes clotted plasma proteins and selectively homes to tumors, covalently tethered to a rigid inert surface in aqueous solution. In the first, we examined the segregation of the two species in the monolayers using different long-chain:short-chain ratios, whereas in the second, we explored the conformational space of CREKA and ions distribution considering densities of peptides per nm(2) ranging from 0.03 to 1.67. Results indicate a spontaneous segregation in alkyl thiol monolayers, which enhances when the concentration of longest chains increases. However, the whole conformational profile of CREKA depends on the number of molecules tethered to the surface pointing out the large influence of molecular density on the intermolecular interactions, even though the bioactive conformation was found as the most stable in all cases.

Exploring Multidimensional Free Energy Landscapes Using Time-Dependent Biases on Collective Variables

A new implementation of the adaptive biasing force (ABF) method is described. This implementation supports a wide range of collective variables and can be applied to the computation of multidimensional energy profiles. It is provided to the community as part of a code that implements several analogous methods, including metadynamics. ABF and metadynamics have not previously been tested side by side on identical systems. Here, numerical tests are carried out on processes including conformational changes in model peptides and translocation of a halide ion across a lipid membrane through a peptide nanotube. On the basis of these examples, we discuss similarities and differences between the ABF and metadynamics schemes. Both approaches provide enhanced sampling and free energy profiles in quantitative agreement with each other in different applications. The method of choice depends on the dimension of the reaction coordinate space, the height of the barriers, and the relaxation times of degrees of freedom in the orthogonal space, which are not explicitly described by the chosen collective variables.

Application of evolutionary algorithm methods to polypeptide folding: comparison with experimental results for unsolvated Ac-(Ala-Gly-Gly)5-LysH+

We present an evolutionary method for finding the low-energy conformations of polypeptides. The application, called FOLDAWAY,is based on a generic framework and uses several evolutionary operators as well as local optimization to navigate the complex energy landscape of polypeptides. It maintains two complementary representations of the structures and uses the CHARMM force field for evaluating the energies. The method is applied to unsolvated Met-enkephalin and Ac-(Ala-Gly-Gly)(5)-Lys(+)H(+). Unsolvated Ac-(Ala-Gly-Gly)(5)-Lys(+)H(+) has been the object of recent experimental studies using ion mobility measurements. It has a flat energy landscape where helical and globular conformations have similar energies. FOLDAWAY locates several large groups of structures not found in previous molecular dynamics simulations for this peptide, including compact globular conformations, which are probably present in the experiments. However, the relative energies of the different conformations found by FOLDAWAY do not accurately match the relative energies expected from the experimental observations.

Performance of SuSi: A method for generating atomistic models of amorphous polymers based on a random search of energy minima

The performance of a recently developed method to generate representative atomistic models of amorphous polymers has been investigated. This method, which is denoted SuSi, can be defined as a random generator of energy minima. The effects produced by different parameters used to define the size of the system and the characteristics of the generation algorithm have been examined. Calculations have been performed on poly(L,D-lactic) acid (rho = 1.25 g/cm3) and nylon 6 (rho = 1.084 g/cm(3)), which are important commercial polymers.

New conformationally homogeneous beta-turn antagonists of the human B2 kinin receptor

We have designed and synthesized a conformationally homogeneous series of cyclic pentapeptides of the general structure c[Pro-aa(i)-D-Tic-Oic-aa(i + 3)] which adopt a type-II' beta-turn conformation believed important for high affinity antagonism of the bradykinin (BK) B2 receptor. We incorporated D-Tic and octahydroindole-2-carboxylic acid (Oic) residues (present in known active antagonists) in a cyclic pentapeptide that would place the D-aa in the i + 1 position of the beta-turn and a proline as a bridge between the C- and N-termini sides of the turn. In positions i and i + 3 alkyl, aromatic, polar or charged amino acids could be introduced without dramatically changing the overall structure. Ten analogues were studied using 1H nuclear magnetic resonance (NMR) and evaluated for their binding affinity for the human B2 receptor. The NMR data in dimethylsulfoxide (DMSO) confirmed the structural homogeneity within the class and, on the basis of this, one representative member of the series was chosen for a detailed structure determination using NMR data in sodium dodecylsulphate (SDS) micelles and molecular dynamics calculations. Despite the structural similarity, the binding affinity of the ten analogues was strongly influenced by the nature of the side-chains in positions i and i + 3, with the doubly charged analogue 49 (pKi = 6.2) proving best. This compound may serve as the starting point for the discovery of new non-peptide bradykinin B2 receptor antagonists.

Molecular determinants of non-specific recognition of delta, mu, and kappa opioid receptors

Identification of the molecular determinants of recognition common to all three opioid receptors embedded in a single three-dimensional (3D) non-specific recognition pharmacophore has been carried out. The working hypothesis that underlies the computational study reported here is that ligands that bind with significant affinity to all three cloned opioid receptors, delta, mu, and kappa, but with different combinations of activation and inhibition properties at these receptors, could be promising behaviorally selective analgesics with diminished side effects. The study presented here represents the first step towards the rational design of such therapeutic agents. The common 3D pharmacophore developed for recognition of delta, mu, and kappa opioid receptors was based on the receptor affinities determined for 23 different opioid ligands that display no specificity for any of the receptor subtypes. The pharmacophore centers identified are a protonated amine, two hydrophobic groups, and the centroid of an aromatic group in a geometric arrangement common to all 23, non-specific, opioid ligands studied. Using this three-dimensional pharmacophore as a query for searching 3D structural databases, novel compounds potentially involved in non-specific recognition of delta, mu, and kappa opioid receptors were retrieved. These compounds can be valuable candidates for novel behaviorally selective analgesics with diminished or no side effects, and thus with potential therapeutic usefulness.

Conformational analysis of the highly potent bradykinin antagonist Hoe-140 by means of two different computational methods

The AMBER 4.0 force field was used to perform the characterization of the conformational profile of the highly potent bradykinin antagonist Hoe-140 (D-Arg0-Arg1-Pro2-Hyp3-Gly4-Thi5-Ser6-D-++ +Tic7-Oic8-Arg9). The structural features of the peptide were assessed using two different computational methods, both capable to provide a good sampling of the low-energy conformations of the molecule. Specifically, the conformational space of the peptide was explored: i) computing molecular dynamics trajectories in cycles of high (900 K) and low (300 K) temperature and ii) using simulated annealing (SA) in an iterative fashion. Analysis of the structures characterized indicates that most of the low-energy conformations of the peptide exhibit a betaII'-turn motif at its C-terminus, in agreement with previous experimental and theoretical studies. On the other hand, about a 50% of the low-energy conformations characterized also exhibit different beta-turn type motifs at the N-terminus, whereas the rest of the conformations can be described as bends. Finally, in order to get new insights into the structural requirements necessary to design more potent and selective antagonists of bradykinin, present results were compared with those previously reported by this laboratory on the conformational preferences of the native nonapeptide and its DPhe7 analog.

Computational study of the conformational domains of peptide T

The conformational preferences of peptide T (ASTTTNYT) were analysed by means of computational methods. A thorough exploration of the conformational space was carried out within the framework of the molecular mechanics approach, using simulated annealing as a searching strategy. Specifically, in order to obtain a subset of low-energy conformations with energies close to the global minimum as complete as possible, a simulated annealing protocol was repeated several times in a recursive fashion. The results of the search indicate that the peptide exhibits a alpha-helical character although most of the conformations characterized, including the global minimum, can be described as bent conformations. Conformations exhibiting beta-turn motives previously proposed from NMR studies were also characterized, although they are not very predominant in the set of low-energy conformations.

On the bioactive conformation of a small peptide and its set of thermodynamically accessible conformations

In order to investigate the relationship between the bioactive conformation of a peptide and its set of thermodynamically accessible structures in solution, the conformational profile of the tetrapeptide Ac-Pro-Ala-Pro-Tyr-OH was characterized by computational methods. Search of the conformational space was performed within the molecular mechanics frame-work using the AMBER4.0 force field with an effective dielectric constant of 80. Unique structures of the peptide were compared with its bioactive conformation for the protein Streptomyces griseus Protease A, as taken from the crystal structure of the enzyme-peptide complex. The results show that the bound conformation is close to one of the unique conformations characterized in the conformational search of the isolated peptide. Moreover, the lowest energy minimum characterized in the conformational search exhibits large deviations when compared to the bound conformation of the crystal structure.

Conformational study of [Met5]enkephalin-Arg-Phe in the presence of phosphatidylserine vesicles

The interaction of [Met5]enkephalin-Arg.Phe with phosphatidylserine (PtdSer) was studied by circular dichroism (CD), two-dimensional nuclear magnetic resonance spectroscopy, hybrid distance geometry simulated annealing (DG-SA) and molecular dynamics (MD) calculations. The very low solubility of [Met5]enkephalin-Arg-Phe and the instability of the solution containing PtdSer vesicles at low pH values did not allow us to observe the amide proton resonances in the usual two-dimensional NMR work. NOESY cross-peaks of protons of side chains from two-dimensional NMR were converted into distances which were used as restraints for modelling with DG-SA and MD. Our results indicate that, in aqueous solutions at pH 7.68 [Met5]enkephalin-Arg-Phe exists in the absence of PtdSer as a random distribution of conformers, whereas in the presence of PtdSer it adopts conformations containing a common orientation of the bonds of C alpha 2, C alpha 3, C alpha 4, and C alpha 5, although different orientations of the peptide planes are consistent with the results. Two of the reported conformers from MD simulations are characterized by the presence of a 2<--4 gamma and inverse gamma turns centered on Gly3. A gradual decline of order was observed when moving from the central moiety of the peptide to both the N-terminus and C-terminus. Finally, the DG-SA and MD calculations resulted in a structure such that the orientation of the Phe4 and Met5 side chains favours hydrophobic interactions with the apolar portion of the PtdSer vesicle to form a hydrophobic cluster. These data support the hypothesis of a role of lipids to modify the conformation of [Met5]enkephalin-Arg-Phe to permit the interactions with the receptor site.

Characterization of the bioactive form of linear peptide antagonists at the omega-opioid receptor

The stereochemical requirements for omega-opioid receptor binding of a series of linear peptide antagonists with a novel conformationally restricted Phe analogue (Tic) as a second residue were examined by using a variety of computational chemistry methods. The omega-opioid receptor analogues with significant affinity, Tyr-Tic-NH2(TI-NH2), Tyr-Tic-Phe-OH(TIP), Tyr-Tic-Phe-NH2(TIP-NH2), Tyr-Tic-Phe-Phe-OH(TIPP), Tyr-Tic-Phe-Phe-NH2)(TIPP-NH2), and the low affinity omega-opioid peptides Tyr-Pro-Phe-Pro-NH2(morphiceptin) and Tyr-Phe-Phe-Phe-NH2 (TPPP-NH2), were included in this study. The conformational profiles of these peptides were obtained by consecutive cycles of high and low temperature molecular dynamic stimulations, coupled to molecular mechanical energy minimization carried out until no new conformational minima were obtained. Comparing the results for TPPP-NH2 and TIPP-NH2, the presence of the conformationally restricted Tic residue did not greatly reduce the number of unique low energy conformations, but did allow low energy conformers involving cis bonds between the first two residues. The conformational libraries of these peptides were examined for their ability to satisfy the three key ligand components for receptor recognition already identified by previous studies of high affinity cyclic (Tyr1-D-Pen2-Gly3-Phe4-D-Pen5) enkephalin (DPDPE) type agonists: a protonated amine group, an aromatic ring, and a lipophilic moiety in a specific geometric arrangement. Two types of conformations common to the five high omega-opioid affinity L-Tic analogues were found that satisfied these requirements, one with a cis and the other with a trans peptide bond between the Tyr1 and Tic2 residues. Moreover, both the Tic2 and Phe3 residues could mimic the hydrophobic interactions with the receptor of the Phe4 moiety in the cyclic DPDPE type agonists, consistent with the appreciable affinity of both di- and tripeptides. The low omega-opioid receptor affinity of morphiceptin can be understood as the result of conformational preferences that prevent the fulfillment of this pharmacophore for recognition.

Molecular modelling and conformational analysis of a GABAB antagonist

Crystallographic database studies and molecular dynamics simulations in different media have enabled us to sample the conformational space of a GABAB antagonist. As a result, we have defined a pharmacophoric pattern for GABAB antagonists. This study has led us to compare the conformational preferences deduced from database studies and molecular dynamics simulations. The influence of the medium on the conformations has also been investigated.

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 20 degrees, 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 beta-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.

Effects of different continium dielectric models in a molecular dynamics and energy minimization study of the antigenic loop of foot-and-mouth disease virus

This study presents the influence of the dielectric constant on the final structure of the major antigenic loop of the FMDV serotype C. Minimizations have been performed on the nine-residue peptide Ac-TASARGDLA-NHMe, using two expressions for the dielectric constant: a distance-independent (epsilon = kappa), and a distance-dependent (epsilon = kappa *r) forms, and for kappa values from 1 to 10. In addition, kappa = 40 and 80 has also been considered for the constant expression of the dielectrics. The calculations suggest, for both expressions of the dielectrics, that for the sequence under study a value of kappa in the range of 2-6 performs an adequate treatment of the electrostatic interactions. Finally, molecular dynamic simulations at 298 K calculated with an effective dielectric constants of epsilon = 1*r and epsilon = 4*r are compared. The results indicated that a dielectric of epsilon = 4*r provides a peptide fluctuations which are in agreement with recent X-ray diffraction studies.

Strategies for indirect computer-aided drug design

This review is intended to describe some of the methods and procedures used for computer-aided drug design when the structure of the macromolecular target is unknown, as is the case for CNS active drugs. Strategies and methods used in computer-aided design of drugs in such instances must be "indirect," i.e., focusing on the characterization of the ligands themselves. This situation is different from one in which the three-dimensional structure of the macromolecular target for a drug is known, for example, for drugs that are enzyme inhibitors, allowing "direct" characterization of ligand-receptor interactions. Two qualitatively different "indirect" approaches are described here. One, called 2D-QSAR, is briefly reviewed. It is based on delineating regression relationships between a specified biological end point and properties of the compounds eliciting it. The other, based on pharmacophore development, constitutes the main part of this review. Several levels of pharmacophore development are described, which differ in the extent to which they encompass fundamental molecular properties that are determinants of receptor recognition and activation. The strengths and limitations of each procedure are discussed and illustrated by examples. Two methods for obtaining model receptor structures are then briefly described. Both rely on the prior success of the indirect methods in obtaining ligand properties that modulate receptor recognition and activation. These emerging capabilities have the potential to bridge the gap between indirect and direct methods of drug design, since, if successful, the design process can continue in a direct mode using explicit characterization of drug-receptor interactions. Strategies for hypothesis validation and use of hypothesis for drug design and discovery are also briefly reviewed. The final sections of this review describe specific computational tools such as molecular mechanics and quantum mechanical methods used to characterize and identify relevant molecular properties and indicate some areas for future development of computational chemistry methods that could increase its effectiveness in the design of novel drugs.

SCF-MO study of the polyglycine II structure

In order to get insight into the conditions that make polyglycine (PG)II a stable structure, the conformational features of three model molecules closely related to the PGII conformation were investigated. The model molecules selected were glycine dipeptide (AGN), glycine tripeptide (AGGN), and glycine tetrapeptide (AGGGN). Environmental effects were mimicked by means of formaldehyde molecules. The calculations were carried out at the SCF semiempirical level, using the AM1 method. The calculations show that of the three systems considered, only the AGGGN molecule presents a minimum energy conformation which corresponds to a PGII structure. The environmental conditions in which this conformation is found were also analyzed.

Conformational study of a nine residue fragment of the antigenic loop of foot-and-mouth disease virus

The nine-residue peptide Ac-TASARGDLA-NHMe was selected as model peptide in order to understand the conformational features of the antigenic loop of foot-and-mouth disease virus (FMDV). A throughout exploration of the conformational space has been carried out by means of molecular dynamics (MD) and energy minimization. The calculations have been carried out using the AMBER force field. Solvent effects have been included by an effective dielectric constant of epsilon = 4r. The lowest energy conformation presents a secondary structure constituted by an alpha-helix at the N-terminal end followed by two gamma-turns in the central region. The rest of the accessible minima found present also a high tendency to form gamma-turns. Finally, a 100 ps MD trajectory calculation at 298 K suggest a stability of the secondary structure elements of the lowest energy conformation.