Abstract: Computational modeling techniques and computer simulations have become a routine in biological sciences and have gained great attention from researchers. Molecular dynamics simulation is a valuable tool towards an understanding of the complex structure of biological systems, especially in the study of the flexibility of the biological molecules such as peptides or proteins. Peptides play a very important role in human physiology and control many of the processes involved in the immune system response. Designing new and optimal peptide vaccines is one of the hottest challenges of the 21(st) century science and it brings together researchers from different fields. Molecular dynamics simulations have proven to be a helpful tool assisting laboratory work, saving financial sources and opening possibilities for exploring properties of the molecular systems that are hardly accessible by conventional experimental methods. Present review is dedicated to the recent contributions in applications of molecular dynamics simulations in peptide design for immunological purposes, such as B or T cell epitopes.
Abstract: In a recent study C8gamma (complement protein) with Cys40Ala substitution and a C8alpha derived peptide with Cys164Ala substitution were co-crystallized and their binding mode was revealed. Computer modeling and molecular dynamics simulations were performed in order to check the hypothesis that the residues Ala164 of C8alpha and Ala40 of C8gamma occupied the right position if cysteine residues were in their place for disulfide bonding. Substitution of these two alanine residues with cysteine along with disulfide bond creation via molecular modeling and subsequent molecular dynamics simulation of the complex corroborated the hypothesis, which was also confirmed from recent crystallographic data. Average RMSD between backbone atoms of the indel peptide during the MD trajectory in comparison with the corresponding sequence of crystal structure of the C8alpha/gamma complex was found only 0.085 nm.
Abstract: Humanized CAMPATH-1H antibody has been found to have biological applications through the recognition of the CD52 antigen. A peptide mimotope of the CD52 antigen with the sequence T1SSPSAD7 has been co-crystallized with the CAMPATH-1H antibody. A plethora of hydrogen bond interactions were found to mediate antigen recognition. An important feature of peptide’s bound conformation was the type I β-turn found in the S3PSA6 peptide’s fragment. Paradoxically, this fact has been underestimated from other researchers. In order to further investigate the importance of this structural feature and its significance in antibody/antigen binding we have performed molecular dynamics simulations in explicit water of the T1SSPSAD7 peptide in both antibody free and bound states. We have found that the turn structure has been perfectly retained in the bound state but it was eliminated in the free state. This fact implies that the turn structure of the peptide is unstable in aqueous environment and it is induced upon antibody binding. Analysis of the trajectories revealed also several other important features of the antibody/antigen binding mode.
Notes: Get reprints or see more at:
http://stavrakoudis.econ.uoi.gr/index.php
Abstract: Resistance to cationic antimicrobial peptide polymyxin B from Gram-negative bacteria is accomplished by two-component systems (TCSs), protein complexes PmrA/PmrB and PhoP/PhoQ. PmrD is the first protein identified to mediate the connectivity between two TCSs. The 3D structure of PmrD has been recently solved by NMR and its unique fold was revealed. Here, a molecular dynamics study is presented started from the NMR structure. Numerous hydrophobic and electrostatic interactions were identified to contribute to PmrD's 3D stability. Moreover, the mobility of the five loops that connect the protein's six beta-strands has been explored. Solvent-accessible surface area calculation revealed that a Leucine-rich hydrophobic cluster of the protein stabilized the protein's structure.
Abstract: The peptide sequence YMESRADRKLAEVGRVYLFL, derived from 313-332 region of the αIIb, has been identified as a potent inhibitor of platelet aggregation and fibrinogen binding to αIIbβ3. More detailed studies have revealed that the Y313MESRADR320 sequence is the shortest octapeptide with strong inhibitory activity. This work provides insight of the solution conformation of these peptides, by performing extensive molecular dynamics simulations of 100 ns. The 8mer peptide has no stable conformation in water while the 20mer peptide retains a relative conformational stability. Analysis of side chain orientation of the RAD fragment revealed the synplanar arrangement of guanidinium and β-carboxylic groups providing a framework for explaining the similar biological activity of the two peptides, despite their differences in sequence and conformation.
Abstract: BcZBP is an LmbE-like, homohexameric, zinc-dependent deacetylase from the opportunistic pathogen Bacillus cereus with three, thus far uncharacterized, homologues in B. anthracis. Although its specific substrate is still unknown, the enzyme has been shown to preferentially deacetylate N-acetylglucosamine and diacetylchitobiose via an active site based on a zinc-binding motif of the type HXDDXnH. In the crystal structure, the active site is located at a deep and partially blocked cleft formed at the interface between monomers related by the molecular 3-fold axis, although the major, in structural terms, building block of the enzyme is not the trimer, but the intertwined dimer. Here, we report results from a 50 ns molecular dynamics simulation of BcZBP in explicit solvent with full electrostatics and show that (i) the view of the intertwined dimer as the major structural and functional building block of this class of hexameric enzymes is possibly an oversimplification of the rather complex dynamics observed in the simulation, (ii) the most mobile (with respect to their atomic fluctuations) parts of the structure coincide with three surface loops surrounding the active site, and (iii) these mobile loops define the active site’s accessibility, and may be implicated in the determination of the enzyme’s specificity.
Abstract: Arenicin-2 is a 21 residue antimicrobial cyclic peptide, possessing one disulphide bond between residues Cys(3) and Cys(20). NMR and CD studies suggested that the structure of arenicin-2 in water represented a well formed, but highly twisted beta-hairpin. To investigate the spatial arrangement of the peptide side chains and to get a clear view of its possible amphipathic properties we performed molecular dynamics in explicit water. Four independent trajectories, 50 ns in length, were produced, starting from various initial conformations or by applying different simulation conditions. Arenicin-2 retained its beta-hairpin structure during simulations, although the residues close to strand ends were found to escape from the ideal hairpin conformation. The type I' beta-turn connecting the two strands fluctuated between type IV and II' beta-turn. Conversely, the right-handed twist of the beta-hairpin was well conserved with average twist value 203 degrees +/- 19 degrees per eight residues. Several nonbonded interactions, like hydrophobic interactions between aliphatic side chains, cation/pi-aromatic interactions, CH...pi aromatic bond and water bridges, contributed to the hairpin stabilization.
Abstract: Disulfide cyclization is a well-known procedure to impose conformational restriction to peptides undergoing backbone flexibility. Rigid conformations are induced only for small rings with a specific combination of amino acids. In this work, we present a computational search of the backbone and backbone-dependent side-chain orientation of two series of linear and cyclic peptide analogs. The -C[XY]C- scaffold (where X,Y is arginine, aspartic acid or alanine residue) in its open and (S,S) cyclic form was used for the design of the studied analogs. Thirty-six compounds, resulting from the extension with one residue at either the N- or the C-terminus were studied with classical MD. The local backbone conformation and the relative orientation of the X and Y side chains induced by either cyclization and/or the presence of the charged residues are discussed. From the present study it is concluded that cyclization has a great impact on the synplanar orientation of the X and Y side chains in the (S,S)Ac-XCYC-NH(2) series of compounds while charge-charge interaction has only a weak synergic effect. On the contrary, the antiplanar orientation is favored in the case of (S,S)Ac-CXCY-NH(2). Copyright (c) 2008 European Peptide Society and John Wiley & Sons, Ltd.
Abstract: Recently, the synthesis of a molecule has been reported that belongs to a Lysine based, branched cyclic peptide class. This work explores the ability of such molecules to preserve the 3D geometry of the Trypsin's Active Site (TAS) by applying an integrated framework of automated computer procedures. The following four factors a) D/L chirality, b) different amino acids at different positions of the molecular scaffold's cyclic part, c) the application of AMBER and CHARMM force fields and d) different implicit solvation models were evaluated against TAS similarity. It was found that a number of molecules exhibit satisfactory geometric affinity to the TAS during extended Molecular Dynamics runs. We estimated that more than 2000 molecules (belonging to this class) could retain good similarity to the TAS arrangement.
Abstract: Molecular dynamics simulations have been performed for the 104–117 α-helical fragment of apolipoprotein A–I using the CHARMM22 force field and the Namd simulation engine. Simulation (50 ns in explicit water) resulted in significant appearance of π-helix conformation, which was totally diminished when the CMAP correction of the CHARMM force field was applied. This is consistent with other similar studies which suggest that the observation of π-helix in peptide conformation was force field biased rather actually existed. This study suggests that the 104–117 fragment of apolipoprotein A–I has a stable α-helical conformation in water.
Notes: Get reprints or see more at: http://stavrakoudis.econ.uoi.gr/stavrakoudis/publications.php
Abstract: We produce Monte Carlo evidence on the size and power of the RESET, a heteroscedasticity test, and a test for autocorrelation applied to realistic distributed-lag models. We find that the autocorrelation test has the correct size and high power to detect not only autocorrelation (given a correct model), but also the erroneous omission of several lags of an explanatory variable, whereas the RESET and heteroscedasticity tests are oversized in the presence of positive disturbance autocorrelation, especially when the regressors are also positively autocorrelated, and have no power to detect such misspecification errors. In large samples, size distortion may be avoided by using autocorrelation-robust methods.
Abstract: T-cell epitopes are important components of the inappropriate response of the immune system to self-proteins in autoimmune diseases. In this study, the candidate T-cell epitopes of the La/SSB autoantigen, the main target of the autoimmune response in patients with Sjogren's Syndrome (SS), and Systemic Lupus Erythematosus (SLE) were predicted using as a template the HLA-DQ2 and DQ7 molecules, which are genetically linked to patients with SS and SLE. Modeling of DQ2 and DQ7 was based on the crystal structure of HLA-DQ8, an HLA molecule of high risk factor of type I diabetes, which is also an autoimmune disease. The quality and reliability of the modeled DQ2 and DQ7 was confirmed by the Ramachandran plot and the TINKER molecular modeling software. Common and/or similar candidate T-cell epitopes, obtained by comparing three different approaches the Taylor's sequence pattern, the TEPITOPE quantitative matrices, and the MULTIPRED artificial neural network, were subjected to homology modeling with the crystal structure of the insulin-B peptide complexed with HLA-DQ8, and the best superposed candidate epitopes were placed into the modeled HLA-DQ2 and DQ7 binding grooves to perform energy minimization calculations. Six T-cell epitopes were predicted for HLA-DQ7 and nine for HLA-DQ2 covering parts of the amino-terminal and the central regions of the La/SSB autoantigen. Residues corresponding to the P1, P4, and P9 pockets of the HLA-DQ2 and DQ7 binding grooves experience very low SASA because they are less exposed to the microenvironment of the groove. The proposed T-cell epitopes complexed with HLA-DQ2/DQ7 were further evaluated for their binding efficiency according to their potential interaction energy, binding affinity, and IC50 values. Our approach constitutes the ground work for a rapid and reliable experimentation concerning the T-cell epitope mapping of autoantigens, and could lead to the development of T-cell inhibitors as immunotherapeutics in autoimmune diseases.
Abstract: Synthetic carriers play an important role in immunogen presentation, due to their ability of inducing improved and specific responses to conjugated epitopes. Their influence on the bioactive conformation of the epitope, though admittedly crucial for relevant in vitro and in vivo applications, is difficult to evaluate, given the usual lack of information on the complex conformational features determined by the nature of the carrier and the mode of ligation. Using the Herpes simplex virus glycoprotein D-1 epitope (Leu(9)-Lys-Nle-Ala-Asp-Pro-Asn-Arg-Phe-Arg-Gly-Lys-Asp-Leu(22)) as a model, we have performed a detailed conformational analysis on the free epitope peptide in solution and on three constructs in which the epitope was conjugated to sequential oligopeptide carriers {Ac-[Lys-Aib-Gly](4)-OH (SOC(4))} (through either a thioether or an amide bond; Ac: acetyl) and polytuftsin oligomers {H-[Thr-Lys-Pro-Lys-Gly](4)-NH(2) (T20)}, (through a thioether bond). The analysis of the epitope conformation in the parent protein, in carrier-conjugated and free form, suggests that the beta-turn structure of the -Asp(13)-Pro-Asn-Arg(16)- segment is highly conserved and independent of the epitope form. However, small conformational variations were observed at the C-terminal part of the epitope, depending on the nature of the carrier.
Abstract: We propose a new variant of RESET that is appropriate for distributed lag models. Monte Carlo evidence on size and power strongly supports the use of the new variant instead of the traditional RESET.
Abstract: The ability of an integrin to distinguish between the RGD-containing extracellular matrix proteins is thought to be due partially to the variety of RGD conformations. Three criteria have been proposed for the evaluation of the structure-activity relationship of RGD-containing peptides. These include: (i) the distance between the charged centres, (ii) the distance between the Arg Cbeta and Asp Cbeta atoms, and (iii) the pseudo-dihedral angle defining the Arg and Asp side-chain orientation formed by the Arg Czeta, Arg Calpha, Asp Calpha and Asp Cgamma atoms. A comparative conformation-activity study was performed between linear RGD peptides and strongly constrained cyclic (S,S) -CDC- bearing compounds, which cover a wide range of inhibition potency of platelet aggregation. It is concluded that the fulfilment of the -45 degrees < or = pseudo-dihedral angle < or = +45 degrees criterion is a prerequisite for an RGD compound to exhibit inhibitory activity. Once this criterion is accomplished, the longer the distance between the charged centres and/or between the Arg and Asp Cbeta atoms, the higher is the biological activity. In addition, the stronger the ionic interaction between Arg and Asp charged side chains, the lower the anti-aggregatory activity.
Abstract: A proline-rich [60]fullerene peptide was synthesized by use of (i) a 1,3-dipolar cycloaddition of an N-substituted glycine derivative to [60]fullerene, (ii) esterification of the isolated alcohol with the C-terminal aminoacid of the desired peptide sequence, and finally (iii) coupling of the remaining hexapeptide to give the final product 8 as a TFA salt, with oxidized methionine. Product 8 was found to be biologically active against sera from MCTD and SLE patients (ELISA experiment).
Abstract: In a previous work we described the design, synthesis and catalytic activity of a branched cyclic peptide as a serine protease mimic. To maximize its catalytic activity we present now a systematic search of a large number of homologous peptides for potential enzyme activity as revealed by the topological arrangement of the catalytic triad residues. This process is accomplished by applying a combined molecular mechanics and molecular dynamics conformational search of about 200 molecules. Starting from a previously synthesized compound that showed some hydrolytic activity several analogues were modelled by aminoacid substitutions in the main molecular framework using the Insight II molecular modelling environment with some script automation. Also presented is an algorithm that: (a) generates all possible combinations of residue substitutions, (b) scans the conformational space for each molecule via high temperature molecular dynamics, (c) picks the set of molecules the trajectories of which retained, to a considerable degree, the catalytic triad molecular arrangement, (d) subjects the selected molecules to layer solvation and energy minimization and chooses the molecules, the conformations of which could preserve the catalytic triad arrangement. Finally, a modelling with periodic boundary conditions, was performed to further support the reported algorithm. We found that at least one of the analogues could be a potential serine protease mimic, as revealed by the root-mean-square comparison between the catalytic triad in two molecular dynamics trajectories of the peptide and the corresponding residues in the crystal structure of trypsin. The most promising model candidate was synthesized and tested for its catalytic activity.
Abstract: According to general belief, the conformational information on short linear peptides in solution derived at ambient temperature from NMR spectrometry represents a population-weighted average over all members of an ensemble of rapidly interconverting conformations. Usually the search for discrete conformations is concentrated at low temperatures especially when sharp NMR resonances are detected at room temperature. Using the peptide Ac-RGD-NH(2) (Ac-Arg-Gly-Asp-NH(2), Ac: acetyl) as a model system and following a new approach, we have been able to demonstrate that short linear peptides can adopt discrete conformational states in DMSO-d(6) (DMSO: dimethylsulfoxide) which vary in a way critically dependent on the reconstitution conditions used before their dissolution in DMSO-d(6). The conformers are stabilized by intramolecular hydrogen bonds, which persist at high temperatures and undergo a very slow exchange with their extended structures in the NMR chemical shift time scale. The reported findings provide clear evidence for the occurrence of solvent-induced conformational exchange and point to DMSO as a valuable medium for folding studies of short linear peptides.
Abstract: Antagonists of fibrinogen at the GPIIb/IIIa receptor, which is the most abundant membrane protein on the platelet surface, are under active investigation as potential antithrombotics. The critical interaction between GPIIb/IIIa and fibrinogen can be inhibited by either linear or cyclic RGDS-containing peptides, which have been proved as lead compounds in the design of platelet aggregation inhibitors. In this study we present the design and construction of a new class of cyclic (S,S) non-RGD containing peptide sequences, using two Cys as a structural scaffold for the development of antiaggregatory agents. The (S,S)-CDC- sequence was incorporated as a conformational constraint, in molecules bearing at least one positive charge with the general formula (S,S)XCDCZ, where X = Ac-Arg, Pro-Arg, Pro-Ser-Lys, and Pro-Ser-Arg, and Z = -NH(2) and Arg-NH(2). Investigation of the structure-function relationships was performed on the basis of (a) the local conformation induced by the (S,S)-CDC motif, (b) the distance of the positively (R-C(zeta) or K-N(zeta)) and negatively (D-C(gamma)) charged centers, (c) the presence of a second positive or negative charge on the molecule, and (d) the orientation of the basic and acidic side chains defined by the pseudo dihedral angle (Pdo), which is formed by the R-C(zeta), R-C(alpha), D-C(alpha), and D-C(gamma) atoms in the case of (S,S)-RCDC and by the K-N(zeta), K-C(alpha), D-C(alpha), and D-C(gamma) atoms in the case of (S,S)-KCDC.
Abstract: The design, synthesis and catalytic properties of a cyclic branched peptide carrier that possesses the catalytic triad residues of the serine proteases is reported. The synthesis of the peptide model was totally completed on solid support using three different orthogonal amino protecting groups. Hydrolytic activity measurements against Suc-Ala-Ala-Ala-pNA substrate showed that it is hydrolysed by the peptide model to a small extent. Despite this small hydrolytic activity, it is the first time, to our knowledge, that hydrolysis of such a substrate is reported by an enzyme model compound. Contrary, this enzyme model peptide showed considerable activity against the Boc-Ala-pNP substrate (k_cat = 0.414 min^-1 and K_m = 0.228 mm). These results suggest that the designed carrier brings in appropriate contact the catalytic triad residues (Ser, His, Asp) resulting in the obtained hydrolytic activity.
Abstract: Humanized CAMPATH-1H antibody has been found to have biological applications through the recognition of the CD52 antigen. A peptide mimotope of the CD52 antigen with the sequence T1SSPSAD7 has been co-crystallized with the CAMPATH-1H antibody. A plethora of hydrogen bond interactions were found to mediate antigen recognition. An important feature of peptide’s bound conformation was the type I β-turn found in the S3PSA6 peptide’s fragment. Paradoxically, this fact has been underestimated from other researchers. In order to further investigate the importance of this structural feature and its significance in antibody/antigen binding we have performed molecular dynamics simulations in explicit water of the T1SSPSAD7 peptide in both antibody free and bound states. We have found that the turn structure has been perfectly retained in the bound state but it was eliminated in the free state. This fact implies that the turn structure of the peptide is unstable in aqueous environment and it is induced upon antibody binding. Analysis of the trajectories revealed also several other important features of the antibody/antigen binding mode.
