1
|
Investigation of the conformational space of hydrophobic-polar heteropolymers by gyration tensor based parameters. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2021.111372] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
2
|
Ravichandran A, Gu G, Escano J, Lu SE, Smith L. The presence of two cyclase thioesterases expands the conformational freedom of the cyclic Peptide occidiofungin. JOURNAL OF NATURAL PRODUCTS 2013; 76:150-156. [PMID: 23394257 PMCID: PMC4142711 DOI: 10.1021/np3005503] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Occidiofungin is a cyclic nonribosomally synthesized antifungal peptide with submicromolar activity produced by the Gram-negative bacterium Burkholderia contaminans. The biosynthetic gene cluster was confirmed to contain two cyclase thioesterases. NMR analysis revealed that the presence of both thioesterases is used to increase the conformational repertoire of the cyclic peptide. The loss of the OcfN cyclic thioesterase by mutagenesis results in a reduction of conformational variants and an appreciable decrease in bioactivity against Candida species. Presumably, the presence of both asparagine and β-hydroxyasparagine variants coordinates the enzymatic function of both of the cyclase thioesterases. OcfN has presumably evolved to be part of the biosynthetic gene cluster due to its ability to produce structural variants that enhance antifungal activity against some fungi. The enhancement of the antifungal activity from the incorporation of an additional cyclase thioesterase into the biosynthetic gene cluster of occidiofungin supports the need to explore new conformational variants of other therapeutic or potentially therapeutic cyclic peptides.
Collapse
Affiliation(s)
- Akshaya Ravichandran
- Department of Biological Sciences, Texas A&M University, College Station, TX 77843
| | - Ganyu Gu
- Department of Biochemistry, Molecular Biology, Entomology and Plant Pathology, Mississippi State University, 32 Creelman St., Mississippi State, MS 39762
| | - Jerome Escano
- Department of Biological Sciences, Texas A&M University, College Station, TX 77843
| | - Shi-En Lu
- Department of Biochemistry, Molecular Biology, Entomology and Plant Pathology, Mississippi State University, 32 Creelman St., Mississippi State, MS 39762
| | - Leif Smith
- Department of Biological Sciences, Texas A&M University, College Station, TX 77843
| |
Collapse
|
3
|
Abstract
Loop modeling is crucial for high-quality homology model construction outside conserved secondary structure elements. Dozens of loop modeling protocols involving a range of database and ab initio search algorithms and a variety of scoring functions have been proposed. Knowledge-based loop modeling methods are very fast and some can successfully and reliably predict loops up to about eight residues long. Several recent ab initio loop simulation methods can be used to construct accurate models of loops up to 12-13 residues long, albeit at a substantial computational cost. Major current challenges are the simulations of loops longer than 12-13 residues, the modeling of multiple interacting flexible loops, and the sensitivity of the loop predictions to the accuracy of the loop environment.
Collapse
|
4
|
Meirovitch H. Methods for calculating the absolute entropy and free energy of biological systems based on ideas from polymer physics. J Mol Recognit 2010; 23:153-72. [PMID: 19650071 PMCID: PMC2823937 DOI: 10.1002/jmr.973] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
The commonly used simulation techniques, Metropolis Monte Carlo (MC) and molecular dynamics (MD) are of a dynamical type which enables one to sample system configurations i correctly with the Boltzmann probability, P(i)(B), while the value of P(i)(B) is not provided directly; therefore, it is difficult to obtain the absolute entropy, S approximately -ln P(i)(B), and the Helmholtz free energy, F. With a different simulation approach developed in polymer physics, a chain is grown step-by-step with transition probabilities (TPs), and thus their product is the value of the construction probability; therefore, the entropy is known. Because all exact simulation methods are equivalent, i.e. they lead to the same averages and fluctuations of physical properties, one can treat an MC or MD sample as if its members have rather been generated step-by-step. Thus, each configuration i of the sample can be reconstructed (from nothing) by calculating the TPs with which it could have been constructed. This idea applies also to bulk systems such as fluids or magnets. This approach has led earlier to the "local states" (LS) and the "hypothetical scanning" (HS) methods, which are approximate in nature. A recent development is the hypothetical scanning Monte Carlo (HSMC) (or molecular dynamics, HSMD) method which is based on stochastic TPs where all interactions are taken into account. In this respect, HSMC(D) can be viewed as exact and the only approximation involved is due to insufficient MC(MD) sampling for calculating the TPs. The validity of HSMC has been established by applying it first to liquid argon, TIP3P water, self-avoiding walks (SAW), and polyglycine models, where the results for F were found to agree with those obtained by other methods. Subsequently, HSMD was applied to mobile loops of the enzymes porcine pancreatic alpha-amylase and acetylcholinesterase in explicit water, where the difference in F between the bound and free states of the loop was calculated. Currently, HSMD is being extended for calculating the absolute and relative free energies of ligand-enzyme binding. We describe the whole approach and discuss future directions.
Collapse
Affiliation(s)
- Hagai Meirovitch
- Department of Computational Biology, University of Pittsburgh School of Medicine, 3059 BST3, Pittsburgh, PA 15260, USA.
| |
Collapse
|
5
|
Cheluvaraja S, Mihailescu M, Meirovitch H. Entropy and free energy of a mobile protein loop in explicit water. J Phys Chem B 2008; 112:9512-22. [PMID: 18613721 PMCID: PMC2671085 DOI: 10.1021/jp801827f] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Estimation of the energy from a given Boltzmann sample is straightforward since one just has to average the contribution of the individual configurations. On the other hand, calculation of the absolute entropy, S (hence the absolute free energy F) is difficult because it depends on the entire (unknown) ensemble. We have developed a new method called "the hypothetical scanning molecular dynamics" (HSMD) for calculating the absolute S from a given sample (generated by any simulation technique). In other words, S (like the energy) is "written" on the sample configurations, where HSMD provides a prescription of how to "read" it. In practice, each sample conformation, i, is reconstructed with transition probabilities, and their product leads to the probability of i, hence to the entropy. HSMD is an exact method where all interactions are considered, and the only approximation is due to insufficient sampling. In previous studies HSMD (and HS Monte CarloHSMC) has been extended systematically to systems of increasing complexity, where the most recent is the seven-residue mobile loop, 304-310 (Gly-His-Gly-Ala-Gly-Gly-Ser) of the enzyme porcine pancreatic alpha-amylase modeled by the AMBER force field and AMBER with the implicit solvation GB/SA (paper I, Cheluvaraja, S.; Meirovitch, H. J. Chem. Theory Comput. 2008, 4, 192). In the present paper we make a step further and extend HSMD to the same loop capped with TIP3P explicit water at 300 K. As in paper I, we are mainly interested in entropy and free energy differences between the free and bound microstates of the loop, which are obtained from two separate MD samples of these microstates. The contribution of the loop to S and F is calculated by HSMD and that of water by a particular thermodynamic integration procedure. As expected, the free microstate is more stable than the bound microstate by a total free energy difference, Ffree-Fbound=-4.8+/-1, as compared to -25.5 kcal/mol obtained with GB/SA. We find that relatively large systematic errors in the loop entropies, Sfree(loop) and Sbound(loop) are cancelled in their difference which is thus obtained efficiently and with high accuracy, i.e., with a statistical error of 0.1 kcal/mol. This cancellation, which has been observed in previous HSMD studies, is in accord with theoretical arguments given in paper I.
Collapse
Affiliation(s)
- Srinath Cheluvaraja
- Department of Computational Biology, University of Pittsburgh School of Medicine, 3059 BST3, Pittsburgh, PA 15260
| | - Mihail Mihailescu
- Department of Computational Biology, University of Pittsburgh School of Medicine, 3059 BST3, Pittsburgh, PA 15260
| | - Hagai Meirovitch
- Department of Computational Biology, University of Pittsburgh School of Medicine, 3059 BST3, Pittsburgh, PA 15260
| |
Collapse
|
6
|
Cheluvaraja S, Meirovitch H. Stability of the Free and Bound Microstates of a Mobile Loop of α-Amylase Obtained from the Absolute Entropy and Free Energy. J Chem Theory Comput 2007; 4:192-208. [DOI: 10.1021/ct700116n] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Srinath Cheluvaraja
- Department of Computational Biology, University of Pittsburgh School of Medicine, 3059 BST3, Pittsburgh, Pennsylvania 15260
| | - Hagai Meirovitch
- Department of Computational Biology, University of Pittsburgh School of Medicine, 3059 BST3, Pittsburgh, Pennsylvania 15260
| |
Collapse
|
7
|
Rezai T, Bock JE, Zhou MV, Kalyanaraman C, Lokey RS, Jacobson MP. Conformational Flexibility, Internal Hydrogen Bonding, and Passive Membrane Permeability: Successful in Silico Prediction of the Relative Permeabilities of Cyclic Peptides. J Am Chem Soc 2006; 128:14073-80. [PMID: 17061890 DOI: 10.1021/ja063076p] [Citation(s) in RCA: 290] [Impact Index Per Article: 16.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report an atomistic physical model for the passive membrane permeability of cyclic peptides. The computational modeling was performed in advance of the experiments and did not involve the use of "training data". The model explicitly treats the conformational flexibility of the peptides by extensive conformational sampling in low (membrane) and high (water) dielectric environments. The passive membrane permeabilities of 11 cyclic peptides were obtained experimentally using a parallel artificial membrane permeability assay (PAMPA) and showed a linear correlation with the computational results with R(2) = 0.96. In general, the results support the hypothesis, already well established in the literature, that the ability to form internal hydrogen bonds is critical for passive membrane permeability and can be the distinguishing factor among closely related compounds, such as those studied here. However, we have found that the number of internal hydrogen bonds that can form in the membrane and the solvent-exposed polar surface area correlate more poorly with PAMPA permeability than our model, which quantitatively estimates the solvation free energy losses upon moving from high-dielectric water to the low-dielectric interior of a membrane.
Collapse
Affiliation(s)
- Taha Rezai
- Department of Chemistry and Biochemistry, University of California at Santa Cruz, Santa Cruz, CA 95064, USA
| | | | | | | | | | | |
Collapse
|
8
|
Cheluvaraja S, Meirovitch H. Calculation of the entropy and free energy of peptides by molecular dynamics simulations using the hypothetical scanning molecular dynamics method. J Chem Phys 2006; 125:24905. [PMID: 16848609 DOI: 10.1063/1.2208608] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Hypothetical scanning (HS) is a method for calculating the absolute entropy S and free energy F from a sample generated by any simulation technique. With this approach each sample configuration is reconstructed with the help of transition probabilities (TPs) and their product leads to the configuration's probability, hence to the entropy. Recently a new way for calculating the TPs by Monte Carlo (MC) simulations has been suggested, where all system interactions are taken into account. Therefore, this method--called HSMC--is in principle exact where the only approximation is due to insufficient sampling. HSMC has been applied very successfully to liquid argon, TIP3P water, self-avoiding walks on a lattice, and peptides. Because molecular dynamics (MD) is considered to be significantly more efficient than MC for a compact polymer chain, in this paper HSMC is extended to MD simulations as applied to peptides. Like before, we study decaglycine in vacuum but for the first time also a peptide with side chains, (Val)(2)(Gly)(6)(Val)(2). The transition from MC to MD requires implementing essential changes in the reconstruction process of HSMD. Results are calculated for three microstates, helix, extended, and hairpin. HSMD leads to very stable differences in entropy TDeltaS between these microstates with small errors of 0.1-0.2 kcal/mol (T=100 K) for a wide range of calculation parameters with extremely high efficiency. Various aspects of HSMD and plans for future work are discussed.
Collapse
Affiliation(s)
- Srinath Cheluvaraja
- Department of Computational Biology, University of Pittsburgh School of Medicine, 3059 BST3, Pittsburgh, PA 15260, USA
| | | |
Collapse
|
9
|
Cheluvaraja S, Meirovitch H. Calculation of the entropy and free energy by the hypothetical scanning Monte Carlo method: application to peptides. J Chem Phys 2006; 122:54903. [PMID: 15740349 DOI: 10.1063/1.1835911] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new approach, the hypothetical scanning Monte Carlo (HSMC), for calculating the absolute entropy, S, and free energy, F, has been introduced recently and applied first to fluids (argon and water) and later to peptides. In this paper the method is further developed for peptide chains in vacuum. S is calculated from a given MC sample by reconstructing each sample conformation i step-by-step, i.e., calculating transition probabilities (TPs) for the dihedral and bond angles and fixing the related atoms at their positions. At step k of the process the chain's coordinates that have already been determined are kept fixed (the "frozen past") and TP(k) is obtained from a MC simulation of the "future" part of the chain whose TPs as yet have not been determined; when the process is completed the contribution of conformation i to the entropy is, S(i) approximately -ln Pi(k) TP(k). In a recent paper we studied polyglycine chains, modeled by the AMBER force field with constant bond lengths and bond angles (the rigid model). Decaglycine [(Gly)(10)] was studied in the helical, extended, and hairpin microstates, while (Gly)(16) was treated only in the first two microstates. In this paper the samples are increased and restudied, (Gly)(16) is also investigated in the hairpin microstate, and for (Gly)(10) approximations are tested where only part of the future is considered for calculating the TPs. We calculate upper and lower bounds for F and demonstrate that like for fluids, F can be obtained from multiple reconstructions of a single conformation. We also test a more realistic model of (Gly)(10) where the bond angles are allowed to move (the flexible model). Very accurate results for S and F are obtained which are compared to results obtained by the quasiharmonic approximation and the local states method. Thus, differences in entropy and free energy between the three microstates are obtained within errors of 0.1-0.3 kcal/mol. The HSMC method can be applied to a macromolecule with any degree of flexibility, ranging from local fluctuations to a random coil. The present results demonstrate that the difference in stability, DeltaF(mn)=F(m)-F(n) between significantly different microstates m and n, can be obtained from two simulations only without the need to resort to thermodynamic integration. Our long-term goal is to extend this method to any peptide and apply it to a peptide immersed in a box with explicit water.
Collapse
Affiliation(s)
- Srinath Cheluvaraja
- Center for Computational Biology and Bioinformatics and Department of Molecular Genetics and Biochemistry, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | | |
Collapse
|
10
|
Szarecka A, Meirovitch H. Optimization of the GB/SA solvation model for predicting the structure of surface loops in proteins. J Phys Chem B 2006; 110:2869-80. [PMID: 16471897 PMCID: PMC1945207 DOI: 10.1021/jp055771+] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Implicit solvation models are commonly optimized with respect to experimental data or Poisson-Boltzmann (PB) results obtained for small molecules, where the force field is sometimes not considered. In previous studies, we have developed an optimization procedure for cyclic peptides and surface loops in proteins based on the entire system studied and the specific force field used. Thus, the loop has been modeled by the simplified solvation function E(tot) = E(FF) (epsilon = 2r) + Sigma(i) sigma(i)A(i), where E(FF) (epsilon = nr) is the AMBER force field energy with a distance-dependent dielectric function, epsilon = nr, A(i) is the solvent accessible surface area of atom i, and sigma(i) is its atomic solvation parameter. During the optimization process, the loop is free to move while the protein template is held fixed in its X-ray structure. To improve on the results of this model, in the present work we apply our optimization procedure to the physically more rigorous solvation model, the generalized Born with surface area (GB/SA) (together with the all-atom AMBER force field) as suggested by Still and co-workers (J. Phys. Chem. A 1997, 101, 3005). The six parameters of the GB/SA model, namely, P(1)-P(5) and the surface area parameter, sigma (programmed in the TINKER package) are reoptimized for a "training" group of nine loops, and a best-fit set is defined from the individual sets of optimized parameters. The best-fit set and Still's original set of parameters (where Lys, Arg, His, Glu, and Asp are charged or neutralized) were applied to the training group as well as to a "test" group of seven loops, and the energy gaps and the corresponding RMSD values were calculated. These GB/SA results based on the three sets of parameters have been found to be comparable; surprisingly, however, they are somewhat inferior (e.g, of larger energy gaps) to those obtained previously from the simplified model described above. We discuss recent results for loops obtained by other solvation models and potential directions for future studies.
Collapse
Affiliation(s)
- Agnieszka Szarecka
- Department of Computational Biology, University of Pittsburgh School of Medicine, Suite 3064, BST 3, 3501 Fifth Avenue, Pittsburgh, PA 15213
| | - Hagai Meirovitch
- Department of Computational Biology, University of Pittsburgh School of Medicine, Suite 3064, BST 3, 3501 Fifth Avenue, Pittsburgh, PA 15213
| |
Collapse
|
11
|
Cheluvaraja S, Meirovitch H. Calculation of the entropy and free energy from monte carlo simulations of a peptide stretched by an external force. J Phys Chem B 2005; 109:21963-70. [PMID: 16853854 PMCID: PMC1540612 DOI: 10.1021/jp052969l] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Hypothetical scanning Monte Carlo (HSMC) is a method for calculating the absolute entropy, S, and free energy, F, from a trajectory generated by any simulation technique. HSMC was applied initially to fluids (argon and water) and later to peptides and self-avoiding walks on a lattice. In this paper we make a step further and apply it to a model of decaglycine (at T = 300 K) in vacuum with constant bond lengths where external stretching forces are exerted at the end points; the changes in S and F are calculated as the forces are increased. The molecule is placed initially in a helical structure, which is changed to an extended structure after a short simulation time due to the exerted forces. This study has relevance to problems in polymers (e.g., rubber elasticity) and to the analysis of experiments where individual molecules are stretched by atomic force microscopy (AFM), for example. The results for S and F are accurate and are significantly better than those obtained by the quasi-harmonic approximation and the local states method. However, the molecule is quite stiff due to the strong bond angle potentials and the extensions are small even for relatively large forces. Correspondingly, as the force is increased the decrease in the entropy is relatively small while the potential energy is enhanced significantly. Still, differences, TDeltaS, for different forces are obtained with very good accuracy of approximately 0.2 kcal/mol.
Collapse
Affiliation(s)
- Srinath Cheluvaraja
- Department of Computational Biology University of Pittsburgh School of Medicine W1058 BST, Pittsburgh, PA 15261
| | - Hagai Meirovitch
- Department of Computational Biology University of Pittsburgh School of Medicine W1058 BST, Pittsburgh, PA 15261
| |
Collapse
|
12
|
Duran D, Aviyente V, Baysa C. Solvent effect on the synthesis of clarithromycin: a molecular dynamics study. J Comput Aided Mol Des 2005; 18:145-54. [PMID: 15287700 DOI: 10.1023/b:jcam.0000030037.67742.cb] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Clarithromycin (6-O-methylerythromycin A) is a 14-membered macrolide antibiotic which is active in vitro against clinically important gram-positive and gram-negative bacteria. The selectivity of the methylation of the C-6 OH group is studied on erythromycin A derivatives. To understand the effect of the solvent on the methylation process, detailed molecular dynamics (MD) simulations are performed in pure DMSO, pure THF and DMSO:THF (1:1) mixture by using the anions at the C-6, C-11 and C-12 positions of 2',4"-[O-bis(TMS)]erythromycin A 9-[O-(dimethylthexylsilyl)oxime] under the assumption that the anions are stable on the sub-nanosecond time scale. The conformations of the anions are not affected by the presence of the solvent mixture. The radial distribution functions are computed for the distribution of different solvent molecules around the 'O-' of the anions. At distances shorter than 5 A, DMSO molecules are found to cluster around the C-11 anion, whereas the anion at the C-12 position is surrounded by the THF molecules. The anion at the C-6 position is not blocked by the solvent molecules. The results are consistent with the experimental finding that the methylation yield at the latter position is increased in the presence of a DMSO:THF (1:1) solvent mixture. Thus, the effect of the solvent in enhancing the yield during the synthesis is not by changing the conformational properties of the anions, but rather by creating a suitable environment for methylation at the C-6 position.
Collapse
Affiliation(s)
- Dilek Duran
- Chemistry, Department, Faculty of Arts and Sciences, Bogazici University, Bebek 34342 Istanbul, Turkey
| | | | | |
Collapse
|
13
|
Rayan A, Senderowitz H, Goldblum A. Exploring the conformational space of cyclic peptides by a stochastic search method. J Mol Graph Model 2004; 22:319-33. [PMID: 15099829 DOI: 10.1016/j.jmgm.2003.12.012] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
A stochastic search algorithm is applied in order to probe the conformations of cyclic peptides. The search is conducted in two stages. In the first stage, random conformations are generated and evaluated by a penalty function for ring closure ability, following a stepwise construction of each amino acid into the peptide by a random choice of one of its allowed conformations. The allowed conformational ranges of backbone dihedral angles for each amino acid have been extracted from a Data Bank of diverse proteins. Values of dihedral angles that do not contribute favorably to the scoring of ring closure are retained or discarded by a statistical test. Values are discarded up to a point from which all remaining combinations of angles are constructed, scored, sorted, and clustered. In the second stage, side chains have been added and fast optimization was applied to the set of diverse conformations in a "united atoms" approach, with the "Kollman forcefield" of Sybyl 6.8. This iterative stochastic elimination algorithm finds the global minimum and most of the best results, when compared to a full exhaustive search in appropriately sized problems. In larger problems, we compare the results to experimental structures. The root mean square deviation (RMSD) of our best results compared to crystal structures of cyclic peptides with sizes from 4 to 15 amino acids are mostly below 1.0 A up to 8 mers and under 2.0 A for larger cyclic peptides.
Collapse
Affiliation(s)
- Anwar Rayan
- Department of Medicinal Chemistry and Natural Products, David R. Bloom Center for Pharmacy, School of Pharmacy, The Hebrew University of Jerusalem, Jerusalem 91120, Israel.
| | | | | |
Collapse
|
14
|
Cheluvaraja S, Meirovitch H. Simulation method for calculating the entropy and free energy of peptides and proteins. Proc Natl Acad Sci U S A 2004; 101:9241-6. [PMID: 15197271 PMCID: PMC438960 DOI: 10.1073/pnas.0308201101] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
A method called complete hypothetical scanning Monte Carlo has been introduced for calculating the absolute entropy, S, and free energy, F, of fluids. Here, the method is extended to peptide chains in vacuum. Thus, S is calculated from a given sample by reconstructing each conformation step-by-step by using transition probabilities (TPs); at each step, part of the chain coordinates have already been determined (the "frozen past"), and the TP is obtained from a Monte Carlo simulation of the (future) part of the chain whose TPs as yet have not been calculated. Very accurate results for S and F are obtained for the helix, extended, and hairpin microstates of a simplified model of decaglycine (Gly)(10) and (Gly)(16). These results agree well with results obtained by the quasiharmonic approximation and the local states method. The complete HSMC method can be applied to a macromolecule with any degree of flexibility, ranging from local fluctuations to a random coil. Also, the difference in stability, Delta F(mn) = F(m) - F(n) between significantly different microstates m and n can be obtained from two simulations only without the need to resort to thermodynamic integration. Our long-term goal is to extend this method to any peptide and apply it to a peptide immersed in a box with explicit water.
Collapse
Affiliation(s)
- Srinath Cheluvaraja
- Center for Computational Biology and Bioinformatics and Department of Molecular Genetics and Biochemistry, University of Pittsburgh School of Medicine, W1058 BST, PA 15261, USA
| | | |
Collapse
|
15
|
Szarecka A, White RP, Meirovitch H. Absolute entropy and free energy of fluids using the hypothetical scanning method. I. Calculation of transition probabilities from local grand canonical partition functions. J Chem Phys 2003. [DOI: 10.1063/1.1625919] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Yaşar F, Arkin H, Celik T, Berg BA, Meirovitch H. Efficiency of the multicanonical simulation method as applied to peptides of increasing size: the heptapeptide deltorphin. J Comput Chem 2002; 23:1127-34. [PMID: 12116381 DOI: 10.1002/jcc.10113] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The advantage of the multicanonical (MUCA) simulation method of Berg and coworkers over the conventional Metropolis method is in its ability to move a system effectively across energy barriers thereby providing results for a wide range of temperatures. However, a MUCA simulation is based on weights (related to the density of states) that should be determined prior to a production run and their calculation is not straightforward. To overcome this difficulty a procedure has been developed by Berg that calculates the MUCA weights automatically. In a previous article (Yaşar et al. J Comput Chem 2000, 14, 1251-1261) we extended this procedure to continuous systems and applied it successfully to the small pentapeptide Leu-enkephalin. To investigate the performance of the automated MUCA procedure for larger peptides, we apply it here to deltorphin, a linear heptapeptide with bulky side chains (H-Tyr(1)-D-Met(2)-Phe(3)-His(4)-Leu(5)-Met(6)-Asp(7)-NH(2)). As for Leu-enkephalin, deltorphin is modeled in vacuum by the potential energy function ECEPP. MUCA is found to perform well. A weak second peak is seen for the specific heat, which is given a special attention. By minimizing the energy of structures along the trajectory it is found that MUCA provides a good conformational coverage of the low energy region of the molecule. These latter results are compared with conformational coverage obtained by the Monte Carlo minimization method of Li and Scheraga.
Collapse
Affiliation(s)
- Fatih Yaşar
- Department of Physics Engineering, Hacettepe University, 06532, Ankara, Turkey
| | | | | | | | | |
Collapse
|
17
|
Groth M, Malicka J, Rodziewicz- Motowidło S, Czaplewski C, Klaudel L, Wiczk W, Liwo A. Determination of conformational equilibrium of peptides in solution by NMR spectroscopy and theoretical conformational analysis: application to the calibration of mean-field solvation models. Biopolymers 2001; 60:79-95. [PMID: 11455544 DOI: 10.1002/1097-0282(2001)60:2<79::aid-bip1006>3.0.co;2-l] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Peptides occur in solution as ensembles of conformations rather than in a fixed conformation. The existing energy functions are usually inadequate to predict the conformational equilibrium in solution, because of failure to account properly for solvation, if the solvent is not considered explicitly (which is usually prohibitively expensive). NMR data are therefore widely incorporated into theoretical conformational analysis. Because of conformational flexibility, restrained molecular dynamics (with restraints derived from NMR data), which is usually applied to determine protein conformation is of limited use in the case of peptides. Instead, (a) the restraints are averaged within predefined time windows during molecular dynamics (MD) simulations (time averaging), (b) multiple-copy MD simulations are carried out and the restraints are averaged over the copies (ensemble averaging), or (c) a representative ensemble of sterically feasible conformations is generated and the weights of the conformations are then fitted so that the computed average observables match the experimental data (weight fitting). All these approaches are briefly discussed in this article. If an adequate force field is used, conformations with large statistical weights obtained from the weight-fitting procedure should also have low energies, which can be implemented in force field calibration. Such a procedure is particularly attractive regarding the parameterization of the solvation energy in nonaqueous solvents, e.g., dimethyl sulfoxide, for which thermodynamic solvation data are scarce. A method for calibration of solvation parameters in dimethyl sulfoxide, which is based on this principle was recently proposed by C. Baysal and H. Meirovitch (Journal of the American Chemical Society, 1998, Vol. 120, pp. 800--812), in which the energy gap between the conformations compatible with NMR data and the alternative conformations is maximized. In this work we propose an alternative method based on the principle that the best-fitting statistical weights of conformations should match the Boltzmann weights computed with the force field applied. Preliminary results obtained using three test peptides of varying conformational mobility: H-Ser(1)-Pro(2)-Lys(3)-Leu(4)-OH, Ac-Tyr(1)-D-Phe(2)-Ser(3)-Pro(4)-Lys(5)-Leu(6)-NH(2), and cyclo(Tyr(1)-D-Phe(2)-Ser(3)-Pro(4)-Lys(5)-Leu(6)) are presented.
Collapse
Affiliation(s)
- M Groth
- Faculty of Chemistry University of Gdańsk Sobieskiego 18 80-952 Gdańsk Poland
| | | | | | | | | | | | | |
Collapse
|
18
|
Das B, Meirovitch H. Optimization of solvation models for predicting the structure of surface loops in proteins. Proteins 2001; 43:303-14. [PMID: 11288180 DOI: 10.1002/prot.1041] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
A novel procedure for optimizing the atomic solvation parameters (ASPs) sigma(i) developed recently for cyclic peptides is extended to surface loops in proteins. The loop is free to move, whereas the protein template is held fixed in its X-ray structure. The energy is E(tot) = E(FF)(epsilon = nr) + summation operator sigma(i)A(i), where E(FF)(epsilon = nr) is the force-field energy of the loop-loop and loop-template interactions, epsilon = nr is a distance-dependent dielectric constant, and n is an additional parameter to be optimized. A(i) is the solvent-accessible surface area of atom i. The optimal sigma(i) and n are those for which the loop structure with the global minimum of E(tot)(n, sigma(i)) becomes the experimental X-ray structure. Thus, the ASPs depend on the force field and are optimized in the protein environment, unlike commonly used ASPs such as those of Wesson and Eisenberg (Protein Sci 1992;1:227-235). The latter are based on the free energy of transfer of small molecules from the gas phase to water and have been traditionally combined with various force fields without further calibration. We found that for loops the all-atom AMBER force field performed better than OPLS and CHARMM22. Two sets of ASPs [based on AMBER (n = 2)], optimized independently for loops 64-71 and 89-97 of ribonuclease A, were similar and thus enabled the definition of a best-fit set. All these ASPs were negative (hydrophilic), including those for carbon. Very good (i.e., small) root-mean-square-deviation values from the X-ray loop structure were obtained with the three sets of ASPs, suggesting that the best-fit set would be transferable to loops in other proteins as well. The structure of loop 13-24 is relatively stretched and was insensitive to the effect of the ASPs.
Collapse
Affiliation(s)
- B Das
- School of Computational Science and Information Technology, Florida State University, Tallahassee, FL 32306-4052, USA
| | | |
Collapse
|
19
|
Baysal C, Meirovitch H. On the transferability of atomic solvation parameters: Ab initio structural prediction of cyclic heptapeptides in DMSO. Biopolymers 2000; 54:416-28. [PMID: 10951328 DOI: 10.1002/1097-0282(200011)54:6<416::aid-bip60>3.0.co;2-2] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
A statistical mechanics methodology for predicting the solution structures and populations of peptides developed recently is based on a novel method for optimizing implicit solvation models, which was applied initially to a cyclic hexapeptide in DMSO (C. Baysal and H. Meirovitch, Journal of American Chemical Society, 1998, vol. 120, pp. 800-812). Thus, the molecule has been described by the simplified energy function E(tot) = E(GRO) + summation operator(k) sigma(k)A(k), where E(GRO) is the GROMOS force-field energy, sigma(k) and A(k) are the atomic solvation parameter (ASP) and the solvent accessible surface area of atom k, respectively. In a more recent study, these ASPs have been found to be transferable to the cyclic pentapeptide cyclo(D-Pro(1)-Ala(2)-Ala(3)-Ala(4)-Ala(5)) in DMSO (C. Baysal and H. Meirovitch, Biopolymers, 2000, vol. 53, pp. 423-433). In the present paper, our methodology is applied to the cyclic heptapeptides axinastatin 2 [cyclo(Asn(1)-Pro(2)-Phe(3)-Val(4)-Leu(5)-Pro(6)-Val(7))] and axinastatin 3 [cyclo(Asn(1)-Pro(2)-Phe(3)-Ile(4)-Leu(5)-Pro(6)-Val(7))], in DMSO, which were studied by nmr by Mechnich et al. (Helvetica Chimica Acta, 1997, vol. 80, pp. 1338-1354). The calculations for axinastatin 2 show that special ASPs should be optimized for the partially charged side-chain atoms of Asn while the rest of the atoms take their values derived in our previous work; this suggests that similar optimization might be needed for other side chains as well. The solution structures of these peptides are obtained ab initio (i.e., without using experimental restraints) by an extensive conformational search based on E(GRO) alone and E(*)(tot), which consists of the new set of ASPs. For E(*)(tot), the theoretical values of proton-proton distances, (3)J coupling constants, and other properties are found to agree very well with the nmr results, and they are always better than those based on E(GRO).
Collapse
Affiliation(s)
- C Baysal
- Supercomputer Computations Research Institute, Florida State University, Tallahassee, Florida 32306, USA
| | | |
Collapse
|
20
|
Baysal C, Meirovitch H. Ab initio prediction of the solution structures and populations of a cyclic pentapeptide in DMSO based on an implicit solvation model. Biopolymers 2000; 53:423-33. [PMID: 10738203 DOI: 10.1002/(sici)1097-0282(20000415)53:5<423::aid-bip6>3.0.co;2-c] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Using a recently developed statistical mechanics methodology, the solution structures and populations of the cyclic pentapeptide cyclo(D-Pro(1)-Ala(2)-Ala(3)-Ala(4)-Ala(5)) in DMSO are obtained ab initio, i.e., without using experimental restraints. An important ingredient of this methodology is a novel optimization of implicit solvation parameters, which in our previous publication [Baysal, C.; Meirovitch, H. J Am Chem Soc 1998, 120, 800-812] has been applied to a cyclic hexapeptide in DMSO. The molecule has been described by the simplified energy function E(tot) = E(GRO) + summation operator(k) sigma(k)A(k), where E(GRO) is the GROMOS force-field energy, sigma(k) and A(k) are the atomic solvation parameter (ASP) and the solvent accessible surface area of atom k. This methodology, which relies on an extensive conformational search, Monte Carlo simulations, and free energy calculations, is applied here with E(tot) based on the ASPs derived in our previous work, and for comparison also with E(GRO) alone. For both models, entropy effects are found to be significant. For E(tot), the theoretical values of proton-proton distances and (3)J coupling constants agree very well with the NMR results [Mierke, D. F.; Kurz, M.; Kessler, H. J Am Chem Soc 1994, 116, 1042-1049], while the results for E(GRO) are significantly worse. This suggests that our ASPs might be transferrable to other cyclic peptides in DMSO as well, making our methodology a reliable tool for an ab initio structure prediction; obviously, if necessary, parts of this methodology can also be incorporated in a best-fit analysis where experimental restraints are used.
Collapse
Affiliation(s)
- C Baysal
- Supercomputer Computations Research Institute, Florida State University, Tallahassee, Florida 32306, USA
| | | |
Collapse
|
21
|
|
22
|
|