Theoretical conformational analysis of tripeptides. N-acetyl-L-alanyl-L-alanyl-L-alanyl methylamide]

The minimization procedure has been used for calculation of the local minimum conformations of threepeptide--Ac-(L-Ala)3-NHMe without intramolecular H-bonds. The significant energy deviations from additivity found, arising with increase backbone length to three links, can be considered as the evidence for mutual dependence of conformational states of the neighbouring and terminal amino acid residues. It have been shown that stability of alpha-helix form for alanine threepeptide in contrary to corresponding dipeptide is noticeably higher due to stabilizing effect of dispersion interactions. The results of calculations are compared with the data on conformational distrubution of the threepeptide fragments in proteins with known three dimensional structure. The important role of the backbone interaction in protein chain have been marked.     

Critical Main-Chain Length for Conformational Conversion From 3(10)-Helix to α-Helix in Polypeptides

To assess the minimal peptide length required for the stabilization of the alpha-helix relative to the 3(10)-helix in Aib-rich peptides, we have solved the X-ray diffraction structures of the terminally blocked sequential hexa- and octapeptides with the general formula-(Aib-L-Ala)n-(n = 3 and 4, respectively). The hexapeptide molecules are completely 3(10)-helical with four 1----4 intramolecular N-H . . . O = C H-bonds. On the other hand, the octapeptide molecules are essentially alpha-helical with four 1----5 H-bonds; however, the helix is elongated at the N-terminus, with two 1----4 H-bonds, giving these molecules a mixed alpha/3(10)-helical character. In both compounds the right-handed screw sense of the helix is dictated by the presence of the Ala residues of L-configuration. This study represents the first experimental proof for a 3(10)----alpha-helix conversion in the crystal state induced by peptide backbone lengthening only.     

Hydrogen-bonding interactions in adrenaline–water complexes: DFT and QTAIM studies of structures, properties, and topologies

ωB97XD/6-311++G(d,p) calculations were carried out to investigate the hydrogen-bonding interactions between adrenaline (Ad) and water. Six Ad-H(2)O complexes possessing various types of hydrogen bonds (H-bonds) were characterized in terms of their geometries, energies, vibrational frequencies, and electron-density topology. Natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM) analyses were performed to elucidate the nature of the hydrogen-bonding interactions in these complexes. The intramolecular H-bond between the amino and carboxyl oxygen atom of Ad was retained in most of the complexes, and cooperativity between the intra- and intermolecular H-bonds was present in some of the complexes. H-bonds in which hydroxyls of Ad/water acted as proton donors were stronger than other H-bonds. Both hydrogen-bonding interactions and structural deformation play important roles in the relative stabilities of the complexes. The intramolecular H-bond was broken during the formation of the most stable complex, which indicates that Ad tends to break the intramolecular H-bond and form two new intermolecular H-bonds with the first water molecule.     

How hydrogen bonds impact P-glycoprotein transport and permeability

The requirement to cross a biological membrane can be a complex process especially if multidrug transporters such as P-gp must be considered. Drug partitioning into the lipid membrane and efflux by P-gp are tightly coupled processes wherein H-bonding interactions play a key role. All H-bond donors and acceptors are not equal in terms of the strength of the H-bonds that they form, hence it is important to consider their relative strength. Using various examples from literature, we illustrate the benefits of considering the relative strengths of individual H-bonds and introducing intramolecular H-bonds to increase membrane permeability and/or decrease P-gp efflux.     

Conformational studies of sphingolipids by NMR spectroscopy. I. Dihydrosphingomyelin

The conformational features of dihydrosphingomyelin (DHSM), the major phospholipid of human lens membranes, were investigated by 1H and 31P nuclear magnetic resonance spectroscopy. Several postulates emerge from the observed trends: (a) in partially hydrated samples of DHSM in CDCl3 above 13 mM, at which lipid-lipid interactions prevail, the amide proton is mostly involved in intermolecular H-bonds that link neighboring phospholipids through bridging water molecules. In the absence of water, the NH group is involved in an intramolecular H-bond that restricts the mobility of the phosphate group. (b) In the monomeric form of the lipid molecule, the amide proton of the major conformer is bound intramolecularly with one of the anionic and/or ester oxygens of the phosphate group. A minor conformer may also be present in which the NH proton participates in an intramolecular H-bond linking to the OH group of the sphingoid base. (c) Complete hydration leads to an extension of the head group as water molecules bind to the phosphate and NH groups via H-bonds, thus disrupting the intramolecular H-bonds prevalent at low concentrations.     

Observation of the closing of individual hydrogen bonds during TFE-induced helix formation in a peptide

Helix formation of an S-peptide analog, comprising the first 20 residues of Ribonuclease A and two additional N-terminal residues, was studied by measuring hydrogen bond (H-bond) (h3)J(NC') scalar couplings as a function of 2,2,2-trifluoroethanol (TFE) concentration. The (h3)J(NC') couplings give direct evidence for the closing of individual backbone N-H***O = C H-bonds during the TFE-induced formation of secondary structure. Whereas no (h3)J(NC') correlations could be detected without TFE, alpha-helical (i,i +4) H-bond correlations were observed for the amides of residues A5 to M15 in the presence of TFE. The analysis of individual coupling constants indicates that alpha-helix formation starts at the center of the S-peptide around residue E11 and proceeds gradually from there to both peptide ends as the TFE concentration is increased. At 60% to 90% TFE, well-formed alpha-helical H-bonds were observed for the amides hydrogens of residues K9 to Q13, whereas H-bonds of residues T5 to A8, H14, and M15 are affected by fraying. No intramolecular backbone H-bonds are present at and beyond the putative helix stop signal D16. As the (h3)J(NC') constants represent ensemble averages and the dependence of (h3)J(NC') on H-bond lengths is very steep, the size of the individual (h3)J(NC') coupling constants can be used as a measure for the population of a closed H-bond. These individual populations are in agreement with results derived from the Lifson-Roig theory for coil-to-helix transitions. The present work shows that the closing of individual H-bonds during TFE-induced helix formation can be monitored by changes in the size of H-bond scalar couplings.     

Molecular dynamics simulations of a guaiacyl beta-O-4 lignin model compound: examination of intramolecular hydrogen bonding and conformational flexibility

The dynamical conformational behavior of a guaiacyl beta-O-4 lignin model compound has been investigated by molecular simulations. The potential energy surface of the molecule in vacuum has been examined by means of an adiabatic map, showing a large accessible conformational space with multiple energy minima separated by low barriers. Molecular dynamics simulations have been performed in vacuum and with explicit solvent molecules for 10 and 2.1 ns, respectively. Molecular dynamics trajectories recorded in vacuum have shown the molecule to be flexible and to visit a large number of conformations. Many intramolecular H-bonds have been observed, existing for more than 90% of the total simulation time. The presence of explicit solvent molecules induces a significant broadening of some regions of the accessible conformational space and also largely reduces the statistical significance of intramolecular H-bonding. Intramolecular H-bonds observed in vacuum do not persist significantly and are preferentially exchanged with intermolecular H-bonds to the surrounding solvent molecules. The theoretical results are in good agreement with experimental NMR data that do not support the existence of strong and persistent intramolecular H-bonds in solution but instead indicate that H-bonds to solvent predominate. Finally, both molecular modeling and NMR approaches predict the guaiacyl beta-O-4 structure to be flexible and indicate that intramolecular H-bonds are not strong and persistent enough to confer rigidity to the molecule in solution.     

Conformations of proline-containing cyclic peptides. 1H and 13C n.m.r. evidence for a solvent stabilized All-Cis X-Pro conformer of Cyclo-(Pro-Gly-Gly-Pro)2

As inferred from 13C, 1H n.m.r. data, CD measurements and ion-binding experiments, the title molecule can assume two major C2 symmetric conformations. One of these has an all-trans X-Pro peptide backbone and two 1 comes from 4 intramolecular H-bonds and represents the predominant (greater than or equal to 95%) form in D2O and nonpolar (CD3CN) solvents. Stabilized by specific solvent-solute interactions, the other conformer becomes competitive (45%) in DMSO solution. It is shown to possess a four-cis X-Pro skeleton and no intramolecular H-bonds. The Mg++ complex of the cyclic peptide in CD3CN is again C2 symmetric and its formation proceeds via a slow trans leads to cis isomerization of two X-Pro peptide bonds.     

Comparative study of SP[6-11] and its analogs using simulated annealing

In this study we compared the steric structures of the bioactive part of substance P (SP[6-11]) and its analogs (NY3460 and pHOPA-SP5). The molecular dynamics-simulated annealing method was used to explore the conformational space, and the structural differences and similarities of these molecules were identified. For the three peptides, the conformational distributions were represented in Ramachandran density plots. The occurring secondary structural elements of the investigated molecules were identified, namely alpha-Helix, type III beta-Turn, gamma-Turn, and inverse gamma-Turn. For SP[6-11] and its two analogs, different intramolecular interactions (H-bonds between the main-chain atoms, aromatic-aromatic interactions, and amino-aromatic interactions) that can stabilize the various conformations of the three peptides were investigated. Detailed examination of these intramolecular interactions revealed that H-bonds between the main-chain atoms are relevant in the determination and stabilization of the conformer structures of the peptides, while the aromatic-aromatic interactions do not play an important stabilizing role. Furthermore, in the conformers of NY3460 and pHOPA-SP5, different types of amino-aromatic interactions were identified that contribute to the formation of the various structures of these peptides. For all three molecules, the orientations of the side chains were investigated and the rotamer populations were determined.     

New Fourier transform infrared based computational method for peptide secondary structure determination. II. Application to study of peptide fragments reproducing processing site of ocytocin-neurophysin precursor

A new method for the quantitative determination of the percentage of intramolecular H-bonds, based on Fourier transform infrared techniques, is applied to the conformational analysis of a series of synthetic peptides spanning the processing site of the ocytocin and neurophysin precursor. Even though the method uses traditional tools such as Fourier self-deconvolution, the Nth derivative, and curve-fitting procedures for the analysis of the spectra, the assignment of the absorptions due to peptide groups participating into secondary structures is based on the direct observation and quantification of the isotopic effect induced on the groups participating in intramolecular H-bonds in the presence of organic solvents. This permits the quantification of the different populations of molecules containing intramolecular H-bonds involved in beta-turns and alpha-helices. The results are consistent with those previously obtained by NMR techniques in the same solvent systems.     

Alpha-helix folding by Monte Carlo simulated annealing in isolated C-peptide of ribonuclease A

Conformation of the C-peptide fragment of RNase A is calculated by Monte Carlo simulated annealing. We adopt the total potential energy as given by the sum of generic interatomic energies whose parameters are determined separately for each amino acid without referring to the empirical structure of the C-peptide. The simulation is carried out in a completely unrestricted way without imposing any weight towards given final destinations. Starting from completely random initial conformations and minimizing the total potential energy with respect to main-chain dihedral angles and side-chain torsion angles, we have obtained partial alpha-helix structure with a high probability (approximately 40%). The energetically most favourable structure exhibits a 2.5-turn alpha-helix at the location identical with that of the 3-turn alpha-helix in the native enzyme molecule. Classification of conformations obtained in the simulation into clusters of similar structure shows that our simulation indeed predicts the alpha-helix structure for the isolated C-peptide with specific charged residues. The results of simulation with various amino acid substitutions are also found to be consistent with the experimental implication for the importance of intramolecular ionic interactions for alpha-helix stability for this peptide.     

Hydrogen bonding interactions in noradrenaline-DMSO complexes: DFT and QTAIM studies of structure,properties and topology

The hydrogen bonding interactions between noradrenaline (NA) and DMSO were studied with density functional theory (DFT) regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules (QTAIM) and the natural bond orbital (NBO) analyses were employed to elucidate the hydrogen bonding interaction characteristics in noradrenaline-DMSO complexes. The H-bonds involving the hydroxyls hydrogen in NA and the O atom in DMSO are dominant intermolecular H-bonds and are stronger than other H-bonds involving the methyl hydrogen of DMSO as a H-donor. The weak H-bonds also include a π H-bond which involves the benzene ring as a H-donor or H-acceptor. QTAIM identified the weak H-bonds formed between the methyl hydrogen of DMSO and the N atom in NA in some complexes (AB5, AB6 and AB7), which cannot be further confirmed by NBO and other methods, so there are probably no interactions between hydrogen and nitrogen atoms among these complexes. A good linear relationship between logarithmic electron density (lnρ _b ) at the bond critical point (BCP) and structural parameter (δR _H···Y) was found. The formations of new H-bonds in some complexes are helpful to strengthen the original intramolecular H-bond, this is attributed to the cooperativity of H-bonds in complexes and can be learned from the structure results and the NBO and QTAIM analyses. Analysis of various physically meaningful contributions arising from the energy decomposition procedures show that the orbital interactions of H-bond is predominant during the formation of the complex, moreover, both the hydrogen bonding interaction and the structural deformation are responsible for the stability of the complexes.     

Estimating the contribution of engineered surface electrostatic interactions to protein stability by using double-mutant cycles

Coulombic interactions between charges on the surface of proteins contribute to stability. It is difficult, however, to estimate their importance by protein engineering methods because mutation of one residue in an ion pair alters the energetics of many interactions in addition to the coulombic energy between the two components. We have estimated the interaction energy between two charged residues, Asp-12 and Arg-16, in an alpha-helix on the surface of a barnase mutant by invoking a double-mutant cycle involving wild-type enzyme (Asp-12, Thr-16), the single mutants Thr----Arg-16 and Asp----Ala-12, and the double mutant Asp----Ala-12, Thr----Arg-16. The changes in free energy of unfolding of the single mutants are not additive because of the coulombic interaction energy. Additivity is restored at high concentrations of salt that shield electrostatic interactions. The geometry of the ion pair in the mutant was assumed to be the same as that in the highly homologous ribonuclease from Bacillus intermedius, binase, which has Asp-12 and Arg-16 in the native enzyme. The ion pair does not form a hydrogen-bonded salt bridge, but the charges are separated by 5-6 A. The mutant barnase containing the ion pair Asp-12/Arg-16 is more stable than wild type by 0.5 kcal/mol, but only a part of the increased stability is attributable to the electrostatic interaction. We present a formal analysis of how double-mutant cycles can be used to measure the energetics of pairwise interactions.     

Parallel-stranded linear homoduplexes of d(A+-G)n > 10 and d(A-G)n > 10 manifesting the contrasting ionic strength sensitivities of poly(A+.A+) and DNA.

In contrast to shorter homologs which only form a single-stranded nucleic acid alpha-helix in acid solution at [Na+]</=0.02 M Na+, d(A-G)20,30 form in addition a parallel-stranded duplex with (A+.A+) and (G.G) base pairs and interstrand dA+...PO2-ionic and dA+NH2... O=P H-bonds. Under conditions where duplex prevails over alpha-helix, the contribution of the base-backbone interactions to stability varies directly with [H+] and inversely with [Na+], just as in poly(A+.A+). These duplexes are characterized by intense circular dichroism and a large cooperative thermally-induced hyperchromic transition that is dependent on oligomer concentration. Dimethylsulfate reactivity of the dG residues indicates G.G and therefore dA+.dA+rather than dA+.G base pairs. At much higher ionic strength (Na+>/=0.2 M) the protonated base-backbone interactions are so weakened that duplex stability becomes increasingly dependent upon H-bonded base pairing and stacking and almost independent of pH. Between pH 6 and 8 this duplex structure is devoid of protonated dA residues and shows positive dependence of T m on ionic strength similar to that of DNA.     

Conformational studies of sphingolipids by NMR spectroscopy. II. Sphingomyelin

Sphingomyelin (SM) is the most prevalent sphingolipid in the majority of mammalian membranes. Proton and 31P nuclear magnetic resonance spectral data were acquired to establish the nature of intra- and intermolecular H-bonds in the monomeric and aggregated forms of SM and to assess possible differences between this lipid and dihydrosphingomyelin (DHSM), which lacks the double bond between carbons 4 and 5 of the sphingoid base. The spectral trends suggest the formation of an intramolecular H-bond between the OH group of the sphingosine moiety and the phosphate ester oxygen of the head group. The narrower linewidth and the downfield shift of the resonance corresponding to OH proton in SM suggest that this H-bond is stronger in SM than in DHSM. The NH group appears to be involved predominantly in intramolecular H-bonding in the monomer. As the concentration of SM increases and the molecules come in closer proximity, these intramolecular bonds are partially disrupted and the NH group becomes involved in lipid-water interactions. The difference between the SM and DHSM appears to be not in the nature of these interactions but rather in the degree to which these intermolecular interactions prevail. As SM molecules cannot come as close together as DHSM molecules can, both the NH and OH moieties remain, on average, more intramolecularly bonded as compared to DHSM.     

Contacts between Tet repressor and tet operator revealed by new recognition specificities of single amino acid replacement mutants.

We have analyzed the DNA binding properties of Tet-repressor mutants with single amino acid residue replacements at eight positions within the alpha-helix-turn-alpha-helix DNA-binding motif. A saturation mutagenesis of Gln38, Pro39, Thr40, Tyr42, Trp43 and His44 in the second alpha-helix was performed; in addition, several substitutions of Thr27 and Arg28 in the first alpha-helix were constructed. The abilities of these mutant repressors to bind a set of 16 operator variants were determined and revealed 23 new binding specificities. All repressor mutants with DNA-binding activity were inducible by tetracycline, while mutants lacking binding activity were trans-dominant over the wild-type. All mutant proteins were present at the same intracellular steady-state concentrations as the wild-type. These results suggest the structural integrity of the mutant repressors. On the basis of the new recognition specificities, five contacts between a repressor monomer and each operator half-site and the chemical nature of these repressor-operator interactions are proposed. We suggest that Arg28 contacts guanine of the G.C base-pair at operator position 2 with two H-bonds, Gln38 binds adenine of the A.T base-pair at position 3 with two H-bonds, and the methyl group of Thr40 participates in a van der Waals' contact with cytosine of the G.C base-pair at position 6 of tet operator. A previously unrecognized type of interaction is proposed for Pro39, which inserts its side-chain between the methyl groups of the thymines of T.A and A.T base-pairs at positions 4 and 5. Computer modeling of these proposed contacts reveals that they are possible using the canonical structures of the helix-turn-helix motif and B-DNA. These contacts suggest an inverse orientation of the Tet repressor helix-turn-helix with respect to the operator center as compared with non-inducible repressor-operator complexes, and are supported by similar contacts of other repressor-operator complexes.     

Surface binding of alamethicin stabilizes its helical structure: molecular dynamics simulations.

Alamethicin is an amphipathic alpha-helical peptide that forms ion channels. An early event in channel formation is believed to be the binding of alamethicin to the surface of a lipid bilayer. Molecular dynamics simulations are used to compare the structural and dynamic properties of alamethicin in water and alamethicin bound to the surface of a phosphatidylcholine bilayer. The bilayer surface simulation corresponded to a loosely bound alamethicin molecule that interacted with lipid headgroups but did not penetrate the hydrophobic core of the bilayer. Both simulations started with the peptide molecule in an alpha-helical conformation and lasted 2 ns. In water, the helix started to unfold after approximately 300 ps and by the end of the simulation only the N-terminal region of the peptide remained alpha-helical and the molecule had collapsed into a more compact form. At the surface of the bilayer, loss of helicity was restricted to the C-terminal third of the molecule and the rod-shaped structure of the peptide was retained. In the surface simulation about 10% of the peptide/water H-bonds were replaced by peptide/lipid H-bonds. These simulations suggest that some degree of stabilization of an amphipathic alpha-helix occurs at a bilayer surface even without interactions between hydrophobic side chains and the acyl chain core of the bilayer.     

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.     

Exploring the interaction between D-xylose isomerase and D-xylose by free energy calculation

The numerical quadrature thermodynamic integration method is used to investigate enzyme-substrate interaction of D-xylose isomerase. A screening function for the coulombic interaction is introduced into the simulation to correct the effect of finite cutoff radius for the non-bonded interaction. The binding free energy difference for D-xylose with D-xylose isomerase and its N184D mutant has been calculated, and the result 3.9 ± 1.2 kJ/mol agrees well with experimental data of 4.38 kJ/mol. In addition, the structure and dynamics of enzyme-substrate complex were simulated for mutant and wild-type enzyme, respectively. Analysis of the structures and intramolecular interactions of the complexes were found to be valuable for understanding the reaction mechanism of the enzyme D-xylose isomerase. © Wiley-Liss, Inc.     

Mechanism by which the amyloid-like fibrils of a beta 2-microglobulin fragment are induced by fluorine-substituted alcohols

Although the formation of an alpha-helix or partial unfolding of proteins has been suggested to be important for amyloid fibrils to form in alcohols, the exact mechanism involved remains elusive. To obtain further insight into the development of amyloid fibrils, we used a 22-residue peptide, K3, corresponding to Ser20 to Lys41 of intact beta2-microglobulin. Although K3 formed an alpha-helix at high concentrations of 2,2,2-trifluoroethanol (TFE) and 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) in 10 mM HCl (pH approximately 2), the helical content was not high, indicating a low preference to do so. The partly alpha-helical conformation was converted with time into a highly ordered beta-sheet with a fibrillar morphology as revealed by atomic force microscopy. Importantly, the TFE and HFIP-induced fibrillation exhibited a concentration dependence with a maximum at approximately 20 and approximately 10% (v/v), respectively, slightly below the concentrations at which these alcohols form dynamic clusters. Focusing on the similarity of the effects of alcohol on proteins with those of sodium dodecyl sulfate (SDS), we examined the effects of SDS on K3. SDS also induced fibrils to form with a maximum at approximately 4 mM, slightly below the critical micelle concentration. These results indicate that, with an increase in the concentration of hydrophobic cosolvent (TFE, HFIP, or SDS), a delicate balance of decreasing hydrophobic interactions and increasing polar interactions (i.e. H-bonds) in and between peptides leads to the formation of ordered fibrils with a bell-shaped concentration dependence.