Comparison of the solution conformations of a cell-adhesive peptide LBE and its reverse sequence EBL

T-cell adhesion is mediated by an ICAM-1/LFA-1 interaction; this interaction plays a crucial role in T-cell activation during immune response. LBE peptide, which is derived from the beta-subunit of LFA-1, has been shown to inhibit ICAM-1/LFA-1-mediated T-cell adhesion. In this work, we studied the solution conformations of LBE peptide and its reverse sequence (EBL) by NMR, CD and molecular dynamics simulations. Reverse peptides have been used as controls in biological studies. The effect of reversing the sequence of LBE to EBL peptides on their respective conformations is important in understanding their biological properties in vitro or in vivo. The NMR studies for these peptides were carried out in water and in TFE/water solvent systems. In 40% TFE/water, both peptides exhibited helical conformation. CD studies suggested that the LBE exhibits 30% helical conformation, while the EBL exhibits 20% helical conformation. From the NMR and MD simulation studies, it was evident that the peptides exhibited a stable helical conformation; a stable helical structure was found at Leu6 to Leu15 for LBE and at Gly9 to Leu17 for EBL. The helical conformations of LBE and EBL may be in equilibrium with other possible conformers; the other conformers contain loop and turn structures. Both peptides bind to divalent cations because the LBE is derived from the cation-binding region of the LFA-1. This study shows that reversing the peptide sequence did not alter the secondary structure of the corresponding sequence. Hence, caution must be exercised when using reverse peptides as controls in biological studies. This report will improve our ability to design a better inhibitor of ICAM-1/LFA-1 interaction.     

Enhanced sampling of the molecular potential energy surface using mutually orthogonal latin squares: application to peptide structures

The computational identification of the optimal three-dimensional fold of even a small peptide chain from its sequence, without reference to other known structures, is a complex problem. There have been several attempts at solving this by sampling the potential energy surface of the molecule in a systematic manner. Here we present a new method to carry out the sampling, and to identify low energy conformers of the molecule. The method uses mutually orthogonal Latin squares to select (of the order of) n(2) points from the multidimensional conformation space of size m(n), where n is the number of dimensions (i.e., the number of conformational variables), and m specifies the fineness of the search grid. The sampling is accomplished by first calculating the value of the potential energy function at each one of the selected points. This is followed by analysis of these values of the potential energy to obtain the optimal value for each of the n-variables separately. We show that the set of the n-optimal values obtained in this manner specifies a low energy conformation of the molecule. Repeated application of the method identifies other low energy structures. The computational complexity of this algorithm scales as the fourth power of the size of the molecule. We applied this method to several small peptides, such as the neuropeptide enkephalin, and could identify a set of low energy conformations for each. Many of the structures identified by this method have also been previously identified and characterized by experiment and theory. We also compared the best structures obtained for the tripeptide (Ala)(3) by the present method, with those obtained by an exhaustive grid search, and showed that the algorithm is successful in identifying all the low energy conformers of this molecule.     

Chrysophanol anthrone: a theoretical study on the potential energy surface

The potential energy surface (PES) of chrysophanol anthrone, the active component of Goa Powder, was systematically explored and thoroughly scrutinised via density functional theory, in order to gain an understanding of its physicochemical properties. In particular, we focused on the rotations of the two hydroxyl-phenyl dihedral angles. A picture with a four stable rotamers emerged where only the most stable conformer has a planar structure and the less stable conformer has the maximum deviation from planarity. The computed PES shows that the energy barriers for the conformer interconversion are less than 15 kcal/mol. From the analysis of the calculated intramolecular hydrogen bond enthalpy, we conclude that the number of the intramolecular hydrogen bonds governs the conformer stability. Additionally, the conformational equilibrium was pursued by means of an analysis of the energy of OH internal rotation barriers. The total energy changes were decomposed in an electrostatic decomposition scheme in order to gain an insight into the effects governing the torsional barrier and preferred conformations. This analysis shows that the interplay between the repulsive and attractive potentials causes the conformer stability, where the attractive term dominates the conformer stabilisation.     

The molecular mechanics program DUPLEX: Computing structures of carcinogen modified DNA by surveying the potential energy surface

The nucleic acids molecular mechanics program DUPLEX has been designed with useful features for surveying the potential energy surface of polynucleotides, especially ones that are modified by polycyclic aromatic carcinogens. The program features helpful strategies for addressing the multiple minimum problem: (1) the reduced variable domain of torsion angle space; (2) search strategies that emphasize large scale searches for smaller subunits, followed by building to larger units by a variety of strategies; (3) the use of penalty functions to aid the minimizer in locating selected structural types in first stage minimizations; penalty functions are released in terminal minimizations to yield final unrestrained minimum energy conformations. Predictive capability is illustrated by DNA modified by activated benzo[a]pyrenes.     

Using a low denaturant model to explore the conformational features of translocation-active SecA

The SecA molecular nanomachine in bacteria uses energy from ATP hydrolysis to drive post-translational secretion of preproteins through the SecYEG translocon. Cytosolic SecA exists in a dimeric, "closed" state with relatively low ATPase activity. After binding to the translocon, SecA undergoes major conformational rearrangement, leading to a state that is structurally more "open", has elevated ATPase activity, and is active in translocation. The structural details underlying this conformational change in SecA remain incompletely defined. Most SecA crystal structures report on the cytosolic form; only one structure sheds light on a form of SecA that has engaged the translocon. We have used mild destabilization of SecA to trigger conformational changes that mimic those in translocation-active SecA and thus study its structural changes in a simplified, soluble system. Results from circular dichroism, tryptophan fluorescence, and limited proteolysis demonstrate that the SecA conformational reorganization involves disruption of several domain-domain interfaces, partial unfolding of the second nucleotide binding fold (NBF) II, partial dissociation of the helical scaffold domain (HSD) from NBF I and II, and restructuring of the 30 kDa C-terminal region. These changes account for the observed high translocation SecA ATPase activity because they lead to the release of an inhibitory C-terminal segment (called intramolecular regulator of ATPase 1, or IRA1) and of constraints on NBF II (or IRA2) that allow it to stimulate ATPase activity. The observed conformational changes thus position SecA for productive interaction with the SecYEG translocon and for transfer of segments of its passenger protein across the translocon.     

Enthalpy measurement using calorimetry shows a significant difference in potential energy between the active and latent conformations of PAI-1

A central feature of the serpin inhibition mechanism is insertion of the reactive center loop into the central beta-sheet (beta-sheet A). This insertion also occurs when the reactive center loop is cleaved without protease inhibition. Using this effect, we have measured the enthalpy (DeltaH) of loop cleavage and insertion for plasminogen activator inhibitor 1 (PAI-1) as -38 kcal/mol. Because loop insertion can be blocked by incorporating a peptide into the central beta-sheet, it was possible to assign -7 kcal/mol to loop cleavage and -31 kcal/mol to loop insertion. These values are lower than values reported for the serpins alpha 1 -proteinase inhibitor and antithrombin of -53 to -63 kcal/mol, respectively, for loop insertion with negligible enthalpy for loop cleavage. A free energy difference of -9 kcal/mol has been reported between the active and spontaneously loop inserted "latent forms" of PAI-1, which is significantly smaller in magnitude than the -31 kcal/mol of enthalpy we measured for loop insertion. Because the enthalpy should relate closely to those regions of PAI-1 that have moved to lower potential energy, a difference distance matrix is presented that identifies regions of PAI-1 that move during loop insertion.     

Fluorescent inhibitors of a cell surface protease used to locate leukaemia cells in kidney sections

Guanidinobenzoatase is a trypsin-like protease on the surface of cells capable of migration, for example leukaemia cells. We have used a number of fluorescent probes that are competitive inhibitors of guanidinobenzoatase to locate leukaemia cells in resin sections of kidney tissue obtained from leukaemic rats. We have demonstrated how this competitive inhibition system can be used to direct desired molecules (such as cytotoxic drugs) to these cells and to monitor the arrival of such compounds at the active site of guanidinobenzoatase. The principles developed in this study could equally well be applied to other enzymes on other cells provided suitable competitive inhibitors were designed. The presence of an enzyme on the surface of a cell can be used to direct molecules to that cell provided that these molecules contain a functional group that acts as an inhibitor for the chosen enzyme.     

Short peptide amyloid organization: stabilities and conformations of the islet amyloid peptide NFGAIL

Experimentally, short peptides have been shown to form amyloids similar to those of their parent proteins. Consequently, they present useful systems for studies of amyloid conformation. Here we simulate extensively the NFGAIL peptide, derived from the human islet amyloid polypeptide (residues 22-27). We simulate different possible strand/sheet organizations, from dimers to nonamers. Our simulations indicate that the most stable conformation is an antiparallel strand orientation within the sheets and parallel between sheets. Consistent with the alanine mutagenesis, we find that the driving force is the hydrophobic effect. Whereas the NFGAIL forms stable oligomers, the NAGAIL oligomer is unstable, and disintegrates very quickly after the beginning of the simulation. The simulations further identify a minimal seed size. Combined with our previous simulations of the prion-derived AGAAAAGA peptide, AAAAAAAA, and the Alzheimer Abeta fragments 16-22, 24-36, 16-35, and 10-35, and the solid-state NMR data for Abeta fragments 16-22, 10-35, and 1-40, some insight into the length and the sequence matching effects may be obtained.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

Solid state and solution conformations of a helical peptide with a central gly-gly segment

The influence of amino acids with contrasting conformational tendencies on the stereochemistry of oligopeptides has been investigated using an octapeptide Boc-Leu-Aib-Val-Gly-Gly-Leu-Aib-Val-OMe, which contains two helix-promoting Aib residues and a central helix-destabilizing Gly-Gly segment. Single crystal x-ray diffraction studies reveal that a 3 ₁₀-helix is formed up to the penultimate Aib residue, at which point there is a helix reversal in the backbone, reminiscent of a C-terminal 6 → I hydrogen bond. The curious feature in the crystal is the solvation of the possible 6 → 1 bond by a CH₃OH molecule, where the OH is inserted between O(3) and N(8) and participates in hydrogen bonds with both. The cell parameters are as follows: space group P2₁2₁2₁, a = 10.649(4) Å, b = 15.694(5) Å, c = 30.181(8) Å, R = 6.7% for 3427 data (| F₀| > 3σF) observed to 0.9 Å. Nuclear magnetic resonance studies in CDCl₃ using NH group solvent accessibility and nuclear Overhauser effects as probes are consistent with a 3 ₁₀-helical conformation. In contrast, in (CD₃)₂SO, unfolding of the central segment results in a multiple β-turn structure, with β-turn conformations populated at residues 1–2, 3–4, and 6–7. CD studies in methanol-2,2,2-trifluoroethanol (TFE) mixtures also provide evidence for a solvent-dependent structural transition. Helical conformations are populated in TFE, while type II β-turn structures are favored in methanol. © 1996 John Wiley & Sons, Inc.     

Crystal structure, conformation, and potential energy calculations of the chemotactic peptide N-formyl-L-Met-L-Leu-L-Phe-OMe

The tripeptide N-formyl-L-Met-L-Leu-L-Phe-OMe (FMLP-OMe) crystallizes in the orthorhombic system, space group P 2(1)2(1)2(1), with the following unit-cell parameters: a = 21.727, b = 21.836, c = 5.133 A, Z = 4. The structure has been solved and refined to a final R of 0.068 for 1838 independent reflexions with I greater than 2 omega (I). The peptide backbone is folded at the Leu residue (phi L = -67.7, psi L = -49.1 degrees) without intramolecular hydrogen bonds. Considering each peptide plane, the Leu side-chain is oriented on the same side of that of the Phe residue and on the opposite side of that of the Met residue, respectively. The crystal conformation differs from all the other conformations proposed for FMLP-OMe and the anionic form of N-formyl-L-Met-L-Leu-L-Phe-OH (FMLP) in solution accounts for the amphiphilic character of the peptide, giving rise, through intermolecular hydrogen bonds, to a stacking of molecules which could be maintained in the aggregation states experimentally observed in solvents of low polarity. Intramolecular potential energy calculations have been carried out in order to compare the energies of the various backbone conformers.     

Comparison of the low energy conformations of an oncogenic and a non-oncogenic p21 protein, neither of which binds GTP or GDP

Oncogenic p21 protein, encoded by theras-oncogene, that causes malignant transformation of normal cells and many human tumors, is almost identical in sequence to its normal protooncogene-encoded counterpart protein, except for the substitution of arbitrary amino acids for the normally occurring amino acids at critical positions such as Gly 12 and Gin 61. Since p21 is normally activated by the binding of GTP in place of GDP, it has been postulated that oncogenic forms must retain bound GTP for prolonged time periods. However, two multiply substituted p21 proteins have been cloned, neither of which binds GDP or GTP. One of these mutant proteins with Val for Gly 10, Arg for Gly 12, and Thr for Ala 59 causes cell transformation, while the other, similar protein with Gly 10, Arg 12, Val for Gly 13 and Thr 59 does not transform cells. To define the critical conformational changes that occur in the p21 protein that cause it to become oncogenic, we have calculated the low energy conformations of the two multiply substituted mutant p21 proteins using a new adaptation of the electrostatically driven Monte Carlo (EDMC) technique, based on the program ECEPP. We have used this method to explore the conformational space available to both proteins and to compute the average structures for both using statistical mechanical averaging. Comparison of the average structures allows us to detect the major differences in conformation between the two proteins. Starting structures for each protein were calculated using the recently deposited x-ray crystal coordinates for the p21 protein, that was energy-refined using ECEPP, and then perturbed using the EDMC method to compute its average structure. The specific amino acid substitutions for both proteins were then generated into the lowest energy structure generated by this procedure, subjected to energy minimization and then to full EDMC perturbations. We find that both mutant proteins exhibit major differences in conformation in specific regions, viz., residues 35–47, 55–78, 81–93, 96–110, 115–126, and 123–134, compared with the EDMC-refined x-ray structure of the wild-type protein. These regions have been found to be the most flexible in the p21 protein bound to GDP from prior molecular dynamics calculations (Dykeset al., 1993). Comparison of the EDMC-average structure of the transforming mutant with that of the nontransforming mutant reveals major structural differences at residues 10–16, 32–40, and 60–68. These structural differences appear to be the ones that are critical in activation of the p21 protein. Analysis of the correlated motions of the different regions of the two mutant proteins reveals that changes in the conformation of regions in the carboxyl half of the protein are caused by changes in conformation around residues 10–16 and are transmitted by means of residues around Gln 61. The latter region therefore constitutes a molecular switch unit, in agreement with conclusions from prior work.On leave from the Department of Chemistry, University of Gdask, ul. Sobieskiego 18, 80-952 Gdask, Poland.     

Mass spectrometry to characterize the binding of a peptide to a lipid surface.

The binding of an amphipathic alpha-helical peptide to small unilamellar lipid vesicles has been examined using chemical derivitization and mass spectrometry. The peptide is derived from the sequence of human apolipoprotein C-II (apoC-II), the protein activator of lipoprotein lipase (LpL). ApoC-II(19-39) forms approximately 60% alpha-helix upon binding to model egg yolk phosphatidylcholine small unilamellar vesicles. Measurement of the affinity of the peptide for lipid by spectrophotometric methods is complicated by the contribution of scattered light to optical signals. Instead, we characterize the binding event using the differential labeling of lysine residues by the lipid- and aqueous-phase cross-linkers, disuccinimidyl suberate (DSS) and bis(sulfosuccinimidyl) suberate (BS(3)), respectively. In aqueous solution, the three lysine residues of the peptide are accessible to both cross-linkers. In the presence of lipid, the C-terminal lysine residue becomes inaccessible to the lipid-phase cross-linker DSS, but remains accessible to the aqueous-phase cross-linker, BS(3). We use mass spectrometry to characterize this binding event and to derive a dissociation constant for the interaction (K(d) = 5 microM). We also provide evidence for the formation of dimeric cross-linked peptide when high densities of peptide are bound to the lipid surface.     

Structural characterization of the (MeSH)4 potential energy surface

A random walk on the PES for (MeSH)₄ clusters produced 50 structural isomers held together by hydrogen-bonding networks according to calculations performed at the B3LYP/6–311++G** and MP2/6–311++G** levels. The geometric motifs observed are somewhat similar to those encountered for the methanol tetramer, but the interactions responsible for cluster stabilization are quite different in origin. Cluster stabilization is not related to the number of hydrogen bonds. Two distinct, well-defined types of hydrogen bonds scattered over a wide range of distances are predicted.

MARS: a program to find potential restriction sites

The figure shown below should have accompanied the paper detailedhere.     

Sorting nexins 1 and 2a locate mainly to the TGN

The subcellular localization of the sorting nexins (SNXs) in higher plants is a matter of controversy. Previous confocal laser scanning microscopy (CLSM studies on root cells from a transgenic Arabidopsis line expressing SNX1-GFP have suggested that this SNX is present on an endosome having characteristics of both the trans-Golgi network (TGN) and the multivesicular body (MVB). In contrast, SNX2a locates exclusively to the TGN when transiently expressed in tobacco mesophyll protoplasts. By performing immunogold electron microscopy on cryofixed Arabidopsis roots, we have tried to clarify the situation. Both SNX1-GFP and endogenous SNX2a locate principally to the TGN. Labeling of MVBs could not be confirmed with any certainty.     

Rearrangement of partially ordered stacked conformations contributes to the rugged energy landscape of a small RNA hairpin

We have studied the fast relaxation kinetics of a small RNA hairpin tetraloop using time-resolved infrared spectroscopy. A laser-induced temperature jump initiated the relaxation by rapidly perturbing the thermal equilibrium of the sample. We probed the relaxation kinetics at two different wavenumbers, 1574 and 1669 cm (-1). The latter is due to the C6O6 carbonyl stretch of the base guanine and is a direct measure of guanine base pairing. The former is assigned to a ring vibration of guanine and tracks structure by sensing base stacking interactions. Overall, the kinetics at 1574 cm (-1) are faster than those observed at 1669 cm (-1). When relaxation occurs at the melting temperature, the kinetics at both wavenumbers are biexponential. When relaxation occurs at a temperature that is higher than the melting temperature, the data at 1669 cm (-1) are still biexponential while only a single fast phase is resolved in the data at 1574 cm (-1). The fast phases are in the range of microseconds, while the slower phases are in the range of tens of microseconds. At both wavenumbers, a portion of the relaxation is not resolved, indicating the existence of a very fast, sub-100 ns phase. Our results provide additional evidence that small, fast folding hairpin loops are characterized by a rugged energy landscape. Furthermore, our data suggest that single-strand stacking interactions and stacking interactions in the loop contribute significantly to the ruggedness of the energy landscape. This work also demonstrates the utility of time-resolved infrared spectroscopy in studying RNA folding.