Spatial organization of the N-terminal beta-endorphin tridecapeptide

Using conformational analysis spatial structure and conformational properties of the N-terminal tridecapeptide--endorphine molecules were investigated. Calculations were based on the fragmental analysis using non-valent, electrostatic, torsional interactions and hydrogen bonds. It was shown that tridecapeptide could exist in several low-energetical conformational states. Enkephaline fragments structure depends on the most perspective structure of free metioninenkephaline. The results can be used for conformational analysis of endorphine molecules, for structure--function relations study.     

Ligand-induced conformational change of a protein reproduced by a linear combination of displacement vectors obtained from normal mode analysis

The conformational change of a protein upon ligand binding was examined by normal mode analysis (NMA) based on an elastic-network model (ENM) for a full-atom system using dihedral angles as independent variables. Specifically, we investigated the extent to which conformational change vectors of atoms from an apo form to a holo form of a protein can be represented by a linear combination of the displacement vectors of atoms in the apo form calculated for the lowest-frequency m normal modes (m=1, 2,…, 20). In this analysis, the latter vectors were best fitted to the former ones by the least-squares method. Twenty-two paired proteins in the holo and apo forms, including three dimer pairs, were examined. The results showed that, in most cases, the conformational change vectors were reproduced well by a linear combination of the displacement vectors of a small number of low-frequency normal modes. The conformational change around an active site was reproduced as well as the entire conformational change, except for some proteins that only undergo significant conformational changes around active sites. The weighting factors for 20 normal modes optimized by the least-squares fitting characterize the conformational changes upon ligand binding for these proteins. The conformational changes sampled around the apo form of a protein by the linear combination of the displacement vectors obtained by ENM-based NMA may help solve the flexible-docking problem of a protein with another molecule because the results presented herein suggest that they have a relatively high probability of being involved in an actual conformational change.     

Theoretical conformation analysis of beta-casomorphin-5 molecule

The spatial structure of beta-casomorphin-5 molecule--H-Tyr-Pro-Phe-Pro-Gly-OH has been investigated by the theoretical conformational analysis method. The results indicate that the conformational properties of the molecule can be represented by a large number of structures having twelve different backbone forms.     

Understanding biomolecular motion,recognition, and allostery by use of conformational ensembles

We review the role conformational ensembles can play in the analysis of biomolecular dynamics, molecular recognition, and allostery. We introduce currently available methods for generating ensembles of biomolecules and illustrate their application with relevant examples from the literature. We show how, for binding, conformational ensembles provide a way of distinguishing the competing models of induced fit and conformational selection. For allostery we review the classic models and show how conformational ensembles can play a role in unravelling the intricate pathways of communication that enable allostery to occur. Finally, we discuss the limitations of conformational ensembles and highlight some potential applications for the future.     

Trajectory analysis of NMR structure calculations

Summary NMR as well as X-ray crystallography are used to determine the three-dimensional structures of macromolecules at atomic resolution. Structure calculation generates coordinates that are compatible with NMR data from randomly generated initial structures. We analyzed the trajectory taken by structures during NMR structure calculation in conformational space, assuming that the distance between two structures in conformational space is the root-mean-square deviation between the two structures. The coordinates of a structure in conformational space were obtained by applying the metric multidimensional scaling method. As an example, we used a 22-residue peptide, -Conotoxin GIIIA, and a simulated annealing protocol of XPLOR. We found that the three-dimensional solution of the multidimensional scaling analysis is sufficient to describe the overall configuration of the trajectories in conformational space. By comparing the trajectories of the entire calculation with those of the converged calculation, random sampling of conformational space is readily discernible. Trajectory analysis can also be used for optimization of protocols of NMR structure calculation, by examining individual trajectories.Abbreviations MD molecular dynamics - MDS multidimensional scaling - rmsd root-mean-square deviation - armsd angular rmsd - R multiple correlation coefficient - YASAP yet another simulated annealing protocol - PCA principal component analysis     

Characterization of conformational epitope of alginate-derived polymannuronates by surface plasmon resonance

In the present study, we developed a mAb to alginate-derived polymannuronates (ADPM) and examined the antigenic epitopes using surface plasmon resonance (SPR) in conjunction with a large panel of oligomannuronate probes. We found that tetrasaccharide is the minimum-binding unit, and that increases in chain length from the tetrasaccharide to the heptsaccharide further enhance monovalent binding. A sharp increase in affinity was observed when increasing from the octasaccharide to the cosasaccharide, which is due to a further enhancement of the individual antigenic epitope combined with multivalency. Kinetic binding studies further suggested that the conformational epitope is discontinuous and infrequent and that a C6-carboxyl group is important in maintaining the conformational epitope. Moreover, CD analysis revealed there were conformational structures in epitopes. The data support our hypothesis that the conformational epitope for the mAb may be an extended helical segment of ADPM. ADPM exists mainly in linear form, but it can infrequently and spontaneously form extended helices. Although helices are not prevalent in ADPM, the immune system preferentially selects these conformational epitopes because they are unique. Together, our results indicate that the antigenic epitopes in beta-d-mannuronates are conformational and require C6-carboxyl groups.     

Quantitative analysis of cyclic beta-turn models.

The beta-turn is a frequently found structural unit in the conformation of globular proteins. Although the circular dichroism (CD) spectra of the alpha-helix and beta-pleated sheet are well defined, there remains some ambiguity concerning the pure component CD spectra of the different types of beta-turns. Recently, it has been reported (Hollósi, M., Kövér, K.E., Holly, S., Radics, L., & Fasman, G.D., 1987, Biopolymers 26, 1527-1572; Perczel, A., Hollósi, M., Foxman, B.M., & Fasman, G.D., 1991a, J. Am. Chem. Soc. 113, 9772-9784) that some pseudohexapeptides (e.g., the cyclo[(delta)Ava-Gly-Pro-Aaa-Gly] where Aaa = Ser, Ser(OtBu), or Gly) in many solvents adopt a conformational mixture of type I and the type II beta-turns, although the X-ray-determined conformation was an ideal type I beta-turn. In addition to these pseudohexapeptides, conformational analysis was also carried out on three pseudotetrapeptides and three pseudooctapeptides. The target of the conformation analysis reported herein was to determine whether the ring stress of the above beta-turn models has an influence on their conformational properties. Quantitative nuclear Overhauser effect (NOE) measurements yielded interproton distances. The conformational average distances so obtained were interpreted utilizing molecular dynamics (MD) simulations to yield the conformational percentages. These conformational ratios were correlated with the conformational weights obtained by quantitative CD analysis of the same compounds. The pure component CD curves of type I and type II beta-turns were also obtained, using a recently developed algorithm (Perczel, A., Tusnády, G., Hollósi, M., & Fasman, G.D., 1991b, Protein Eng. 4(6), 669-679). For the first time the results of a CD deconvolution, based on the CD spectra of 14 beta-turn models, were assigned by quantitative NOE results. The NOE experiments confirmed the ratios of the component curves found for the two major beta-turns by CD analysis. These results can now be used to enhance the conformational determination of globular proteins on the basis of their CD spectra.     

An analysis of side-chain conformation in proteins.

The crystal structures of a number of globular proteins are currently available. An analysis of the distribution of side-chains among different allowed conformations in these proteins has been carried out. The observed conformations of individual residues are discussed on the basis of well-known stereochemical criteria. The population distribution of side-chains in different allowed regions in conformational space can be explained largely on the basis of simple steric considerations. In addition to examining the conformational behaviour of individual residues, some population distributions of conformational angles of general interest involving groups of residues have also been analyzed.     

Partition decoupling for multi-gene analysis of gene expression profiling data

Background

Molecular dynamics (MD) simulations are powerful tools to investigate the conformational dynamics of proteins that is often a critical element of their function. Identification of functionally relevant conformations is generally done clustering the large ensemble of structures that are generated. Recently, Self-Organising Maps (SOMs) were reported performing more accurately and providing more consistent results than traditional clustering algorithms in various data mining problems. We present a novel strategy to analyse and compare conformational ensembles of protein domains using a two-level approach that combines SOMs and hierarchical clustering.

Results

The conformational dynamics of the α-spectrin SH3 protein domain and six single mutants were analysed by MD simulations. The Cα's Cartesian coordinates of conformations sampled in the essential space were used as input data vectors for SOM training, then complete linkage clustering was performed on the SOM prototype vectors. A specific protocol to optimize a SOM for structural ensembles was proposed: the optimal SOM was selected by means of a Taguchi experimental design plan applied to different data sets, and the optimal sampling rate of the MD trajectory was selected. The proposed two-level approach was applied to single trajectories of the SH3 domain independently as well as to groups of them at the same time. The results demonstrated the potential of this approach in the analysis of large ensembles of molecular structures: the possibility of producing a topological mapping of the conformational space in a simple 2D visualisation, as well as of effectively highlighting differences in the conformational dynamics directly related to biological functions.

Conclusions

The use of a two-level approach combining SOMs and hierarchical clustering for conformational analysis of structural ensembles of proteins was proposed. It can easily be extended to other study cases and to conformational ensembles from other sources.     

A model of threshold behavior reveals rescue mechanisms of bystander proteins in conformational diseases

Conformational diseases result from the failure of a specific protein to fold into its correct functional state. The misfolded proteins can lead to the toxic aggregation of proteins. Protein misfolding in conformational diseases often displays a threshold behavior characterized by a sudden shift between nontoxic to toxic levels of misfolded proteins. In some conformational diseases, evidence suggests that misfolded proteins interact with bystander proteins (unfolded and native folded proteins), eliciting a misfolded phenotype. These bystander isomers would follow their normal physiological pathways in absence of misfolded proteins. In this article, we present a general mechanism of bystander and misfolded protein interaction which we have used to investigate how the threshold behavior in protein misfolding is triggered in conformational diseases. Using a continuous flow reactor model of the endoplasmic reticulum, we found that slight changes in the bystander protein residence time in the endoplasmic reticulum or the ratio of basal misfolded to bystander protein inflow rates can trigger the threshold behavior in protein misfolding. Our analysis reveals three mechanisms to rescue bystander proteins in conformational diseases. The results of our model can now help direct experiments to understand the threshold behavior and develop therapeutic strategies targeting the modulation of conformational diseases.     

Effect of colloidal gold size on the conformational changes of adsorbed cytochrome c: probing by circular dichroism, UV-visible, and infrared spectroscopy

The conformational changes of bovine heart cytochrome c (cyt c) induced by the adsorption on gold nanoparticles with different sizes have been investigated by electronic absorption, circular dichroism (CD), and Fourier transform infrared spectra. The combination of these techniques can give complementary information about adsorption-induced conformational changes. The results show that there are different conformational changes for cyt c adsorbed on gold nanoparticles with different sizes due to the different interaction forces between cyt c and gold nanoparticles. The colloidal gold concentration-dependent conformation distribution curves of cyt c obtained by analysis of CD spectra using the singular value decomposition least-squares method show that the coverage of cyt c on the gold nanoparticles surface also affects the conformational changes of the adsorbed cyt c.     

Calculation of the structure of the serotonin-antibiotic lasalocid A complex by conformation analysis and the Monte Carlo method

The conformational aspects of interaction of the antibiotic X537A at complexation with serotonin, hydration of molecules and their complex were studied by the methods of theoretical conformational analysis and Monte-Carlo.     

Tetracycline and its analogues as inhibitors of amyloid fibrils: searching for a geometrical pharmacophore by theoretical investigation of their conformational behavior in aqueous solution

Tetracycline (TC) and its derivatives have recently been proposed as a new class of antagonists in prion diseases as they prevent the aggregation of prion protein peptides and their acquisition of protease resistance in vitro and in vivo. Looking for relationships between conformational flexibility and biological activity, we searched for a geometrical pharmacophore by investigating, in aqueous solution, the conformational behavior of 15 TCs in both the zwitterionic and the anionic forms. For TC similar conformational flexibility was found for the two forms and two main conformational families were detected, an extended and a folded conformation characterized by different intramolecular hydrogen-bond networks. On comparing the Molecular Mechanics results with the ab initio ones and the experimental evidence, it can be seen that the conformational behavior of TC is reasonably well predicted by the MM2 force field, whereas the conformational energies provided by the Amber force field are unreliable. The conformational analysis of the other TC derivatives was then performed by the MM2 force field. As a result, their conformational behavior was similar to that observed for TC itself. Despite the hydronaphthacene moiety's conformational flexibility, no geometrical pharmacophore was found among the TCs, i.e. properties other than geometrical ones should play a crucial role in determining their anti-fibrillogenic ability.     

Molecular mechanisms of activation in CDK2

Cyclin-dependent kinases (CDKs) are enzymes involved in crucial cellular processes. Their biological activity is directly linked to their high conformational variability, which involves large protein conformational rearrangements. We present here the application of an enhancing sampling technique to the study of conformational transitions between the open and closed state of CDKs. The analysis of the conformational intermediates supports the idea that the process is regulated by two important protein regions, which sequentially rearrange in order to allow the protein to reach its final conformation. Furthermore, the two paths involve additional (minor) protein rearrangements which are specific to the paths. Our results show that our procedure can provide reasonable transition pathways between the two protein forms at a very reduced computational cost. The robustness and the simplicity of our approach make it of general application to describe virtually any macromolecular conformational transitions.     

Effect of an induced conformational change on the physical properties of two chemotactic receptor molecules.

The physical properties and conformational dynamics of the Salmonella typhimurium ribose and galactose receptors have been examined. Studies involving circular dichroism, fluorescence, absorption spectroscopy, and sedimentation analysis show that the two receptor proteins have different morphologies and exhibit diverse responses to sugar binding. The ribose receptor lacks both tryptophan and disulfide residues, and the galactose receptor lacks disulfides and has only a single tryptophan residue. By virtue of these fortuitous properties, the conformational changes induced in these proteins by sugar binding can be dissected by utilizing a variety of physical probes. A ligand-induced conformational change in the ribose receptor is shown by circular dichroism and fluorescence spectroscopy, which reveal spectral changes assignable to tyrosine, phenylalanine, and methionine residues. A conformational change in the galactose receptor has been demonstrated by fluorescence spectroscopy involving the distant reporter group method, which shows changes assignable to tryptophan and methionine sites and which is corroborated by sedimentation analysis. It is clear that there are extensive conformational changes in the two receptor proteins and that the different physical methods provide complementary information on the nature of these changes.     

Structure-function organization of neurokinin A and neurokinin B molecules. III. A conformational study of glycine-monosubstituted analogues of neurokinins A and B

The conformational features of some glycine-monosubstituted analogues of neurokinins A and B were investigated by the method of theoretical conformational analysis. The calculated geometry and energy parameters permitted one to determine the structural role of each of these substituted amino acids in the mechanism of folding of the low-energy conformational states of neuropeptides. On the basis of the calculated data and the results of biological tests of these analogues, the structure-function relationships of neurokinins A and B were discussed.     

Conformational properties of the CalliFMRF-amide series neuropeptides

Conformational properties of five neuropeptides belonging to the calliFMRF-amide series with the Xaa-Pro-Yaa-Gln-Asp-Phe-Met-Arg-Phe-NH2 homologous sequences were studied by the method of theoretical conformational analysis. Three members of these group (1) (Xaa = Thr, Yaa = Gln), (2) (Xaa = Thr, Yaa = Ser), and (3) (Xaa = Yaa = Ser) can stimulate the saliva secretion from the separated salivary gland of the Calliphora vomitoria fly, whereas two other calliFMRF-amides (4) (Xaa = Lys, Yaa = Asn) and (5) (Xaa = Ala, Yaa = Gly) are inactive in this biological test. Low-energy spatial structures of the studied compounds were determined by a conformational analysis. A comparison of the stable structures of the biologically active and inactive neuropeptides revealed a similarity in their conformational properties and allowed determination of the role of separate residues in the peptide folding. The calculations demonstrated that the C-terminal hexapeptide fragment identical in all the five peptides tends to form alpha-helical structure, whereas the variable N-terminal tripeptide regions of CalliFMRF-amides (1)-(5) form more conformationally flexible structures.     

Passive transfer of collagen arthritis: heterogeneity of anti-collagen IgG

Using a combination of affinity chromatography procedures, anticollagen IgG was fractionated into three distinct populations. One population reacted only to conformational determinants, another population reacted only to structural determinants, and the third population reacted to both conformation and structural determinants. When these populations were tested for their arthritogenicity, only those fractions that reacted to conformational and to conformational and structural determinants were active in inducing clinical arthritis. Immunofluorescence analysis of the hind paw of recipient rats indicated that antibodies directed only to conformational and to both conformational and structural determinants bound to articular cartilage and activated the complement system. Antibodies directed strictly to structural determinants did not bind to articular cartilage and were nonarthritogenic.     

Empirical Bayes Methods Enable Advanced Population-Level Analyses of Single-Molecule FRET Experiments

Many single-molecule experiments aim to characterize biomolecular processes in terms of kinetic models that specify the rates of transition between conformational states of the biomolecule. Estimation of these rates often requires analysis of a population of molecules, in which the conformational trajectory of each molecule is represented by a noisy, time-dependent signal trajectory. Although hidden Markov models (HMMs) may be used to infer the conformational trajectories of individual molecules, estimating a consensus kinetic model from the population of inferred conformational trajectories remains a statistically difficult task, as inferred parameters vary widely within a population. Here, we demonstrate how a recently developed empirical Bayesian method for HMMs can be extended to enable a more automated and statistically principled approach to two widely occurring tasks in the analysis of single-molecule fluorescence resonance energy transfer (smFRET) experiments: 1), the characterization of changes in rates across a series of experiments performed under variable conditions; and 2), the detection of degenerate states that exhibit the same FRET efficiency but differ in their rates of transition. We apply this newly developed methodology to two studies of the bacterial ribosome, each exemplary of one of these two analysis tasks. We conclude with a discussion of model-selection techniques for determination of the appropriate number of conformational states. The code used to perform this analysis and a basic graphical user interface front end are available as open source software.