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 271 records · Page 15Linked to original sources

Conformational preferences of proline oligopeptides.

A conformational study on the terminally blocked proline oligopeptides, Ac-(Pro)(n)()-NMe(2) (n = 2-5), is carried out using the ab initio Hartree-Fock level of theory with the self-consistent reaction field method in the gas phase and in solutions (chloroform, 1-propanol, and water) to explore the preference and transition between polyproline II (PPII) and polyproline I (PPI) conformations depending on the chain length, the puckering, and the solvent. The mean differences in the free energy per proline of the up-puckered conformations relative to the down-puckered conformations for both diproline and triproline increases for the PPII-like conformations and decreases for the PPI-like conformations as the solvent polarity increases. These calculated results indicate that the PPII-like structures have preferentially all-down puckerings in solutions, whereas the PPI-like structures have partially mixed puckerings. The free energy difference per proline residue between the PPII- and PPI-like structures decreases as the proline chain becomes longer in the gas phase but increases as the proline chain becomes longer in solutions and the solvent polarity increases. In particular, our calculated results indicate that each of the proline oligopeptides can exist as an ensemble of conformations with the trans and cis peptide bonds in solutions, although the PPII-like structure with all-trans peptide bonds is dominantly preferred, which is reasonably consistent with the previously observed results. In diproline Ac-(Pro)(2)-NMe(2), the rotational barrier to the cis-to-trans isomerization for the first prolyl peptide bond increases as the solvent polarity increases, whereas the rotational barrier for the second prolyl peptide bond does not show the monotonic increase as the solvent polarity increases. When the rotational barriers for these two prolyl peptide bonds were compared, it could be deduced that the conformational transition from PPI with the cis peptide bond to PPII with the trans peptide bond is initiated at the C-terminus and proceeds to the N-terminus in water. This is consistent with the results from NMR experiments on polyproline in D(2)O but opposite to the results from enzymatic hydrolysis kinetics experiments on polyproline.

Chemistry, Physical↗

Discretized torsional dynamics and the folding of an RNA chain.

The aim of this work is to implement a discrete coarse codification of local torsional states of the RNA chain backbone in order to explore the long-time limit dynamics and ultimately obtain a coarse solution to the RNA folding problem. A discrete representation of the soft-mode dynamics is turned into an algorithm for a rough structure prediction. The algorithm itself is inherently parallel, as it evaluates concurrent folding possibilities by pattern recognition, but it may be implemented in a personal computer as a chain of perturbation-translation-renormalization cycles performed on a binary matrix of local topological constraints. This requires suitable representational tools and a periodic quenching of the dynamics for system renormalization. A binary coding of local topological constraints associated with each structural motif is introduced, with each local topological constraint corresponding to a local torsional state. This treatment enables us to adopt a computation time step far larger than hydrodynamic drag time scales. Accordingly, the solvent is no longer treated as a hydrodynamic drag medium. Instead we incorporate its capacity for forming local conformation-dependent dielectric domains. Each translation of the matrix of local topological constraints (LTM's) depends on the conformation-dependent local dielectric created by a confined solvent. Folding pathways are resolved as transitions between patterns of locally encoded structural signals which change within the 1 ns-100 ms time scale range. These coarse folding pathways are generated by a search at regular intervals for structural patterns in the LTM. Each pattern is recorded as a base-pairing pattern (BPP) matrix, a consensus-evaluation operation subject to a renormalization feedback loop. Since several mutually conflicting consensus evaluations might occur at a given time, the need arises for a probabilistic approach appropriate for an ensemble of RNA molecules. Thus, a statistical dynamics of consensus formation is determined by the time evolution of the base pairing probability matrix. These dynamics are generated for a functional RNA molecule, a representative of the so-called group I ribozymes, in order to test the model. The resulting ensemble of conformations is sharply peaked and the most probable structure features the predominance of all phylogenetically conserved intrachain helices tantamount to ribozyme function. Furthermore, the magnesium-aided cooperativity that leads to the shaping of the catalytic core is elucidated. Once the predictive folding algorithm has been implemented, the validity of the so-called "adiabatic approximation" is tested. This approximation requires that conformational microstates be lumped up into BPP's which are treated as quasiequilibrium states, while folding pathways are coarsely represented as sequences of BPP transitions. To test the validity of this adiabatic ansatz, a computation of the coarse Shannon information entropy sigma associated to the specific partition of conformation space into BPP's is performed taking into account the LTM evolution and contrasted with the adiabatic computation. The results reveal a subordination of torsional microstate dynamics to BPP transitions within time scales relevant to folding. This adiabatic entrainment in the long-time limit is thus identified as responsible for the expediency of the folding process.

Algorithms↗

A statistical mechanical model for hydrogen exchange in globular proteins.

We develop a statistical mechanical theory for the mechanism of hydrogen exchange in globular proteins. Using the HP lattice model, we explore how the solvent accessibilities of chain monomers vary as proteins fluctuate from their stable native conformations. The model explains why hydrogen exchange appears to involve two mechanisms under different conditions of protein stability: (1) a "global unfolding" mechanism by which all protons exchange at a similar rate, approaching that of the denatured protein, and (2) a "stable-state" mechanism by which protons exchange at rates that can differ by many orders of magnitude. There has been some controversy about the stable-state mechanism: does exchange take place inside the protein by solvent penetration, or outside the protein by the local unfolding of a subregion? The present model indicates that the stable-state mechanism of exchange occurs through an ensemble of conformations, some of which may bear very little resemblance to the native structure. Although most fluctuations are small-amplitude motions involving solvent penetration or local unfolding, other fluctuations (the conformational distant relatives) can involve much larger transient excursions to completely different chain folds.

Drug Stability↗

Structure of the C-terminal fragment 300-320 of the rat angiotensin II AT1A receptor and its relevance with respect to G-protein coupling.

Angiotensin II AT1A receptor is coupled to G-protein, and the molecular mechanism of signal transduction is still unclear. The solution conformation of a synthetic peptide corresponding to residues 300-320 of the rat AT1A receptor, located in the C-terminal cytoplasmic tail and indicated by mutagenesis work to be critical for the G-protein coupling, has been investigated by circular dichroism (CD), nuclear magnetic resonance (NMR) and restrained molecular dynamics calculations. The CD data indicate that, in acidic water, at concentration below 0.8 mM, the peptide exists in a predominantly coil structure while at higher concentration it can form helical aggregates; addition of small amounts of trifluoroethanol induces a secondary structure, mostly due to the presence of helical elements. Using NMR-derived constraints, an ensemble of conformers for the peptide has been determined by restrained molecular dynamics calculations. Analysis of the converged three-dimensional structures indicates that a significant population of them adopts an amphipathic alpha-helical conformation that, depending upon experimental conditions, presents a variable extension in the stretch Leu6-Tyr20. An equilibrium with nonhelical structured conformers is also observed. We suggest that the capability of the peptide to modulate its secondary structure as a function of the medium dielectric constant, as well as its ability to form helical aggregates by means of intermolecular hydrophobic interactions, can play a significant role for G-protein activation.

Animals↗

Cross-validation tests of time-averaged molecular dynamics refinements for determination of protein structures by X-ray crystallography.

Time-averaged structure-factor restraints have been used in two molecular dynamics refinement schemes to define ensembles of conformations for myoglobin that fit the experimentally measured Bragg scattering from P6 crystals. The geometries of the structures have been maintained to the same currently acceptable limits in all cases. Free R value analysis was used to assess the validity of the two approaches. In the first scheme, where atoms have no B values, the decrease in R value was found to be spurious as judged by a concomitant increase in the free R value. The other scheme, however, which retains individual B values, was found to yield both low R values and low free R values; thus, here the additional variables introduced by modeling the protein in terms of an evolving ensemble of states do not overfit the data. For comparison, refinements were also carried out on the system using several other techniques for isotropic and anisotropic crystallographic refinement. The time-averaged refinements with B values compare quite favorably with the standard methods, but yield additional information about substates of the system. Hence, correctly applied time-averaged refinements can yield accurate models for protein molecules; moreover, by essentially relaxing the harmonic approximation from the refinement process, these refinements allow a more detailed description of the motions of complex molecules, such as proteins, to be determined from X-ray crystallographic data.

Journal Article↗

Application of 2D and 3D NMR experiments to the conformational study of a diantennary oligosaccharide.

By the application of homonuclear 3D NOE-HOHAHA and heteronuclear 3D HMQC-NOE experiments in studies of complex oligosaccharides, NOEs can be investigated which are hidden in conventional 2D NOE spectra. In the 3D NOE-HOHAHA spectrum omega 3 cross sections were considered to be the most suitable for assignment of NOEs. Alternatively, these cross sections could be measured separately in selective 2D HOHAHA-NOE spectroscopy. The advantages and limitations of the 2D alternative are compared with those of the 3D NOE-HOHAHA approach. In 3D HMQC-NOE spectroscopy the larger chemical shift displacement of the carbon spectrum with respect to the proton spectrum can be used to unmask NOEs hidden in the bulk region. If the extra proton dimension is not needed, 2D HMQC-NOE is a good alternative. The suitability of 2D and 3D NOE-HOHAHA and HMQC-NOE experiments for the estimation of proton-proton distances is demonstrated by comparing the results of these experiments on a diantennary asparagine-linked oligosaccharide with those of a conventional 2D NOE experiment. NOEs identified in the 2D and 3D NOE-HOHAHA as well as HMQC-NOE experiments, so far not identified or not quantified in 2D NOE experiments, are discussed in relation to each glycosidic linkage. The flexibility of the Man alpha (1-3)Man linkage is demonstrated, confirming the existence of an ensemble of conformations for this linkage.

Carbohydrate Conformation↗

Functional significance of the conformational dynamics of the N-terminal segment of secreted phospholipase A2 at the interface.

The kinetic and fluorescence properties of several pig pancreatic phospholipase A2 (PLA2) with substitutions and deletions in the N-terminal region and of tyrosines 52 and 73 are characterized. The substitutions Ala-1-D-Ala or -Gly, Trp-3-Phe, Gln-4-Nle, Arg-6-Glu, Tyr-52-Phe, and Tyr-73-Phe had at the most only a modest effect on the interfacial catalytic activity on the anionic interface to which they bind with high affinity. The observed rate of hydrolysis in the scooting mode by deletion mutants lacking one or more successive residues from the N-terminal region was lower by 50-95%. Detailed kinetic analysis of the deletion mutant lacking Ala-1 (des-1-AMPA) showed that the 50% decrease in the rate is due to a 5-fold increase in the interfacial Michaelis-Menten parameter, KM*, without a significant change in kcat. These results and direct measurements show that the primary effect of Ala-1 deletion is to lower the affinity for the active site directed ligands. Although the affinity of these mutants for anionic interface remains the same as for the wild type, the affinity for zwitterionic neutral diluents is considerably lower. Significant differences in the fluorescence quantum yields and the heterogeneity in the frequency-domain fluorescence intensity decays of these enzymes suggest that both in solution and at the interface the N-terminal region is an ensemble of conformations rather than a discrete state. Additional results suggest that the interfacial microenvironment of Trp-3 in des-1-AMPA is more polar and Trp-3 is more accessible to quenching by acrylamide.(ABSTRACT TRUNCATED AT 250 WORDS)

Animals↗

Analysis of ligand binding and protein dynamics of human retinoid X receptor alpha ligand-binding domain by nuclear magnetic resonance.

Retinoid X receptors (RXRs) are nuclear receptors that can activate transcription as homodimers or as obligate heterodimeric partners of other nuclear receptors. While the crystal structures of the RXR ligand-binding domains (LBD) have been previously determined, the dynamics of activation is less well characterized at an atomic level. To probe the effect of ligand binding on RXR LBD dynamics, we initiated nuclear magnetic resonance studies of recombinant human RXRalpha LBD (T223-T462) with and without bound 9-cis-retinoic acid (9cRA). The 1HN, 15N, 13C(alpha), 13CO, and 13C(beta) resonance assignments were established for 164 of 240 residues in apo-RXRalpha LBD. Resonances corresponding to an additional 47 residues emerged upon 9cRA binding. These additional residues included those located in the vicinity of the ligand-binding pocket (helices H3, H5, and strands S1, S2), as well as residues located at the dimerization interface (helices H9 and H10). Thus 9cRA binding stabilized the ligand-binding pocket and had allosteric effects on the dimerization interface. Ligand-induced chemical shift perturbations outside the binding cavity were mapped to helix H3 and the AF-2 helix H12, indicating conformational changes in these regions. However, helix H11, a component of the tetramerization interface, and a large part of helix H10, a component of the dimerization interface, remained undetectable even after 9cRA binding. Although apo- and holo-hRXRalpha LBD existed predominantly as homodimers in solution, exchange between monomeric, dimeric, and tetrameric forms of the protein could have contributed to line broadening of cross-peaks corresponding to helices H10 and H11. 15N T1, T2, and steady-state {1H}-15N NOE data collected at 500 and 700 MHz static magnetic fields showed that the internal motions for the residues in the H1-H3 loop (K245-D263) were much less restricted than those in the protein core for both apo- and holo-forms. Significant exchange R(ex) contributions to the transverse relaxation rate were detected for most of the residues measured in both apo- and holo-RXRalpha LBDs by transverse relaxation optimized spectroscopy-Carr-Purcell-Meiboom-Gill (CPMG) experiments at two B1 field strengths. Taken together these results suggest that the RXRalpha LBD exists as a dynamic ensemble of conformations, even after binding its cognate ligand. Such dynamic characteristics may allow RXRalpha to partner with multiple nuclear receptors.

Alitretinoin↗

Efficient molecular docking of NMR structures: application to HIV-1 protease.

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.

Algorithms↗

Direct demonstration of the flexibility of the glycosylated proline-threonine linker in the Cellulomonas fimi Xylanase Cex through NMR spectroscopic analysis.

The modular xylanase Cex (or CfXyn10A) from Cellulomonas fimi consists of an N-terminal catalytic domain and a C-terminal cellulose-binding domain, joined by a glycosylated proline-threonine (PT) linker. To characterize the conformation and dynamics of the Cex linker and the consequences of its modification, we have used NMR spectroscopy to study full-length Cex in its nonglycosylated ( approximately 47 kDa) and glycosylated ( approximately 51 kDa) forms. The PT linker lacks any predominant structure in either form as indicated by random coil amide chemical shifts. Furthermore, heteronuclear (1)H-(15)N nuclear Overhauser effect relaxation measurements demonstrate that the linker is flexible on the ns-to-ps time scale and that glycosylation partially dampens this flexibility. The catalytic and cellulose-binding domains also exhibit identical amide chemical shifts whether in isolation or in the context of either unmodified or glycosylated full-length Cex. Therefore, there are no noncovalent interactions between the two domains of Cex or between either domain and the linker. This conclusion is supported by the distinct (15)N relaxation properties of the two domains, as well as their differential alignment within a magnetic field by Pf1 phage particles. These data demonstrate that the PT linker is a flexible tether, joining the structurally independent catalytic and cellulose-binding domains of Cex in an ensemble of conformations; however, more extended forms may predominate because of restrictions imparted by the alternating proline residues. This supports the postulate that the binding-domain anchors Cex to the surface of cellulose, whereas the linker provides flexibility for the catalytic domain to hydrolyze nearby hemicellulose (xylan) chains.

Catalytic Domain↗

Structural factors controlling ligand binding to myoglobin: a kinetic hole-burning study.

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as "kinetic hole burning." The study of these correlations gives direct insights into the structure-function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Animals↗

Predicting the equilibrium protein folding pathway: structure-based analysis of staphylococcal nuclease.

The equilibrium folding pathway of staphylococcal nuclease (SNase) has been approximated using a statistical thermodynamic formalism that utilizes the high-resolution structure of the native state as a template to generate a large ensemble of partially folded states. Close to 400,000 different states ranging from the native to the completely unfolded states were included in the analysis. The probability of each state was estimated using an empirical structural parametrization of the folding energetics. It is shown that this formalism predicts accurately the stability of the protein, the cooperativity of the folding/unfolding transition observed by differential scanning calorimetry (DSC) or urea denaturation and the thermodynamic parameters for unfolding. More importantly, this formalism provides a quantitative account of the experimental hydrogen exchange protection factors measured under native conditions for SNase. These results suggest that the computer-generated distribution of states approximates well the ensemble of conformations existing in solution. Furthermore, this formalism represents the first model capable of quantitatively predicting within a unified framework the probability distribution of states seen under native conditions and its change upon unfolding.

Enzyme Stability↗

1H- and natural abundance 15N-NMR studies of a derivative of a rabies glycoprotein fragment.

Using a combination of one- and two-dimensional methods, 1H- and 15N-nmr spectroscopy has been employed to perform the complete assignment and the structural determination of the immunogenic undecapeptide CTTTNSRGTTT in DMSO solution. Nuclear Overhauser enhancement spectroscopy experiments indicated the presence of secondary structures, mainly turn-like structures, which only represent a family, albeit a dominant one, of an ensemble of conformations available to the peptide. Since reverse turns may play an important role as intermediates in protein folding, the experimental observations described here may link the immunological and theoretical approaches to protein folding.

Amino Acid Sequence↗

Development of an extended simulated annealing method: application to the modeling of complementary determining regions of immunoglobulins.

An extended simulated annealing process (ESAP) has been developed in order to obtain an ensemble of conformations of a peptide segment from a protein fluctuating at a given temperature. The annealing process was performed with a fast Monte Carlo method using the scaled collective variables developed by Noguti and Go. The system was divided into two parts: one consists of one or more peptide segments and is flexible around the main-chain and side-chain torsional angles; the other represents the rest of the molecule and was maintained fixed at the atomic positions determined by x-ray experiments. The target function included the nonbonding atomic interactions and a distance function to anchor the N and C terminal ends of each segment to the fixed part. Three systems of complementary determining regions (CDR) of antibodies were tested and compared to x-ray data: L2 loop (7 residues) of the light chain of lambda-type Bence-Jones protein, H1 and the H2 loops (14 residues) of McPC603, and H1 and H2 loops (12 residues) of HyHEL-5. Each state of CDR conformations was characterized at room temperature by the average of their coordinates (average conformation) and the internal energy. With a limited number of annealing processes (10), starting from the extended conformation, we have obtained states with conformations close to the observed x-ray structures, from 1.1 to 1.7 A root mean square deviation (rmsd) of main-chain atoms depending on the system. These states were identical or within 0.25 A rmsd of those of lowest internal energy. For unknown CDR structures the criteria of lowest internal energies from ESAP can be used to predict hypervariable loop structures in antibodies with an accuracy comparable to other methods.

Amino Acid Sequence↗

Stability issues of covalently and noncovalently bonded peptide subunits.

The present study focuses on important questions associated with modeling of peptide and protein stability. Computing at different levels of theory (RHF, B3LYP) for a representative ensemble of conformers of di- and tripeptides of alanine, we found that the Gibbs Free Energy values correlate significantly with the total electronic energy of the molecules (0.922 < or = R2). For noncovalently attached but interacting peptide subunits, such as [For-NH2]2 or [For-L-Ala-NH2]2, we have found, as expected, that the basis set superimposition error (BSSE) is large in magnitude for small basis set but significantly smaller when larger basis sets [e.g., B3LYP/6-311++G(d,p)] are used. Stability of the two hydrogen bonds of antiparallel beta-pleated sheets were quantitatively determined as a function of the molecular structure, S10 and S14, computed as 4.0 +/- 0.5 and 8.1 +/- 1.1 kcal/mol, respectively. Finally, a suitable thermoneutral isodesmic reaction was introduced to scale both covalently and noncovalently attached peptide units onto a common stability scale. We found that a suitable isodesmic reaction can result in the total electronic energy as well as the Gibbs free energy of a molecule, from its "noninteracting" fragments, as accurate as a few tenths of a kcal per mol. The latter observation seems to hold for peptides regardless of their length (1 < or = n < or = 8) or the level of theory applied.

Alanine↗

Lack of coupling between secondary structure formation and collapse in a model polypeptide that mimics early folding intermediates, the F2 fragment of the Escherichia coli tryptophan-synthase beta chain.

The isolated, 101-residue long C-terminal (so called F2) fragment of the beta chain from Escherichia coli tryptophan synthase was shown previously to fold into an ensemble of conformations that are condensed, to contain large amounts of highly dynamic secondary structures, and to behave as a good model of structured intermediates that form at the very early stages of protein folding. Here, solvent perturbations were used to investigate the forces that are involved in stabilizing the secondary structure (monitored by far-UV CD) and the condensation of the polypeptide chain (monitored by dynamic light scattering) in isolated F2. It was observed that neither the ionic strength, nor the pH (between 7 and 10), nor salts of the Hofmeister series affected the global secondary structure contents of F2, whereas some of these salts affected the collapse slightly. Addition of trifluoroethanol resulted in a large increase in both the amount of secondary structure and the Stokes radius of F2. Conversely, F2 became more condensed upon raising the temperature from 4 to 60 degrees C, whereas in this temperature range, the secondary structure undergoes significant melting. These observations lead to the conclusion that, in isolated F2, there is no coupling between the hydrophobic collapse and the secondary structure. This finding will be discussed in terms of early events in protein folding.

Chemical Phenomena↗

Observation of residual dipolar couplings in short peptides.

Residual dipolar couplings provide information on the orientation of individual bond vectors with respect to a unique set of molecular axes. We report that short peptides from 2 to 15 amino acids in length of arbitrary sequence exhibit a modest range of residual dipolar couplings when aligned in either strained polyacrylamide gels or alkyl-PEG bicelles. The absence of significant line broadening in gels suggests peptides align predominantly through steric interactions with the polyacrylamide matrix. However, broadening of NMR lines for a subset of residues aligned in bicelles indicates some peptides bind weakly to these lipid disks, yet a weak negative correlation between the couplings measured in gels and bicelles is consistent with steric hindrance playing a role in both media. The observation of dipolar couplings for peptides of length 10-15 suggests the statistical segment lengths of polypeptide chains must often be >10-15 residues, with data from denatured proteins indicating even larger values. Presumably, local side-chain backbone interactions severely restrict chain flexibility, with the cumulative effect of many such restrictions giving rise to biases in chain direction that may persist for the entire length of a protein chain. Comparison of experimental dipolar couplings for peptides with couplings calculated for ensembles of conformations generated by molecular dynamics should permit evaluation of the accuracy of molecular mechanics potentials in reproducing sequence-specific preferences for phi and psi angles.

Amino Acid Sequence↗

Flexible structures of SIBLING proteins, bone sialoprotein, and osteopontin.

Bone sialoprotein (BSP) and osteopontin (OPN) are two members of the SIBLING (Small Integrin-Binding LIgand, N-linked Glycoprotein) family of genetically related proteins that are clustered on human chromosome 4. We present evidence that this entire family is the result of duplication and subsequent divergent evolution of a single ancient gene. The solution structures of these two post-translationally modified recombinant proteins were solved by one dimensional proton NMR and transverse relaxation times. The polypeptide backbones of both free BSP and OPN rapidly sample an ensemble of conformations consistent with them both being completely unstructured in solution. This flexibility appears to enable these relatively small glycoproteins to rapidly associate with a number of different binding partners including other proteins as well as the mineral phase of bones and teeth. These proteins often function by bridging two proteins of fixed structures into a biologically active complex.

Amino Acid Sequence↗