Hsc70-interacting HPD loop of the J domain of polyomavirus T antigens fluctuates in ps to ns and micros to ms

The backbone dynamics of the J domain from polyomavirus T antigens have been investigated using 15N NMR relaxation and molecular dynamics simulation. Model-free relaxation analysis revealed picosecond to nanosecond motions in the N terminus, the I-II loop, the C-terminal end of helix II through the HPD loop to the beginning of helix III, and the C-terminal end of helix III to the C terminus. The backbone dynamics of the HPD loop and termini are dominated by motions with moderately large amplitudes and correlation times of the order of a nanosecond or longer. Conformational exchange on the microsecond to millisecond timescale was identified in the HPD loop, the N and C termini, and the I-II loop. A 9.7ns MD trajectory manifested concerted swings of the HPD loop. Transitions between major and minor conformations of the HPD loop featured distinct patterns of change in backbone dihedral angles and hydrogen bonds. Fraying of the C-terminal end of helix II and the N-terminal end of helix III correlated with displacements of the HPD loop. Correlation of crankshaft motions of Gly46 and Gly47 with the collective motions of the HPD loop suggested an important role of the two glycine residues in the mobility of the loop. Fluctuations of the HPD loop correlated with relative reorientation of side-chains of Lys35 and Asp44 that interact with Hsc70.     

Crystal structures of a pantothenate synthetase from M. tuberculosis and its complexes with substrates and a reaction intermediate

Pantothenate biosynthesis is essential for the virulence of Mycobacterium tuberculosis, and this pathway thus presents potential drug targets against tuberculosis. We determined the crystal structure of pantothenate synthetase (PS) from M. tuberculosis, and its complexes with AMPCPP, pantoate, and a reaction intermediate, pantoyl adenylate, with resolutions from 1.6 to 2 A. PS catalyzes the ATP-dependent condensation of pantoate and beta-alanine to form pantothenate. Its structure reveals a dimer, and each subunit has two domains with tight association between domains. The active-site cavity is on the N-terminal domain, partially covered by the C-terminal domain. One wall of the active site cavity is flexible, which allows the bulky AMPCPP to diffuse into the active site to nearly full occupancy when crystals are soaked in solutions containing AMPCPP. Crystal structures of the complexes with AMPCPP and pantoate indicate that the enzyme binds ATP and pantoate tightly in the active site, and brings the carboxyl oxygen of pantoate near the alpha-phosphorus atom of ATP for an in-line nucleophilic attack. When crystals were soaked with, or grown in the presence of, both ATP and pantoate, a reaction intermediate, pantoyl adenylate, is found in the active site. The flexible wall of the active site cavity becomes ordered when the intermediate is in the active site, thus protecting it from being hydrolyzed. Binding of beta-alanine can occur only after pantoyl adenylate is formed inside the active site cavity. The tight binding of the intermediate pantoyl adenylate suggests that nonreactive analogs of pantoyl adenylate may be inhibitors of the PS enzyme with high affinity and specificity.     

Investigation of Domain Motions in Bacteriophage T4 Lysozyme

Abstract

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge loci were found in the simulation. One locus comprises residues 8–14 near the C-terminal of the A helix; the other site, residues 77–83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19–24) acted as a mobile ‘flap’ and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.     

Dynamics of Klebsiella pneumoniae OmpA transmembrane domain: The four extracellular loops display restricted motion behavior in micelles and in lipid bilayers

The transmembrane domain of Klebsiella pneumoniae OmpA (KpOmpA) possesses four long extracellular loops that exhibit substantial sequence variability throughout OmpA homologs in Enterobacteria, in comparison with the highly conserved membrane-embedded β-barrel core. These loops are responsible for the immunological properties of the protein, including cellular and humoral recognition. In addition to key features revealed by structural elucidation of the KpOmpA transmembrane domain in detergent micelles, studies of protein dynamics provide insight into its function and/or mechanism of action. We have investigated the dynamics of KpOmpA in a lipid bilayer, using magic angle spinning solid-state NMR. The dynamics of the β-barrel and loop regions were probed by the spin-lattice relaxation times of the C(α) and C(β) atoms of the serine and threonine residues, and by cross-polarization dynamics. The β-barrel core of the protein is rigid; the C-terminal halves of two of the four extracellular loops (L1 and L3), which are particularly long in KpOmpA, are highly mobile. The other two loops (L2 and L4), which are very similar to their homologs in Escherichia coli OmpA, and the N-terminal halves of L1 and L3 exhibit more restricted motions. We suggest a correlation between the sequence variability and the dynamics of certain loop regions, which accounts for their respective contributions to the structural and immunological properties of the protein.     

Conformational dynamics of the major yeast phosphatidylinositol transfer protein sec14p: insight into the mechanisms of phospholipid exchange and diseases of sec14p-like protein deficiencies

Molecular dynamics simulations coupled with functional analyses of the major yeast phosphatidylinositol/phosphatidylcholine transfer protein Sec14p identify structural elements involved in regulating the ability of Sec14p to execute phospholipid exchange. The molecular dynamics simulations suggest large rigid body motions within the Sec14p molecule accompany closing and opening of an A(10)/T(4)/A(11) helical gate, and that "state-of-closure" of this helical gate determines access to the Sec14p phospholipid binding cavity. The data also project that conformational dynamics of the helical gate are controlled by a hinge unit (residues F(212), Y(213), K(239), I(240), and I(242)) that links to the N- and C-terminal ends of the helical gate, and by a novel gating module (composed of the B(1)LB(2) and A(12)LT(5) substructures) through which conformational information is transduced to the hinge. The (114)TDKDGR(119) motif of B(1)LB(2) plays an important role in that transduction process. These simulations offer new mechanistic possibilities for an important half-reaction of the Sec14p phospholipid exchange cycle that occurs on membrane surfaces after Sec14p has ejected bound ligand, and is reloading with another phospholipid molecule. These conformational transitions further suggest structural rationales for known disease missense mutations that functionally compromise mammalian members of the Sec14-protein superfamily.     

Molecular dynamics simulations of peptides and proteins with amplified collective motions

We present a novel method that uses the collective modes obtained with a coarse-grained model/anisotropic network model to guide the atomic-level simulations. Based on this model, local collective modes can be calculated according to a single configuration in the conformational space of the protein. In the molecular dynamics simulations, the motions along the slowest few modes are coupled to a higher temperature by the weak coupling method to amplify the collective motions. This amplified-collective-motion (ACM) method is applied to two test systems. One is an S-peptide analog. We realized the refolding of the denatured peptide in eight simulations out of 10 using the method. The other system is bacteriophage T4 lysozyme. Much more extensive domain motions between the N-terminal and C-terminal domain of T4 lysozyme are observed in the ACM simulation compared to a conventional simulation. The ACM method allows for extensive sampling in conformational space while still restricting the sampled configurations within low energy areas. The method can be applied in both explicit and implicit solvent simulations, and may be further applied to important biological problems, such as long timescale functional motions, protein folding/unfolding, and structure prediction.     

Conformational dynamics of the ligand-binding domain of inward rectifier K channels as revealed by molecular dynamics simulations: toward an understanding of Kir channel gating

Inward rectifier (Kir) potassium channels are characterized by two transmembrane helices per subunit, plus an intracellular C-terminal domain that controls channel gating in response to changes in concentration of various ligands. Based on the crystal structure of the tetrameric C-terminal domain of Kir3.1, it is possible to build a homology model of the ATP-binding C-terminal domain of Kir6.2. Molecular dynamics simulations have been used to probe the dynamics of Kir C-terminal domains and to explore the relationship between their dynamics and possible mechanisms of channel gating. Multiple simulations, each of 10 ns duration, have been performed for Kir3.1 (crystal structure) and Kir6.2 (homology model), in both their monomeric and tetrameric forms. The Kir6.2 simulations were performed with and without bound ATP. The results of the simulations reveal comparable conformational stability for the crystal structure and the homology model. There is some decrease in conformational flexibility when comparing the monomers with the tetramers, corresponding mainly to the subunit interfaces in the tetramer. The beta-phosphate of ATP interacts with the side chain of K185 in the Kir6.2 model and simulations. The flexibility of the Kir6.2 tetramer is not changed greatly by the presence of bound ATP, other than in two loop regions. Principal components analysis of the simulated dynamics suggests loss of symmetry in both the Kir3.1 and Kir6.2 tetramers, consistent with "dimer-of-dimers" motion of subunits in C-terminal domains of the corresponding Kir channels. This is suggestive of a gating model in which a transition between exact tetrameric symmetry and dimer-of-dimers symmetry is associated with a change in transmembrane helix packing coupled to gating of the channel. Dimer-of-dimers motion of the C-terminal domain tetramer is also supported by coarse-grained (anisotropic network model) calculations. It is of interest that loss of exact rotational symmetry has also been suggested to play a role in gating in the bacterial Kir homolog, KirBac1.1, and in the nicotinic acetylcholine receptor channel.     

The crystal structure of E. coli pantothenate synthetase confirms it as a member of the cytidylyltransferase superfamily

BACKGROUND: Pantothenate synthetase (EC 6.3.2.1) is the last enzyme of the pathway of pantothenate (vitamin B(5)) synthesis. It catalyzes the condensation of pantoate with beta-alanine in an ATP-dependent reaction. RESULTS: We describe the overexpression, purification, and crystal structure of recombinant pantothenate synthetase from E. coli. The structure was solved by a selenomethionine multiwavelength anomalous dispersion experiment and refined against native data to a final R(cryst) of 22.6% (R(free) = 24.9%) at 1.7 A resolution. The enzyme is dimeric, with two well-defined domains per protomer: the N-terminal domain, a Rossmann fold, contains the active site cavity, with the C-terminal domain forming a hinged lid. CONCLUSIONS: The N-terminal domain is structurally very similar to class I aminoacyl-tRNA synthetases and is thus a member of the cytidylyltransferase superfamily. This relationship has been used to suggest the location of the ATP and pantoate binding sites and the nature of hinge bending that leads to the ternary enzyme-pantoate-ATP complex.     

The C-Terminal Variable Region Specifies the Dynamic Properties of Arabidopsis Microtubule-Associated Protein MAP65 Isotypes

The microtubule-associated protein, MAP65, is a member of a family of divergent microtubule-associated proteins from different organisms generally involved in maintaining the integrity of the central spindle in mitosis. The dicotyledon Arabidopsis thaliana and the monocotyledon rice (Oryza sativa) genomes contain 9 and 11 MAP65 genes, respectively. In this work, we show that the majority of these proteins fall into five phylogenetic clades, with the greatest variation between clades being in the C-terminal random coil domain. At least one Arabidopsis and one rice isotype is within each clade, indicating a functional specification for the C terminus. In At MAP65-1, the C-terminal domain is a microtubule binding region (MTB2) harboring the phosphorylation sites that control its activity. The At MAP65 isotypes show differential localization to microtubule arrays and promote microtubule polymerization with variable efficiency in a MTB2-dependent manner. In vivo studies demonstrate that the dynamics of the association and dissociation of different MAP65 isotypes with microtubules can vary up to 10-fold and that this correlates with their ability to promote microtubule polymerization. Our data demonstrate that the C-terminal variable region, MTB2, determines the dynamic properties of individual isotypes and suggest that slower turnover is conditional for more efficient microtubule polymerization.     

A novel isoform of pantothenate synthetase in the Archaea

The linear biosynthetic pathway leading from alpha-ketoisovalerate to pantothenate (vitamin B5) and on to CoA comprises eight steps in the Bacteria and Eukaryota. Genes for up to six steps of this pathway can be identified by sequence homology in individual archaeal genomes. However, there are no archaeal homologs to known isoforms of pantothenate synthetase (PS) or pantothenate kinase. Using comparative genomics, we previously identified two conserved archaeal protein families as the best candidates for the missing steps. Here we report the characterization of the predicted PS gene from Methanosarcina mazei, which encodes a hypothetical protein (MM2281) with no obvious homologs outside its own family. When expressed in Escherichia coli, MM2281 partially complemented an auxotrophic mutant without PS activity. Purified recombinant MM2281 showed no PS activity on its own, but the enzyme enabled substantial synthesis of [14C]4'-phosphopantothenate from [14C]beta-alanine, pantoate and ATP when coupled with E. coli pantothenate kinase. ADP, but not AMP, was detected as a coproduct of the coupled reaction. MM2281 also transferred the 14C-label from [14C]beta-alanine to pantothenate in the presence of pantoate and ADP, presumably through isotope exchange. No exchange took place when pantoate was removed or ADP replaced with AMP. Our results indicate that MM2281 represents a novel type of PS that forms ADP and is strongly inhibited by its product pantothenate. These properties differ substantially from those of bacterial PS, and may explain why PS genes, in contrast to other pantothenate biosynthetic genes, were not exchanged horizontally between the Bacteria and Archaea.     

Presenilins interact with armadillo proteins including neural-specific plakophilin-related protein and beta-catenin

Abstract : Missense substitutions in the presenilin 1 (PS1) and presenilin 2 (PS2) proteins are associated with early-onset familial Alzheimer's disease. We have used yeast-two-hybrid and coimmunoprecipitation methods to show that the large cytoplasmic loop domains of PS1 and PS2 interact specifically with three members of the armadillo protein family, including β-catenin, p0071, and a novel neuronal-specific armadillo protein—neural plakophilin-related armadillo protein (NPRAP). The PS1 : NPRAP interaction occurs between the arm repeats of NPRAP and residues 372-399 at the C-terminal end of the large cytoplasmic loop of PS1. The latter residues contain a single arm -like domain and are highly conserved in the presenilins, suggesting that they form a functional armadillo protein binding site for the presenilins.     

Conformational regulation of the EGFR kinase core by the juxtamembrane and C-terminal tail: a molecular dynamics study

The catalytic domain of epidermal growth factor receptor (EGFR) is activated by dimerization, which requires allosteric coupling between distal dimerization and catalytic sites. Although crystal structures of EGFR kinases, solved in various conformational states, have provided important insights into EGFR activation by dimerization, the atomic details of how dimerization signals are dynamically coupled to catalytic regions of the kinase core are not fully understood. In this study, we have performed unrestrained and targeted molecular dynamics simulations on the active and inactive states of EGFR, followed by principal component analysis on the simulated trajectories, to identify correlated motions in the EGFR kinase domain upon dimerization. Our analysis reveals that the conformational changes associated with the catalytic functions of the kinase core are highly correlated with motions in the juxtamembrane (JM) and C-terminal tail, two flexible structural elements that play an active role in EGFR kinase activation and dimerization. In particular, the opening and closing of the ATP binding lobe relative to the substrate binding lobe is highly correlated with motions in the JM and C-terminal tail, suggesting that ATP and substrate binding can be coordinated with dimerization through conformational changes in the JM and C-terminal tail. Our study pinpoints key residues involved in this conformational coupling, and provides new insights into the role of the JM and C-terminal tail segments in EGFR kinase functions.     

Gating at both ends and breathing in the middle: conformational dynamics of TolC

Drug extrusion via efflux through a tripartite complex (an inner membrane pump, an outer membrane protein, and a periplasmic protein) is a widely used mechanism in Gram-negative bacteria. The outer membrane protein (TolC in Escherichia coli; OprM in Pseudomonas aeruginosa) forms a tunnel-like pore through the periplasmic space and the outer membrane. Molecular dynamics simulations of TolC have been performed, and are compared to simulations of Y362F/R367S mutant, and to simulations of its homolog OprM. The results reveal a complex pattern of conformation dynamics in the TolC protein. Two putative gate regions, located at either end of the protein, can be distinguished. These regions are the extracellular loops and the mouth of the periplasmic domain, respectively. The periplasmic gate has been implicated in the conformational changes leading from the closed x-ray structure to a proposed open state of TolC. Between the two gates, a peristaltic motion of the periplasmic domain is observed, which may facilitate transport of the solutes from one end of the tunnel to the other. The motions observed in the atomistic simulations are also seen in coarse-grained simulations in which the protein tertiary structure is represented by an elastic network model.     

Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.     

Computer simulations of the OmpF porin from the outer membrane of Escherichia coli.

Molecular dynamics simulations were used to study the structure and dynamics of the Escherichia coli OmpF porin, which is composed of three identical 16-stranded beta-barrels. Simulations of the full trimer in the absence of water and the membrane led to significant contraction of the channel in the interior of each beta-barrel. With very weak harmonic constraints (0.005 kcal/mol A2/atom) applied to the main-chain C alpha atoms of the beta-barrel, the structure was stabilized without alteration of the average fluctuations. The resulting distribution of the fluctuations (small for beta-strands, large for loops and turns) is in good agreement with the x-ray B factors. Dynamic cross-correlation functions showed the importance of coupling between the loop motions and barrel flexibility. This was confirmed by the application of constraints corresponding to the observed temperature factors to the barrel C alpha atoms. With these constraints, the beta-barrel fluctuations were much smaller than the experimental values because of the intrinsic restrictions on the atomic motions, and the loop motions were reduced significantly. This result indicates that considerable care is required in introducing constraints to keep proteins close to the experimental structure during simulations, as has been done in several recent studies. Loop 3, which is thought to be important in gating the pore, undergoes a displacement that shifts it away from the x-ray structure. Analysis shows that this arises from the breakdown of a hydrogen bond network, which appears to result more from the absence of solvent that from the use of standard ionization states for the side chains of certain beta-barrel residues.     

Comparison of simulated and experimentally determined dynamics for a variant of the Lacl DNA-binding domain, Nlac-P.

Recent advances in the experimentally determined structures and dynamics of the domains within LacI provide a rare context for evaluating dynamics calculations. A 1500-ps trajectory was simulated for a variant of the LacI DNA-binding domain, which consists of the first three helices in LacI and the hinge helix of the homologous PurR. Order parameters derived from dynamics simulations are compared to those obtained for the LacI DNA-binding domain with 15N relaxation NMR spectroscopy (Slijper et al., 1997. Biochemistry. 36:249-254). The MD simulations suggest that the unstructured loop between helices II and III does not exist in a discrete state under the conditions of no salt and neutral pH, but occupies a continuum of states between the DNA-bound and free structures. Simulations also indicate that the unstructured region between helix III and the hinge helix is very mobile, rendering motions of the hinge helix essentially independent of the rest of the protein. Finally, the alpha-helical hydrogen bonds in the hinge helix are broken after 1250 ps, perhaps as a prelude to helix unfolding.     

Structure,Dynamics, and Allosteric Potential of Ionotropic Glutamate Receptor N-Terminal Domains

Ionotropic glutamate receptors (iGluRs) are tetrameric cation channels that mediate synaptic transmission and plasticity. They have a unique modular architecture with four domains: the intracellular C-terminal domain (CTD) that is involved in synaptic targeting, the transmembrane domain (TMD) that forms the ion channel, the membrane-proximal ligand-binding domain (LBD) that binds agonists such as L-glutamate, and the distal N-terminal domain (NTD), whose function is the least clear. The extracellular portion, comprised of the LBD and NTD, is loosely arranged, mediating complex allosteric regulation and providing a rich target for drug development. Here, we briefly review recent work on iGluR NTD structure and dynamics, and further explore the allosteric potential for the NTD in AMPA-type iGluRs using coarse-grained simulations. We also investigate mechanisms underlying the established NTD allostery in NMDA-type iGluRs, as well as the fold-related metabotropic glutamate and GABA_B receptors. We show that the clamshell motions intrinsically favored by the NTD bilobate fold are coupled to dimeric and higher-order rearrangements that impact the iGluR LBD and ultimately the TMD. Finally, we explore the dynamics of intact iGluRs and describe how it might affect receptor operation in a synaptic environment.     

Correlated motions in the U1 snRNA stem/loop 2:U1A RBD1 complex

The complex formed by U1A RBD1 and the U1 snRNA stem/loop II is noted for its high affinity and exquisite specificity. Here, that complex is investigated by 5 ns molecular dynamics simulations and analyzed by reorientational eigenmode dynamics to determine the dynamic properties of the RNA:protein interface that could contribute to the binding mechanism. The analysis shows that there is extensive correlation between motions of the RNA and protein, involving 7 of the 10 RNA loop nucleotides, the protein beta-sheet surface, two of its loops, and its C-terminal tripeptide sequence. Order parameters of these regions of the complex are uniformly high, indicating restricted motion. However, several regions of both RNA and protein retain local flexibility, notably three nucleotides of the RNA loop and one loop of RBD1 that does not contact RNA. The highly correlated motions involving both molecules reflect the intricate network of interactions that characterize this complex and could account in part for the thermodynamic coupling observed for complex formation.