Helix Propensity Research Articles

Relationship Between Biophysical Properties of Antimicrobial Peptides (AMPs) and their Associated Drug Efficacies

Antibiotic resistance is a growing global health threat. One consequence is that patients with cystic fibrosis (CF) are prone to developing antibiotic resistant lung infections caused by multiple strains of bacteria, including Pseudomonas aeruginosa. Due to the limited number of treatment options for patients with chronic antibiotic resistant infections, there is a need for finding new antibiotics that allow for effective eradication of bacterial infections, such as those in the CF lung. Many antimicrobial peptides (AMPs) have been annotated in databases and are considered as potential alternatives for current antibiotics. However, in many instances, the suitability of AMPs as drug molecules has not been extensively explored. Here, we propose that certain molecular properties of AMPs favor high antibiotic efficacy. Using information from AMP databases, we combined statistical analyses and machine learning techniques to identify relationships between various biophysical properties of AMPs and their drug efficacies. Analyses from classification and regression trees (CART) and random forests suggest that net charge and maximum average hydrophobic moment are the most important properties in determining if a peptide is useful against P. aeruginosa infections in CF patients. Maximum average hydrophobic residue, average alpha helix propensity score, hydrophobic proportion, and peptide length still contribute to this determination but to lesser degrees. Cation-pi interactions, on the other hand, do not appear to factor into this decision at all. Based on these properties, our current work is focused on designing and experimentally testing new peptides that may have activity against P. aeruginosa infections.

Effects of single-point mutations on translocation of non-interacting heteropolymer protein chains

We analyse the role of single-point mutations on the translocation of heteropolymer protein chains with helical folding. We propose a statistical mechanics model that computes the equilibrium partition function of a given protein by incorporating experimental helix propensities of amino acids. From that we evaluate the free energy barrier for translocation and the average translocation time for single-point protein mutants and compared with wild-type proteins used as reference. In our calculations, single-point mutations that increase energy of the folded state in increase the energy barrier for translocation up to , , in consequence, the translocation time undergoes variations between and with the largest variations associated to substitutions of amino acids glycine and proline by alanine. In fact, substitutions of amino acids with the lowest helix propensities by amino acids with the largest helix propensities induced the largest mutational effects. These results indicate that effects of single-point mutations become evident when mutations increase energy of the folded state at least in . Our model predicts translocation time delays for proline to alanine single-point mutations which agree with deficiencies of protein translocation observed recently in experiments.

Induction and destruction of folding motifs in cyclotides by cyclic cystine knots

Multivariate roles of cyclic cystine knots (CCK) present in cyclotides are well known. More than often this marvellous combination of disulphide linkages embedded in a macrocyclic backbone is seen to induce and stabilize otherwise unfavourable secondary structural motifs in peptides. Partial or complete reduction of the disulphide bonds leads to significant overhauling of the peptide structure and associated dynamics with disappearance of some of the CCK-stabilized local motifs. In this work, we explore structural and dynamical intricacies in two prototypical bracelet and a Möbius cyclotide, respectively cycloviolacin O1 (cyO1), cycloviolacin O2 (cyO2) and kalata B1 (kB1), either in their native (N, with CCK), partially reduced (with either one or two S-S bonds) or completely reduced (D, without CCK) forms using molecular dynamics simulations. The S-S linkage(s) primarily responsible for sustaining a given organized motif in each case is identified from the simulation of the partially reduced forms. Correlation of helix propensity, strand propensity and hydropathicity indices of the amino acids conserved across families of different bracelet/Möbius cyclotides give insight to the natural inclination of the CCK stabilized backbones to exhibit a given motif. Dihedral principal component analysis (dPCA) and normal mode analysis (NMA) help scan through the effect of structural reorganization on peptide dynamics. The relatively restricted dynamical modes in CCK-intact native peptides are replaced with large amplitude dynamical fluctuations in the reduced peptides. Concomitant structural variance is observed across the clusters in the space spanned by the dominant principal components (PCs).

Effects of Replacing a Central Glycine Residue in GLP-1 on Receptor Affinity and Signaling Profile.

Agonists of the glucagon-like peptide-1 receptor (GLP-1R) are used to treat diabetes and obesity. Cryo-EM structures indicate that GLP-1 is completely α-helical when bound to the GLP-1R. The mature form of this hormone, GLP-1(7-36), contains a glycine residue near the center (Gly22). Since glycine has the second-lowest α-helix propensity among the proteinogenic α-amino acid residues, and Gly22 does not appear to make direct contact with the receptor, we were motivated to explore the impact on agonist activity of altering the α-helix propensity at this position. We examined GLP-1 analogues in which Gly22 was replaced with L-Ala, D-Ala, or β-amino acid residues with varying helix propensities. The results suggest that the receptor is reasonably tolerant of variations in helix propensity, and that the functional receptor-agonist complex may comprise a conformational spectrum rather than a single fixed structure.

Hydrophobic mismatch and sequence specificity compete when transmembrane helix-helix interactions are measured with the TOXCAT assay.

Genetic assays capable of measuring the propensity of transmembrane helices to oligomerize within the cytoplasmic membrane of the bacterium E. coli are frequently used when sequence-specificity in transmembrane helix-helix interactions is investigated. In the present study, dimerization of the well-investigated wild-type and G83I-mutated transmembrane helix of the human glycophorin A protein was studied. Gradual prolongation of the transmembrane helix at the C-terminus with Leu residues lead to pronounced changes in the dimerization propensity when measured with the TOXCAT assay. Thus, besides sequence specificity, hydrophobic mismatch between the hydrophobic core of a studied transmembrane helix and the E. coli membrane can impact the oligomerization propensity of a transmembrane helix. This suggests that the results of genetic assays aiming at determining interactions of heterologous transmembrane helices within the E. coli membrane do not necessarily solely reflect sequence specificity in transmembrane helix-helix interactions, but might be additionally modulated by topological and structural effects caused by hydrophobic mismatch.

Resemblance-Ranking Peptide Library to Screen for Binders to Antibodies on a Peptidomic Scale.

A novel resemblance-ranking peptide library with 160,000 10-meric peptides was designed to search for selective binders to antibodies. The resemblance-ranking principle enabled the selection of sequences that are most similar to the human peptidome. The library was synthesized with ultra-high-density peptide arrays. As proof of principle, screens for selective binders were performed for the therapeutic anti-CD20 antibody rituximab. Several features in the amino acid composition of antibody-binding peptides were identified. The selective affinity of rituximab increased with an increase in the number of hydrophobic amino acids in a peptide, mainly tryptophan and phenylalanine, while a total charge of the peptide remained relatively small. Peptides with a higher affinity exhibited a lower sum helix propensity. For the 30 strongest peptide binders, a substitutional analysis was performed to determine dissociation constants and the invariant amino acids for binding to rituximab. The strongest selective peptides had a dissociation constant in the hundreds of the nano-molar range. The substitutional analysis revealed a specific hydrophobic epitope for rituximab. To show that conformational binders can, in principle, be detected in array format, cyclic peptide substitutions that are similar to the target of rituximab were investigated. Since the specific binders selected via the resemblance-ranking peptide library were based on the hydrophobic interactions that are widespread in the world of biomolecules, the library can be used to screen for potential linear epitopes that may provide information about the cross-reactivity of antibodies.

Additivity in the helix-coil transition of alpha-helical peptides with varying length

The mechanical unfolding of alanine-rich helices with various lengths and sequences has been characterized using Adaptive Steered Molecular Dynamics (ASMD). We have found that the potential of mean force (PMF) and hydrogen bonding profiles for these various helices exhibit a correspondence between different chains, and the required energetics for the overall unfolding process is additive with respect to the chain length. [Biophys. J, 120, 2009 (2021)] Polyalanine peptides were used to guide the hypotheses because of their high helix propensity for helical formation, and the low computational cost needed to simulate them.

The role of alpha-helix on the structure-targeting drug design of amyloidogenic proteins

The most accredited hypothesis links the toxicity of amyloid proteins to their harmful effects on membrane integrity through the formation of prefibrillar-transient oligomers able to disrupt cell membranes. However, damage mechanisms necessarily assume a first step in which the amyloidogenic protein transfers from the aqueous phase to the membrane hydrophobic core. This determinant step is still poorly understood. However, according to our lipid-chaperon hypothesis, free lipids in solution play a crucial role in facilitating this footfall. Free phospholipid concentration in the aqueous phase acts as a switch between ion channel-like pore and fibril formation, so that high free lipid concentration in solution promotes pore and repress fibril formation. Conversely, low free lipids in the solution favor fibril and repress pore formation. This behavior is due to the formation of stable lipid-protein complexes. Here, we hypothesize that the helix propensity is a fundamental requirement to fulfill the lipid-chaperon model. The alpha-helix region seems to be responsible for the binding with amphiphilic molecules fostering the proposed mechanism. Indeed, our results show the dependency of protein-lipid binding from the helical structure presence. When the helix content is substantially lower than the wild type, the contact probability decreases. Instead, if the helix is broadening, the contact probability increases. Our findings open a new perspective for in silico screening of secondary structure-targeting drugs of amyloidogenic proteins.

Physicochemical and structural properties of lunasin revealed by spectroscopic, chromatographic and molecular dynamics approaches

Lunasin is a 43-amino acid peptide from seeds and grains with bioavailability in humans and potent chemotherapeutic action against several cancer cell lines. Here, we investigate new information about the physicochemical and structural properties of lunasin using circular dichroism (CD), fluorescence spectroscopy, electrospray ionization-ion mobility spectrometry-mass spectrometry (ESI-IMS-MS), size exclusion chromatography (SEC), molecular dynamics (MD), and bioinformatics. CD analysis and disorder prediction obtained by PONDR indicate that lunasin has a mostly unordered structure. Double wavelength [θ]222nm x [θ]200nm plot data suggests that lunasin is an intrinsically disordered peptide (IDP) in a pre-molten globule-like (PMG-like) state, while CD spectrum deconvolution and MD simulation indicate small β-strand content. The presence of residual structure was supported by loss of CD signal at 222 nm after treatment with urea and by increasing fluorescence emission upon bis-ANS binding. Lunasin also demonstrated stability to heating up to the temperature of 100 °C, as verified by CD. MD and CD analyses in the presence of TFE and MoRFpred prediction indicated the helix propensity of lunasin. ESI-IMS-MS data revealed that lunasin shows a propensity to form disulfide bonds at the conditions used. MD data also indicated that disulfide bond formation affects the adopted structure, showing a possible role of aspartyl-end in structure stabilization and compaction. In conclusion, our data support a characterization of lunasin as a peptide with an intrinsic disorder in a PMG-like state and reveal new aspects about its structural stability and plasticity, as well as the effects of disulfide bond formation and electrostatic attractions.

Secondary structure effects on internal proton transfer in poly-peptides.

A pump–probe approach was designed to determine the internal proton transfer (PT) rate in a series of poly-peptide radical cations containing both histidine and tryptophan. The proton transfer is driven by the gas-phase basicity difference between residues. The fragmentation scheme indicates that the gas-phase basicity of histidine is lower than that of radical tryptophan so that histidine is always pulling the proton away from tryptophan. However, the proton transfer requires the two basic sites to be in close proximity, which is rate limited by the peptide conformational dynamics. PT rate measurements were used to probe and explore the peptide conformational dynamics in several poly-glycines/prolines/alanines. For small and unstructured peptides, the PT rate decreases with the size, as expected from a statistical point of view in a flat conformational space. Conversely, if structured conformations are accessible, the structural flexibility of the peptide is decreased. This slows down the occurrence of conformations favorable to proton transfer. A dramatic decrease in the PT rates was observed for peptides HAnW, when n changes from 5 to 6. This is attributed to the onset of a stable helix for n = 6. No such discontinuity is observed for poly-glycines or poly-prolines. In HAnW, the gas-phase basicity and helix propensity compete for the position of the charge. Interestingly, in this competition between PT and helix formation in HA6W, the energy gain associated with helix formation is large enough to slow down the PT beyond experimental time but does not ultimately prevail over the proton preference for histidine.

Structural preferences of cysteine sulfinic acid: The sulfinate engages in multiple local interactions with the peptide backbone

Cysteine sulfinic acid (Cys-SO2–) is a non-enzymatic oxidative post-translational modification (PTM) that has been identified in hundreds of proteins. However, the effects of cysteine sulfination are in most cases poorly understood. Cys-SO2– is structurally distinctive, with long sulfur-carbon and sulfur-oxygen bonds, and with tetrahedral geometry around sulfur due to its lone pair. Cys-SO2– thus has a unique range of potential interactions with the protein backbone which could facilitate protein structural changes. Herein, the structural effects of cysteine oxidation to the sulfinic acid were investigated in model peptides and folded proteins using NMR spectroscopy, circular dichroism, bioinformatics, and computational studies. In the PDB, Cys-SO2– shows a greater preference for α-helix than Cys. In addition, Cys-SO2– is more commonly found in structures with φ > 0, including in multiple types of β-turn. Sulfinate oxygens engage in hydrogen bonds with adjacent (i or i + 1) amide hydrogens. Over half of sulfinates have at least one hydrogen bond with an adjacent amide, and several structures have hydrogen bonds with both adjacent amides. Alternately, sulfur or either oxygen can act as an electron donor for n→π* interactions with the backbone carbonyl of the same residue, as indicated by frequent S⋯CO or O⋯CO distances below the sums of their van der Waals radii in protein structures. In peptides, Cys-SO2– favored α-helical structure at the N-terminus, consistent with helix dipole effects and backbone hydrogen bonds with the sulfinate promoting α-helix. Cys-SO2– has only modestly greater polyproline II helix propensity than Cys-SH, likely due to competition from multiple side chain-backbone interactions. Cys-SO2– stabilizes the i+1 position of a β-turn relative to Cys-SH. Within proteins, the range of side chain-main chain interactions available to Cys-SO2– compared to Cys-SH provides a basis for potential changes in protein structure and function due to cysteine oxidation to the sulfinic acid.

On predicting foldability of a protein from its sequence.

Several properties of amino acid sequences corresponding to proteins that are known to fold are compared to those of randomly generated sequences and to sequences of intrinsically disordered proteins in order to find properties that distinguish folding sequences from the rest. The properties studied included helix and sheet propensities from secondary structure prediction, adjacency correlations, directionality correlations, as well as propensities of all possible triplets and quadruplets. Small differences between known folded and random sequences were observed for the adjacency and directional correlations, and significant differences were seen on the triplet and especially on the quadruplet propensities. Based on the differences in the adjacency, triplet or quadruplet propensities folding scores were defined and used to test the accuracy of foldability prediction based on these statistics. The best predictions were obtained from the quadruplet propensities.

A Bone Morphogenetic Protein (BMP)-derived Peptide Based on the Type I Receptor-binding Site Modifies Cell-type Dependent BMP Signalling

Bone morphogenetic proteins (BMPs) are multifunctional cytokines of the transforming growth factor β (TGFβ) superfamily with potential therapeutic applications due to their broad biological functionality. Designing BMP mimetics with specific activity will contribute to the translational potential of BMP-based therapies. Here, we report a BMP9 peptide mimetic, P3, designed from the type I receptor binding site, which showed millimolar binding affinities for the type I receptor activin receptor like kinase 1 (ALK1), ALK2 and ALK3. Although showing no baseline activity, P3 significantly enhanced BMP9-induced Smad1/5 phosphorylation as well as ID1, BMPR2, HEY1 and HEY2 gene expression in pulmonary artery endothelial cells (hPAECs), and this activity is dependent on its alpha helix propensity. However, in human dermal microvascular endothelial cells, P3 did not affect BMP9-induced Smad1/5 phosphorylation, but potently inhibited ALK3-dependent BMP4-induced Smad1/5 phosphorylation and gene expression. In C2C12 mouse myoblast cells, P3 had no effect on BMP9-induced osteogenic signalling, which is primarily mediated by ALK2. Interestingly, a previously published peptide from the knuckle region of BMP9 was found to inhibit BMP4-induced Smad1/5 phosphorylation. Together, our data identify a BMP9-derived peptide that can selectively enhance ALK1-mediated BMP9 signalling in hPAECs and modulate BMP9 and BMP4 signalling in a cell type-specific manner.

Effects of Arginine Deimination and Citrulline Side-Chain Length on Peptide Secondary Structure Formation.

Post-translational modifications expand the chemical functionality of peptides and proteins beyond that originating from the encoded amino acids, but studies on the structural effects of these modifications have been limited. Arginine undergoes deimination to give citrulline (Cit), converting the positively charged guanidinium moiety into a neutral urea group. Herein, we report the effect of Arg deimination on secondary structure formation. To understand the reason for the number of methylene units in Cit, the effect of Cit side-chain length on secondary structure formation was also studied. Ala-based peptides and β-hairpin peptides were used to study α-helix and β-sheet formation, respectively. Peptides containing Cit analogues were prepared by an orthogonal protecting group strategy coupled with solid-phase carbamylation. The CD data for the Ala-based peptides were analyzed by using modified Lifson-Roig theory, showing that the helix propensity of Arg decreased upon deimination and that either shortening or lengthening Cit also decreased the helix propensity. The β-hairpin peptides were analyzed by NMR methods, showing minimal change in strand formation energetics upon Arg deimination. Altering the Cit side-chain length did not affect strand formation energetics either. These results should be useful for the preparation of urea-bearing systems and the design of peptides incorporating urea-bearing residues with varying side-chain length.

Directed evolution of α-amylase from Bacillus licheniformis to enhance its acid-stable performance

α-Amylases (1,4-α-D-glucanohydrolases) are widely used in starch liquefaction, but their acid stability needs to be continuously explored to reduce the costs of raw material and operation. In this study, to better meet the industrial requirements, the acid stability of Bacillus licheniformis α-amylase (BLA) was further improved by directed evolution using error prone polymerase chain reaction (PCR). The mutant BLA (G81R) was selected with the improved acid stability based on a high-throughput activity assay. After incubating at pH 4.5 for 40 min, G81R still retained 10% of its initial activity, but the wild-type (WT) was already inactive. The kcat/Km value of G81R at pH 4.5 was 1.4-fold higher than that of WT. Combined with the three-dimensional structural modeling analysis, the improved stability of G81R under low pH condition might be due to the interactions of electrostatic, hydrophilicity, and helix propensity. Therefore, these findings would be beneficial for developing BLA with properties suitable for applications in industrial starch processing.

Thermodynamics of Helix formation in small peptides of varying lengthin vacuo, implicit solvent and explicit solvent: Comparison between AMBER force fields

Helix formation is of great significance in protein folding. The helix-forming tendencies of amino acids are accumulated along the sequence to determine the helix-forming tendency of peptides. Computer simulation can be used to model this process in atomic details and give structural insights. In the current work, we employ equilibrate-state free energy simulation to systematically study the folding/unfolding thermodynamics of a series of mutated peptides. Two AMBER force fields including AMBER99SB and AMBER14SB are compared. The new 14SB force field uses refitted torsion parameters compared with 99SB and they share the same atomic charge scheme. We find that in vacuo the helix formation is mutation dependent, which reflects the different helix propensities of different amino acids. In general, there are helix formers, helix indifferent groups and helix breakers. The helical structure becomes more favored when the length of the sequence becomes longer, which arises from the formation of additional backbone hydrogen bonds in the lengthened sequence. Therefore, the helix indifferent groups and helix breakers will become helix formers in long sequences. Also, protonation-dependent helix formation is observed for ionizable groups. In 14SB, the helical structures are more stable than in 99SB and differences can be observed in their grouping schemes, especially in the helix indifferent group. In solvents, all mutations are helix indifferent due to protein–solvent interactions. The decrease in the number of backbone hydrogen bonds is the same with the increase in the number of protein–water hydrogen bonds. The 14SB in explicit solvent is able to capture the free energy minima in the helical state while 14SB in implicit solvent, 99SB in explicit solvent and 99SB in implicit solvent cannot. The helix propensities calculated under 14SB agree with the corresponding experimental values, while the 99SB results obviously deviate from the references. Hence, implicit solvent models are unable to correctly describe the thermodynamics even for the simple helix formation in isolated peptides. Well-developed force fields and explicit solvents are needed to correctly describe the protein dynamics. Aside from the free energy, differences in conformational ensemble under different force fields in different solvent models are observed. The numbers of hydrogen bonds formed under different force fields agree and they are mostly determined by the solvent model.

Structural determinants of coiled coil mechanics

The natural abundance of coiled coil (CC) motifs in the cytoskeleton and the extracellular matrix suggests that CCs play a crucial role in the bidirectional mechanobiochemical signaling between cells and the matrix. Their functional importance and structural simplicity has allowed the development of numerous applications, such as protein-origami structures, drug delivery systems and biomaterials. With the goal of establishing CCs as nanomechanical building blocks, we investigated the importance of helix propensity and hydrophobic core packing on the mechanical stability of 4-heptad CC heterodimers. Using single-molecule force spectroscopy, we show that both parameters determine the force-induced dissociation in shear loading geometry; however, with different effects on the energy landscape. Decreasing the helix propensity lowers the transition barrier height, leading to a concomitant decrease in the distance to the transition state. In contrast, a less tightly packed hydrophobic core increases the distance to the transition state. We propose that this originates from a larger side chain dynamics, possible water intrusion at the interface as well as differences in solvation of the hydrophobic amino acids at the transition state. In conclusion, the different contributions of helix propensity and hydrophobic core packing need to be considered when tuning the mechanical properties of CCs for applications.

Thermodynamics of helix formation in small peptides of varying length in vacuo, in implicit solvent, and in explicit solvent.

Helix formation plays a critical role in the functionality of proteins. Various biological processes such as protein-ligand, protein-protein, and protein-nucleotide interactions are strongly influenced by the conformational stability of peptides. Small peptides represent ideal models to study the process of helix formation due to their small size and simple structures. As is widely known in the biophysical community, different amino acids have different propensities for helix formation. In this work, we used model peptides to study the peptide-length dependence and the sequence dependence of the thermodynamic profiles along the helix formation pathway in vacuo and in water. The peptides showed interesting unfolding behavior in vacuo, as this process was found to be protonation-state dependent, sequence-length dependent, and to change when residues in the peptide were mutated. Ionizable group protonation alters helix propensity by changing the degree of charge-charge repulsion in the peptide. Longer sequences tend to form more hydrogen bonds and are more likely to form helical structures. Mutation of an alanine residue to another naturally occurring amino acid changes the interactions within the peptide, which in turn alters its helix-forming tendency. By contrast, solvation leads to the generation of protein-solvent hydrogen bonds, which compensate for the energy penalty associated with the deformation of backbone hydrogen bonds in the extended conformation. Flat free-energy profiles were obtained for all mutated peptides of various lengths. Although several implicit-solvent models were able to reproduce free energies of helix formation in explicit solvent, the conformational ensembles varied with the solvent model applied, particularly in terms of the populations of secondary structures and the formation of backbone hydrogen bonds. Also, due to the absence of explicit descriptions of solvent molecules and thus protein-solvent hydrogen bonds in implicit-solvent models, they did not yield accurate hydrogen-bonding profiles.

Thermodynamic Scale of β-Amino Acid Residue Propensities for an α-Helix-like Conformation.

A thiol-thioester exchange system has been used to measure the propensities of diverse β-amino acid residues to participate in an α-helix-like conformation. These measurements depend on formation of a parallel coiled-coil tertiary structure when two peptide segments become linked by thioester formation. One peptide segment contains a "guest" site that accommodates diverse β residues and is distal to the coiled-coil interface. We find that helix propensity is influenced by side chain placement within the β residue [β3 (side chain adjacent to nitrogen) slightly favored relative to β2 (side chain adjacent to carbonyl)]. The previously recognized helix stabilization resulting from five-membered ring incorporation is quantified. These results are significant because so few quantitative thermodynamic measurements have been reported for α/β-peptide folding.

Hidden Structural Codes in Protein Intrinsic Disorder

Protein intrinsic disorder is a major structural category in biology yet its definition is often limited to the absence of folding. The explosion of information in the genomic era showed that it may account for over 30% of coding regions across life domains, and it is particularly overrepresented in viruses. Papillomaviruses are an unparalleled case for sequence to structure correlation analysis because of the existence of hundreds of anciently evolved and stable virus types which are divergent enough in sequence, but conserving the function of each protein. E7, the main transforming oncoprotein from human papillomaviruses, is a paradigmatic example of an intrinsically disordered protein with pathological moonlighting activities evolved for hijacking cell cycle control. Despite of being intrinsically disordered, the N-terminal domain shows more conserved residues than the globular C-terminal domain. Mutation of five hyper conserved residues precisely distributed along the sequence lead to a marked increase in both α-helix and ß-sheet structural content, reflected by drastic effects on equilibrium propensities and oligomerization kinetics. These results strongly suggest the existence of local nuclei, yet to be defined in structural terms, that oppose to canonical folding as expected for globular proteins. Moreover, direct coupling analysis pinpoint interactions surprisingly coincident with the proposed nuclei. In addition, these match with regions of helix propensities found in TFE/NMR experiments.