首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
    
β‐Peptides possess the ability to fold into secondary structure elements, and this property, together with resistance to biodegradation, makes these compounds interesting for pharmaceutical applications. Recently, a novel class of β‐peptides containing methylidene moieties was described. The GROMOS 53A6 force field was used to simulate the folding equilibrium of a β3‐hexapeptide with methylidene (CH2?) groups at all six CA‐atoms. Due to the rotational barriers induced by these methylidene groups, the helical secondary‐structure elements, normally found in β3‐peptides, are disfavored in this molecule. Simulations, started from fully extended and 314‐helical conformations, showed that the molecule adopts a complete 28‐helix for ca. 5% of the time and partial 28‐helical conformations for ca. 20% of the time. Yet, as suggested by experiments, the folding equilibrium is dominated by unfolded conformations.  相似文献   

2.
Molecular dynamics simulations were performed on four large ring cyclodextrins (LR-CDs), CD14, CD21, CD26 and CD28, both in gas phase and in water solution. Four different force fields (parm94, parm99, glycam2000a and MM3*) were used to examine their dependence with the results. The differences were not significant when parm99 and glycam2000a were used. Parm94 and MM3* results differ considerably due to inadequate bending parameters (parm94) or due to a strong stabilization by the intramolecular hydrogen bonds between hydroxyl groups. Simulations in the gas phase show that the CDs are elongated and twisted, and that cavities typical for the native CDs (CDn, n = 6, 7, 8) are not present. Simulations in water solution produced average structures that do not correspond to the conformations in the crystalline state. The LR-CDs are highly flexible and this could favor the formation of inclusion complexes through the intermediation of some other molecules that could fit to the smaller and more specific cavities they have.On sabbatical leave from the Institute of Organic Chemistry with Center of Phytochemistry, Bulgarian Academy of Science, ul. Acad. G. Bonchev, bloc 9, 1113 Sofia, Bulgaria  相似文献   

3.
Cyclic peptides are therapeutically attractive due to their high bioavailability, potential selectivity, and scaffold novelty. Furthermore, the presence of D-residues induces conformational preferences not followed by peptides consisting of naturally abundant L-residues. Therefore, comprehending how amino acids induce turns in peptides, subsequently facilitating cyclization, is significant in peptide design. Here, we performed 20-ns explicit-solvent molecular dynamics simulations for three diastereomeric peptides with stereochemistries: LLLLL, LLLDL, and LDLDL. Experimentally LLLLL and LDLDL readily cyclize, whereas LLLDL cyclizes in low yield. Simulations at 310 K produced conformations with inter-terminal hydrogen bonds that correlated qualitatively with the experimental cyclization trend. Energies obtained for representative structures from quantum chemical (B3LYP/PCM/cc-pVTZ//HF/6-31G*) calculations predicted pseudo-cyclic and extended conformations as the most stable for LLLLL and LLLDL, respectively, in agreement with the experimental data. In contrast, the most stable conformer predicted for peptide LDLDL was not a pseudo-cyclic structure. Moreover, D-residues preferred the experimentally less populated αL rotamers even when simulations were performed at a higher temperature and with strategically selected starting conformations. Energies calculated with molecular mechanics were consistent only with peptide LLLLL. Thus, the conformational preferences obtained for the all L-amino acid peptide were in agreement with the experimental observations. Moreover, refinement of the force field is expected to provide far-reaching conformational sampling of peptides containing D-residues to further develop force field-based conformational-searching methods. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.  相似文献   

4.
    
Glucokinase (GK) plays a critical role in maintaining glucose homeostasis in the human liver and pancreas. In the liver, the activity of GK is modulated by the glucokinase regulatory protein (GKRP) which functions as a competitive inhibitor of glucose to bind to GK. Moreover, the inhibitory intensity of GKRP–GK is suppressed by fructose 1-phosphate (F1P), and reinforced by fructose 6-phosphate (F6P). Here, we employed a series of computational techniques to explore the interactions of fructose phosphates with GKRP. Calculation results reveal that F1P and F6P can bind to the same active site of GKRP with different binding modes, and electrostatic interaction provides a major driving force for the ligand binding. The presence of fructose phosphate severely influences the motions of protein and the conformational space, and the structural change of sugar phosphate influences its interactions with GKRP, leading to a large conformational rearrangement of loop2 in the SIS2 domain. In particular, the binding of F6P to GKRP facilitates the protruding loop2 contacting with GK to form the stable GK–GKRP complex. The conserved residues 179–184 of GKRP play a major role in the binding of phosphate group and maintaining the stability of GKRP. These results may provide deep insight into the regulatory mechanism of GKRP to the activity of GK.  相似文献   

5.
    
A molecular‐dynamics (MD) simulation study of two heptapeptides containing α‐ and β‐amino acid residues is presented. According to NMR experiments, the two peptides differ in dominant fold when solvated in MeOH: peptide 3 adopts predominantly β‐hairpin‐like conformations, while peptide 8 adopts a 14/15‐helical fold. The MD simulations largely reproduce the experimental data. Application of NOE atom? atom distance restraining improves the agreement with experimental data, but reduces the conformational sampling. Peptide 3 shows a variety of conformations, while still agreeing with the NOE and 3J‐coupling data, whereas the conformational ensemble of peptide 8 is dominated by one helical conformation. The results confirm the suitability of the GROMOS 54A7 force field for simulation or structure refinement of mixed α/β‐peptides in MeOH.  相似文献   

6.
    
The ability to design properly folded β‐peptides with specific biological activities requires detailed insight into the relationship between the amino acid sequence and the secondary and/or tertiary structure of the peptide. One of the most frequently used spectroscopic techniques for resolving the structure of a biomolecule is NMR spectroscopy. Because only signal intensities and frequencies are recorded in the experiment, a conformational interpretation of the recorded data is not straightforward, especially for flexible molecules. The occurrence of conformational and/or time averaging, and the limited amount and accuracy of experimental data hamper the precise conformational determination of a biomolecule. In addition, the relation between experimental observables with the underlying conformational ensemble is often only approximately known, thereby aggravating the difficulty of structure determination of biomolecules. The problematic aspects of structure refinement based on NMR nuclear Overhauser effect (NOE) intensities and 3J‐coupling data are illustrated by simulating a β‐octapeptide in explicit MeOH and H2O as solvents using three different force fields. NMR Data indicated that this peptide would fold into a 314‐helix in MeOH and into a hairpin in H2O. Our analysis focused on the conformational space visited by the peptide, on structural properties of the peptide, and on agreement of the MD trajectories with available NMR data. We conclude that 1) although the 314‐helical structure is present when the peptide is solvated in MeOH, it is not the only relevant conformation, and that 2) the NMR data set available for the peptide, when solvated in H2O, does not provide sufficient information to derive a single secondary structure, but rather a multitude of folds that fulfill the NOE data set.  相似文献   

7.
Classical molecular dynamics (MD) and non-equilibrium steered molecular dynamics (SMD) simulations were performed on the molecular structure of the potassium channel KcsA using the GROMOS 87 force fields. Our simulations focused on mechanistic and dynamic properties of the permeation of potassium ions through the selectivity filter of the channel. According to the SMD simulations a concerted movement of ions inside the selectivity filter from the cavity to extracellular side depends on the conformation of the peptide linkage between Val76 and Gly77 residues in one subunit of the channel. In SMD simulations, if the carbonyl oxygen of Val76 is positioned toward the ion bound at the S3 site (gate-opened conformation) the net flux of ions through the filter is observed. When the carbonyl oxygen leaped out from the filter (gate-closed conformation), ions were blocked at the S3 site and no flux occurred. A reorientation of the Thr75-Val76 linkage indicated by the CHARMM-based MD simulations performed Berneche and Roux [(2005) Structure 13:591–600; (2000) Biophys J 78:2900–2917] as a concomitant process of the Val76-Gly77 conformational interconversion was not observed in our GROMOS-based MD simulations.  相似文献   

8.
An original procedure approach taking into account the implementation of parameters determined, using calculations based on density functional theory, for the amidocyano-pyridinium methylide in the MM2 augmented harmonic potential function has been proposed. A good agreement between theoretical force field calculation and X-ray diffraction data has been observed. Thus an empirical force field for cycloimmonium ylides has been established. It provides good quality geometries for cycloimmonium ylide molecules by energy minimization. In this study we proposed a new MM2 augmented atom-type for the ylidic carbon atom. To our knowledge no attempt has been done in this way for such organic systems. Thus, we have shown that parameterization established by the DFT method is able to reproduce or to predict with good accuracy the structures of the cycloimmonium ylide compounds. This study also includes a full conformational analysis.  相似文献   

9.
Giordano Lesma 《Tetrahedron》2007,63(25):5567-5578
New peptidomimetics containing the Tic moiety were synthesized in enantiomerically pure form and their conformational features were studied by NMR, IR, and molecular modeling techniques. The presence of a reverse turn conformation was observed in all the structures, suggesting the key role of the scaffold as reverse turn inducer. In particular, all the analyses led to the conclusion that a β-turn conformation is mostly stabilized in tetrapeptide mimetic 4b and in hexapeptide mimetics 5a,b. In the case of 5a,b, the C1 stereochemistry plays a central role in determining stable conformations, supporting the formation of a β-hairpin arrangement with a 14-membered intramolecular hydrogen bond ring only in 5b.  相似文献   

10.
高岭石-水体系中水分子结构的分子动力学模拟   总被引:1,自引:0,他引:1  
以Hendricks模型为初始结构, 利用CLAYFF力场对高岭石-水体系进行无晶体学限制的分子动力学模拟. 结果表明, 层间水有三种类型: I型类似于Costanzo提出的“洞水”分子, 其HH矢量(水分子中从一个氢原子位置指向另一个氢原子位置的方向矢量)平行于(001)平面, 而C2轴稍微倾斜于(001)面法线; II型类似于“连接水”, 一个氢氧键指向临近的层间四面体氧形成氢键, 另一个氢氧键与(001)面近似平行; III型水分子在层间近似保持为竖直状, 一个氢与层间四面体氧形成氢键, 而另一个氢与对面层的羟基氧形成氢键. 高岭石羟基氢沿(001)晶面法线的浓度曲线显示一部分羟基指向变为近似平行于(001)面, 羟基氧因此能够暴露出来与层间水分子氢形成氢键. 此外, 模拟中还观察到部分II型水分子氧偏离于层间的平均位置而更靠近四面体层, 这和Costanzo的实验结果一致, 可能是X射线谱图中(002)弱衍射峰出现的原因.  相似文献   

11.
During recent decades it has become feasible to simulate the dynamics of molecular systems on a computer. The method of molecular dynamics (MD) solves Newton's equations of motion for a molecular system, which results in trajectories for all atoms in the system. From these atomic trajectories a variety of properties can be calculated. The aim of computer simulations of molecular systems is to compute macroscopic behavior from microscopic interactions. The main contributions a microscopic consideration can offer are (1) the understanding and (2) interpretation of experimental results, (3) semiquantitative estimates of experimental results, and (4) the capability to interpolate or extrapolate experimental data into regions that are only difficultly accessible in the laboratory. One of the two basic problems in the field of molecular modeling and simulation is how to efficiently search the vast configuration space which is spanned by all possible molecular conformations for the global low (free) energy regions which will be populated by a molecular system in thermal equilibrium. The other basic problem is the derivation of a sufficiently accurate interaction energy function or force field for the molecular system of interest. An important part of the art of computer simulation is to choose the unavoidable assumptions, approximations and simplifications of the molecular model and computational procedure such that their contributions to the overall inaccuracy are of comparable size, without affecting significantly the property of interest. Methodology and some practical applications of computer simulation in the field of (bio)chemistry will be reviewed.  相似文献   

12.
13.
    
The correlation between β2‐, β3‐, and β2,3‐amino acid‐residue configuration and stability of helix and hairpin‐turn secondary structures of peptides consisting of homologated proteinogenic amino acids is analyzed (Figs. 1–3). To test the power of Zn2+ ions in fortifying and/or enforcing secondary structures of β‐peptides, a β‐decapeptide, 1 , four β‐octapeptides, 2 – 5 , and a β‐hexadecapeptide, 10 , have been devised and synthesized. The design was such that the peptides would a) fold to a 14‐helix ( 1 and 3 ) or a hairpin turn ( 2 and 4 ), or form neither of these two secondary structures (i.e., 5 ), and b) carry the side chains of cysteine and histidine in positions, which will allow Zn2+ ions to use their extraordinary affinity for RS? and the imidazole N‐atoms for stabilizing or destabilizing the intrinsic secondary structures of the peptides. The β‐hexadecapeptide 10 was designed to a) fold to a turn, to which a 14‐helical structure is attached through a β‐dipeptide spacer, and b) contain two cysteine and two histidine side chains for Zn complexation, in order to possibly mimic a Zn‐finger motif. While CD spectra (Figs. 6–8 and 17) and ESI mass spectra (Figs. 9 and 18) are compatible with the expected effects of Zn2+ ions in all cases, it was shown by detailed NMR analyses of three of the peptides, i.e., 2, 3, 5 , in the absence and presence of ZnCl2, that i) β‐peptide 2 forms a hairpin turn in H2O, even without Zn complexation to the terminal β3hHis and β3hCys side chains (Fig. 11), ii) β‐peptide 3 , which is present as a 14‐helix in MeOH, is forced to a hairpin‐turn structure by Zn complexation in H2O (Fig. 12), and iii) β‐peptide 5 is poorly ordered in CD3OH (Fig. 13) and in H2O (Fig. 14), with far‐remote β3hCys and β3hHis residues, and has a distorted turn structure in the presence of Zn2+ ions in H2O, with proximate terminal Cys and His side chains (Fig. 15).  相似文献   

14.
15.
    
Free‐energy differences govern the equilibrium between bound and unbound states of a host and its guest molecules. The understanding of the underlying entropic and enthalpic contributions, and their complex interplay are crucial for the design of new drugs and inhibitors. In this study, molecular dynamics (MD) simulations were performed with inclusion complexes of α‐cyclodextrin (αCD) and three monosubstituted benzene derivatives to investigate host–guest binding. αCD Complexes are an ideal model system, which is experimentally and computationally well‐known. Thermodynamic integration (TI) simulations were carried out under various conditions for the free ligands in solution and bound to αCD. The two possible orientations of the ligand inside the cavity were investigated. Agreement with experimental data was only found for the more stable orientation, where the substituent resides inside the cavity. The better stability of this conformation results from stronger Van der Waals interactions and a favorable antiparallel host–guest dipole–dipole alignment. To estimate the entropic contributions, simulations were performed at three different temperatures (250, 300, and 350 K) and using positional restraints for the host. The system was found to be insensitive to both factors, due to the large and symmetric cavity of αCD, and the nondirectional nature of the host–guest interactions.  相似文献   

16.
To further study the preference of the antiperiplanar (ap) conformation in α‐fluoro‐amide groups, two β‐peptides, 1 and 2 , containing a (2‐F)‐β3hAla and a (2‐F)‐β2hPhe residue, have been synthesized. Their NMR‐solution structures in CD3OH were determined and compared with those of non‐F‐substituted analogs, 3 and 4a . While we have found in a previous investigation (Helv. Chim. Acta 2005 , 88, 266) that a stereospecifically introduced F‐substituent in the central position of a βheptapeptide is capable of ‘breaking’ the 314‐helical structure by enforcing the F? C? C?O ap‐conformation, we could now demonstrate that the same procedure leads to a structure with the unfavorable ca. 90° F? C? C?O dihedral angle, enforced by the 314‐helical folding in a βtridecapeptide (cf. 1 ; Fig. 4). This is interpreted as a consequence of cooperative folding in the longer β‐peptide. A F‐substituent placed in the turn section of a β‐peptidic hairpin turn was shown to be in an ap‐arrangement with respect to the neighboring C?O bond (cf. 2 ; Fig. 7). Analysis of the non‐F‐substituted β‐tetrapeptides (with helix‐preventing configurations of the two central β2/β3‐amino acid residues) provides unusually tight hairpin structural clusters (cf. 3 and 4a ; Figs. 8 and 9). The skeleton of the β‐tetrapeptide H‐(R)β3hVal‐(R)β2hVal‐(R)β3hAla‐(S)β3hPhe‐OH ( 4a ) is proposed as a novel, very simple backbone structure for mimicking α‐peptidic hairpin turns.  相似文献   

17.
以量子化学计算作为起点, 为最简单的糖类分子——乙醇醛开发了两套分子力学力场参数: 基于肽类的力场和基于醛类的力场. 分子动力学模拟结果表明, 所开发的类醛力场参数能够较好地描述乙醇醛分子在水中的结构以及水分子在其周围的分布. 通过瞬时简正模式分析, 得到了3N-6个模式的瞬时振动频率和振动跃迁偶极矩等振动光谱参数的统计分布及其相关性. 结合量子化学计算和分子动力学模拟对生物分子体系的多元振动光谱参数进行预测和评估, 为从化学键水平出发模拟宽带飞秒二维红外相关光谱提供了一个新方法.  相似文献   

18.
    
4‐Hydroxyphenylpyruvate dioxygenase is a relevant target in both pharmaceutical and agricultural research. We report on molecular dynamics simulations and free energy calculations on this enzyme, in complex with 12 inhibitors for which experimental affinities were determined. We applied the thermodynamic integration approach and the more efficient one‐step perturbation. Even though simulations seem well converged and both methods show excellent agreement between them, the correlation with the experimental values remains poor. We investigate the effect of slight modifications on the charge distribution of these highly conjugated systems and find that accurate models can be obtained when using improved force field parameters. This study gives insight into the applicability of free energy methods and current limitations in force field parameterization. © 2011 Wiley Periodicals, Inc. J Comput Chem 2011  相似文献   

19.
In order to study the conformational space of new natural products derivatives using molecular mechanics (MM) and molecular dynamics (MD), the Gromos force field has been expanded to include epoxidic systems. The parameterization and validation of Gromos were done to simulate 22, 23 epoxides in brassinosteroids analogs. The parameters were derived with an emphasis on the dependence between energy and dihedral angle due to its relevance in the conformational analysis. Molecular dynamics simulations of two model systems similar to those of interest were performed to validate the force field with the proposed parameters. Excellent agreement has been obtained between the MD simulation and the results of a potential energy surface (PES) calculated at B3LYP/6-31G** level.  相似文献   

20.
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号