Normal mode based flexible fitting of high-resolution structure into low-resolution experimental data from cryo-EM

A new method for the flexible fitting of high-resolution structures into low-resolution maps of macromolecular complexes from electron microscopy has been recently described in applications to simulated electron density maps. This method uses a linear combination of low-frequency normal modes in an iterative manner to deform the structure optimally to conform to the low-resolution electron density map. Gradient-following techniques in the coordinate space of collective normal modes are used to optimize the overall correlation coefficient between computed and measured electron densities. With this approach, multi-scale flexible fitting can be performed using all-atoms or Calpha atoms. In this paper, illustrative studies of normal mode based flexible fitting to experimental cryo-EM maps are presented for three different systems. Large, functionally relevant conformational changes for elongation factor G bound to the ribosome, Escherichia coli RNA polymerase and cowpea chlorotic mottle virus are elucidated as the result of the application of NMFF from high-resolution structures to cryo-electron microscopy maps.     

Keap1 in adhesion complexes

Cell adhesion complexes are sensors that interact with the extracellular environment and allow for the transmission of signals found outside the cell across the plasma membrane to the cell interior. Keap1 is a newly identified component of cell adhesion complexes. We investigated Keap1's association with these complexes in diverse tissues and cell types. Keap1 is present in focal adhesion (FA)-like assemblies in kidney proximal tubule cells where it colocates with actin. In liver, Keap1 is found in the adherens junctions (AJ) and at the base of the bile canaliculi. To study Keap1's involvement in both the integrin-based FA and the cadherin-based AJ, we induced formation of these complexes in fibroblasts, using a serum starvation followed by a serum supplementation method. When compared with vinculin, a component of all FA, we found that Keap1 assembles only in the peripheral FA. Within the peripheral FA, Keap1 was present in distinct foci along the length of the FA and these foci were different from vinculin, talin, paxillin, and phospho-tyrosine rich regions of the FA. Unlike most FA components, Keap1 was also recruited to the newly formed AJ. As Keap1 homologues are actin-bundling proteins, we hypothesize that Keap1's function is to bundle F-actin within these diverse types of cell adhesion components.     

Cryo-Cooling Effect on DHFR Crystal Studied by Replica-Exchange Molecular Dynamics Simulations

Cryo-cooling is routinely performed before x-ray diffraction image collection to reduce the damage to crystals due to ionizing radiation. It has been suggested that although backbone structures are usually very similar between room temperature and cryo-temperature, cryo-cooling may hamper biologically relevant dynamics. In this study, the crystal of Escherichia coli dihydrofolate reductase is studied with replica-exchange molecular dynamics simulation, and the results are compared with the crystal structure determined at cryo-temperature and room temperature with the time-averaged ensemble method. Although temperature dependence of unit cell compaction and root mean-square fluctuation of Cα is found in accord with experiment, it is found that the protein structure at low temperature can be more heterogeneous than the ensemble of structures reported by using the time-averaged ensemble method, encouraging further development of the time-averaged ensemble method and indicating that data should be examined carefully to avoid overinterpretation of one average structure.     

Replica Exchange Molecular Dynamics Simulations Provide Insight into Substrate Recognition by Small Heat Shock Proteins

The small heat shock proteins (sHSPs) are a virtually ubiquitous and diverse group of molecular chaperones that can bind and protect unfolding proteins from irreversible aggregation. It has been suggested that intrinsic disorder of the N-terminal arm (NTA) of sHSPs is important for substrate recognition. To investigate conformations of the NTA that could recognize substrates we performed replica exchange molecular dynamics simulations. Behavior at normal and stress temperatures of the dimeric building blocks of dodecameric HSPs from wheat (Ta16.9) and pea (Ps18.1) were compared because they display high sequence similarity, but Ps18.1 is more efficient in binding specific substrates. In our simulations, the NTAs of the dimer are flexible and dynamic; however, rather than exhibiting highly extended conformations they retain considerable α-helical character and contacts with the conserved α-crystallin domain (ACD). Network analysis and clustering methods reveal that there are two major conformational forms designated either “open” or “closed” based on the relative position of the two NTAs and their hydrophobic solvent accessible surface area. The equilibrium constant for the closed to open transition is significantly different for Ta16.9 and Ps18.1, with the latter showing more open conformations at elevated temperature correlated with its more effective chaperone activity. In addition, the Ps18.1 NTAs have more hydrophobic solvent accessible surface than those of Ta16.9. NTA hydrophobic patches are comparable in size to the area buried in many protein-protein interactions, which would enable sHSPs to bind early unfolding intermediates. Reduced interactions of the Ps18.1 NTAs with each other and with the ACD contribute to the differences in dynamics and hydrophobic surface area of the two sHSPs. These data support a major role for the conformational equilibrium of the NTA in substrate binding and indicate features of the NTA that contribute to sHSP chaperone efficiency.     

Simulations of substrate transport in the multidrug transporter EmrD

EmrD is a multidrug resistance (MDR) transporter from Escherichia coli, which is involved in the efflux of amphipathic compounds from the cytoplasm, and the first MDR member of the major facilitator superfamily to be crystallized. Molecular dynamics simulation of EmrD in a phospholipid bilayer was used to characterize the conformational dynamics of the protein. Motions that support a previously proposed lateral diffusion pathway for substrate from the cytoplasmic membrane leaflet into the EmrD central cavity were observed. In addition, the translocation pathway of meta-chloro carbonylcyanide phenylhydrazone (CCCP) was probed using both standard and steered molecular dynamics simulation. In particular, interactions of a few specific residues with CCCP have been identified. Finally, a large motion of two residues, Val 45 and Leu 233, was observed with the passage of CCCP into the periplasmic space, placing a lower bound on the extent of opening required at this end of the protein for substrate transport. Overall, our simulations probe details of the transport pathway, motions of EmrD at an atomic level of detail, and offer new insights into the functioning of MDR transporters.     

Twelve Transmembrane Helices Form the Functional Core of Mammalian MATE1 (Multidrug and Toxin Extruder 1) Protein

The x-ray structure of the prototypic MATE family member, NorM from Vibrio cholerae, reveals a protein fold composed of 12 transmembrane helices (TMHs), confirming hydropathy analyses of the majority of (prokaryotic and plant) MATE transporters. However, the mammalian MATEs are generally predicted to have a 13(th) TMH and an extracellular C terminus. Here we affirm this prediction, showing that the C termini of epitope-tagged, full-length human, rabbit, and mouse MATE1 were accessible to antibodies from the extracellular face of the membrane. Truncation of these proteins at or near the predicted junction between the 13(th) TMH and the long cytoplasmic loop that precedes it resulted in proteins that (i) trafficked to the membrane and (ii) interacted with antibodies only after permeabilization of the plasma membrane. CHO cells expressing rbMate1 truncated at residue Gly-545 supported levels of pH-sensitive transport similar to that of cells expressing the full-length protein. Although the high transport rate of the Gly-545 truncation mutant was associated with higher levels of membrane expression (than full-length MATE1), suggesting the 13(th) TMH may influence substrate translocation, the selectivity profile of the mutant indicated that TMH13 has little impact on ligand binding. We conclude that the functional core of MATE1 consists of 12 (not 13) TMHs. Therefore, we used the x-ray structure of NorM to develop a homology model of the first 12 TMHs of MATE1. The model proved to be stable in molecular dynamic simulations and agreed with topology evident from preliminary cysteine scanning of intracellular versus extracellular loops.     

Allosteric Regulation of DNA Cleavage and Sequence-Specificity through Run-On Oligomerization

1. Download : Download high-res image (402KB)
2. Download : Download full-size image

     

Electrostatic properties of cowpea chlorotic mottle virus and cucumber mosaic virus capsids

Electrostatic properties of cowpea chlorotic mottle virus (CCMV) and cucumber mosaic virus (CMV) were investigated using numerical solutions to the Poisson-Boltzmann equation. Experimentally, it has been shown that CCMV particles swell in the absence of divalent cations when the pH is raised from 5 to 7. CMV, although structurally homologous, does not undergo this transition. An analysis of the calculated electrostatic potential confirms that a strong electrostatic repulsion at the calcium-binding sites in the CCMV capsid is most likely the driving force for the capsid swelling process during the release of calcium. The binding interaction between the encapsulated genome material (RNA) inside of the capsid and the inner capsid shell is weakened during the swelling transition. This probably aids in the RNA release process, but it is unlikely that the RNA is released through capsid openings due to unfavorable electrostatic interaction between the RNA and capsid inner shell residues at these openings. Calculations of the calcium binding energies show that Ca(2+) can bind both to the native and swollen forms of the CCMV virion. Favorable binding to the swollen form suggests that Ca(2+) ions can induce the capsid contraction and stabilize the native form.     

A Myo6 mutation destroys coordination between the myosin heads, revealing new functions of myosin VI in the stereocilia of mammalian inner ear hair cells

Myosin VI, found in organisms from Caenorhabditis elegans to humans, is essential for auditory and vestibular function in mammals, since genetic mutations lead to hearing impairment and vestibular dysfunction in both humans and mice. Here, we show that a missense mutation in this molecular motor in an ENU-generated mouse model, Tailchaser, disrupts myosin VI function. Structural changes in the Tailchaser hair bundles include mislocalization of the kinocilia and branching of stereocilia. Transfection of GFP-labeled myosin VI into epithelial cells and delivery of endocytic vesicles to the early endosome revealed that the mutant phenotype displays disrupted motor function. The actin-activated ATPase rates measured for the D179Y mutation are decreased, and indicate loss of coordination of the myosin VI heads or ‘gating’ in the dimer form. Proper coordination is required for walking processively along, or anchoring to, actin filaments, and is apparently destroyed by the proximity of the mutation to the nucleotide-binding pocket. This loss of myosin VI function may not allow myosin VI to transport its cargoes appropriately at the base and within the stereocilia, or to anchor the membrane of stereocilia to actin filaments via its cargos, both of which lead to structural changes in the stereocilia of myosin VI–impaired hair cells, and ultimately leading to deafness.