Role of spin state on the geometry and nuclear quadrupole resonance parameters in hemin complex

Theoretical calculations of structural parameters, 57Fe, 14N and 17 O electric field gradient (EFG) tensors for full size-hemin group have been carried out using density functional theory. These calculations are intended to shed light on the difference between the geometry parameters, nuclear quadrupole coupling constants (QCC), and asymmetry parameters (eta Q) found in three spin states of hemin; doublet, quartet and sextet. The optimization results reveal a significant change for propionic groups and porphyrin plane in different spin states. It is found that all principal components of EFG tensor at the iron site are sensitive to electronic and geometry structures. A relationship between the EFG tensor at the 14N and 17 O sites and the spin state of hemin complex is also detected.     

Theoretical investigation of hydrogen bonding effects on oxygen, nitrogen, and hydrogen chemical shielding and electric field gradient tensors of chitosan/HI salt

A density functional theory study has been carried out to calculate the (17)O, (15)N, (13)C, and (1)H chemical shielding as well as (17)O, (14)N, and (2)H electric field gradient tensors of chitosan/HI type I salt. These calculations were performed using the B3LYP functional and 6-311++G (d,p) and 6-31++G (d,p) basis sets. Calculated EFG and chemical shielding tensors were used to evaluate the (17)O, (14)N, and (2)H nuclear quadruple resonance, NQR, and (17)O, (15)N, (13)C, and (1)H nuclear magnetic resonance, NMR, parameters in the cluster model, which are in good agreement with the available experimental data. The difference in the isotropic shielding (sigma(iso)) and quadrupole coupling constant (C(Q)) between monomer and target molecule in the cluster was analyzed in detail. It was shown that both EFG and CS tensors are sensitive to hydrogen-bonding interactions, and calculating both tensors is an advantage. A different influence of various hydrogen bond types, N-Hcdots, three dots, centeredI, O-Hcdots, three dots, centeredI, and N-Hcdots, three dots, centeredO was observed on the calculated CS and EFG tensors. On the basis of this study, nitrogen and O-6 are the most important nuclei to confirm crystalline structure of chitosan/HI. These nuclei have large change in their CS and EFG tensors because of forming intermolecular hydrogen bonds. Moreover, the quantum chemical calculations indicated that the intermolecular hydrogen-bonding interactions play an essential role in determining the relative orientation of CS and EFG tensors of O-6 and nitrogen atoms in the molecular frame axes.     

A density functional study of (17)O, (14)N and (2)H electric field gradient tensors in the real crystalline structure of alpha-glycine

A density functional theory (DFT) study was carried out to calculate (17)O, (14)N and (2)H electric field gradient (EFG) tensors in accurate neutron diffraction structures of alpha-glycine at 288 and 427 K. B3LYP is the used method and 6-311+G(*) and 6-311++G(**) are the basis sets in the calculations of EFG tensors at the sites of (17)O, (14)N and (2)H nuclei in the monomer and the octameric cluster of alpha-glycine at two temperatures. Quadrupole coupling constants and asymmetry parameters are the converted parameters of calculated EFG tensors to experimentally measurable ones. The calculated results of monomer and the target molecule in octameric cluster reveal that hydrogen-bonding interactions play an important role in the crystalline structure of alpha-glycine where the results of the target molecule in octameric cluster are in good agreement with the experiments.     

Investigation of C-H...O=C and N-H...OC hydrogen-bonding interactions in crystalline thymine by DFT calculations of O-17, N-14 and H-2 NQR parameters

A computational study at the level of density functional theory (DFT) was carried out to investigate C-H...O=C and N-H...O=C hydrogen-bonding interactions (HBs) in the real crystalline cluster of thymine by O-17, N-14 and H-2 calculated nuclear quadrupole resonance (NQR) parameters. To perform the calculations, a hydrogen-bonded pentameric cluster of thymine was created using X-ray coordinates where the hydrogen atoms positions are optimized and the electric field gradient (EFG) tensors were calculated for the target molecule. Additional EFG calculations were also performed for crystalline monomer and an optimized isolated gas-phase thymine. The calculated EFG tensors at the level of B3LYP and B3PW91 DFT methods and 6-311++G**and CC-pVTZ basis sets were converted to those experimentally measurable NQR parameters, quadrupole coupling constants and asymmetry parameters. The results reveal that because of strong contribution to N-H...O=C HBs, NQR parameters of O2, N1 and N3 undergo significant changes from monomer to the target molecule in cluster. Furthermore, the NQR parameters of O2 also undergo some changes because of non-classical C-H...O=C HBs.     

<Emphasis Type="Italic">Ab initio</Emphasis> DFT study of bisphosphonate derivatives as a drug for inhibition of cancer: NMR and NQR parameters

DFT computations were carried out to characterize the ¹⁷Oand ²H electric field gradient, EFG, in various bisphosphonate derivatives. The computations were performed at the B3LYP level with 6-311++G (d,P) standard basis set. Calculated EFG tensors were used to determine the ¹⁷O and ²H nuclear quadrupole coupling constant, χ and asymmetry parameter, η. For better understanding of the bonding and electronic structure of bisphosphonates, isotropic and anisotropic NMR chemical shieldings were calculated for the ¹³C, ¹⁷O and ³¹P nuclei using GIAO method for the optimized structure of intermediate bisphosphonates at B3LYP level of theory using 6-311++G (d, p) basis set. The results showed that various substituents have a strong effect on the nuclear quadrupole resonance (NQR) parameters (χ, η) of ¹⁷O in contrast with ²H NQR parameters. The NMR and NQR parameters were studied in order to find the correlation between electronic structure and the activity of the desired bisphosphonates. In addition, the effect of substitutions on the bisphosphonates polarity was investigated. Molecular polarity was determined via the DFT calculated dipole moment vectors and the results showed that substitution of bromine atom on the ring would increase the activity of bisphosphonates.     

Investigation of intermolecular hydrogen bond interactions in crystalline l-Cysteine by DFT calculations of the oxygen-17, nitrogen-14, and hydrogen-2 EFG tensors and AIM analysis

A systematic computational study is carried out to investigate hydrogen bond (HB) interactions in the real crystalline structures of l-Cysteine at 30 and 298 K by density functional theory (DFT) calculations of electric field gradient (EFG) tensors at the sites of O-17, N-14, and H-2 nuclei. One-molecule (monomer) and nine-molecule (cluster) models of l-Cysteine are created by available crystal coordinates at both temperatures and the EFG tensors are calculated for both models to indicate the effect of HB interactions on the tensors. The calculated EFG tensors at the level of B3LYP and B3PW91 DFT methods and 6-311++G?? and cc-pVTZ basis sets are converted to those experimentally measurable nuclear quadrupole resonance (NQR) parameters i.e. quadrupole coupling constants (qcc) and asymmetry parameters (η_Q). The evaluated NQR parameters reveal that the EFG tensors of ¹⁷O, ¹⁴N, and ²H are influenced and show particular trends from monomer to the target molecule in the cluster due to the contribution of target molecule to classic N–H…O, and non-classic S–H…O and S–H…S types of HB interactions. On the other hand, atoms in molecules (AIM) analyses confirm the presence of HB interactions and rationalize the observed EFG trends. The results indicate different contribution of various nuclei to HB interactions in the cluster where O2 and N1 have major contributions. The EFG tensors as well as AIM analysis at the H6 site show that the N1-H6…O2 HB undergoes a significant change from 30 to 298 K where changes in other N–H…O interactions are almost negligible. There is a good agreement between the calculated ¹⁴N NQR parameters and reported experimental data.     

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.     

All-atom molecular dynamics simulations using orientational constraints from anisotropic NMR samples

Orientational constraints obtained from solid state NMR experiments on anisotropic samples are used here in molecular dynamics (MD) simulations for determining the structure and dynamics of several different membrane-bound molecules. The new MD technique is based on the inclusion of orientation dependent pseudo-forces in the COSMOS-NMR force field. These forces drive molecular rotations and re-orientations in the simulation, such that the motional time-averages of the tensorial NMR properties approach the experimentally measured parameters. The orientational-constraint-driven MD simulations are universally applicable to all NMR interaction tensors, such as chemical shifts, dipolar couplings and quadrupolar interactions. The strategy does not depend on the initial choice of coordinates, and is in principle suitable for any flexible molecule. To test the method on three systems of increasing complexity, we used as constraints some deuterium quadrupolar couplings from the literature on pyrene, cholesterol and an antimicrobial peptide embedded in oriented lipid bilayers. The MD simulations were able to reproduce the NMR parameters within experimental error. The alignment of the three membrane-bound molecules and some aspects of their conformation were thus derived from the NMR data, in good agreement with previous analyses. Furthermore, the new approach yielded for the first time the distribution of segmental orientations with respect to the membrane and the order parameter tensors of all three systems.     

Influence of N-H...O and O-H...O hydrogen bonds on the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen: a density functional theory study

A computational investigation was carried out to characterize the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen. We found that N-H...O and O-H...O hydrogen bonds around the acetaminophen molecule in the crystal lattice have different influences on the calculated (17)O, (15)N and (13)C chemical shielding eigenvalues and their orientations in the molecular frame of axes. The calculations were performed with the B3LYP method and 6-311++G(d, p) and 6-311+G(d) standard basis sets using the Gaussian 98 suite of programs. Calculated chemical shielding tensors were used to evaluate the (17)O, (15)N, and (13)C NMR chemical shift tensors in crystalline acetaminophen, which are in reasonable agreement with available experimental data. The difference between the calculated NMR parameters of the monomer and molecular clusters shows how much hydrogen-bonding interactions affect the chemical shielding tensors of each nucleus. The computed (17)O chemical shielding tensor on O(1), which is involved in two intermolecular hydrogen bonds, shows remarkable sensitivity toward the choice of the cluster model, whereas the (17)O chemical shielding tensor on O(2) involved in one N-H...O hydrogen bond, shows smaller improvement toward the hydrogen-bonding interactions. Also, a reasonably good agreement between the experimentally obtained solid-state (15)N and (13)C NMR chemical shifts and B3LYP/6-311++G(d, p) calculations is achievable only in molecular cluster model where a complete hydrogen-bonding network is considered. Moreover, at the B3LYP/6-311++G(d, p) level of theory, the calculated (17)O, (15)N and (13)C chemical shielding tensor orientations are able to reproduce the experimental values to a reasonably good degree of accuracy.     

Equilibrium unfolding of the poly(glutamic acid)20 helix

The equilibrium structural ensemble of a 20-residue polyglutamic acid peptide (E(20)) was studied with FRET, circular dichroism, and molecular dynamics (MD) simulations. A FRET donor, o-aminobenzamide, and acceptor, 3-nitrotyrosine, were introduced at the N- and C-termini, respectively. Circular dichroism, steady state FRET, and time-resolved FRET measurements were employed to characterize the fraction helix and end-to-end distance under different pH conditions: pH 4 (60% alpha-helix), pH 6 (0% alpha-helix), and pH 9 (0% alpha-helix). At pH 4, the end-to-end distance was measured at 24 A and determined to be considerably less than the 31 A predicted for an alpha-helix of the same length. At pH 6 and 9, the end-to-end distance was measured at > 31 and 39 A respectively, both which are determined to be considerably greater than the 27 A predicted for a freely jointed random coil of the same length. To better understand the physical forces underlying the unusual helix-coil transition in this peptide, three theoretical MD models of E(20) were constructed: (1) a pure alpha-helix, (2) an alpha-helix with equivalent attractive intramolecular contacts, and (3) a weak alpha-helix with termini-weighted intramolecular contacts ("sticky ends"). Using MD simulations, the bent helix structure calculated from Model 3 was found to be the closest in agreement with the experimental data.     

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.     

Dynamics of Ca2+-saturated calmodulin D129N mutant studied by multiple molecular dynamics simulations

Fifteen independent 1-nsec MD simulations of fully solvated Ca(2+) saturated calmodulin (CaM) mutant D129N were performed from different initial conditions to provide a sufficient statistical basis to gauge the significance of observed dynamical properties. In all MD simulations the four Ca(2+) ions remained in their binding sites, and retained a single water ligand as observed in the crystal structure. The coordination of Ca(2+) ions in EF-hands I, II, and III was sevenfold. In EF-hand IV, which was perturbed by the mutation of a highly conserved Asp129, an anomalous eightfold Ca(2+) coordination was observed. The Ca(2+) binding loop in EF-hand II was observed to dynamically sample conformations related to the Ca(2+)-free form. Repeated MD simulations implicate two well-defined conformations of Ca(2+) binding loop II, whereas similar effect was not observed for loops I, III, and IV. In 8 out of 15 MD simulations Ca(2+) binding loop II adopted an alternative conformation in which the Thr62 >C=O group was displaced from the Ca(2+) coordination by a water molecule, resulting in the Ca(2+) ion ligated by two water molecules. The alternative conformation of the Ca(2+) binding loop II appears related to the "closed" state involved in conformational exchange previously detected by NMR in the N-terminal domain fragment of CaM and the C-terminal domain fragment of the mutant E140Q. MD simulations suggest that conformations involved in microsecond exchange exist partially preformed on the nanosecond time scale.     

Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.     

Molecular dynamics study of kaliotoxin in water

Kaliotoxin (KTX), a potassium channel blocker found in the venom of the scorpion Androctonous Mauretanicus is a 38 residue polypeptide with a well defined structure consisting of a alpha-helix and a three strand antiparallel beta-sheet interconnected by three disulfide bonds. Although the 3D structure has been determined by NMR, there is a number of features, mainly concerning the conformation and flexibility of the side chains, but also the long range order in the peptide and its fluctuations, that may have escaped the experimental study. These questions are now being addressed using molecular dynamics (MD) simulations. Accordingly, the present work reports the analysis of a 430 ps molecular dynamics trajectory of the polypeptide soaked with 4700 TIP3 water molecules inside a 56 A box. MD calculations were performed with periodic boundary conditions. Analysis of the conformational space sampled by each of the residues along the trajectory, suggests a special behavior of Pro17 and Lys19 both located on the helix. Furthermore, analysis of the relative movements of the secondary structure elements indicates that the alpha-helix and beta-sheets fluctuate in a correlated motion, preserving the tertiary structure of the polypeptide along the trajectory. Finally, analysis of the charge distribution was also examined. The direction of the dipole moment, computed from the center of masses appears to be an interesting feature of the structure probably related to the biological function of the molecule.     

Conformation and aggregation of M13 coat protein studied by molecular dynamics.

Molecular dynamics (MD) simulations are performed on M13 coat protein, a small membrane protein for which both alpha- and beta-structures have been suggested. The simulations are started from initial conformations that are either monomers or dimers of alpha-helices or U-shaped beta-sheets. The lipid bilayer is represented by a hydrophobic potential. The results are analyzed in terms of stability, energy and secondary structure. The U-shaped beta-structure changes from a planar to a twisted form with larger twist for the monomer than the dimer. The beta-sheet is much more flexible than the alpha-helix as monitored by the root mean square (rms) fluctuations of the C alpha atoms. A comparison of the energies after 100 ps MD simulation shows that of the monomers, the alpha-helix has the lowest energy. The energy difference between alpha- and beta-structures decreases from 266 kJ/mol to 148 kJ/mol, when going from monomers to dimers. It is expected that this difference will decrease with higher aggregation numbers.     

Understanding interactions of gastric inhibitory polypeptide (GIP) with its G-protein coupled receptor through NMR and molecular modeling.

Gastric inhibitory polypeptide (GIP, or glucose-dependent insulinotropic polypeptide) is a 42-amino acid incretin hormone moderating glucose-induced insulin secretion. Antidiabetic therapy based on GIP holds great promise because of the fact that its insulinotropic action is highly dependent on the level of glucose, overcoming the sideeffects of hypoglycemia associated with the current therapy of Type 2 diabetes. The truncated peptide, GIP(1-30)NH2, has the same activity as the full length native peptide. We have studied the structure of GIP(1-30)NH2 and built a model of its G-protein coupled receptor (GPCR). The structure of GIP(1-30)NH2 in DMSO-d6 and H2O has been studied using 2D NMR (total correlation spectroscopy (TOCSY), nuclear overhauser effect spectroscopy (NOESY), double quantum filtered-COSY (DQF-COSY), 13C-heteronuclear single quantum correlation (HSQC) experiments, and its conformation built by MD simulations with the NMR data as constraints. The peptide in DMSO-d6 exhibits an alpha-helix between residues Ile12 and Lys30 with a discontinuity at residues Gln19 and Gln20. In H2O, the alpha-helix starts at Ile7, breaks off at Gln19, and then continues right through to Lys30. GIP(1-30)NH2 has all the structural features of peptides belonging to family B1 GPCRs, which are characterized by a coil at the N-terminal and a long C-terminal alpha-helix with or without a break. A model of the seven transmembrane (TM) helices of the GIP receptor (GIPR) has been built on the principles of comparative protein modeling, using the crystal structure of bovine rhodopsin as a template. The N-terminal domain of GIPR has been constructed from the NMR structure of the N-terminal of corticoptropin releasing factor receptor (CRFR), a family B1 GCPR. The intra and extra cellular loops and the C-terminal have been modeled from fragments retrieved from the PDB. On the basis of the experimental data available for some members of family B1 GPCRs, four pairs of constraints between GIP(1-30)NH2 and its receptor were used in the FTDOCK program, to build the complete model of the GIP(1-30)NH2:GIPR complex. The model can rationalize the various experimental observations including the potency of the truncated GIP peptide. This work is the first complete model at the atomic level of GIP(1-30)NH2 and of the complex with its GPCR.     

Secondary conformations and temperature effect on structural transformation of amyloid beta (1-28), (1-40) and (1-42) peptides

Secondary structure of three amyloid b-peptides [A beta(1-28), A beta(1-40) and A beta(1-42)] in the solid state was respectively determined by Fourier transform infrared (FT-IR) microspectroscopy. Their thermal-dependent structural transformation were also investigated by FT-IR microspectroscopy equipped with a thermal analyzer. The present result demonstrates that the solid-state A beta(1-28), A beta(1-40) and A beta(1-42) peptides showed a significant IR spectral difference in the amide I and II bands. The secondary conformation of A beta(1-28) peptide was the combination of major beta-sheet and minor alpha-helix with little random coil structures, but A beta(1-40) peptide showed the co-existence of major beta-sheet and minor random coil with little alpha-helix structures. A beta(1-42) peptide mainly consisted of the predominant b-sheet structure. Although the intact A beta(1-28), A beta(1-40) or A beta(1-42) peptide exhibits a different secondary structure, a similar beta-conformation may form after thermal treatment. A thermal-dependent transition was found for solid A beta(1-28) and A beta(1-40) peptides near 40 degrees C and 45 degrees C, respectively. There was no transition temperature for solid A beta(1-42) peptide, however, due to only a very little level of alpha-helix and random coil structure containing in the solid A beta(1-42) peptide. The thermal denaturation plays an important role in the structural transformation from alpha-helix/random coil to beta-sheet.     

Loop propensity of the sequence YKGQP from staphylococcal nuclease: implications for the folding of nuclease.

Recently we performed molecular dynamics (MD) simulations on the folding of the hairpin peptide DTVKLMYKGQPMTFR from staphylococcal nuclease in explicit water. We found that the peptide folds into a hairpin conformation with native and nonnative hydrogen-bonding patterns. In all the folding events observed in the folding of the hairpin peptide, loop formation involving the region YKGQP was an important event. In order to trace the origins of the loop propensity of the sequence YKGQP, we performed MD simulations on the sequence starting from extended, polyproline II and native type I' turn conformations for a total simulation length of 300 ns, using the GROMOS96 force field under constant volume and temperature (NVT) conditions. The free-energy landscape of the peptide YKGQP shows minima corresponding to loop conformation with Tyr and Pro side-chain association, turn and extended conformational forms, with modest free-energy barriers separating the minima. To elucidate the role of Gly in facilitating loop formation, we also performed MD simulations of the mutated peptide YKAQP (Gly --> Ala mutation) under similar conditions starting from polyproline II conformation for 100 ns. Two minima corresponding to bend/turn and extended conformations were observed in the free-energy landscape for the peptide YKAQP. The free-energy barrier between the minima in the free-energy landscape of the peptide YKAQP was also modest. Loop conformation is largely sampled by the YKGQP peptide, while extended conformation is largely sampled by the YKAQP peptide. We also explain why the YKGQP sequence samples type II turn conformation in these simulations, whereas the sequence as part of the hairpin peptide DTVKLMYKGQPMTFR samples type I' turn conformation both in the X-ray crystal structure and in our earlier simulations on the folding of the hairpin peptide. We discuss the implications of our results to the folding of the staphylococcal nuclease.     

Effect of lipid composition on buforin II structure and membrane entry

Buforin II is a 21-amino acid polycationic antimicrobial peptide derived from a peptide originally isolated from the stomach tissue of the Asian toad Bufo bufo gargarizans. It is hypothesized to target a wide range of bacteria by translocating into cells without membrane permeabilization and binding to nucleic acids. Previous research found that the structure and membrane interactions of buforin II are related to lipid composition. In this study, we used molecular dynamics (MD) simulations along with lipid vesicle experiments to gain insight into how buforin II interacts differently with phosphatidylcholine (PC), phosphatidylglycerol (PG), and phosphatidylethanolamine (PE) lipids. Fluorescent spectroscopic measurements agreed with the previous assertion that buforin II does not interact with pure PC vesicles. Nonetheless, the reduced entry of the peptide into anionic PG membranes versus neutral PC membranes during simulations correlates with the experimentally observed reduction in BF2 translocation through pure PG membranes. Simulations showing membrane entry into PC also provide insight into how buforin II may initially penetrate cell membranes. Our MD simulations also allowed us to consider how neutral PE lipids affect the peptide differently than PC. In particular, the peptide had a more helical secondary structure in simulations with PE lipids. A change in structure was also apparent in circular dichroism measurements. PE also reduced membrane entry in simulations, which correlates with decreased translocation in the presence of PE observed in previous studies. Together, these results provide molecular-level insight into how lipid composition can affect buforin II structure and function and will be useful in efforts to design peptides with desired antimicrobial and cell-penetrating properties.     

NMR conformational analysis of proadrenomedullin N-terminal 20 peptide, a proangiogenic factor involved in tumor growth

The preferred conformation of Proadrenomedullin N-Terminal 20 Peptide (PAMP; ARLDVASEFRKKWNKWALSR-amide) has been determined using 1H and 13C two-dimensional nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. PAMP is a peptide that has various physiological functions, including its role as a proangiogenic factor in facilitating tumor growth and its inhibitory effect on catecholamine secretion at nicotinic receptors. The preferred conformation of PAMP was determined in a helix-inducing trifluoroethanol and water (TFE/H2O) solution, and in a membrane-mimetic sodium dodecylsulfate-d25 (SDS) micellar solution. The secondary structure consists of an alpha-helix for residues Arg2 to Arg20 in TFE/H2O solution and an alpha-helix for residues Arg2 to Ala17 in SDS solution. We postulate that the polar charged residues Arg2, Lys12, and Arg20 are responsible for the initial interaction of the peptide with the micelle, and that this is followed by the binding of the hydrophobic residues Leu3, Val5, Phe9, Trp13, and Trp16 to the micellar core. The three C-terminal amino acid residues adopt an extended structure in SDS, suggesting that they are important in receptor recognition and binding. This is supported by truncation studies done by Mahata et al. (Hypertension, 1998, Vol. 32, pp. 907-916), which show the importance of the C-terminal in physiological activity. Furthermore, Belloni et al. (Hypertension, 1999, Vol. 33, pp. 1185-1189), and Martinez et al. (Cancer Research, 2004, Vol. 64, pp. 6489-6494) suggested that the N-terminal was also important in PAMP activity. However, no differences in conformational preference of the N-terminal were observed between the two solvent systems.