Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.

Chemical modification of Methionines of Ribonuclease a with O-Benzoquinone

The accessibility of methionines in RNAase A to reaction with OBQ has been studied at highly acidic pH. The differences between the rate constants of reactions of the methionine and methionines of RNAase A with OBQ is a reflection on the limited accessibility of methionines in the protein conformation. Nevertheless, at sufficiently high OBQ concentration, all the four methionines of the enzyme can be modified. At lower concentration of OBQ, a derivative may be prepared in which a specific methionine is modified. The introduced chromophore ionizes at around pH 3 in this derivative. The derivative has partial activity towards RNA which is enhanced on addition of S-protein.

Ligand dependent intra and inter subunit communication in human tryptophanyl tRNA synthetase as deduced from the dynamics of structure networks

Homodimeric protein tryptophanyl tRNA synthetase (TrpRS) has a Rossmann fold domain and belongs to the 1c subclass of aminoacyl tRNA synthetases. This enzyme performs the function of acylating the cognate tRNA. This process involves a number of molecules (2 protein subunits, 2 tRNAs and 2 activated Trps) and thus it is difficult to follow the complex steps in this process. Structures of human TrpRS complexed with certain ligands are available. Based on structural and biochemical data, mechanism of activation of Trp has been speculated. However, no structure has yet been solved in the presence of both the tRNA(Trp) and the activated Trp (TrpAMP). In this study, we have modeled the structure of human TrpRS bound to the activated ligand and the cognate tRNA. In addition, we have performed molecular dynamics (MD) simulations on these models as well as other complexes to capture the dynamical process of ligand induced conformational changes. We have analyzed both the local and global changes in the protein conformation from the protein structure network (PSN) of MD snapshots, by a method which was recently developed in our laboratory in the context of the functionally monomeric protein, methionyl tRNA synthetase. From these investigations, we obtain important information such as the ligand induced correlation between different residues of this protein, asymmetric binding of the ligands to the two subunits of the protein as seen in the crystal structure analysis, and the path of communication between the anticodon region and the aminoacylation site. Here we are able to elucidate the role of dimer interface at a level of detail, which has not been captured so far.

Ricin-membrane interaction: membrane penetration depth by fluorescence quenching and resonance energy transfer

The entry of the plant toxin ricin and its A- and B-subunits in model membranes in the presence as well as absence of monosialoganglioside (GM(1)) has been studied. Dioleoylphosphatidylcholine and 5-, 10-, and 12-doxyl- or 9,10-dibromophosphatidylcholines serve as quenchers of intrinsic tryptophan fluorescence of the proteins. The parallax method of Chattopadhyay and London [(1987) Biochemistry 26, 39-45] has been employed to measure the average membrane penetration depth of tryptophans of ricin and its B-chain and the actual depth of the sole Trp 211 in the A-chain. The results indicate that both of the chains as well as intact ricin penetrate the membrane deeply and the C-terminal end of the A-chain is well inside the bilayer, especially at pH 4.5. An extrinsic probe N-(iodoacetyl)-N'-(5-sulfo-1-naphthyl) ethylenediamine (I-AEDANS) has been attached to Cys 259 of the A-chain, and the kinetics of penetration has been followed by monitoring the increase in AEDANS fluorescence at 480 nm. The insertion follows first-order kinetics, and the rate constant is higher at a lower pH. The energy transfer distance analysis between Trp 211 and AEDANS points out that the conformation of the A-chain changes as it inserts into the membrane. CD studies indicate that the helicity of the proteins increases after penetration, which implies that some of the unordered structure in the native protein is converted to the ordered form during this process. Hydrophobic forces seem to be responsible for stabilizing a particular protein conformation inside the membrane.

Preservation of native conformation during aluminium-induced aggregation of tau protein

ALUMINIUM exposure has been shown to result in aggregation of microtubule-associated protein tau in vitro. In the light of recent observations that the native random structure of tau protein is maintained in its monomeric and dimeric states as well as in the paired helical filaments characteristic of Alzheimer's disease, it is likely that factors playing a causative role in neurofibrillary pathology would not drastically alter the native conformation of tau protein. We have studied the interaction of tau protein with aluminium using circular dichroism (CD) and 27(Al) NMR spectroscopy. The CD studies revealed a five-fold increase in the observed ellipticity of the tau-aluminium assembly. The increase in elipticity was not associated with a change in the general conformation of the protein and was most likely due to an aggregation of the tau protein induced by aluminium. Al-27 NMR spectroscopy confirmed the binding of aluminium to tau protein. Hyperphosphorylation of tau in Alzheimer's disease is known to be associated with defective microtubule assembly in this condition. Abnormally phosphorylated tau exists in a polymerized form in the paired helical filaments (PHF) which constitute the neurofibrillary tangles found in Alzheimer's disease. While it is hypothesized that its altered biophysical characteristics render abnormally phosphorylated tau resistant to proteolysis, causing the formation of stable deposits,the sequence of events resulting in the polymerization of tau are little understood, as are the additional factors or modifications required for tills process. Based on the results of our spectroscopic studies, a model for the sequence of events occurring in neurofibrillary pathology is proposed.

A hypothesis on a specific sequence-dependent conformation of DNA and its relation to the binding of the lac-repressor protein

The structure and properties of the double-helical form of the alternating copolymer poly(dA-dT) are considered. Different lines of evidence are interpreted in terms of a structure in which every second phosphate-diester linkage has a conformation different from that of the normal B form. A rationale for this “alternating-B” structure is given which provides an explanation for the effects of chemical modifications of the T residues on the binding of the poly(dA-dT)· poly(dA-dT) to the lac repressor of Escherichia coli.

Conformation-activity correlations for chemotactic tripeptide analogs incorporating dialkyl residues with linear and cyclic alkyl sidechains at position 2

Five stereochemically constrained analogs of the chemotactic tripeptide incorporating 1-aminocycloalkane-1-carboxylic acid (Ac(n)c) and alpha,alpha-dialkylglycines (Deg, diethylglycine; Dpg, n,n-dipropylglycine and Dbg, n,n-dibutylglycine) at position 2 have been synthesized. NMR studies of peptides For-Met-Xxx-Phe-OMe (Xxx = Ac(7)c, I; Ac(8)c, II; Deg, III; Dpg, IV and Dbg, V; For, formyl) establish that peptides with cycloalkyl residues, I and II, adopt folded beta-turn conformations in CDCl3 and (CD3)(2)SO. In contrast, analogs with linear alkyl sidechains, III-V, favour fully extended (C-5) conformations in solution. Peptides I-V exhibit high activity in inducing beta-glucosaminidase release from rabbit neutrophils, with ED(50) values ranging from 1.4-8.0 x 10(-11)M. In human neutrophils the Dxg peptides III-V have ED(50) values ranging from 2.3 x 10(-8) to 5.9 x 10(-10) M, with the activity order being V > IV > III. While peptides I-IV are less active than the parent. For-Met-Leu-Phe-OH, in stimulating histamine release from human basophils, the Dbg peptide V is appreciably more potent, suggesting its potential utility as a probe for formyl peptide receptors.

Effect of Crowding Agents, Signal Peptide, and Chaperone SecB on the Folding and Aggregation of E. coli Maltose Binding Protein

SecB, a soluble cytosolic chaperone component of the Secexport pathway, binds to newly synthesized precursor proteins and prevents their premature aggregation and folding and subsequently targets them to the translocation machinery on the membrane. PreMBP, the precursor form of maltose binding protein, has a 26-residue signal sequence attached to the N-terminus of MBP and is a physiological substrate of SecB. We examine the effect of macromolecular crowding and SecB on the stability and refolding of denatured preMBP and MBP. PreMBP was less stable than MBP (ΔTm =7( 0.5 K) in both crowded and uncrowded solutions. Crowding did not cause any substantial changes in the thermal stability ofMBP(ΔTm=1(0.4 K) or preMBP (ΔTm=0(0.6 K), as observed in spectroscopically monitored thermal unfolding experiments. However, both MBP and preMBP were prone to aggregation while refolding under crowded conditions. In contrast to MBP aggregates, which were amorphous, preMBP aggregates form amyloid fibrils.Under uncrowded conditions, a molar excess of SecB was able to completely prevent aggregation and promote disaggregation of preformed aggregates of MBP. When a complex of the denatured protein and SecB was preformed, SecB could completely prevent aggregation and promote folding of MBP and preMBP even in crowded solution. Thus, in addition to maintaining substrates in an unfolded, export-competent conformation, SecB also suppresses the aggregation of its substrates in the crowded intracellular environment. SecB is also able to promote passive disaggregation of macroscopic aggregates of MBP in the absence of an energy source such as ATP or additional cofactors. These experiments also demonstrate that signal peptide can reatly influence protein stability and aggregation propensity.

Sequence-dependent conformation of an A-DNA double helix: The crystal structure of the octamer d(G-G-T-A-T-A-C-C)

The crystal structures of the synthetic self-complementary octamer d(G-G-T-A-T-A-C-C) and its 5-bromouracil-containing analogue have been refined to R values of 20% and 14% at resolutions of 1·8 and 2·25 Å, respectively. The molecules adopt an A-DNA type double-helical conformation, which is minimally affected by crystal forces. A detailed analysis of the structure shows a considerable influence of the nucleotide sequence on the base-pair stacking patterns. In particular, the electrostatic stacking interactions between adjacent guanine and thymine bases produce symmetric bending of the double helix and a major-groove widening. The sugar-phosphate backbone appears to be only slightly affected by the base sequence. The local variations in the base-pair orientation are brought about by correlated adjustments in the backbone torsion angles and the glycosidic orientation. Sequence-dependent conformational variations of the type observed here may contribute to the specificity of certain protein-DNA interactions.

Solid state and solution conformation of Boc-L-Met-Aib-L-Phe-OMe. Beta-turn conformation of a sequence related to an active chemotactic peptide analog

The conformation of the peptide Boc-L-Met-Aib-L-Phe-OMe has been studied in the solid state and solution by X-ray diffraction and 1H n.m.r., respectively. The peptide differs only in the N-terminal protecting group from the biologically active chemotactic peptide analog formyl-L-Met-Aib-L-Phe-OMe. The molecules adopt a type-II beta-turn in the solid state with Met and Aib as the corner residues (phi Met = -51.8 degrees, psi Met = 139.5 degrees, phi Aib = 58.1 degrees, psi Aib = 37.0 degrees). A single, weak 4----1 intramolecular hydrogen bond is observed between the Boc CO and Phe NH groups (N---O 3.25 A, N-H---O 128.4 degrees). 1H n.m.r. studies, using solvent and temperature dependencies of NH chemical shifts and paramagnetic radical induced line broadening of NH resonances, suggest that the Phe NH is solvent shielded in CDCl3 and (CD3)2SO. Nuclear Overhauser effects observed between Met C alpha H and Aib NH protons provide evidence of the occurrence of Met-Aib type-II beta-turns in these solvents.

Conformation of the LL and LD hairpin bends with internal hydrogen bonds in proteins and peptides

The conformation of three linked peptide units having an internal 4 → 1 type of hydrogen bond has been studied in detail, and the low energy conformations are listed. These conformations all lead to the reversal of the chain direction, and may therefore be called as “hairpin bends” or “U-bends”. Since this bend can occur at the end of two chains hydrogen-bonded in the antiparallel β-conformation, it is also known as the “β-bend”. Two types of conformation are possible when the residues at the second and third Cα atoms are both of type L (the LL bend), while only one type is possible for the LD and the DL bend. The LL bend can also accommodate the sequences LG, GL, GG (G = glycine), while the LD bend can accommodate the sequences LG, GD and GG. The conformations for the sequences DD and DL are exact inverses (or mirror images) of those for the sequences LL and LD, respectively, and have dihedral angles (phi2, ψ2), (phi3, ψ3) of the same magnitudes, but of opposite signs as those for the former types, which are listed, along with the characteristics (length, angle and energy) of the hydrogen bonds. A comparison of the theoretical predictions with experimental data (from X-ray diffraction and NMR studies) on proteins and peptides, show reasonably good agreement. However, a systematic trend is observable in the experimental data, slightly deviating from theory, which indicates that some deformations occur in the shapes of the peptide units forming the bend, differing from that of the standard planar peptide unit.

Conformation of polypeptide-chains containing both l-residues and D-residues .2. Double-helical structures of poly-ld-peptides

Polypeptides with alternating L- and D-amino acid residues can take up stereochemically satisfactory coaxial double-helical structures, both antiparallel and parallel, which are stabilized by systematic interchain NH O hydrogen bonds. Semiempirical energy calculations over allowed regions of conformational space have yielded the characteristics of these double-helices. There are four possible types of antiparallel double-helices - A3, A4, A5 and A6, with n, the number of LD peptide units per turn, around 2.8, 3.6, 4.5 and 5.5 respectively, while for the parallel double-helices there are two types, P3 and P4, having similar helical parameters as in A3 and A4. The hydrogen-bonding scheme restricts the pitch in all the models to the narrow range of 10.0 to 11.5 Å. All these helices have large central cores whose radii increase proportionately with n. In this respect, A3 and A4 are suitable models for the structure of gramicidin A. In terms of their relative energies, antiparallel double-helices are marginally more stable than those with parallel strands. Our results indicate that the energy differences amongst the members in the antiparallel family are not significant and thus provide an explanation for the polymorphism reported for poly(γ-benzyl-LD-glutamate).

Pi-turns in proteins and peptides: Classification, conformation, occurrence, hydration and sequence.

The i + 5-->i hydrogen bonded turn conformation (pi-turn) with the fifth residue adopting alpha L conformation is frequently found at the C-terminus of helices in proteins and hence is speculated to be a "helix termination signal." An analysis of the occurrence of i + 5-->i hydrogen bonded turn conformation at any general position in proteins (not specifically at the helix C-terminus), using coordinates of 228 protein crystal structures determined by X-ray crystallography to better than 2.5 A resolution is reported in this paper. Of 486 detected pi-turn conformations, 367 have the (i + 4)th residue in alpha L conformation, generally occurring at the C-terminus of alpha-helices, consistent with previous observations. However, a significant number (111) of pi-turn conformations occur with (i + 4)th residue in alpha R conformation also, generally occurring in alpha-helices as distortions either at the terminii or at the middle, a novel finding. These two sets of pi-turn conformations are referred to by the names pi alpha L and pi alpha R-turns, respectively, depending upon whether the (i + 4)th residue adopts alpha L or alpha R conformations. Four pi-turns, named pi alpha L'-turns, were noticed to be mirror images of pi alpha L-turns, and four more pi-turns, which have the (i + 4)th residue in beta conformation and denoted as pi beta-turns, occur as a part of hairpin bend connecting twisted beta-strands. Consecutive pi-turns occur, but only with pi alpha R-turns. The preference for amino acid residues is different in pi alpha L and pi alpha R-turns. However, both show a preference for Pro after the C-termini. Hydrophilic residues are preferred at positions i + 1, i + 2, and i + 3 of pi alpha L-turns, whereas positions i and i + 5 prefer hydrophobic residues. Residue i + 4 in pi alpha L-turns is mainly Gly and less often Asn. Although pi alpha R-turns generally occur as distortions in helices, their amino acid preference is different from that of helices. Poor helix formers, such as His, Tyr, and Asn, also were found to be preferred for pi alpha R-turns, whereas good helix former Ala is not preferred. pi-Turns in peptides provide a picture of the pi-turn at atomic resolution. Only nine peptide-based pi-turns are reported so far, and all of them belong to pi alpha L-turn type with an achiral residue in position i + 4. The results are of importance for structure prediction, modeling, and de novo design of proteins.

Hyperexpression of chicken riboflavin carrier protein: Antibodies to the recombinant protein curtail pregnancy in rodents

Chicken riboflavin carrier protein (RCP) is a phosphoglycoprotein present in the egg white and yolk of egg-laying animals and in the sera of laying hens and of estrogenized chicks. The RCP cDNA, encoding a protein of predictedMr27,000, has been cloned into a T7 polymerase-driven vector, and high-level expression was observed on induction with IPTG inEscherichia coli.The protein was largely localized in inclusion bodies when expressed at 37°C but was present in the cytosolic fraction when induced at 22°C. At 37°C, two major bands were detected in whole-cell lysates of the strain expressing the protein. N-terminal sequence analysis indicated that the two proteins represented translated products with and without the pelB leader sequence encoded in the pET20b vector, but both included an additional 10 amino acids generated during cloning procedures. The inclusion body obtained at 37°C, on extraction with detergent, led to preferential solubilization of the protein without the pelB signal sequence. The solubilized recombinant RCP was recognized by polyclonal antisera to native RCP but radioimmunoassay revealed quantitative differences in the epitopes exhibited by the recombinant protein. Thus, sequence-specific monoclonal antibodies to chicken RCP also cross-reacted with the recombinant protein with almost equal efficiency, but antibodies which recognize conformation-dependent epitopes showed relatively reduced cross-reactivity with the recombinant protein. Polyclonal antibodies to recombinant RCP were able to recognize both the native and the denatured RCP. Administration of recombinant RCP antisera to pregnant mice led to embryonic resorption leading to early pregnancy termination. These findings reveal that the recombinant protein will be useful for investigations related to the mechanism of pregnancy termination on immunoneutralization of RCP in mammals, as well as in unraveling folding properties of RCP in terms of its ligand binding and antigenetic determinants exposed at its surface.