Interaction Energy Based Protein Structure Networks

The three dimensional structure of a protein is formed and maintained by the noncovalent interactions among the amino acid residues of the polypeptide chain These interactions can be represented collectively in the form of a network So far such networks have been investigated by considering the connections based on distances between the amino acid residues Here we present a method of constructing the structure network based on interaction energies among the amino acid residues in the protein We have investigated the properties of such protein energy based networks (PENs) and have shown correlations to protein structural features such as the clusters of residues involved in stability formation of secondary and super secondary structural units Further we demonstrate that the analysis of PENs in terms of parameters such as hubs and shortest paths can provide a variety of biologically important information such as the residues crucial for stabilizing the folded units and the paths of communication between distal residues in the protein Finally the energy regimes for different levels of stabilization in the protein structure have clearly emerged from the PEN analysis

Interactions of tris buffer with nucleotides: the crystal structure of tris(hydroxymethyl)methylammonium adenosine 5'-diphosphate dihydrate

The crystal and molecular structures of the Tris salt of adenosine 5'-diphosphate were determined from X-ray diffraction data. The crystals are monoclinic, space P21, and Z = 2 with a=9.198 (2) A, b=6.894 (1) A, c=18.440 (4) A, and beta = 92.55 (2) degrees. Intensity data were collected on an automated diffractometer. The structure was solved by the heavy-atom technique and refined by least squares to R = 0.047. The ADP molecule adopts a folded conformation. The conformation about the glycosidic bond is anti. The conformation of the ribose ring is close to a perfect C(2')-endo-C-(3')-exo puckering. The conformation about C(4')-C(5') is gauche-gauche, similar to other nucleotide structures. The pyrophosphate chain displays a nearly eclipsed geometry when viewed down the P-P vector, unlike the staggered conformation observed in crystal structures of other pyrophosphates. The less favorable eclipsed conformation probably results from the observed association of Tris molecules with the polar diphosphate chain through electrostatic interactions and hydrogen bonds. Such interactions may play an important role in Tris-buffered aqueous solutions of nucleotides and metal ions.

Temperature dependent free energy surface of polymer folding from equilibrium and quench studies

Langevin dynamics simulation studies have been employed to calculate the temperature dependent free energy surface and folding characteristics of a 500 monomer long linear alkane (polyethylene) chain with a realistic interaction potential. Both equilibrium and temperature quench simulation studies have been carried out. Using the shape anisotropy parameter (S) of the folded molecule as the order parameter, we find a weakly first order phase transition between the high-temperature molten globule and low-temperature rodlike crystalline states separated by a small barrier of the order of k(B)T. Near the melting temperature (580 K), we observe an intriguing intermittent fluctuation with pronounced ``1/f noise characteristics'' between these two states with large difference in shape and structure. We have also studied the possibilities of different pathways of folding to states much below the melting point. At 300 K starting from the all-trans linear configuration, the chain folds stepwise into a very regular fourfold crystallite with very high shape anisotropy. Whereas, when quenched from a high temperature (900 K) random coil regime, we identify a two step transition from the random coiled state to a molten globulelike state and, further, to a anisotropic rodlike state. The trajectory reveals an interesting coupling between the two order parameters, namely, radius of gyration (R-g) and the shape anisotropy parameter (S). The rodlike final state of the quench trajectory is characterized by lower shape anisotropy parameter and significantly larger number of gauche defects as compared to the final state obtained through equilibrium simulation starting from all-trans linear chain. The quench study shows indication of a nucleationlike pathway from the molten globule to the rodlike state involving an underlying rugged energy landscape. (C) 2010 American Institute of Physics. doi:10.1063/1.3509398]

Conformational Polymorphism in Telomeric Structures: Loop Orientation and Interloop Pairing in d(G4TnG4),

Sequence repeats constituting the telomeric regions of chromosomes are known to adopt a variety of unusual structures, consisting of a G tetraplex stem and short stretches of thymines or thymines and adenines forming loops over the stem. Detailed model building and molecular mechanics studies have been carried out for these telomeric sequences to elucidate different types of loop orientations and possible conformations of thymines in the loop. The model building studies indicate that a minimum of two thymines have to be interspersed between guanine stretches to form folded-back structures with loops across adjacent strands in a G tetraplex (both over the small as well as large groove), while the minimum number of thymines required to build a loop across the diagonal strands in a G tetraplex is three. For two repeat sequences, these hairpins, resulting from different types of folding, can dimerize in three distinct ways-i.e., with loops across adjacent strands and on same side, with loops across adjacent strands and on opposite sides, and with loops across diagonal strands and on opposite sides-to form hairpin dimer structures. Energy minimization studies indicate that all possible hairpin dimers have very similar total energy values, though different structures are stabilized by different types of interactions. When the two loops are on the same side, in the hairpin dimer structures of d(G(4)T(n)G(4)), the thymines form favorably stacked tetrads in the loop region and there is interloop hydrogen bonding involving two hydrogen bonds for each thymine-thymine pair. Our molecular mechanics calculations on various folded-back as well as parallel tetraplex structures of these telomeric sequences provide a theoretical rationale for the experimentally observed feature that the presence of intervening thymine stretches stabilizes folded-back structures, while isolated stretches of guanines adopt a parallel tetraplex structure

Conformational Polymorphism in G-tetraplex Structures: Strand Reversal by Base Flipover or Sugar Flipover

Guanine rich sequences adopt a variety of four stranded structures, which differ in strand orientation and conformation about the glycosidic bond even though they are all stabilised by Hoogsteen hydrogen bonded guanine tetrads. Detailed model building and molecular mechanics calculations have been carried out to investigate various possible conformations of guanines along a strand and different possible orientations of guanine strands in a G-tetraplex structure. It is found that for an oligo G stretch per se, a parallel four stranded structure with all guanines in anti conformation is favoured over other possible tetraplex structures. Hence an alternating syn-anti arrangement of guanines along a strand is likely to occur only in folded back tetraplex structures with antiparallel G strands. Our study provides a theoretical rationale for the observed alternation of glycosidic conformation and the inverted stacking arrangement arising from base flipover, in antiparallel G-tetraplex structures and also highlights the various structural features arising due to different types of strand orientations. The molecular mechanics calculations help in elucidating the various interactions which stabilize different G-tetraplex structures and indicate that screening of phosphate charge by counterions could have a dramatic effect on groove width in these four stranded structures.

Conformational polymorphism in G-tetraplex structures: strand reversal by base flipover or sugar flipover

Guanlne rich sequences adopt a variety of four stranded structures, which differ in strand orientation and conformation about the glycosldic bond even though they are all stabilised by Hoogsteen hydrogen bonded guanlne tetrads. Detailed model building and molecular mechanics calculations have been carried out to investigate various possible conformations of guanlnes along a strand and different possible orientations of guanlne strands In a G-tetraplex structure. It is found that for an ollgo G stretch per se, a parallel four stranded structure with all guanines In anti conformation is favoured over other possible tetraplex structures. Hence an alternating syn-anti arrangement of guanlnes along a strand is likely to occur only in folded back tetraplex structures with antiparallel G strands. Our study provides a theoretical rationale for the observed alternation of glycosldic conformation and the inverted stacking arrangement arising from base filpover, In antlparallel G-tetraplex structures and also highlights the various structural features arising due to different types of strand orientations. The molecular mechanics calculations help in elucidating the various interactions which stabilize different G-tetraplex structures and indicate that screening of phosphate charge by counterions could have a dramatic effect on groove width in these four stranded structures.

Folding of a Donor-Containing Ionene by Intercalation with an Acceptor

Cationic ionenes that bear electron-rich 1,5-dialkoxynaphthalene (DAN) units within the alkylene segment were allowed to interact with different types of electron-deficient, acceptor-containing molecules in an effort to realize intercalation-induced folding of the ionenes; the collapse of the chains was expected to occur in such a way that the donor and acceptor units become arranged in an alternating fashion. Several acceptor-bearing molecules were prepared by the derivatization of pyromellitic dianhydride and naphthalene tetracarboxylic dianhydride with two different oligoethylene glycol monomethyl ether monoamines. This yielded acceptor molecules with different water solubility and allowed the examination of solvophobic effects in the folding process. UV/Vis spectroscopic studies were carried out by using a 1:1 mixture of the DAN-ionenes and different acceptor molecules in water/DMSO solvent mixtures. The intensity of the charge-transfer (CT) band was seen to increase with the water content in the solvent mixture, thereby suggesting that the intercalation is indeed aided by solvophobic effects. The naphthalene diimide (NDI) bearing acceptor molecules consistently formed significantly stronger CT complexes when compared to the pyromellitic diimide (PDI) bearing acceptor molecules, which is a reflection of the stronger pi-stacking tendency of the former. AFM studies of drop-cast films of different ionene-acceptor combinations revealed that compact folded structures are formed most effectively under conditions in which the strongest CT complex is formed.

Hybrid Polypeptides: Gabapentin as a Stereochemically Constrained gamma-Amino Acid Residue

The design of folded structures in peptides containing the higher homologues of alpha-amino acid residues requires the restriction of the range of local conformational choices In alpha-amino acids stereochemically constrained residues like alpha,alpha-dialkylated residue, aminoisobutyric acid (Aib), and D-Proline ((D)Pro) have proved extremely useful in the design of helices and hairpins in short peptides Extending this approach, backbone substitution and cyclization are anticipated to bc useful in generating conformationally constrained beta- and gamma-residues This brief review provides a survey of work on hybrid peptide sequences concerning the conformationally constrained gamma-amino acid residue 1-aminomethyl cyclohexane acetic acid, gabapentin (Gpn) This achiral, beta,beta-disubstituted, gamma-residue strongly favors gauche-gauche conformations about the C-alpha-C-beta (0(2)) and C-alpha-C-gamma (0(1)) bonds, facilitating local folding The Gpn residue can adopt both C-7 (NH1 -> CO1) and C-9 (CO1 (I)<- NH1+I) hydrogen bonds which are analogous to the C-5 and C7 (gamma-turn) conformations at alpha-residues In conjunction with adjacent residues, Gpn may be used in ay and gamma alpha segments to generate C-12 hydrogen bonded conformations which may be considered as expanded analogs of conventional beta-turns The structural characterization of C-12 helices, C-12/C-10 helices with mixed hydrogen bond directionalities and beta-hairpins incorporating Gpn residues at the turn segment is illustrated (C) 2010 Wiley Periodicals, Inc Biopolymers (Pept Sci) 94 733-741 2010

Structure of N-tert-butyloxycarbonyl-[alpha]-aminoisobutyryl-DL-pipecolyl-[alpha]-aminoisobutyric acid methyl ester

C20H35N3O6 (Boc-Aib-DL-Pip-Aib-OMe, Boc = tert-butyloxycarbonyl, Aib = alpha-aminoisobutyric acid, Pip = pipecolic acid, OMe = methoxy), M(r) = 413.5, monoclinic, P2(1)/c, a = 18.055 (3), b = 15.048 (3), c = 17.173 (3) angstrom, beta = 91.7 (1)-degrees, V = 4663.8 (9) angstrom3, Z = 8, D(m) = 1.16, D(x) = 1.178 Mg m-3, lambda(Mo Kalpha) = 0.71069 angstrom, mu = 0.081 mm-1, F(000) = 1792, T = 297 K. The final R value for 4925 [I greater-than-or-equal-to 3sigma(I)] reflections is 0.065 (wR = 0.067). The peptide backbone of the two independent molecules in the asymmetric unit is folded at the -Aib-Pip- sequence to form a type-I (I') beta-bend stabilized by a 1 <-- 4 intramolecular N-H...O=C hydrogen bond between the Aib(3) peptide N-H and Boc urethane C=O groups.

Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG(3)) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

Crystal and molecular structure of Boc-Phe-Val-OMe; comparison of the peptide conformation with its dehydro analogue

The crystal structure of the peptide Boc-Phe-Val-OMe determined by X-ray diffraction methods is reported in this paper. The crystals grown from aqueous methanol are orthorhombic, space group P2(1)2(1)2(1), a = 11.843(2), b = 21.493(4), c = 26.676(4)Angstrom and V = 6790 Angstrom(3). Data were collected on a CAD4 diffractometer using MoK2 radiation (lambda = 0.7107 Angstrom) up to Bragg angle theta = 26 degrees. The structure was solved by direct methods and refined by a least-squares procedure to an R value of 6.8% for 3288 observed reflections. There are three crystallographically independent peptide molecules in the asymmetric unit. All the three molecules exhibit extended conformation. The sidechain of the Val(2) residue shows two different conformations. The conformation of the peptide Boc-Phe-Val-OMe is compared with the conformation of Ac-Delta Phe-Val-OH. It is observed that while Boc-Phe-Val-OMe exhibits an extended conformation, Ac-Delta Phe-Val-OH shows a folded conformation. The results of this comparison highlight the conformation constraining property of the Delta Phe residue. Interestingly, even though Boc-Phe-Val-OMe and Ac-Delta Phe-Val-OH are conformationally different, they exhibit similar packing patterns in the solid state. (C) Munksgaard 1995.

Refolding of riboflavin carrier protein as probed by biochemical and immunological parameters

The unfolding of the chicken egg white riboflavin carrier protein by disulfide reduction with dithiothreitol led to aggregation with concomitant loss of ligand binding characteristics and the capacity to interact with six monoclonal antibodies directed against surface-exposed discontinuous epitopes. The reduced protein could, however, bind to a monoclonal antibody recognizing sequential epitope. Under optimal conditions of protein refolding, the vitamin carrier protein regained its folded structure with high efficiency with simultaneous complete restoration of hydrophobic flavin binding site as well as the epitopic conformations exposed at the surface in a manner comparable to its native form.

The Folding Mechanism of Barstar: Evidence for Multiple Pathways and Multiple Intermediates

The mechanism of folding of the small protein barstar in the pre-transition zone at pH 7, 25 degrees C has been characterized using rapid mixing techniques. Earlier studies had established the validity of the three-state U-S reversible arrow U-F reversible arrow N mechanism for folding and unfolding in the presence of guanidine hydrochloride (GdnHCl) at concentrations greater than 2.0 M, where U-S and U-F are the slow-refolding and fast-refolding unfolded forms, respectively, and N is the fully folded form. It is now shown that early intermediates, I-S1 and I-S2 as well as a late native-like intermediate, I-N, are present on the folding pathways of U-S, and an early intermediate I-F1 on the folding pathway of U-F, when bars tar is refolded in concentrations of GdnHCl below 2.0 M. The rates of formation and disappearance of I-N, and the rates of formation of N at three different concentrations of GdnHCl in the pre-transition zone have been measured. The data indicate that in 1.5 M GdnHCl, I-N is not fully populated on the U-S --> I-S1 --> I-N --> N pathway because the rate of its formation is so slow that the U-S reversible arrow U-F reversible arrow N pathway can effectively compete with that pathway. In 1.0 M GdnHCl, the U-S --> I-S1 --> I-N transition is so fast that I-N is fully populated. In 0.6 M GdnHCl, I-N appears not to be fully populated because an alternative folding pathway, U-S --> I-S2 --> N, becomes available for the folding of U-S, in addition to the U-S --> I-S1 --> I-N --> N pathway Measurement of the binding of the hydrophobic dye 1-anilino-8-naphthalenesulphonate (ANS) during folding indicates that ANS binds to two distinct intermediates, I-M1 and I-M2, that form within 2 ms on the U-S --> I-M1 --> I-S1 --> I-N --> N and U-S --> I-M2 --> I-S2 --> N pathways. There is no evidence for the accumulation of intermediates that can bind ANS on the folding pathway of U-F.

Evaluation of strain-rate effects in transitional round jets using direct numerical simulation

The coherent flame model uses the strain rate to predict reaction rate per unit flame surface area and some procedure that solves for the dynamics of flame surfaces to predict species distributions. The strainrate formula for the reaction rate is obtained from the analytical solution for a flame in a laminar, plane stagnation point flow. Here, the formula's effectiveness is examined by comparisons with data from a direct numerical simulation (DNS) of a round jetlike flow that undergoes transition to turbulence. Significant differences due to general flow features can be understood qualitatively: Model predictions are good in the braids between vortex rings, which are present in the near field of round jets, as the strain rate is extensional and reaction surfaces are isolated. In several other regions, the strain rate is compressive or flame surfaces are folded close together. There, the predictions are poor as the local flow no longer resembles the model flow. Quantitative comparisons showed some discrepancies. A modified, consistent application of the strain-rate solution did not show significant changes in the prediction of mean reaction rate distributions.

Structural studies of protein unfolding

Methods for macromolecular structure determination (NMR and crystallography) are now being used to get structural information on partially folded and unfolded states of proteins. These techniques, in combination with proton hydrogen exchange studies are powerful tools to extract information on non-native states of proteins. This review discusses progress In this area of protein folding.