Peptide T revisited: conformational mimicry of epitopes of anti-HIV proteins.

Peptide T (ASTTTNYT), a fragment corresponding to residues 185-192 of gp120, the coat protein of HIV, is endowed with several biological properties in vitro, notably inhibition of the binding of both isolated gp120 and HIV-1 to the CD4 receptor, and chemotactic activity. Based on previous nuclear magnetic resonance (NMR) studies performed in our laboratory, which were consistent with a regular conformation of the C-terminal pentapeptide, and SAR studies showing that the C-terminal pentapeptide retains most of the biological properties, we designed eight hexapeptides containing in the central part either the TNYT or the TTNY sequence, and charged residues (D/E/R) at the two ends. Conformational analysis based on NMR and torsion angle dynamics showed that all peptides assume folded conformations. albeit with different geometries and stabilities. In particular, peptides carrying an acidic residue at the N-terminus and a basic residue at the C-terminus are characterized by stable helical structures and retain full chemotactic activity. The solution conformation of peptide ETNYTR displays strong structural similarity to the region 19-26 of both bovine pancreatic and bovine seminal ribonuclease, which are endowed with anti-HIV activity. Moreover, the frequent occurrence, in many viral proteins, of TNYT and TTNY, the two core sequences employed in the design of the hexapeptides studied in the present work, hints that the sequence of the C-terminal pentapeptide TTNYT is probably representative of a widespread viral recognition motif.     

Molecular dynamics studies of a DNA-binding protein: 1. A comparison of the trp repressor and trp aporepressor aqueous simulations.

The results of two 30-ps molecular dynamics simulations of the trp repressor and trp aporepressor proteins are presented in this paper. The simulations were obtained using the AMBER molecular mechanical force field and in both simulations a 6-A shell of TIP3P waters surrounded the proteins. The trp repressor protein is a DNA-binding regulatory protein and it utilizes a helix-turn-helix (D helix-turn-E helix) motif to interact with DNA. The trp aporepressor, lacking two molecules of the L-tryptophan corepressor, cannot bind specifically to DNA. Our simulations show that the N- and C-termini and the residues in and near the helix-turn-helix motifs are the most mobile regions of the proteins, in agreement with the X-ray crystallographic studies. Our simulations also find increased mobility of the residues in the turn-D helix-turn regions of the proteins. We find the average distance separating the DNA-binding motifs to be larger in the repressor as compared to the aporepressor. In addition to examining the protein residue fluctuations and deviations with respect to X-ray structures, we have also focused on backbone dihedral angles and corepressor hydrogen-bonding patterns in this paper.     

基于温度副本交换分子动力学模拟方法研究温度致H1小肽结构变化特征

采用GROMOS43A1力场,用温度副本交换分子动力学模拟方法研究水溶液中H1小肽在4个不同温度下的结构特征.选择H1小肽的初始构象分别为α螺旋和β折叠片,完成了两组独立的36个温度副本交换的分子动力学模拟,一组从α螺旋出发的模拟用来对该小肽的结构特征进行研究,另一组从β折叠片出发的模拟用于验证构象采样的收敛性,每个副本的模拟时间为300ns,共计模拟时间长达21.6μs.在此基础上,研究了H1小肽在温度300K、330K、350K和370K下的结构特征,分析了其主链二面角分布、天然氢键数、β转角的形成概率以及不同温度下偏好采样构象的变化特征等.模拟结果表明,在4个不同温度下,均能够采样到同β折叠片结构的Cα原子均方根偏差最小为0.05nm的构象类,该构象类在4个不同温度300K、330K、350K和370K下分别包含了全部构象的39%、23%、13%和11%.GROMOS43A1力场在刻画小肽的结构方面具有一定的精度,但是在描述氢键方面仍需要加强,H1小肽在不同温度下结构特征的比较能够为分子力场的优化提供重要的帮助.     

Peptide Free Energy Landscapes Calibrated by Molecular Orbital Calculations

Free energy landscapes of peptide conformations werecalibrated by ab initiomolecular orbital calculations, after enhancedconformational sampling using the multicanonical molecular dynamicssimulations. Three different potentials of mean force for an isolateddipeptide were individually obtained using the conventional force fields,AMBER parm94, AMBER parm96, and CHARMm22. Each potential ofmean force was calibrated based on the umbrella sampling algorithm fromthe adiabatic energy map that was calculated separately by the abinitiomolecular orbital method. All the calibrated potentials of mean forcecoincided well. The calibration was applied to a peptide in explicit water,and the calibrated free energy landscapes did not depend on the force fieldused in conformational sampling, as far as the conformational space waswell sampled.     

Conformational studies of human islet amyloid peptide using molecular dynamics and simulated annealing methods

Molecular dynamics simulations and simulated annealing in vacuum, model aqueous solution, and simulated membrane were used to analyze the conformational preferences of a segment spanning 20–29 residues of human islet amyloid polypeptide, [referred to as IAPP^H(20–29)]. Molecular dynamics simulations were conducted at 300 K on IAPP^H(20–29). The minimum energy conformers obtained in model aqueous solution and vacuum exhibited similar structures. Even in the absence of any constraints on peptide bonds, trans conformation was preferred consistently by all the peptide bonds. Analysis of the minimum energy conformers indicated that IAPP^H(20–29) showed a strong preference for turn structures in all the environments. These turn structures were stabilized by the formation of hydrogen bonds between the backbone amide and carbonyl groups. A good agreement was found between the results obtained from the molecular dynamics simulation and solid-state nmr experimental studies. © 1998 John Wiley & Sons, Inc. Biopoly 45: 9–20, 1998     

Improved model of hydrated calcium ion for molecular dynamics simulations using classical biomolecular force fields

Calcium ions (Ca²⁺) play key roles in various fundamental biological processes such as cell signaling and brain function. Molecular dynamics (MD) simulations have been used to study such interactions, however, the accuracy of the Ca²⁺ models provided by the standard MD force fields has not been rigorously tested. Here, we assess the performance of the Ca²⁺ models from the most popular classical force fields AMBER and CHARMM by computing the osmotic pressure of model compounds and the free energy of DNA–DNA interactions. In the simulations performed using the two standard models, Ca²⁺ ions are seen to form artificial clusters with chloride, acetate, and phosphate species; the osmotic pressure of CaAc₂ and CaCl₂ solutions is a small fraction of the experimental values for both force fields. Using the standard parameterization of Ca²⁺ ions in the simulations of Ca²⁺‐mediated DNA–DNA interactions leads to qualitatively wrong outcomes: both AMBER and CHARMM simulations suggest strong inter‐DNA attraction whereas, in experiment, DNA molecules repel one another. The artificial attraction of Ca²⁺ to DNA phosphate is strong enough to affect the direction of the electric field‐driven translocation of DNA through a solid‐state nanopore. To address these shortcomings of the standard Ca²⁺ model, we introduce a custom model of a hydrated Ca²⁺ ion and show that using our model brings the results of the above MD simulations in quantitative agreement with experiment. Our improved model of Ca²⁺ can be readily applied to MD simulations of various biomolecular systems, including nucleic acids, proteins and lipid bilayer membranes. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 752–763, 2016.     

Conformational Flexibility of Two RNA Trimers Explored by Computational Tools and Database Search

Abstract

Two RNA sequences, AAA and AUG, were studied by the conformational search program CICADA and by molecular dynamics (MD) in the framework of the AMBER force field, and also via thorough PDB database search. CICADA was used to provide detailed information about conformers and conformational interconversions on the energy surfaces of the above molecules. Several conformational families were found for both sequences. Analysis of the results shows differences, especially between the energy of the single families, and also in flexibility and concerted conformational movement. Therefore, several MD trajectories (altogether 16 ns) were run to obtain more details about both the stability of conformers belonging to different conformational families and about the dynamics of the two systems. Results show that the trajectories strongly depend on the starting structure. When the MD start from the global minimum found by CICADA, they provide a stable run, while MD starting from another conformational family generates a trajectory where several different conformational families are visited. The results obtained by theoretical methods are compared with the thorough database search data. It is concluded that all except for the highest energy conformational families found in theoretical result also appear in experimental data.

Registry numbers:

adenylyl-(3′ →5′)-adenylyl-(3′ →5′)-adenosine [917-44-2]

adenylyl-(3′ →5′)-uridylyl-(3′ →5′)-guanosine [3494-35-7]     

Conformation-activity relationship of peptide T and new pseudocyclic hexapeptide analogs.

Peptide T (ASTTTNYT), a segment corresponding to residues 185-192 of gp120, the coat protein of HIV, has several important biological properties in vitro that have stimulated the search for simpler and possibly more active analogs. We have previously shown that pseudocyclic hexapeptide analogs containing the central residues of peptide T retain considerable chemotactic activity. We have now extended the design of this type of analogs to peptides containing different aromatic residues and/or Ser in lieu of Thr. The complex conformation-activity relationship of these analogs called for a reexamination of the basic conformational tendencies of peptide T itself. Here, we present an exhaustive NMR conformational study of peptide T in different media. Peptide T assumes a gamma-turn in aqueous mixtures of ethylene glycol, a type-IV beta-turn conformation in aqueous mixtures of DMF, and a type-II beta-turn conformation in aqueous mixtures of DMSO. The preferred conformations for the analogs were derived from modeling, starting from the preferred conformations of peptide T. The best models derived from the gamma-turn conformation of peptide T are those of peptides XII (DSNYSR), XIII (ETNYTK) and XVI (ESNYSR). The best models derived from the type-IV beta-turn conformation of peptide T are those of peptides XIV (KTTNYE) and XV (DSSNYR). No low-energy models could be derived starting from the type-II beta-turn conformation of peptide T. The analogs with the most favored conformations are also the most active in the chemotactic test.     

Conformational Studies of Tachykinin Peptides Using NMR Spectroscopy

A conformational analysis in water and DMSO of two tachykinin family peptides (scyliorhinin I (ScyI) and scyliorhinin II (ScyII)) was carried out by 1D and 2D NMR (DQF-COSY, TOCSY, HMQC, HMBC, NOESY and ROESY) and molecular dynamics calculation methods. In DMSO, two groups of conformations (major and minor) were obtained for both peptides based on the experimental data. The conformations proposed for ScyI represent a folded structure, which shows certain similarities to the structures reported for other NK-1 and NK-2 tachykinin agonists. In water ScyII displays a flexible, extended structure, whereas in DMSO the structure is more compact and, in the fragment from the centre to the C-terminus, several -turns may be present.     

Conformational Studies of Tachykinin Peptides Using NMR Spectroscopy

Summary A conformational analysis in water and DMSO of two tachykinin family peptides (scyliorhinin I (ScyI) and scyliorhinin II (ScyII)) was carried out by 1D and 2D NMR (DQF-COSY, TOCSY, HMQC, HMBC, NOESY and ROESY) and molecular dynamics calculation methods. In DMSO, two groups of conformations (major and minor) were obtained for both peptides based on the experimental data. The conformations proposed for ScyI represent a folded structure, which shows certain similarities to the structures reported for other NK-1 and NK-2 tachykinin agonists. In water ScyII displays a flexible, extended structure, whereas in DMSO the structure is more compact and, in the fragment from the centre to the C-terminus, several β-turns may be present.     

Structural insight into epothilones antitumor activity based on the conformational preferences and tubulin binding modes of epothilones A and B obtained from molecular dynamics simulations

Molecular dynamics simulations were employed to analyze the conformational preferences and binding modes of epothilones A and B as a source of structural information regarding the antitumor properties of these species. Our results suggest that the conformation of free and tubulin-bound epothilones is strongly influenced by the presence of a methyl group at C12 and that epothilones A and B exploit the binding cavity in a unique and different way. The binding sites of epothilones A and B share a common region of association (Leu215, Leu217, His227, Leu228, Ala231, Phe270, Gly360, and Leu361), but lead to different ligand–residue interactions. Average interaction energies predict a larger stabilization for the epothilone B–tubulin complex, which is mainly driven by the enhancement of the electrostatic component of ligand–residue interactions compared to the epothilone A–tubulin complex. These structural and energetic results can be useful to account for the activity difference between epothilones A and B, and to design more active and potent analogs that resemble the mechanism of action of epothilones against cancer cells.     

Conformational Flexibility and Peptide Interaction of the Translocation ATPase SecA

The SecA ATPase forms a functional complex with the protein-conducting SecY channel to translocate polypeptides across the bacterial cell membrane. SecA recognizes the translocation substrate and catalyzes its unidirectional movement through the SecY channel. The recent crystal structure of the Thermotoga maritima SecA-SecYEG complex shows the ATPase in a conformation where the nucleotide-binding domains (NBDs) have closed around a bound ADP-BeFx complex and SecA's polypeptide-binding clamp is shut. Here, we present the crystal structure of T. maritima SecA in isolation, determined in its ADP-bound form at 3.1 Å resolution. SecA alone has a drastically different conformation in which the nucleotide-binding pocket between NBD1 and NBD2 is open and the preprotein cross-linking domain has rotated away from both NBDs, thereby opening the polypeptide-binding clamp. To investigate how this clamp binds polypeptide substrates, we also determined a structure of Bacillus subtilis SecA in complex with a peptide at 2.5 Å resolution. This structure shows that the peptide augments the highly conserved β-sheet at the back of the clamp. Taken together, these structures suggest a mechanism by which ATP hydrolysis can lead to polypeptide translocation.     

谷氨酰胺活性肽营养液对大鼠小肠营养作用的研究

用氨甲蝶呤诱发小肠炎后 ,SD大鼠分别饲喂Gln含量不等的氨基酸营养液 ,结果表明 ,72h内 ,第一组 (零剂量组 )大鼠的死亡率为 6 6 .7% ,第二组 (含Gln质量分数 2 % )为 16 .7% ,第三组 (活性肽组 ,含Gln质量分数 2 0 % )和第四组 (结晶氨基酸液 ,含Gln质量分数 2 0 % )为 0 ,并且第三组和第四组的体重、蔗糖酶活、氨肽酶活、DNA和RNA的含量增加 ,说明Gln活性肽营养液对小肠具有明显的营养作用。     

Computation techniques in the conformational analysis of carbohydrates

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

Molecular modeling of formate dehydrogenase: the formation of the Michaelis complex

The formation of the reactive enzyme–substrate complex of formate dehydrogenase has been investigated by molecular dynamics techniques accounting for different conformational states of the enzyme. Simulations revealed that the transport of substrate to the active site through the substrate channel proceeds in the open conformation of enzyme due to the crucial role of the Arg284 residue acting as a vehicle. However, formate binding in the active site of the open conformation leads to the formation of a nonproductive enzyme–substrate complex. The productive Michaelis complex is formed only in the closed enzyme conformation after the substrate and coenzyme have bound, when required rigidity of the binding site and reactive formate orientation due to interactions with Arg284, Asn146, Ile122, and His332 residues is attained. Then, the high occupancy (up to 75%) of the reactive substrate–coenzyme conformation is reached, which was demonstrated by hybrid quantum mechanics/molecular mechanics simulations using various semiempirical Hamiltonians.     

Benchmarking of copper(II) LFMM parameters for studying amyloid-β peptides

Ligand field molecular mechanics (LFMM) parameters have been benchmarked for copper (II) bound to the amyloid-β_1–16 peptide fragment. Several density functional theory (DFT) optimised small test models, representative of different possible copper coordination modes, have been used to test the accuracy of the LFMM copper bond lengths and angles, resulting in errors typically less than 0.1 Å and 5°. Ligand field molecular dynamics (LFMD) simulations have been carried out on the copper bound amyloid-β_1–16 peptide and snapshots extracted from the subsequent trajectory. Snapshots have been optimised using DFT and the semi-empirical PM7 method resulting in good agreement against the LFMM calculated geometry. Analysis of substructures within snapshots shows that the larger contribution of geometrical difference, as measured by RMSD, lies within the peptide backbone, arising from differences in DFT and AMBER, and the copper coordination sphere is reproduced well by LFMM. PM7 performs excellently against LFMM with an average RMSD of 0.2 Å over 21 tested snapshots. Further analysis of the LFMD trajectory shows that copper bond lengths and angles have only small deviations from average values, with the exception of a carbonyl moiety from the N-terminus, which can act as a weakly bound fifth ligand.     

Vibrations-determined properties of green fluorescent protein

The physicochemical characteristics of the green fluorescent protein (GFP), including the thermodynamic properties (entropy, enthalpy, Gibbs' free energy, heat capacity), normal mode vibrations, and atomic fluctuations, were investigated. The Gaussian 03 computational chemistry program was employed for normal mode analysis using the AMBER force field. The thermodynamic parameters and atomic fluctuations were then calculated from the vibrational eigenvalues (frequencies) and eigenvectors. The regions of highest rigidity were shown to be the beta-sheet barrel with the central alpha-helix, which bears the chromophore. The most flexible parts of the GFP molecule were the outlying loops that cover the top and bottom of the beta-barrel. This way, the balance between rigidity and flexibility is maintained, which is the optimal relationship for protein stability in terms of Gibbs' free energy. This dual-schemed structure satisfies the requirements for GFP function. In this sense, the structure of GFP resembles a nanoscale drum: a stiff cylinder with flexible vibrating end(s).     

An online parameter and property database for the TraPPE force field

The transferable potentials for phase equilibria (TraPPE) force field aims to be accurate, computationally efficient and applicable to a wide range of chemical compounds, state points and thermophysical properties. When new users wish to implement TraPPE models into their chosen simulation program, they face several obstacles: the TraPPE models are dispersed over many separate publications and misinterpretations of the primary literature are possible; the TraPPE force field makes specific choices for standard conventions that may require non-trivial code modifications for some simulation software. Therefore, the TraPPE developers report here a resource website and online searchable parameter and property database designed to provide new and experienced users with tools for successful implementation and validation (http://www.chem.umn.edu/groups/siepmann/trappe/).     

Computational studies and peptidomimetic design for the human p53-MDM2 complex

The interaction between human p53 and MDM2 is a key event in controlling cell growth. Many studies have suggested that a p53 mimic would be sufficient to inhibit MDM2 to reduce cell growth in cancerous tissue. In order to design a potent p53 mimic, molecular dynamics (MD) simulations were used to examine the binding interface and the effect of mutating key residues in the human p53-MDM2 complex. The Generalized Born surface area (GBSA) method was used to estimate free energies of binding, and a computational alanine-scanning approach was used to calculate the relative effects in the free energy of binding for key mutations. Our calculations determine the free energy of binding for a model p53-MDM2 complex to be -7.4 kcal/mol, which is in very good agreement with the experimentally determined values (-6.6--8.8 kcal/mol). The alanine-scanning results are in good agreement with experimental data and calculations by other groups. We have used the information from our studies of human p53-MDM2 to design a beta-peptide mimic of p53. MD simulations of the mimic bound to MDM2 estimate a free energy of binding of -8.8 kcal/mol. We have also applied alanine scanning to the mimic-MDM2 complex and reveal which mutations are most likely to alter the binding affinity, possibly giving rise to escape mutants. The mimic was compared to nutlins, a new class of inhibitors that block the formation of the p53-MDM2 complex. There are interesting similarities between the nutlins and our mimic, and the differences point to ways that both inhibitors may be improved. Finally, an additional hydrophobic pocket is noted in the interior of MDM2. It may be possible to design new inhibitors to take advantage of that pocket.     

Conformational basis of N-glycosylation of proteins: conformational analysis of Ac-Asn-Ala-Thr-NH2

The structural requirements of the Asn-X-Thr(Ser) sequence for the N-glycosylation of proteins has been traced to a local conformation acting as a signal for the enzymatic process. The conformational space of the smallest in vivoN-glycosylation substrate, Ac-Asn-Ala-Thr-NH₂, has been thoroughly explored using energy calculations. All the lowest energy conformers have been characterized as bended structures.