Molecular modeling of an antigenic complex between a viral peptide and a class I major histocompatibility glycoprotein.

Computer simulation of the conformations of short antigenic peptides (5-10 residues) either free or bound to their receptor, the major histocompatibility complex (MHC)-encoded glycoprotein H-2 Ld, was employed to explain experimentally determined differences in the antigenic activities within a set of related peptides. Starting for each sequence from the most probable conformations disclosed by a pattern-recognition technique, several energy-minimized structures were subjected to molecular dynamics simulations (MD) either in vacuo or solvated by water molecules. Notably, antigenic potencies were found to correlate to the peptides propensity to form and maintain an overall alpha-helical conformation through regular i,i + 4 hydrogen bonds. Accordingly, less active or inactive peptides showed a strong tendency to form i,i + 3 hydrogen bonds at their N-terminal end. Experimental data documented that the C-terminal residue is critical for interaction of the peptide with H-2 Ld. This finding could be satisfactorily explained by a 3-D Q.S.A.R. analysis postulating interactions between ligand and receptor by hydrophobic forces. A 3-D model is proposed for the complex between a high-affinity nonapeptide and the H-2 Ld receptor. First, the H-2 Ld molecule was built from X-ray coordinates of two homologous proteins: HLA-A2 and HLA-Aw68, energy-minimized and studied by MD simulations. With HLA-A2 as template, the only realistic simulation was achieved for a solvated model with minor deviations of the MD mean structure from the X-ray conformation. Water simulation of the H-2 Ld protein in complex with the antigenic nonapeptide was then achieved with the template-derived optimal parameters. The bound peptide retains mainly its alpha-helical conformation and binds to hydrophobic residues of H-2 Ld that correspond to highly polymorphic positions of MHC proteins. The orientation of the nonapeptide in the binding cleft is in accordance with the experimentally determined distribution of its MHC receptor-binding residues (agretope residues). Thus, computer simulation was successfully employed to explain functional data and predicts alpha-helical conformation for the bound peptide.     

A probabilistic approach for estimating water permeability in pressure-driven membranes

A probabilistic approach is proposed to estimate water permeability in a cellulose triacetate (CTA) membrane. Water transport across the membrane is simulated in reverse osmosis mode by means of non-equilibrium molecular dynamics (MD) simulations. Different membrane configurations obtained by an annealing MD simulation are considered and simulation results are analyzed by using a hierarchical Bayesian model to obtain the permeability of the different membranes. The estimated membrane permeability is used to predict full-scale water flux by means of a process-level Monte Carlo simulation. Based on the results, the parameters of the model are observed to converge within 5-ns total simulation time. The results also indicate that the use of unique structural configurations in MD simulations is essential to capture realistic membrane properties at the molecular scale. Furthermore, the predicted full-scale water flux based on the estimated permeability is within the same order of magnitude of bench-scale experimental measurement of 1.72×10^?5 m/s.     

磷酸化对多肽中Xaa-Pro片段肽脯酰胺键顺反异构的影响研究

多肽和蛋白质中Xaa-Pro片段肽脯酰胺键顺反异构对其构象与功能有重要影响.设计合成了一系列模型多肽及其磷酸化多肽,并采用核磁共振实验和分子动力学模拟的方法,研究了所合成多肽中肽脯酰胺键的顺反异构化.结果表明,对脯氨酸之前的Xaa残基进行侧链O-磷酸化会极大地影响该顺反异构化过程,进而调节肽链构象.此外,磷酸化使得多肽顺式构象比例增加,且当磷酸基团不带负电荷时顺式构象所占比例最大.同时,分子动力学模拟所得结果与核磁共振实验相一致,包括最稳定构象和顺反构象统计分布.磷酸基团所带电荷及其空间位阻可能是影响这类磷酸化多肽构象变化的主要因素.     

Molecular dynamics simulation of adrenocorticotropin (1-10) peptide in a solvated dodecylphosphocholine micelle

Adrenocorticotropin (ACTH) (1-10), an adrenocorticotropin hormone fragment, has been studied by molecular dynamics (MD) simulation in an NPT ensemble in an explicit dodecylphosphocholine (DPC) micelle. Two starting configurations of the peptide/micelle system, corresponding to the insertion and surface-binding modes, were used. A common equilibrated configuration, in which the peptide lies parallel to the micellar surface, was reached from both simulations. In the initial part of the simulations, distance restraints derived from NMR nuclear Overhauser enhancements were incorporated before the peptide reached an equilibrium configuration with respect to the micelle. Analyses of the trajectories from the subsequent free (unrestrained) MD simulation showed that ACTH (1-10) does not conform strictly to a helical structure. The loss of the helical structure is due to decreased intramolecular hydrogen bonding accompanied by an increase of hydrogen bonding between the amide protons of the peptide and the micellar head groups. However, the extent of the latter interaction is less pronounced than in the negatively charged SDS micelle. The final structure enhances the amphipathic nature of the peptide, facilitating better interactions at the water-hydrophobic interface. The primary hydrophobic interactions with the micelle came from the side chains of Met4, Phe7, and Trp9. All peptide bonds were either hydrated or were involved in intramolecular hydrogen bonding. The interactions with the DPC micelle, the conformation of the bound peptide, and the dynamics of the peptide, as revealed by the time correlation functions of the N-H bonds, were compared with those of the ACTH (1-10)/SDS system studied previously by MD simulations.     

Long time dynamics of Met-enkephalin: comparison of explicit and implicit solvent models

Met-enkephalin is one of the smallest opiate peptides. Yet, its dynamical structure and receptor docking mechanism are still not well understood. The conformational dynamics of this neuron peptide in liquid water are studied here by using all-atom molecular dynamics (MD) and implicit water Langevin dynamics (LD) simulations with AMBER potential functions and the three-site transferable intermolecular potential (TIP3P) model for water. To achieve the same simulation length in physical time, the full MD simulations require 200 times as much CPU time as the implicit water LD simulations. The solvent hydrophobicity and dielectric behavior are treated in the implicit solvent LD simulations by using a macroscopic solvation potential, a single dielectric constant, and atomic friction coefficients computed using the accessible surface area method with the TIP3P model water viscosity as determined here from MD simulations for pure TIP3P water. Both the local and the global dynamics obtained from the implicit solvent LD simulations agree very well with those from the explicit solvent MD simulations. The simulations provide insights into the conformational restrictions that are associated with the bioactivity of the opiate peptide dermorphin for the delta-receptor.     

Synthesis and ex vivo profiling of chemically modified cytomegalovirus CMVpp65 epitopes.

The effect of substituting unnatural hydrophobic amino acids into the critical MHC binding residues of an HLA-A*0201-restricted cytomegalovirus CMVpp65 epitope, NLVPMVATV, has been investigated. A new set of peptides containing the amino acids tert-butyl glycine (Tgl), cyclohexyl glycine (Chg), neo-pentyl glycine (Npg), cyclohexyl alanine (Cha) and cyclo leucine (Cyl), at either position 2, to mimic Leu, or position 9, to mimic Val, have been synthesised. Immunological profiling using class I MHC stabilisation assays to assess MHC binding affinity, and enzyme-linked immunospot (ELISPOT) assays to assess the ability of the modified peptides to re-stimulate a specific cytotoxic T-lymphocyte (CTL) response, compared to the native epitope, have been performed. It was found that the majority of the unnatural substitutions resulted in a decrease in either HLA-A*0201 binding affinity or cytotoxic T-cell activity. However, the HLA-A*0201 binding affinity was unrelated to the ability to re-stimulate a T-cell response. Minimisation and molecular dynamics studies proved helpful in dissecting the ELISPOT responses. Two principal peptide binding modes were found by minimisation, designated kinked and straight. Peptides that bound in a kinked conformation were poor at re-stimulating a T-cell response. Of the peptides that bound in a straight conformation, molecular dynamics (MD) simulations revealed that those capable of re-stimulating the strongest responses had the greatest degree of flexibility (as determined by RMSD values across the MD simulation) around the P6 residue, one of the residues important for T-cell receptor recognition.     

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.     

Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.

Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.     

A structure for the trimeric MHC class II-associated invariant chain transmembrane domain

The major histocompatibility complex (MHC)-associated invariant chain (Ii) contains a single transmembrane domain that forms trimers. Ii is involved in the assembly of the MHC and antigen presentation, and is thus central to the function of the immune system. Here, we show by attenuated total reflectance, Fourier transform infrared (ATR-FTIR) spectroscopy that the transmembrane domain is alpha-helical and we provide a structural model of the transmembrane domain obtained by a combination of site-specific infrared dichroism and molecular dynamics (MD) simulations. This work resolves the backbone structure of a transmembrane peptide by multiple (13)C=(18)O labelling at ten different residues. A second purely computational approach, based on MD simulations of Ii transmembrane homologous sequences, yields a similar structure that is consistent with our experimental results. The structure presented forms a left-handed coiled coil with an average helix tilt of 13(+/-6) degrees; the residue Gln47 implicated in trimer formation forms strong interhelical contacts, Thr50 points to the inside of the trimeric coil and forms a network of hydrogen bonds.     

Self-assembling of peptide/membrane complexes by atomistic molecular dynamics simulations

Model biological membranes consisting of peptide/lipid-bilayer complexes can nowadays be studied by classical molecular dynamics (MD) simulations at atomic detail. In most cases, the simulation starts with an assumed state of a peptide in a preformed bilayer, from which equilibrium configurations are difficult to obtain due to a relatively slow molecular diffusion. As an alternative, we propose an extension of reported work on the self-organization of unordered lipids into bilayers, consisting of including a peptide molecule in the initial random configuration to obtain a membrane-bound peptide simultaneous to the formation of the lipid bilayer. This strategy takes advantage of the fast reorganization of lipids, among themselves and around the peptide, in an aqueous environment. Model peptides of different hydrophobicity, CH3-CO-W2L18W2-NH2 (WL22) and CH3-CO-W2A18W2-NH2 (WA22), in dipalmitoyl-phosphatidylcholine (DPPC) are used as test cases. In the equilibrium states of the peptide/membrane complexes, achieved in time ranges of 50-100 ns, the two peptides behave as expected from experimental and theoretical studies. The strongly hydrophobic WL22 is inserted in a transmembrane configuration and the marginally apolar, alanine-based WA22 is found in two alternative states: transmembrane inserted or parallel to the membrane plane, embedded close to the bilayer interface, with similar stability. This shows that the spontaneous assembly of peptides and lipids is an unbiased and reliable strategy to produce and study models of equilibrated peptide/lipid complexes of unknown membrane-binding mode and topology.     

Epitope-Based Peptides Prediction from Proteome of Nipah Virus

The identification of MHC class II epitope-based peptides are urgently needed for appropriate vaccination against Nipah virus (NiV) because there are currently no approved vaccines for human NiV infection. In the present study, prediction and modeling of T cell epitopes of NiV antigenic proteins nucleocapsid, phosphoprotein, matrix, fusion, glycoprotein, L protein, W protein, V protein and C protein followed by the binding simulation studies of predicted highest binding scores with their corresponding MHC class II alleles were done. Immunoinformatic tool ProPred was used to predict the promiscuous MHC class II epitopes of viral antigenic proteins. PEPstr server did the 3D structure models of the epitopes and Modeller 9.10 did alleles. We docked epitope with allele structure using the AutoDock 4.2 Tool. The docked peptide–allele complex structure was optimized using molecular dynamics simulation for 5 ps with the CHARMM-22 force field using NAnoscale Molecular Dynamics program incorporated in visual molecular dynamics (VMD 1.9.2) and then evaluating the stability of complex structure by calculating RMSD values. Epitope MKLQFSLGS of Matrix protein has considerable binding energy and score with DRBI*0421 MHC class II allele. This predicted peptide has potential to induce T cell-mediated immune response and is expected to useful in designing epitope-based vaccines against NiV after further testing by wet lab studies.     

Interpretation of H-NMR Experiments on the Orientation of the Transmembrane Helix WALP23 by Computer Simulations

Orientation, dynamics, and packing of transmembrane helical peptides are important determinants of membrane protein structure, dynamics, and function. Because it is difficult to investigate these aspects by studying real membrane proteins, model transmembrane helical peptides are widely used. NMR experiments provide information on both orientation and dynamics of peptides, but they require that motional models be interpreted. Different motional models yield different interpretations of quadrupolar splittings (QS) in terms of helix orientation and dynamics. Here, we use coarse-grained (CG) molecular dynamics (MD) simulations to investigate the behavior of a well-known model transmembrane peptide, WALP23, under different hydrophobic matching/mismatching conditions. We compare experimental ²H-NMR QS (directly measured in experiments), as well as helix tilt angle and azimuthal rotation (not directly measured), with CG MD simulation results. For QS, the agreement is significantly better than previously obtained with atomistic simulations, indicating that equilibrium sampling is more important than atomistic details for reproducing experimental QS. Calculations of helix orientation confirm that the interpretation of QS depends on the motional model used. Our simulations suggest that WALP23 can form dimers, which are more stable in an antiparallel arrangement. The origin of the preference for the antiparallel orientation lies not only in electrostatic interactions but also in better surface complementarity. In most cases, a mixture of monomers and antiparallel dimers provides better agreement with NMR data compared to the monomer and the parallel dimer. CG MD simulations allow predictions of helix orientation and dynamics and interpretation of QS data without requiring any assumption about the motional model.     

Brownian and essential dynamics studies of the HIV-1 integrase catalytic domain

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computational approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.     

Structure of Met-enkephalin in explicit aqueous solution using replica exchange molecular dynamics.

Replica exchange molecular dynamics (MD) simulations of Met-enkephalin in explicit solvent reveal helical and nonhelical structures. Four predominant structures of Met-enkephalin are sampled with comparable probabilities (two helical and two nonhelical). The energy barriers between these configurations are low, suggesting that Met-enkephalin switches easily between configurations. This is consistent with the requirement that Met-enkephalin be sufficiently flexible to bind to several different receptors. Replica exchange simulations of 32 ns are shown to sample approximately five times more configurational space than constant temperature MD simulations of the same duration. The energy landscape for the replica exchange simulation is presented. A detailed study of replica trajectories demonstrates that the significant increases in temperature provided by the replica exchange technique enable transitions from nonhelical to helical structures that would otherwise be prevented by kinetic trapping. Met-enkephalin (Type Entrez Proteins; Value A61445; Service Entrez Proteins).     

Generalized ordering parameter S2 for N-H peptide bonding as measure of the conformational flexibility of proteins: comparison of algorithms of S2 calculation from molecular dynamics simulation data

NMR spectroscopy in combination with molecular dynamics (MD) simulation technique is a powerful tool for investigation of biomolecule conformations. In the present study we compare the two algorithms of calculation of the generalized order parameter S2 for the peptide N-H bond from MD simulation data. The program based on this algorithm has been written and the effect of an adjustible parameter on results of calculation was studied. HIV-1 protease MD simulation was performed during 7.36 ns using GROMACS 3.1.4. software package. The calculated order parameter values are in good agreement with those obtained by other authors using NMR.     

Solution conformation of a cyclic neurokinin antagonist: A nmr and molecular dynamics study

The solution structure of a hexapeptide, cyclo(Gln-Trp-Phe-Gly-Leu-Met), which is a selective NK-2 antagonist, has been studied by a combination of two-dimensional nmr and molecular dynamics (MD) techniques. The simulation based on nmr and MD data resulted in the convergence to a family of structures. Free molecular dynamics for 50 ps in the presence of DMSO solvent molecules shows that the structure is energetically stable. One intramolecular hydrogen bond between the amide proton of Gin and the carbonyl oxygen of Gly was revealed. This result is consistent with the results from the measurement of the temperature coefficient of the amide protons. The extent of intermolecular hydrogen bonding between the amide protons of the peptide and DMSO was also revealed by the free MD simulation. The resulting structure of the cyclic peptide contains a variation type I′ β-turn in the Gly-Leu-Met-Gln segment. Comparison of the structure of this peptide with that of other NK-2 antagonist cyclic hexapeptides was made, and the activity of cyclic antagonists appears to be inversely related to the conformational rigidity of the cyclic peptides. © 1994 John Wiley & Sons, Inc.     

Molecular dynamics of water in the neighborhood of aquaporins

Present knowledge obtained by molecular dynamics (MD) simulation studies regarding the dynamics of water, both in the vicinity of biological membranes and within the proteinaceous water channels, also known as aquaporins (AQPs), is reviewed. A brief general summary of the water models most extensively employed in MD simulations (SPC, SPC/E, TIP3P, TIP4P), indicating their most relevant pros and cons, is likewise provided. Structural considerations of water are also discussed, based on different order parameters, which can be extracted from MD simulations as well as from experiments. Secondly, the behaviour of water in the neighbourhood of membranes by means of molecular dynamics simulations is addressed. Consequently, the comparison with previous experimental evidence is pointed out. In living cells, water is transported across the plasma membrane through the lipid bilayer and the aforementioned AQPs, which motivates this review to focus mostly on MD simulation studies of water within AQPs. Relevant contributions explaining peculiar properties of these channels are discussed, such as selectivity and gating. Water models used in these studies are also summarised. Finally, based on the information presented here, further MD studies are encouraged.     

Brownian and Essential Dynamics Studies of the HIV-1 Integrase Catalytic Domain

Abstract

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computional approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a “gating”-like motion with respect to the active site.     

Computer simulations of the flexibility of a series of synthetic cyclic peptide analogues

The flexibility of a series of cyclic peptides derived from the epitope of a snake toxin is investigated using computer simulation techniques. Molecular dynamics (MD) simulations and vibrational analyses are performed on chemically constrained peptides modeled on the parent protein loop. In the 50 ps MD simulations, large variations in the atomic fluctuations are seen between the peptides, and can be related to the nature of the chemical constraints present in the molecules. Normal mode analyses are performed on energy-minimized configurations derived from the dynamics trajectories. The atomic fluctuations calculated from the normal modes are about 30% of those of the molecular dynamics for the more flexible peptides and 70% for the more constrained molecules. The calculated differences in flexibility between the molecules are much less significant in the harmonic approximation. © 1993 John Wiley & Sons, Inc.     

Cutoff size need not strongly influence molecular dynamics results for solvated polypeptides

The correct treatment of van der Waals and electrostatic nonbonded interactions in molecular force fields is essential for performing realistic molecular dynamics (MD) simulations of solvated polypeptides. The most computationally tractable treatment of nonbonded interactions in MD utilizes a spherical distance cutoff (typically, 8-12 A) to reduce the number of pairwise interactions. In this work, we assess three spherical atom-based cutoff approaches for use with all-atom explicit solvent MD: abrupt truncation, a CHARMM-style electrostatic shift truncation, and our own force-shifted truncation. The chosen system for this study is an end-capped 17-residue alanine-based alpha-helical peptide, selected because of its use in previous computational and experimental studies. We compare the time-averaged helical content calculated from these MD trajectories with experiment. We also examine the effect of varying the cutoff treatment and distance on energy conservation. We find that the abrupt truncation approach is pathological in its inability to conserve energy. The CHARMM-style shift truncation performs quite well but suffers from energetic instability. On the other hand, the force-shifted spherical cutoff method conserves energy, correctly predicts the experimental helical content, and shows convergence in simulation statistics as the cutoff is increased. This work demonstrates that by using proper and rigorous techniques, it is possible to correctly model polypeptide dynamics in solution with a spherical cutoff. The inherent computational advantage of spherical cutoffs over Ewald summation (and related) techniques is essential in accessing longer MD time scales.