Lipid Conformational Nomenclature: A General Method

Abstract

We propose a conformational nomenclature for amphiphilic lipid molecules that is general and compatible with the stereospecific numbering scheme, in contrast to earlier methods in which discrepancies with the sn-scheme lead to contradictory assignments of the absolute configuration of the system. The present method can be rationally extended to different classes of lipids, both natural and synthetic. It is simple and provides a convenient framework for conformational studies on widely varying classes of lipids.     

Conformational Conversion May Precede or Follow Aggregate Elongation on Alternative Pathways of Amyloid Protofibril Formation

A major goal in the study of protein aggregation is to understand how the conformational heterogeneity characteristic of the process leads to structurally distinct amyloid fibrils. The small protein barstar is known to form amyloid protofibrils in multiple steps at low pH: a small oligomer, the A-form, first transforms into a larger spherical higher oligomeric intermediate (HOI), which then self-associates to form the elongated protofibril. To determine how the conformational conversion reaction during aggregation is coupled to the process of protofibril formation, cysteine-scanning mutagenesis was first used to identify specific residue positions in the protein sequence, which are important in defining the nature of the aggregation process. Two classes of mutant proteins, which are distinguished by their kinetics of aggregation at high protein concentration, have been identified: Class I mutant proteins undergo conformational conversion, as measured by an increase in thioflavin T binding ability and an increase in circular dichroism at 216 nm, significantly faster than Class II mutant proteins. At low protein concentration, the rates of conformational conversion are, however, identical for both classes of mutant proteins. At high protein concentration, the two classes of mutant proteins can be further distinguished on the basis of their rates of protofibril growth, as determined from dynamic light-scattering measurements. For Class I mutant proteins, protofibril elongation occurs at the same, or slightly faster, rate than conformational conversion. For Class II mutant proteins, protofibril elongation is significantly slower than conformational conversion. Dynamic light scattering measurements and atomic force microscopy imaging indicate that for the Class I mutant proteins, conformational conversion occurs concurrently with the self-association of prefibrillar HOIs into protofibrils. On the other hand, for the Class II mutant proteins, the prefibrillar HOI first undergoes conformational conversion, and the conformationally converted HOIs then self-associate to form protofibrils. The two classes of mutant proteins appear, therefore, to use structurally distinct pathways to form amyloid protofibrils. On one pathway, conformational conversion occurs along with, or after, elongation of the oligomers; on the other pathway, conformational conversion precedes elongation of the oligomers. Single mutations in the protein can cause aggregation to switch from one pathway to the other. Importantly, the protofibrils formed by the two classes of mutant proteins have significantly different diameters and different internal structures.     

Protein local conformations arise from a mixture of Gaussian distributions

The classical approaches for protein structure prediction rely either on homology of the protein sequence with a template structure or on ab initio calculations for energy minimization. These methods suffer from disadvantages such as the lack of availability of homologous template structures or intractably large conformational search space, respectively. The recently proposed fragment library based approaches first predict the local structures,which can be used in conjunction with the classical approaches of protein structure prediction. The accuracy of the predictions is dependent on the quality of the fragment library. In this work, we have constructed a library of local conformation classes purely based on geometric similarity. The local conformations are represented using Geometric Invariants, properties that remain unchanged under transformations such as translation and rotation, followed by dimension reduction via principal component analysis. The local conformations are then modeled as a mixture of Gaussian probability distribution functions (PDF). Each one of the Gaussian PDF's corresponds to a conformational class with the centroid representing the average structure of that class. We find 46 classes when we use an octapeptide as a unit of local conformation. The protein 3-D structure can now be described as a sequence of local conformational classes. Further, it was of interest to see whether the local conformations can be predicted from the amino acid sequences. To that end,we have analyzed the correlation between sequence features and the conformational classes.     

Protein local conformations arise from a mixture of Gaussian distributions

The classical approaches for protein structure prediction rely either on homology of the protein sequence with a template structure or on ab initio calculations for energy minimization. These methods suffer from disadvantages such as the lack of availability of homologous template structures or intractably large conformational search space, respectively. The recently proposed fragment library based approaches first predict the local structures, which can be used in conjunction with the classical approaches of protein structure prediction. The accuracy of the predictions is dependent on the quality of the fragment library. In this work, we have constructed a library of local conformation classes purely based on geometric similarity. The local conformations are represented using Geometric Invariants, properties that remain unchanged under transformations such as translation and rotation, followed by dimension reduction via principal component analysis. The local conformations are then modeled as a mixture of Gaussian probability distribution functions (PDF). Each one of the Gaussian PDF’s corresponds to a conformational class with the centroid representing the average structure of that class. We find 46 classes when we use an octapeptide as a unit of local conformation. The protein 3-D structure can now be described as a sequence of local conformational classes. Further, it was of interest to see whether the local conformations can be predicted from the amino acid sequences. To that end, we have analyzed the correlation between sequence features and the conformational classes.     

Angular disorder in actin: is it consistent with general principles of protein structure?

Harold Erickson has recently provided a useful analysis of helical structures having one class versus two classes of intersubunit bonds. His analysis is based upon an assumption that the subunits themselves are essentially unchanged upon bond formation (polymerization). He shows that such a structure having two classes of bonds (i.e. one in which each subunit interacts with four of its neighbors rather than two) can explain some of the features of actin. While he acknowledges that for actin there could be a conformational change and that, in principle, it could explain such features, he argues that the allowed magnitude of such a conformational change is inadequate. Since kinetics and thermodynamics cannot distinguish between the energy derived from the formation of a bond from that due to a conformational change, the question of whether the features of F-actin are derived from a conformational change or a system of two classes of bonds or both must be answered with high-resolution structural information. Recent studies by K. C. Holmes and others suggest that the second possibility might be closest to the truth. The heart of our disagreement is not whether Erickson's thermodynamic analysis is correct, given rigid subunits, but whether all protein polymers are characterized by rigid subunits with rigid intersubunit contacts. Erickson maintains that the observation of an angular disorder of 12 degrees per subunit within the actin filament conflicts with his formalism of rigid subunit interfaces and must therefore result from the erroneous interpretation of measurements. He presents an alternative model to explain the observations. His model, however, does not account for the observations and we will argue that, ultimately, like the resolution of the matter of the number of classes of bonds and the extent of their contact, the amount of angular disorder will require higher-resolution structural studies.     

MOUSE-III: Learning rules of conformational analysis from X-ray data

MOUSE-III is learning program that finds rules of conformational analysis from raw cristallographie data. The program perceives molecular features, finds conformational classes in the data and then learns rules that link features to classes. The rules that MOUSE learns are capable of correctly assigning conformations to ring systems that were not used for training with greater than 95% accuracy, when MOUSE was presented with sufficient data. The rules also show a compression of as much as 99% when compared to the raw data. This is accomplished through abstraction and generalization. The algorithm is presented along with a carefully worked example. An example of a learned rule is also presented and analyzed. Some conclusions about the scope and limitations of the learning process are presented.     

Side-chain conformational changes upon Protein-Protein Association

Conformational changes upon protein-protein association are the key element of the binding mechanism. The study presents a systematic large-scale analysis of such conformational changes in the side chains. The results indicate that short and long side chains have different propensities for the conformational changes. Long side chains with three or more dihedral angles are often subject to large conformational transition. Shorter residues with one or two dihedral angles typically undergo local conformational changes not leading to a conformational transition. A relationship between the local readjustments and the equilibrium fluctuations of a side chain around its unbound conformation is suggested. Most of the side chains undergo larger changes in the dihedral angle most distant from the backbone. The frequencies of the core-to-surface interface transitions of six nonpolar residues and Tyr are larger than the frequencies of the opposite surface-to-core transitions. The binding increases both polar and nonpolar interface areas. However, the increase of the nonpolar area is larger for all considered classes of protein complexes, suggesting that the protein association perturbs the unbound interfaces to increase the hydrophobic contribution to the binding free energy. To test modeling approaches to side-chain flexibility in protein docking, conformational changes in the X-ray set were compared with those in the docking decoy sets. The results lead to a better understanding of the conformational changes in proteins and suggest directions for efficient conformational sampling in docking protocols.     

三肽和四肽构象空间的可视化方法

研究蛋白质寡肽构象在构象空间中的分布情况,对提取寡肽模式并构建短肽库具有重要意义。通过构建一个保距映射,将以主链原子均方根距离(root mean square distance,RMSD)为距离测度的三肽构象空间变换为一维直线上的欧氏距离空间,从而直观地展现三肽构象的聚集情况,表明三肽主链构象可以用单一变量编码。应用该特性对四肽的构象空间加以分析,将四肽构象映射到三维空间中,从而以可视的方式描述四肽构象空间的聚集情况。对短肽构象空间的初步分析表明,短肽的聚集性和二级结构有着密切的联系。在四肽构象空间中存在有自然边界的离散区域(与螺旋等结构相关),也有一些区域(与折叠等结构有关)难以进一步划分。这种方法也为以可视方式分析高维空间中肽段的聚集性给出了一种可能的方案。     

Virus membrane-fusion proteins: more than one way to make a hairpin

Structure-function studies have defined two classes of viral membrane-fusion proteins that have radically different architectures but adopt a similar overall 'hairpin' conformation to induce fusion of the viral and cellular membranes and therefore initiate infection. In both classes, the hairpin conformation is achieved after a conformational change is triggered by interaction with the target cell. This review will focus in particular on the properties of the more recently described class II proteins.     

SIMS: A Hybrid Method for Rapid Conformational Analysis

Proteins are at the root of many biological functions, often performing complex tasks as the result of large changes in their structure. Describing the exact details of these conformational changes, however, remains a central challenge for computational biology due the enormous computational requirements of the problem. This has engendered the development of a rich variety of useful methods designed to answer specific questions at different levels of spatial, temporal, and energetic resolution. These methods fall largely into two classes: physically accurate, but computationally demanding methods and fast, approximate methods. We introduce here a new hybrid modeling tool, the Structured Intuitive Move Selector (sims), designed to bridge the divide between these two classes, while allowing the benefits of both to be seamlessly integrated into a single framework. This is achieved by applying a modern motion planning algorithm, borrowed from the field of robotics, in tandem with a well-established protein modeling library. sims can combine precise energy calculations with approximate or specialized conformational sampling routines to produce rapid, yet accurate, analysis of the large-scale conformational variability of protein systems. Several key advancements are shown, including the abstract use of generically defined moves (conformational sampling methods) and an expansive probabilistic conformational exploration. We present three example problems that sims is applied to and demonstrate a rapid solution for each. These include the automatic determination of “active” residues for the hinge-based system Cyanovirin-N, exploring conformational changes involving long-range coordinated motion between non-sequential residues in Ribose-Binding Protein, and the rapid discovery of a transient conformational state of Maltose-Binding Protein, previously only determined by Molecular Dynamics. For all cases we provide energetic validations using well-established energy fields, demonstrating this framework as a fast and accurate tool for the analysis of a wide range of protein flexibility problems.     

Additivity and non-additivity effects in the conformational changes of 1,2-ethanediol diformatein vacuo and in solution.Ab-initio and molecular mechanics descriptions

The performances of MM2 andab-initio SCF STO-3G calculations to describe conformational changes in (CH₂OCHO)₂ in vacuo and in solution are examined. We present and justify a simple procedure to add solvent effects to the MM2 conformational energies. The analysis is focused on the detection of non-additive effects in simultaneous conformational changes. We have found that the -CH₂-CH₂-group is generally effective in decoupling simultaneous conformational changes. The non-additive effects, actually present in specific classes of conformational changes, are reduced by a solvent contribution when full SCF calculations are employed, and emphasized when MM2 + solvent calculations are employed. A rationale of this finding is presented, and some suggestions are made for the correction of this artifact, due to the use of rigid charges.     

HIV coreceptors: role of structure,posttranslational modifications,and internalization in viral-cell fusion and as targets for entry inhibitors

The human immunodeficiency virus (HIV) envelope glycoprotein forms trimers on the virion surface, with each monomer consisting of two subunits, gp120 and gp41. The gp120 envelope component binds to CD4 on target cells and undergoes conformational changes that allow gp120 to interact with certain G-protein-coupled receptors (GPCRs) on the same target membranes. The GPCRs that function as HIV coreceptors were found to be chemokine receptors. The primary coreceptors are CCR5 and CXCR4, but several other chemokine receptors were identified as "minor coreceptors", indicating their ability support entry of some HIV strains in tissue cultures. Formation of the tri-molecular complexes stabilizes virus binding and triggers a series of conformational changes in gp41 that facilitate membrane fusion and viral cell entry. Concerted efforts are underway to decipher the specific interactions between gp120/CD4, gp120/coreceptors, and their contributions to the subsequent membrane fusion process. It is hoped that some of the transient conformational intermediates in gp120 and gp41 would serve as targets for entry inhibitors. In addition, the CD4 and coreceptors are primary targets for several classes of inhibitors currently under testing. Our review summarizes the current knowledge on the interactions of HIV gp120 with its receptor and coreceptors, and the important properties of the chemokine receptors and their regulation in primary target cells. We also summarize the classes of coreceptor inhibitors under development.     

Molecular orbital studies of the conformation around the 5 position in angiotensin II

Analogs of angiotensin are more potent when the side chain at position 5 is branched rather than unbranched. We have performed molecular orbital calculations of conformational preferences of l-valine (branched side chain) and l-a-aminobutyric acid (unbranched side chain) as amino acid residues. The results of these calculations illustrate the differing conformational preferences of the two classes of amino acids and support a stereochemical role for position 5 of angiotensin.     

Effects of calcium binding on the side-chain methyl dynamics of calbindin D9k: a 2H NMR relaxation study

The effects of Ca(2+) binding on the side-chain methyl dynamics of calbindin D(9k) have been characterized by (2)H NMR relaxation rate measurements. Longitudinal, transverse in-phase, quadrupolar order, transverse anti-phase and double quantum relaxation rates are reported for both the apo and Ca(2+)-loaded states of the protein at two magnetic field strengths. The relatively large size of the data set allows for a detailed analysis of the underlying conformational dynamics by spectral density mapping and model-free fitting procedures. The results reveal a correlation between a methyl group's distance from the Ca(2+) binding sites and its conformational dynamics. Several methyl groups segregate into two limiting classes, one proximal and the other distal to the binding sites. Methyl groups in these two classes respond differently to Ca(2+) binding, both in terms of the timescale and amplitude of their fluctuations. Ca(2+) binding elicits a partial immobilization among methyl groups in the proximal class, which is consistent with previous studies of calbindin's backbone dynamics. The distal class, however, exhibits a trend that could not be inferred from the backbone data in that its mobility actually increases with Ca(2+) binding. We have introduced the term polar dynamics to describe this type of organization across the molecule. The trend may represent an important mechanism by which calbindin D(9k) achieves high affinity binding while minimizing the corresponding loss of conformational entropy.     

Conformational changes in DNA‐binding proteins: Relationships with precomplex features and contributions to specificity and stability

Both Proteins and DNA undergo conformational changes in order to form functional complexes and also to facilitate interactions with other molecules. These changes have direct implications for the stability and specificity of the complex, as well as the cooperativity of interactions between multiple entities. In this work, we have extensively analyzed conformational changes in DNA‐binding proteins by superimposing DNA‐bound and unbound pairs of protein structures in a curated database of 90 proteins. We manually examined each of these pairs, unified the authors' annotations, and summarized our observations by classifying conformational changes into six structural categories. We explored a relationship between conformational changes and functional classes, binding motifs, target specificity, biophysical features of unbound proteins, and stability of the complex. In addition, we have also investigated the degree to which the intrinsic flexibility can explain conformational changes in a subset of 52 proteins with high quality coordinate data. Our results indicate that conformational changes in DNA‐binding proteins contribute significantly to both the stability of the complex and the specificity of targets recognized by them. We also conclude that most conformational changes occur in proteins interacting with specific DNA targets, even though unbound protein structures may have sufficient information to interact with DNA in a nonspecific manner. Proteins 2014; 82:841–857. © 2013 Wiley Periodicals, Inc.     

Conformational analysis and clustering of short and medium size loops connecting regular secondary structures: a database for modeling and prediction.

Loops are regions of nonrepetitive conformation connecting regular secondary structures. We identified 2,024 loops of one to eight residues in length, with acceptable main-chain bond lengths and peptide bond angles, from a database of 223 protein and protein-domain structures. Each loop is characterized by its sequence, main-chain conformation, and relative disposition of its bounding secondary structures as described by the separation between the tips of their axes and the angle between them. Loops, grouped according to their length and type of their bounding secondary structures, were superposed and clustered into 161 conformational classes, corresponding to 63% of all loops. Of these, 109 (51% of the loops) were populated by at least four nonhomologous loops or four loops sharing a low sequence identity. Another 52 classes, including 12% of the loops, were populated by at least three loops of low sequence similarity from three or fewer nonhomologous groups. Loop class suprafamilies resulting from variations in the termini of secondary structures are discussed in this article. Most previously described loop conformations were found among the classes. New classes included a 2:4 type IV hairpin, a helix-capping loop, and a loop that mediates dinucleotide-binding. The relative disposition of bounding secondary structures varies among loop classes, with some classes such as beta-hairpins being very restrictive. For each class, sequence preferences as key residues were identified; those most frequently at these conserved positions than in proteins were Gly, Asp, Pro, Phe, and Cys. Most of these residues are involved in stabilizing loop conformation, often through a positive phi conformation or secondary structure capping. Identification of helix-capping residues and beta-breakers among the highly conserved positions supported our decision to group loops according to their bounding secondary structures. Several of the identified loop classes were associated with specific functions, and all of the member loops had the same function; key residues were conserved for this purpose, as is the case for the parvalbumin-like calcium-binding loops. A significant number, but not all, of the member loops of other loop classes had the same function, as is the case for the helix-turn-helix DNA-binding loops. This article provides a systematic and coherent conformational classification of loops, covering a broad range of lengths and all four combinations of bounding secondary structure types, and supplies a useful basis for modelling of loop conformations where the bounding secondary structures are known or reliably predicted.     

Clustering of protein structural fragments reveals modular building block approach of nature

Structures of peptide fragments drawn from a protein can potentially occupy a vast conformational continuum. We co-ordinatize this conformational space with the help of geometric invariants and demonstrate that the peptide conformations of the currently available protein structures are heavily biased in favor of a finite number of conformational types or structural building blocks. This is achieved by representing a peptides' backbone structure with geometric invariants and then clustering peptides based on closeness of the geometric invariants. This results in 12,903 clusters, of which 2207 are made up of peptides drawn from functionally and/or structurally related proteins. These are termed "functional" clusters and provide clues about potential functional sites. The rest of the clusters, including the largest few, are made up of peptides drawn from unrelated proteins and are termed "structural" clusters. The largest clusters are of regular secondary structures such as helices and beta strands as well as of beta hairpins. Several categories of helices and strands are discovered based on geometric differences. In addition to the known classes of loops, we discover several new classes, which will be useful in protein structure modeling. Our algorithm does not require assignment of secondary structure and, therefore, overcomes the limitations in loop classification due to ambiguity in secondary structure assignment at loop boundaries.     

HIV coreceptors: role of structure, posttranslational modifications, and internalization in viral-cell fusion and as targets for entry inhibitors

The human immunodeficiency virus (HIV) envelope glycoprotein forms trimers on the virion surface, with each monomer consisting of two subunits, gp120 and gp41. The gp120 envelope component binds to CD4 on target cells and undergoes conformational changes that allow gp120 to interact with certain G-protein-coupled receptors (GPCRs) on the same target membranes. The GPCRs that function as HIV coreceptors were found to be chemokine receptors. The primary coreceptors are CCR5 and CXCR4, but several other chemokine receptors were identified as “minor coreceptors”, indicating their ability support entry of some HIV strains in tissue cultures. Formation of the tri-molecular complexes stabilizes virus binding and triggers a series of conformational changes in gp41 that facilitate membrane fusion and viral cell entry. Concerted efforts are underway to decipher the specific interactions between gp120/CD4, gp120/coreceptors, and their contributions to the subsequent membrane fusion process. It is hoped that some of the transient conformational intermediates in gp120 and gp41 would serve as targets for entry inhibitors. In addition, the CD4 and coreceptors are primary targets for several classes of inhibitors currently under testing. Our review summarizes the current knowledge on the interactions of HIV gp120 with its receptor and coreceptors, and the important properties of the chemokine receptors and their regulation in primary target cells. We also summarize the classes of coreceptor inhibitors under development.