A molecular dynamics study of the C-terminal fragment of the L7/L12 ribosomal protein

The crystallographic dimer of the C-terminal fragment (CTF) of the L7/L12 ribosomal protein has been subjected to molecular dynamics (MD) simulations. A 90 picosecond (ps) trajectory for the protein dimer, 19 water molecules and two counter ions has been calculated at constant temperature. Effects of intermolecular interactions on the structure and dynamics have been studied. The exact crystallographic symmetry is lost and the atomic fluctuations differ from one monomer to the other. The average MD structure is more stable than the X-ray one, as judged by accessible surface area and energy calculations. Crystal (non-dimeric) interactions have been simulated in another 40 ps trajectory by using harmonic restraints to represent intermolecular hydrogen bonds. The conformational changes with respect ot the X-ray structure are then virtually suppressed.The unrestrained dimer trajectory has been scanned for cooperative motions involving secondary structure elements. The intrinsic collective motions of the monomer are transmitted via intermolecular contacts to the dimer structure.The existence of a stable dimeric form of CTF, resembling the crystallographic one, has been documented. At the cost of fairly small energy expenditure the dimer has considerable conformational flexibility. This flexibility may endow the dimer with some functional potential as an energy transducer.     

1H-NMR study of protein L7/L12 from bacterial ribosomes

The 500 MHz 1H-NMR spectra of dimeric protein L12 from ribosomes shows a limited number of unusually sharp signals at room temperature. This is interpreted as evidence for substantial segmental flexibility of the region in the protein molecule. We have analysed the extent of the flexible region and also the size of the organized structures of the molecule. Thus residues 37-50 were found to be highly mobile whereas the N-terminal and C-terminal region are organized into folded domains.     

A molecular dynamics study of the C-terminal fragment of the L7/L12 ribosomal protein. Secondary structure motion in a 150 picosecond trajectory

A 150 picosecond molecular dynamics computer simulation of the C-terminal fragment of the L7/L12 ribosomal protein from Escherichia coli is reported. The molecular dynamics results are compared with the available high-resolution X-ray data in terms of atomic positions, distances and positional fluctuations. Good agreement is found between the molecular dynamics results and the X-ray data. The form and parameters of the interaction potential energy function and the procedures for deriving it are discussed. Some current misunderstandings concerning the ways of evaluating the efficiency of molecular dynamics algorithms and of application of bond-length constraints in protein simulations are cleared up. The 150 picosecond trajectory has been scanned in a search for correlated motions within and between secondary structure elements. The beta-strands have diffusional stretching modes, and uncorrelated transversal displacements. The dynamic analysis of alpha-helices shows a variety of features. The atomic fluctuations differ between the helix ends; this effect reflects long time-scale motions. Two alpha-helices, alpha A and alpha C, show diffusive longitudinal stretching modes. The third helix, alpha B, has a correlated asymmetric longitudinal stretching; the N-terminal part dominates this behaviour. Furthermore, alpha B presents a librational motion with respect to the other parts of the molecule with a frequency of approximately 5 cm-1. This motion is coupled to helix stretching. Interestingly, the regions of highly conserved residues contain the most mobile parts of the molecule.     

Structure of the C-terminal domain of the ribosomal protein L7/L12 from Escherichia coli at 1.7 A

The structure of a C-terminal fragment of the ribosomal protein L7/L12 from Escherichia coli has been refined using crystallographic data to 1.7 A resolution. The R-value is 17.4%. Six residues at the N terminus are too disordered in the structure to be localized. These residues are probably part of a hinge in the complete L7/L12 molecule. The possibility that a 2-fold crystallographic axis is a molecular 2-fold axis is discussed. A patch of invariant residues on the surface of the dimer is probably involved in functional interactions with elongation factors.     

Copper binding sites in the C-terminal domain of mouse prion protein: A hybrid (QM/MM) molecular dynamics study

We present a hybrid QM/MM Car-Parrinello molecular dynamics study of the copper-loaded C-terminal domain of the mouse prion protein. By means of a statistical analysis of copper coordination in known protein structures, we localized the protein regions with the highest propensity for copper ion binding. The identified candidate structures were subsequently refined via QM/MM simulations. Their EPR characteristics were computed to make contact with the experimental data and to probe the sensitivity to structural and chemical changes. Overall best agreement with the experimental EPR data (Van Doorslaer et al., J Phys Chem B 2001; 105: 1631-1639) and the information currently available in the literature is observed for a binding site involving H187. Moreover, a reinterpretation of the experimental proton hyperfine couplings was possible in the light of the present computational findings.     

Computed spatial homology between the L12 protein of chloroplast ribosome and 1.7 A structure of Escherichia coli L12 domain

A computer-graphic model of the tertiary structure of a functional domain in an organelle ribosomal protein was generated using the amino acid sequence of chloroplast ribosomal protein L12 from spinach (Bartsch, Kimura and Subramanian, Proc. Natl. Acad. Sci. USA 79, 6871-6875, 1982) and 1.7 A resolution coordinates of the E. coli L12 C-terminal fragment crystal (Leijonmarck, Eriksson and Liljas, Nature 286, 824-826, 1980). A comparison between the model and the experimentally derived structure shows that although 40% of the primary structure of this part of the two proteins has undergone amino acid replacements, the gross spatial structure of the domain is maintained and the character of the surfaces of possible functional importance are not significantly altered.     

From structure and dynamics of protein L7/L12 to molecular switching in ribosome

Based on the (1)H-(15)N NMR spectroscopy data, the three-dimensional structure and internal dynamic properties of ribosomal protein L7 from Escherichia coli were derived. The structure of L7 dimer in solution can be described as a set of three distinct domains, tumbling rather independently and linked via flexible hinge regions. The dimeric N-terminal domain (residues 1-32) consists of two antiparallel alpha-alpha-hairpins forming a symmetrical four-helical bundle, whereas the two identical C-terminal domains (residues 52-120) adopt a compact alpha/beta-fold. There is an indirect evidence of the existence of transitory helical structures at least in the first part (residues 33-43) of the hinge region. Combining structural data for the ribosomal protein L7/L12 from NMR spectroscopy and x-ray crystallography, it was suggested that its hinge region acts as a molecular switch, initiating "ratchet-like" motions of the L7/L12 stalk with respect to the ribosomal surface in response to elongation factor binding and GTP hydrolysis. This hypothesis allows an explanation of events observed during the translation cycle and provides useful insights into the role of protein L7/L12 in the functioning of the ribosome.     

Studies on the structure of protein L7/L12 from Escherichia coli ribosomes.

Circular dichroism, infrared and proton magnetic resonance spectroscopy as well as microcalorimetry methods were used to investigate the intact proteins L7/L12 in solution and their different derivatives (L7 with oxidized residues of methionine, fragments 27--120, 1--73 and 74--120)- On the basis of the data obtained the following conclusions have been drawn: (a) there is no beta structure in the protein L7, (B) the N-terminal region of L7 forms a long alpha helix (c) the Phe-30 residue within the N-terminal region of L7 takes part in the dimerization, (d) the C-terminal of L7 is globular and (e) the Phe-54 residue is included in the hydrophobic core of the globular C-terminal region.     

Deletion of C-terminal residues of Escherichia coli ribosomal protein L10 causes the loss of binding of one L7/L12 dimer: ribosomes with one L7/L12 dimer are active

Escherichia coli ribosomal protein L10 binds the two L7/L12 dimers and thereby anchors them to the large ribosomal subunit. C-Terminal deletion variants (Delta10, Delta20, and Delta33 amino acids) of ribosomal protein L10 were constructed in order to define the binding sites for the two L7/L12 dimers and then to make and test ribosomal particles that contain only one of the two dimers. None of the deletions interfered with binding of L10 variants to ribosomal core particles. Deletion of 20 or 33 amino acids led to the inability of the proteins to bind both dimers of protein L7/L12. The L10 variant with deletion of 10 amino acids bound one L7/L12 dimer in solution and when reconstituted into ribosomes promoted the binding of only one L7/L12 dimer to the ribosome. The ribosomes that contained a single L7/L12 dimer were homogeneous by gel electrophoresis where they had a mobility between wild-type 50S subunits and cores completely lacking L7/L12. The single-dimer ribosomal particles supported elongation factor G dependent GTP hydrolysis and protein synthesis in vitro with the same activity as that of two-dimer particles. The results suggest that amino acids 145-154 in protein L10 determine the binding site ("internal-site") for one L7/L12 dimer (the one reported here), and residues 155-164 ("C-terminal-site") are involved in the interaction with the second L7/L12 dimer. Homogeneous ribosomal particles containing a single L7/L12 dimer in each of the distinct sites present an ideal system for studying the location, conformation, dynamics, and function of each of the dimers individually.     

Ca(2+) dissociation from the C-terminal EF-hand pair in calmodulin: a steered molecular dynamics study

We used steered molecular dynamics (SMD) to simulate the process of Ca²⁺ dissociation from the EF-hand motifs of the C-terminal lobe of calmodulin. Based on an analysis of the pulling forces, the dissociation sequences and the structural changes, we show that the Ca²⁺-coordinating residues lose their binding to Ca²⁺ in a stepwise fashion. The two Ca²⁺ ions dissociate from the two EF-hands simultaneously, with two distinct groups among the five Ca²⁺-coordinating residues affecting the EF-hand conformational changes differently. These results provide new insights into the effects of Ca²⁺ on calmodulin conformation, from which a novel sequential mechanism of Ca²⁺-calmodulin dissociation is proposed.     

Defining a pathway of communication from the C-terminal peptide binding domain to the N-terminal ATPase domain in a AAA protein

AAA proteins remodel other proteins to affect a multitude of biological processes. Their power to remodel substrates must lie in their capacity to couple substrate binding to conformational changes via cycles of nucleotide binding and hydrolysis, but these relationships have not yet been deciphered for any member. We report that when one AAA protein, Hsp104, engages polypeptide at the C-terminal peptide-binding region, the ATPase cycle of the C-terminal nucleotide-binding domain (NBD2) drives a conformational change in the middle region. This, in turn, drives ATP hydrolysis in the N-terminal ATPase domain (NBD1). This interdomain communication pathway can be blocked by mutation in the middle region or bypassed by antibodies that bind there, demonstrating the crucial role this region plays in transducing signals from one end of the molecule to the other.     

Interaction of the G' domain of elongation factor G and the C-terminal domain of ribosomal protein L7/L12 during translocation as revealed by cryo-EM

During tRNA translocation on the ribosome, an arc-like connection (ALC) is formed between the G' domain of elongation factor G (EF-G) and the L7/L12-stalk base of the large ribosomal subunit in the GDP state. To delineate the boundary of EF-G within the ALC, we tagged an amino acid residue near the tip of the G' domain of EF-G with undecagold, which was then visualized with three-dimensional cryo-electron microscopy (cryo-EM). Two distinct positions for the undecagold, observed in the GTP-state and GDP-state cryo-EM maps of the ribosome bound EF-G, allowed us to determine the movement of the labeled amino acid. Molecular analyses of the cryo-EM maps show: (1) that three structural components, the N-terminal domain of ribosomal protein L11, the C-terminal domain of ribosomal protein L7/L12, and the G' domain of EF-G, participate in formation of the ALC; and (2) that both EF-G and the ribosomal protein L7/L12 undergo large conformational changes to form the ALC.     

A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

The proximity of the C-terminal domain of Escherichia coli ribosomal protein L7/L12 to L10 determined by cysteine site-directed mutagenesis and protein-protein cross-linking.

Oligonucleotide-directed mutagenesis was used to produce a serine 89 to cysteine 89 substitution in the C-terminal globular domain of Escherichia coli ribosomal protein L7/L12. Cys-89 represented the only cysteine residue in the protein. L7/L12Cys89 was overproduced in E. coli and purified. An allele replacement strain was also constructed. Growth of this strain was indistinguishable from that of wild type. Ribosomes from the allele replacement strain were used to determine the location of the C-terminal domains of L7/L12 by disulfide cross-linking. A new homobifunctional cysteine-specific cross-linking reagent, 1,4-di[3'-(2'-pyridyldithio)-propionamido]butane, and diagonal gel electrophoresis were used to identify ribosomal proteins cross-linked to L7/L12Cys89. A cross-link between L7/L12 and the single cysteine in L10 was found, in addition to L7/L12 dimers. The L7/L12Cys89-L10 cross-link locates the C-terminal domain of at least one L7/L12 dimer on the body of the large subunit and supports our previous model (Olson, H. M., Sommer, A., Tewari, D. S., Traut, R. R., and Glitz, D. G. (1986) J. Biol. Chem. 261, 6924-6932) that depicts one of the two dimers of L7/L12 on the surface of the body of the 50 S subunit in a bent conformation with the C-terminal domain in close proximity to the N-terminal domain at the base of the stalk.     

Crystal structure of the C-terminal domain of a flagellar hook-capping protein from Xanthomonas campestris

The crystal structure of the C-terminal domain of a hook-capping protein FlgD from the plant pathogen Xanthomonas campestris (Xc) has been determined to a resolution of ca 2.5 Å using X-ray crystallography. The monomer of whole FlgD comprises 221 amino acids with a molecular mass of 22.7 kDa, but the flexible N-terminus is cleaved for up to 75 residues during crystallization. The final structure of the C-terminal domain reveals a novel hybrid comprising a tudor-like domain interdigitated with a fibronectin type III domain. The C-terminal domain of XcFlgD forms three types of dimers in the crystal. In agreement with this, analytical ultracentrifugation and gel filtration experiments reveal that they form a stable dimer in solution. From these results, we propose that the Xc flagellar hook cap protein FlgD comprises two individual domains, a flexible N-terminal domain that cannot be detected in the current study and a stable C-terminal domain that forms a stable dimer.     

Trypsinogen-trypsin transition: a molecular dynamics study of induced conformational change in the activation domain

The trypsinogen to trypsin transition has been investigated by a stochastic boundary molecular dynamics simulation that included a major portion of the trypsin molecule and the surrounding solvent. Attention focused on the "activation domain", which crystallographic studies have shown to be ordered in trypsin and disordered in its zymogen, trypsinogen. The chain segments that form the activation domain were found to exhibit large fluctuations during the simulation of trypsin. To model a difference between trypsin and trypsinogen, the N-terminal residues Ile-16 and Val-17 were removed in the former and replaced by water molecules. As a result of the perturbation, a structural drift of 1-2 A occurred that is limited to the activation domain. Glycine residues are found to act as hinges for the displaced chain segments.