PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “conformational ensemble”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9Linked to original sources

Minireview: vitamin D receptor: new assignments for an already busy receptor.

With its discovery in 1920, the molecule vitamin D achieved prominence as a nutritionally essential vitamin important for calcium homeostasis, particularly in the intestine and bone. Then in 1932, the elucidation of vitamin D's chemical structure revealed that this vitamin was in fact a steroid. But it was not until the late 1960s that it was appreciated that the steroid vitamin D was a precursor of a new steroid hormone, 1alpha,25(OH)2-vitamin D3 [1alpha,25(OH)2D3], that is produced by the kidney acting as an endocrine gland. The discovery in 1969 of the nuclear vitamin D receptor (VDR) for 1alpha,25(OH)2D3 initiated a two-decade-long proliferation of reports that collectively described the broad sphere of influence of the vitamin D endocrine system that is defined by the presence of the VDR in over 30 tissue/organs of man. The new genomic frontiers defined by the cellular presence of the VDR include the immune system's B and T lymphocytes, hair follicle, muscle, adipose tissue, bone marrow, and cancer cells. Unexpectedly in the mid 1980s, a new world of 1alpha,25(OH)2D3-mediated rapid responses (RR) was discovered. These were responses that occurred too rapidly (minutes to an hour) to be explained as the simple consequence of the nuclear VDR regulating gene transcription. Some RR examples include the rapid intestinal absorption of calcium (transcaltachia), secretion of insulin by pancreatic beta-cells, opening of voltage-gated Ca2+ and Cl- channels in osteoblasts, and the rapid migration of endothelial cells. The question then arose as to whether there was a second receptor, apart from the nuclear VDR, which responded to the presence of 1alpha,25(OH)2D3 to generate RR? After some false starts, it now appears that the classic VDR, long known to reside in the cell nucleus, in some cells is also associated with caveolae present in the plasma membrane. Furthermore, the chemical properties of the conformationally flexible 1alpha,25(OH)2D3 allow it to generate different shaped ligands for the VDR that are selective either for genomic or for RR. This minireview summarizes a proposed conformational ensemble model of the VDR that provides insight into how different ligand shapes of 1alpha,25(OH)2D3 acting through the VDR in different cellular locations can selectively mediate both genomic and RR.

Animals↗

Structural evidence that propofol stabilizes different GABA(A) receptor states at potentiating and activating concentrations.

The GABA(A) receptor is a target of many general anesthetics, such as propofol. General anesthetic binding sites are distinct from the GABA binding sites. At low concentrations, the anesthetics potentiate the currents induced by submaximal GABA concentrations. At higher concentrations the anesthetics directly activate GABA(A) receptors. In contrast, benzodiazepines, such as diazepam, only potentiate currents induced by submaximal GABA concentrations. Channel kinetic studies suggest that these drugs stabilize different receptor states. We previously showed that the accessibility of the anionic sulfhydryl reagent p-chloromercuribenzenesulfonate (pCMBS(-)) applied extracellularly to cysteines substituted for residues in the GABA(A) alpha1 subunit M3 membrane-spanning segment was state-dependent. The subset of pCMBS(-)-accessible, M3 segment cysteine mutants acts as a reporter for receptor conformation. Here we show that pCMBS(-), applied in the presence of a potentiating concentration of propofol, reacts with a subset of alpha1 subunit, M3 segment, cysteine-substitution mutants (Y294C, V297C, I302C, F304C). In the presence of a directly activating concentration of propofol pCMBS(-) reacts with a different subset of the M3 cysteine-substitution mutants (Y294C, S299C, I302C, E303C, A305C). These subsets are distinct from the subsets of M3 cysteine-substitution mutants that are reactive with pCMBS(-) in the absence and presence of GABA and in the presence of diazepam. We hypothesize that distinct subsets of reactive residues represent distinct conformations or ensembles of conformations of the receptor. These results provide structural evidence for at least five distinct receptor states, three nonconducting states, resting, diazepam-bound and potentiating propofol-bound, and two conducting-desensitized states, the activating propofol-bound and GABA-bound states.

4-Chloromercuribenzenesulfonate↗

Interpreting NMR data for beta-peptides using molecular dynamics simulations.

NMR is one of the most used techniques to resolve structure of proteins and peptides in solution. However, inconsistencies may occur due to the fact that a polypeptide may adopt more than one conformation. Since the NOE distance bounds and (3)J-values used in such structure determination represent a nonlinear average over the total ensemble of conformers, imposition of NOE or (3)J-value restraints to obtain one unique conformation is not an appropriate procedure in such cases. Here, we show that unrestrained MD simulation of a solute in solution using a high-quality force field yields a conformational ensemble that is largely compatible with the experimental NMR data on the solute. Four 100 ns MD simulations of two forms of a nine-residue beta-peptide in methanol at two temperatures produced conformational ensembles that were used to interpret the NMR data on this molecule and resolve inconsistencies between the experimental NOEs. The protected and unprotected forms of the beta-peptide adopt predominantly a 12/10-helix in agreement with the qualitative interpretation of the NMR data. However, a particular NOE was not compatible with this helix indicating the presence of other conformations. The simulations showed that 3(14)()-helical structures were present in the ensemble of the unprotected form and that their presence correlates with the fulfillment of the particular NOE. Additionally, all inter-hydrogen distances were calculated to compare NOEs predicted by the simulations to the ones observed experimentally. The MD conformational ensembles allowed for a detailed and consistent interpretation of the experimental data and showed the small but specific conformational differences between the protected and unprotected forms of the peptide.

Computer Simulation↗

Towards in silico lead optimization: scores from ensembles of protein/ligand conformations reliably correlate with biological activity.

Accurately ranking protein/ligand interactions and distinguishing subtle differences between homologous compounds in a virtual focused library in silico is essential in a structure-based drug discovery program. In order to establish a predictive model to design novel inhibitors of dihydrofolate reductase (DHFR) from the parasitic protozoa, Cryptosporidium hominis, we docked a series of 30 DHFR inhibitors with measured inhibition constants against the crystal structure of the protein. By including protein flexibility and averaging the energies of the 25 lowest protein/ligand conformers we obtained more accurate total nonbonded energies from which we calculated a predicted biological activity. The calculated and measured biological activities showed reliable correlations of 72.9%. Additionally, visual analysis of the ensemble of protein/ligand conformations revealed alternative ligand binding pockets in the active site. Using the same principles we then created a homology model of DHFR from Toxoplasma gondii and docked 11 inhibitors. A correlation of 50.2% between docking score and activity validates both the method and the model. The correlations presented here are particularly compelling considering the high structural similarity of the ligands and the fact that we have used structures derived from crystallographic data and homology modeling. These docking principles may be useful in any lead optimization study where accurate ranking of similar compounds is desired.

Algorithms↗

Determination of the populations and structures of multiple conformers in an ensemble from NMR data: multiple-copy refinement of nucleic acid structures using floating weights.

A new algorithm is presented for determination of structural conformers and their populations based on NMR data. Restrained Metropolis Monte Carlo simulations or restrained energy minimizations are performed for several copies of a molecule simultaneously. The calculations are restrained with dipolar relaxation rates derived from measured NOE intensities via complete relaxation matrix analysis. The novel feature of the algorithm is that the weights of individual conformers are determined in every refinement step, by the quadratic programming algorithm, in such a way that the restraint energy is minimized. Its design ensures that the calculated populations of the individual conformers are based only on experimental restraints. Presence of internally inconsistent restraints is the driving force for determination of distinct multiple conformers. The method is applied to various simulated test systems. Conformational calculations on nucleic acids are carried out using generalized helical parameters with the program DNAminiCarlo. From different mixtures of A- and B-DNA, minor fractions as low as 10% could be determined with restrained energy minimization. For B-DNA with three local conformers (C2'-endo, O4'-exo, C3'-endo), the minor O4'-exo conformer could not be reliably determined using NOE data typically measured for DNA. The other two conformers, C2'-endo and C3'-endo, could be reproduced by Metropolis Monte Carlo simulated annealing. The behavior of the algorithm in various situations is analyzed, and a number of refinement protocols are discussed. Prior to application of this algorithm to each experimental system, it is suggested that the presence of internal inconsistencies in experimental data be ascertained. In addition, because the performance of the algorithm depends on the type of conformers involved and experimental data available, it is advisable to carry out test calculations with simulated data modeling each experimental system studied.

Algorithms↗

Folding funnels and binding mechanisms.

The long-held views on lock-and-key versus induced fit in binding arose from the notion that a protein exists in a single, most stable conformation, dictated by its sequence. However, in solution proteins exist in a range of conformations, which may be described by statistical mechanical laws and their populations follow statistical distributions. Upon binding, the equilibrium will shift in favor of the bound conformation from the ensemble of conformations around the bottom of the folding funnel. Hence here we extend the implications and the usefulness of the folding funnel concept to explain fundamental binding mechanisms.

Peptides↗

Comparative evaluation of 11 scoring functions for molecular docking.

Eleven popular scoring functions have been tested on 100 protein-ligand complexes to evaluate their abilities to reproduce experimentally determined structures and binding affinities. They include four scoring functions implemented in the LigFit module in Cerius2 (LigScore, PLP, PMF, and LUDI), four scoring functions implemented in the CScore module in SYBYL (F-Score, G-Score, D-Score, and ChemScore), the scoring function implemented in the AutoDock program, and two stand-alone scoring functions (DrugScore and X-Score). These scoring functions are not tested in the context of a particular docking program. Instead, conformational sampling and scoring are separated into two consecutive steps. First, an exhaustive conformational sampling is performed by using the AutoDock program to generate an ensemble of docked conformations for each ligand molecule. This conformational ensemble is required to cover the entire conformational space as much as possible rather than to focus on a few energy minima. Then, each scoring function is applied to score this conformational ensemble to see if it can identify the experimentally observed conformation from all of the other decoys. Among all of the scoring functions under test, six of them, i.e., PLP, F-Score, LigScore, DrugScore, LUDI, and X-Score, yield success rates higher than the AutoDock scoring function. The success rates of these six scoring functions range from 66% to 76% if using root-mean-square deviation < or =2.0 A as the criterion. Combining any two or three of these six scoring functions into a consensus scoring scheme further improves the success rate to nearly 80% or even higher. However, when applied to reproduce the experimentally determined binding affinities of the 100 protein-ligand complexes, only X-Score, PLP, DrugScore, and G-Score are able to give correlation coefficients over 0.50. All of the 11 scoring functions are further inspected by their abilities to construct a descriptive, funnel-shaped energy surface for protein-ligand complexation. The results indicate that X-Score and DrugScore perform better than the other ones at this aspect.

Algorithms↗

Reconciling observations of sequence-specific conformational propensities with the generic polymeric behavior of denatured proteins.

Radii of gyration of denatured proteins vary with chain length and are insensitive to details of amino acid sequence. Observations of sequence independence in polymeric properties conflict with results from spectroscopic experiments, which suggest the presence of sequence-specific residual structure in denatured states. Can we reconcile the two apparently conflicting sets of observations? To answer this question, we need knowledge of the ensemble of conformations accessible to proteins in good solvents. The excluded-volume limit provides an ideal mimic of polymers in good solvents. Therefore, we attempt to solve the "reconciliation problem" by simulating conformational ensembles accessible to peptides and proteins in the excluded-volume limit. Analysis of these ensembles for a variety of polypeptide sequences leads to results that are consistent with experimental observations of sequence-specific conformational preferences in short peptides and the scaling behavior of polymeric quantities for denatured proteins. Reconciliation in the excluded-volume limit comes about due to a tug of war between two factors, namely, minimization of steric overlap and the competing effects of conformational entropy. Minimization of steric overlap promotes chain stretching and leads to experimentally observed sequence-dependent preferences for locally extended segments such as polyproline II helices, beta-strands, and very short stretches of alpha-helix. Conformational entropy opposes chain stretching, and the calculated persistence length for sequence-dependent conformational preferences is less than five amino acids. This estimate does not vary with amino acid sequence. The short persistence lengths lead directly to experimental observations of generic sequence-independent behavior of radii of gyration for denatured proteins.

Amino Acid Sequence↗

Polyproline II conformation is one of many local conformational states and is not an overall conformation of unfolded peptides and proteins.

The alanine-based peptide Ac-XX(A)7OO-NH2, referred to as XAO (where X, A, and O denote diaminobutyric acid, alanine, and ornithine, respectively), has recently been proposed to possess a well defined polyproline II (P(II)) conformation at low temperatures. Based on the results of extensive NMR and CD investigations combined with theoretical calculations, reported here, we present evidence that, on the contrary, this peptide does not have any significant amount of organized P(II) structure but exists in an ensemble of conformations with a distorted bend in the N- and C-terminal regions. The conformational ensemble was obtained by molecular dynamics/simulated annealing calculations using the amber suite of programs with time-averaged distance and dihedral-angle restraints obtained from rotating-frame nuclear Overhauser effect (ROE) volumes and vicinal coupling constants 3J(HN Eta alpha), respectively. The computed ensemble-averaged radius of gyration Rg (7.4 +/- 1.0) A is in excellent agreement with that measured by small-angle x-ray scattering (SAXS) whereas, if the XAO peptide were in the P(II) conformation, Rg would be 11.6 A. Depending on the pH, peptide concentration, and temperature, the CD spectra of XAO do or do not possess the maximum with positive ellipticity in the 217-nm region, which is characteristic of the P(II) structure, reflecting a shifting conformational equilibrium rather than an all-or-none transition. The "P(II) conformation" should, therefore, be considered as one of the accessible conformational states of individual amino acid residues in peptides and proteins rather than as a structure of most of the chain in the early stage of folding.

Alanine↗

Conformational filtering in polypeptides and proteins.

We present a method for assigning an ensemble of conformational states to each amino acid residue of a sequence. The states are defined as regions in the (phi, psi) map. The procedure is based on the use of conformational filters. In each filter we use a different set of approximations to estimate the probability of conformational states, and retain only the ones whose probability exceeds an acceptance probability. The resulting state assignment is not necessarily unique, but provides information that can be further exploited in searches for the tertiary structure. This conformational filtering approach to the de novo analysis of a sequence has a number of advantages over traditional structure prediction. First, it is possible to select acceptance probabilities such that the true conformational state is retained for up to 87% of residues, while substantially reducing the number of potential conformations. Second, in solution most linear peptides are present as ensembles of rapidly interconverting conformers, and such ensembles can be well predicted by filtering. Third, we can use Markov chains instead of a statistical mechanical (Ising) treatment, and avoid the need for estimating statistical weight matrices valid for the molecule as a whole. Markov models can use local transition matrices that are assumed to be independent of the rest of the chain, and are directly calculated from pairwise data. We show here that the locally identifiable transition matrices are transferable from the crystal structures of proteins to the solution structures of short peptides, and the ensembles of filtered conformations are in good agreement with nuclear magnetic resonance data. When applied to proteins, the filters retain several conformational states for most residues, and provide a measure of conformational variability. Small variability means that the segment is well defined by local interactions alone, and hence is likely to preserve its structure when isolated from the rest of the chain. Conversely, the structure of a segment with above-average conformational variability is likely to be significantly affected by its protein environment.

Amino Acid Sequence↗

Relationship between ion pair geometries and electrostatic strengths in proteins.

The electrostatic free energy contribution of an ion pair in a protein depends on two factors, geometrical orientation of the side-chain charged groups with respect to each other and the structural context of the ion pair in the protein. Conformers in NMR ensembles enable studies of the relationship between geometry and electrostatic strengths of ion pairs, because the protein structural contexts are highly similar across different conformers. We have studied this relationship using a dataset of 22 unique ion pairs in 14 NMR conformer ensembles for 11 nonhomologous proteins. In different NMR conformers, the ion pairs are classified as salt bridges, nitrogen-oxygen (N-O) bridges and longer-range ion pairs on the basis of geometrical criteria. In salt bridges, centroids of the side-chain charged groups and at least a pair of side-chain nitrogen and oxygen atoms of the ion-pairing residues are within a 4 A distance. In N-O bridges, at least a pair of the side-chain nitrogen and oxygen atoms of the ion-pairing residues are within 4 A distance, but the distance between the side-chain charged group centroids is greater than 4 A. In the longer-range ion pairs, the side-chain charged group centroids as well as the side-chain nitrogen and oxygen atoms are more than 4 A apart. Continuum electrostatic calculations indicate that most of the ion pairs have stabilizing electrostatic contributions when their side-chain charged group centroids are within 5 A distance. Hence, most (approximately 92%) of the salt bridges and a majority (68%) of the N-O bridges are stabilizing. Most (approximately 89%) of the destabilizing ion pairs are the longer-range ion pairs. In the NMR conformer ensembles, the electrostatic interaction between side-chain charged groups of the ion-pairing residues is the strongest for salt bridges, considerably weaker for N-O bridges, and the weakest for longer-range ion pairs. These results suggest empirical rules for stabilizing electrostatic interactions in proteins.

Biophysical Phenomena↗

Ab initio construction of polypeptide fragments: efficient generation of accurate, representative ensembles.

We describe a novel method to generate ensembles of conformations of the main-chain atoms [N, C(alpha), C, O, Cbeta] for a sequence of amino acids within the context of a fixed protein framework. Each conformation satisfies fundamental stereo-chemical restraints such as idealized geometry, favorable phi/psi angles, and excluded volume. The ensembles include conformations both near and far from the native structure. Algorithms for effective conformational sampling and constant time overlap detection permit the generation of thousands of distinct conformations in minutes. Unlike previous approaches, our method samples dihedral angles from fine-grained phi/psi state sets, which we demonstrate is superior to exhaustive enumeration from coarse phi/psi sets. Applied to a large set of loop structures, our method samples consistently near-native conformations, averaging 0.4, 1.1, and 2.2 A main-chain root-mean-square deviations for four, eight, and twelve residue long loops, respectively. The ensembles make ideal decoy sets to assess the discriminatory power of a selection method. Using these decoy sets, we conclude that quality of anchor geometry cannot reliably identify near-native conformations, though the selection results are comparable to previous loop prediction methods. In a subsequent study (de Bakker et al.: Proteins 2003;51:21-40), we demonstrate that the AMBER forcefield with the Generalized Born solvation model identifies near-native conformations significantly better than previous methods.

Algorithms↗

NMR investigations of protein-carbohydrate interactions: insights into the topology of the bound conformation of a lactose isomer and beta-galactosyl xyloses to mistletoe lectin and galectin-1.

A hallmark of oligosaccharides is their often limited spatial flexibility, allowing them to access a distinct set of conformers in solution. Viewing each individual or even the complete ensemble of conformations as potential binding partner(s) for lectins in protein-carbohydrate interactions, it is pertinent to address the question on the characteristics of bound state conformation(s) in solution. Also, it is possible that entering the lectin's binding site distorts the low-energy topology of a glycosidic linkage. As a step to delineate the strategy of ligand selection for galactosides, a common physiological docking point, we have performed a NMR study on two non-homologous lectins showing identical monosaccharide specificity. Thus, the conformation of lactose analogues bound to bovine heart galectin-1 and to mistletoe lectin in solution has been determined by transferred nuclear Overhauser effect measurements. It is demonstrated that the lectins select the syn conformation of lactose and various structural analogues (Galbeta(1-->4)Xyl, Galbeta(1-->3)Xyl, Galbeta(1-->2)Xyl, and Galbeta(1-->3)Glc) from the ensemble of presented conformations. No evidence for conformational distortion was obtained. Docking of the analogues to the modeled binding sites furnishes explanations, in structural terms, for exclusive recognition of the syn conformer despite the non-homologous design of the binding sites.

Binding Sites↗

Conformational propensities of protein folding intermediates: distribution of species in the 1S, 2S, and 3S ensembles of the [C40A,C95A] mutant of bovine pancreatic ribonuclease A.

A key problem in experimental protein folding is that of characterizing the conformational ensemble of denatured proteins under folding conditions. We address this problem by studying the conformational propensities of reductively unfolded RNase A under folding conditions, since earlier work has indicated that the equilibrium conformational ensemble of fully reduced RNase A resembles the transient conformational ensemble of a burst-phase folding intermediate of disulfide-intact RNase A. To assess these propensities, the relative disulfide-bond populations of the 1S, 2S, and 3S ensembles of the [C40A,C95A] mutant of RNase A were measured. Thirteen of the fifteen possible disulfide bonds are observed, consistent with earlier results and with the rapid reshuffling and lack of stable tertiary structure in these ensembles. This broad distribution contradicts recent observations by another group, but rigorous cross-checks show unambiguously that our data are self-consistent whereas their data are not. The distributions of disulfide bonds in the wild-type and mutant proteins show a power-law dependence on loop length, with an exponent that is significantly smaller than the exponents of either ideal or excluded-volume polymers. The 65-72 disulfide bond is much more strongly favored than would be predicted by this power law, consistent with earlier peptide studies and the disulfide-bond distributions of the 1S and 2S ensembles in wild-type RNase A. Experimental evidence suggests that this preference results from conformational biases in the backbone, rather than from differing accessibilities or reactivities of the two cysteine residues. In general, the other disulfide species do not deviate significantly from the power-law dependence, indicating that the conformational biases are relatively weak.

Alanine↗

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.

Animals↗

Analysis of N-terminal capping using carbonyl-carbon chemical shift measurements.

The model peptide XAAAAEAAARAAAARamide is used to examine the contributions of an N-terminal capping interaction to the conformation and stability of a helical ensemble. The reference peptide has an alanine residue at position X while the capping peptide has a serine residue at this position. The helical ensemble was characterized using circular dichroism measurements and carbonyl-carbon chemical shift measurements of selectively enriched residues. The distribution of helicity within the ensemble of the reference peptide at pH 11 and 0 degrees C appears symmetrical, having a uniform central helix and frayed ends. This distribution is truncated at pH 6 by the repulsive electrostatic interaction between the positively charged alpha-amino group and the positively charged end of the helical macrodipole. The capping peptide forms a side-chain/ main-chain hydrogen bond involving the serine residue and amide of alanine 4. The presence of this hydrogen bond generates a unique motif in the chemical shift profile of its helical ensemble. The conformational stabilization contributed by this hydrogen bond, although cooperatively distributed throughout the helical ensemble, is preferentially focused within the first helical turn. The stabilization provided by this hydrogen bond is able to offset the truncation of the helical ensemble generated by the repulsive electrostatic interaction observed at pH 6.

Amino Acid Sequence↗

OLDERADO: on-line database of ensemble representatives and domains. On Line Database of Ensemble Representatives And DOmains.

In cases where the structure of a single protein is represented by an ensemble of conformations, there is often a need to determine the common features and to choose a "representative" conformation. This occurs, for example, with structures determined by NMR spectroscopy, analysis of the trajectory from a molecular dynamics simulation, or an ensemble of structures produced by comparative modeling. We reported previously automatic methods for (1) defining the atoms with low spatial variance across an ensemble (i.e., the "core" atoms) and the domains in which these atoms lie, and (2) clustering an ensemble into conformationally related subfamilies. To extend the utility of these methods, we have developed a freely available server on the World Wide Web at http:/(/)neon.chem.le.ac.uk/olderado/. This (1) contains an automatically generated database of representative structures, core atoms, and domains determined for 449 ensembles of NMR-derived protein structures in the Protein Data Bank (PDB) in May 1997, and (2) allows the user to upload a PDB-formatted file containing the coordinates of an ensemble of structures. The server returns in real time: (1) information on the residues constituting domains: (2) the structures that constitute each conformational subfamily; and (3) an interactive java-based three-dimensional viewer to visualise the domains and clusters. Such information is useful, for example, when selecting conformations to be used in comparative modeling and when choosing parts of structures to be used in molecular replacement. Here we describe the OLDERADO server.

Computer Communication Networks↗

Photomodulation of conformational states. II. Mono- and bicyclic peptides with (4-aminomethyl)phenylazobenzoic acid as backbone constituent.

It has been reported that backbone cyclization of octapeptides with the photoresponsive (4-aminomethyl)phenylazobenzoic acid imparts sufficient restraints to induce and stabilize ordered conformations of the peptide backbone in both the cis- and trans-azo-isomers (L. Ulysse, J. Cubillos, and J. Chmielewski, Journal of the American Chemical Society, 1995, Vol. 117, pp. 8466-8467). Correspondingly, the active-site octapeptide fragment H-Ala-Cys-Ala-Thr-Cys-Asp-Gly-Phe-OH [134-141] of thioredoxin reductase, with its high preference for a 3(10)-helix turn conformation centered on the Thr-Cys sequence, was backbone cyclized with this azobenzene moiety in the attempt to design a photoresponsive system where the conformational states of the peptide backbone are dictated by the configuration of the azobenzene and can be further modulated by the disulfide bridge. Nuclear magnetic resonance conformational analysis of the monocyclic compound clearly revealed the presence of two conformational families in both the cis- and trans-azo configuration. Of the higher populated conformational families, the structure of the trans-isomer seems like a pretzel-like folding, while the cis-isomer relaxes into a significantly less defined conformational state that does not exhibit any regular structural elements. Further restrictions imparted by disulfide bridging of the peptide moiety leads to an even better defined conformation for the trans-azo-isomer, whereas the cis-isomer can be described as a frustrated system without pronounced energy minima and thus with little conformational preferences. Our findings would suggest that this photoresponsive peptide template may not be of general usefulness for light-induced conformational transitions between two well-defined conformational states at least under the experimental conditions employed, even in the bicyclic form. However, trans --> cis isomerization of the bicyclic peptide is accompanied by a switch from a well-defined conformation to an ensemble of possible conformations.

Aminobenzoates↗