ATP-binding cassette (ABC) transporters are ubiquitous in all three kingdoms of life and are implicated in many clinically relevant physiological processes. They couple the energy released by ATP hydrolysis to facilitate substrate translocation across cell membranes. The crystal structures of type II ABC importers have revealed their unique transmembrane domain architecture consisting of 10 transmembrane helices and their structurally conserved nucleotide-binding domains among all ABC transporters. However, molecular details of the interactions between the importers and their substrate remain largely elusive. Taking vitamin B12 importer BtuCD as an exemplar of type II importers, we investigated the dynamics of its occluded state and the detailed protein–substrate interactions using molecular dynamics simulation. Our trajectories show that the importer accommodates the substrate through a nonspecific binding mode as the substrate undergoes evident vertical and tilt motions inside the translocation cavity. Extensive hydrogen bond and hydrophobic interactions were observed between the substrate and the importer; however, most of these interactions are weak, with <38% occurrence. The presence of substrate leads to enlargement of the translocation cavity, especially at its cytoplasmic end, which may activate cytoplasmic regions and probably facilitate the transportation. The perturbations caused by periplasmic binding protein and nucleotides were also investigated. The study provides deeper insight into the translocation mechanism of BtuCD.

Maltose transporter MalFGK2 is a type-I importer in the ATP-binding cassette (ABC) transporter superfamily. Upon the binding of its periplasmic binding protein, MalE, the ATPase activity of MalFGK2 can be greatly enhanced. Crystal structures of the MalFGK2–MalE–maltose complex in a so-called “pretranslocation” (“pre-T”) state with a partially closed conformation suggest that the formation of this MalE-stabilized intermediate state is a key step leading to the outward-facing catalytic state. On the contrary, crosslinking and fluorescence studies suggest that ATP binding alone is sufficient to promote the outward-facing catalytic state, thereby doubting the role of MalE binding. To clarify the role of MalE binding and to gain deeper understanding of the molecular mechanisms of MalFGK2, we calculated the free energy surfaces (FESs) related to the lateral motion in the presence and absence of MalE using atomistic metadynamics simulations. The results showed that, in the absence of MalE, laterally closing motion was energetically forbidden but, upon MalE binding, more closed conformations similar to the pre-T state become more stable. The significant effect of MalE binding on the free energy landscapes was in agreement with crystallographic studies and confirmed the important role of MalE in stabilizing the pre-T state. Our simulations also revealed that the allosteric effect of MalE stimulation originates from the MalE-binding-promoted vertical motion between MalF and MalG cores, which was further supported by MD simulation of the MalE-independent mutant MalF500.

AbstractProtein interactions involving intrinsically disordered proteins (IDPs) comprise a variety of binding modes, from the well-characterized folding upon binding to dynamic fuzzy complexes. To date, most studies concern the binding of an IDP to a structured protein, while the interaction between two IDPs is poorly understood. In this study, NMR, smFRET, and molecular dynamics (MD) simulation are combined to characterize the interaction between two IDPs, the C-terminal domain (CTD) of protein 4.1G and the nuclear mitotic apparatus (NuMA) protein. It is revealed that CTD and NuMA form a fuzzy complex with remaining structural disorder. Multiple binding sites on both proteins were identified by molecular dynamics and mutagenesis studies. This study provides an atomic scenario in which two IDPs bearing multiple binding sites interact with each other in dynamic equilibrium. The combined approach employed here could be widely applicable for investigating IDPs and their dynamic interactions.

Peptide modification of nanoparticles with high efficiency is critical in determining the properties and bioapplications of nanoparticles, but the methodology remains a challenging task. Here, by using the phosphorylated linear and cyclic peptide with the arginine–glycine–aspartic acid (RGD) targeting motifs as typical examples, the peptides binding efficiency for the inorganic metal compound nanoparticles was increased significantly after the phosphorylation treatment, and the modification allowed for improving the selectivity and signal-to-noise ratio for cancer targeting and reduced the toxicity derived from nonspecific interactions of nanoparticles with cells owing to the higher amount of phosphopeptide binding. In addition, molecular dynamics (MD) simulations of various peptides on inorganic metal compound surfaces revealed that the peptide adsorption on the surface is mainly driven by electrostatic interactions between phosphate oxygen and the polarized interfacial water layer, consistent with the experimental observation of the strong binding propensity of phosphorylated peptides. Significantly, with the RGD phosphopeptide surface modification, these nanoparticles provide a versatile tool for tuning material–cell interactions to achieve the desired level of autophagy and may prove useful for various diagnostic and therapeutic applications.Keywords: autophagy; cancer target; lanthanide upconversion nanoparticles; peptide; phosphorylation

The resistance-nodulation-cell division transporter AcrB is responsible for energy transduction and substrate recognition in the tripartite AcrAB-TolC efflux system in Escherichia coli. Despite a broad substrate specificity, only a few compounds have been cocrystallized with AcrB inside the distal binding pocket (DBP), including doxorubicin (DOX) and D13-9001. D13-9001 is a promising efflux pump inhibitor that potentiates the efficacy of a wide variety of antibiotics. To understand its inhibition effect under the framework of functional rotating mechanism, we performed targeted and steered molecular dynamics simulations to compare the binding and extrusion processes of this inhibitor and the substrate DOX in AcrB. The results demonstrate that, with respect to DOX, the interaction of D13-9001 with the hydrophobic trap results in delayed disassociation from the DBP. Notably, the detachment of D13-9001 is tightly correlated with the side-chain reorientation of Phe628 and large-scale displacement of Tyr327. Furthermore, the inhibitor induces much more significant conformational changes at the exit gate than DOX does, thereby causing higher energy cost for extrusion and contributing to the inhibitory effect in addition to the tight binding at DBP.

The conformation and subcellular localization of R-SNARE protein Ykt6 are regulated by the lipidation state of its C-terminal CCAIM motif. Biochemical and crystallography studies showed that lipid molecules binding at a hydrophobic pocket at the interface between the longin domain and the SNARE core can lock Ykt6 at a closed conformation and mimic the farnesylated state of Ykt6. In this study, we performed in silico farnesylation of Ykt6 and explored the conformational dynamics of Ykt6 using conventional and steered MD simulations. We found that the farnesylated Ykt6 model structure is stable during the 2 μs simulation and the farnesyl group adopts conformations similar to those of the DPC molecule bound to Ykt6. Both DPC binding and farnesylation were found to reduce the conformational flexibility of Ykt6 and hinder the dissociation of SNARE core from the longin domain. The dissociation of the αF-αG segment is the rate-limiting step during the putative closed-to-open conformational transition of Ykt6, and the key residues involved in this process are consistent with the experimental mutagenesis study.

The oxidative dehydrogenation (ODH) of propane to propene over a vanadium-based catalyst suffers from side reactions of further and complete oxidations of propene and other intermediates, which limit the yield of propene. These further oxidation reactions are also referred to as deep oxidation reactions of ODH. In this paper, we present a comprehensive study of the deep oxidation reactions in the ODH of propane over the V2O5(001) surface using the periodic density functional theory method. It is shown that the main source of deep oxidation byproducts originates from the dehydrogenation reaction of the surface intermediate isopropoxide, leading to acetone and the following deep oxidation reactions of acetone and propene. Thorough oxidation of acetone is more difficult than that of propene. Noticeably, formation of acetone and deep oxidation of acetone and propene are only feasible on the terminal oxygen site O(1) of the V2O5(001) surface. The bridging site O(2) has similar reactivity for propene formation but is inert for the side reactions, showing its superiority for selectivity of propane ODH.

The membrane fusion protein (MFP) AcrA is proposed to link the inner membrane transporter AcrB and outer membrane protein TolC, forming the tripartite AcrAB–TolC efflux pump, and was shown to be functionally indispensible. Structural and EPR studies showed that AcrA has high conformational flexibility and exhibited pH-induced conformational change. In this study, we built the complete structure of AcrA through homology modeling and performed atomistic simulations of AcrA at different pH values. It was shown that the conformational flexibility of AcrA originates from the motions of α-hairpin and MP domains. The conformational dynamics of AcrA is sensitive to specific point mutations and pH values. In agreement with the EPR experiments, the interdomain motions were restrained upon lowering pH from 7.0 to 5.0 in the simulations. It was found that the protonation/deprotonation of His285 underlies the pH-regulated conformational dynamics of AcrA by disturbing the local hydrogen bond interactions, suggesting that the changes of pH in the periplasm accompanying the drug efflux could act as a signal to trigger the action of AcrA, which undergoes reversible conformational rearrangement.

Dynein light chain LC8 is a highly conserved, dimeric protein involved in a variety of essential cellular events. Phosphorylation at Ser88 was found to promote mammalian cell survival and regulate the dimer to monomer transition at physiological pH. Combining molecular dynamics (MD) simulation and free energy calculation methods, we explored the atomistic mechanism of the phosphorylation-induced dimer dissociation. The MD simulation revealed that phosphorylation/phosphomimetic mutation at Ser88 opens an entrance into the dimer interface for water molecules, which disturb the hydrogen bond network around His55 and is expected to raise the pKa value and protonation ratio of His55 as well. The free energy calculations showed that the S88E mutation destabilized the dimer by 6.6 kcal/mol, in good agreement with the experimental value of 8.1 kcal/mol. The calculated destabilization upon phosphorylation is 50.8 kcal/mol, showing that phosphorylation definitely prevents dimer formation under physiological conditions. Further analysis of the calculated free energy changes demonstrated that the electrostatic contribution dominates the impact of phosphorylation on dimer dissociation.

The first report on the transferable, plasmid-mediated quinolone-resistance determinant qnrA1 was in 1998. Since then, qnr alleles have been discovered worldwide in clinical strains of Gram-negative bacilli. Qnr proteins confer quinolone resistance, and belong to the pentapeptide repeat protein (PRP) family. Several PRP crystal structures have been solved, but little is known about the functional significance of their structural arrangement.We conducted random and site-directed mutagenesis on qnrA1 and on qnrC, a newly identified quinolone-resistance gene from Proteus mirabilis. Many of the Qnr mutants lost their quinolone resistance function. The highly conserved hydrophobic Leu or Phe residues at the center of the pentapeptide repeats are known as i sites, and loss-of-function mutations included replacement of the i site hydrophobic residues with charged residues, replacing the i-2 site, N-terminal to the i residues, with bulky side-chain residues, introducing Pro into the β-helix coil, deletion of the N- and C-termini, and excision of a central coil. Molecular dynamics simulations and homology modeling demonstrated that QnrC overall adopts a stable β-helix fold and shares more similarities with MfpA than with other PRP structures. Based on homology modeling and molecular dynamics simulation, the dysfunctional point mutations introduced structural deformations into the quadrilateral β-helix structure of PRPs. Of the pentapeptides of QnrC, two-thirds adopted a type II β-turn, while the rest adopted type IV turns. A gap exists between coil 2 and coil 3 in the QnrC model structure, introducing a structural flexibility that is similar to that seen in MfpA.The hydrophobic core and the β-helix backbone conformation are important for maintaining the quinolone resistance property of Qnr proteins. QnrC may share structural similarity with MfpA.

The two N-terminal PDZ domains of postsynaptic density protein-95 (PDS-95 PDZ1 and PDZ2) are closely connected in tandem by a conserved peptide linker of five amino acids. The interdomain orientation between PDZ1 and PDZ2 of the ligand-free PDZ12 tandem is restrained, and this conformational arrangement facilitates the synergistic binding of PDZ12 to multimeric targets.(1) The interdomain orientation of the target-bound state of PDZ12 is not known. Here, we have solved the structure of PDZ12 in complex with its binding domain from cypin. Both chemical shift data and residual dipolar coupling measurements showed that the restrained interdomain orientation disappeared upon cypin peptide binding. NMR-based relaxation experiments revealed slow interdomain motions in the PDZ12/cypin peptide complex. Molecular dynamics simulations also showed that the PDZ12/cypin complex has larger conformational flexibility than the ligand-free PDZ12. This dramatic change of protein dynamics provides extra conformational entropy upon ligand binding, thus enhancing the ligand binding affinity of the PDZ12 tandem. Modulation of ligand binding affinity through concerted interdomain structural and dynamic rearrangements may represent a general property of multidomain scaffold proteins.

The parallel π–π stacking of benzene molecules was studied using Density Functional Theory (DFT) and second-order Moller–Plesset perturbation theory (MP2) methods. The DFT methods are proved to be inadequate in prediction of π–π stacking conformation and interaction energy. Cluster model calculations at the MP2/6-311 + G** level predicted an optimized conformation which is very close to the structure of the parallel-displaced benzene dimer. The calculated inter-plane distance of 3.3 Å is in good agreement with the observation in organic molecular crystals. The interaction energy predicted at MP2 level revealed that the pairwise interaction energy increases with the number of the parallel-stacked benzene molecules.

Asymmetric cell division requires the establishment of cortical cell polarity and the orientation of the mitotic spindle along the axis of cell polarity. Evidence from invertebrates demonstrates that the Par3/Par6/aPKC and NuMA/LGN/Gαi complexes, which are thought to be physically linked by the adaptor protein mInscuteable (mInsc), play indispensable roles in this process. However, the molecular basis for the binding of LGN to NuMA and mInsc is poorly understood. The high-resolution structures of the LGN/NuMA and LGN/mInsc complexes presented here provide mechanistic insights into the distinct and highly specific interactions of the LGN TPRs with mInsc and NuMA. Structural comparisons, together with biochemical and cell biology studies, demonstrate that the interactions of NuMA and mInsc with LGN are mutually exclusive, with mInsc binding preferentially. Our results suggest that the Par3/mInsc/LGN and NuMA/LGN/Gαi complexes play sequential and partially overlapping roles in asymmetric cell division.Graphical AbstractDownload high-res image (240KB)Download full-size imageHighlights► Structures of LGN/NuMA and LGN/mInsc complexes are solved ► mInsc and NuMA compete for LGN binding ► mInsc can displace NuMA from LGN ► mInsc cannot couple the Par3/Par6/aPKC and LGN/NuMA/Gαi complexes together physically

Putative metal-chelate-type ABC transporter HI1470/1 is homologous with vitamin B12 importer BtuCD but exhibits a distinct inward-facing conformation in contrast to the outward-facing conformation of BtuCD. Normal-mode analysis of HI1470/1 reveals the intrinsic asymmetric conformational flexibility in this transporter and demonstrates that the transition from the inward-facing to the outward-facing conformation is realized through the asymmetric motion of individual subunits of the transporter. This analysis suggests that the asymmetric arrangement of the BtuC dimer in the crystal structure of the BtuCD-F complex represents an intermediate state relating HI1470/1 and BtuCD. Furthermore, a twisting motion between transmembrane domains and nucleotide-binding domains encoded in the lowest-frequency normal mode of this type of importer is found to contribute to the conformational transitions during the whole cycle of substrate transportation. A more complete translocation mechanism of the BtuCD type importer is proposed.

•Two Goloco motifs connected in tandem form a TPR repeat binding unit•The structure of the autoinhibited conformation of LGN is solved•GL motifs bind to TPR repeats with a mode distinct from that in GL/Gα complexes•LGN responds to G protein in a receptor-independent mannerLGN plays essential roles in asymmetric cell divisions via its N-terminal TPR-motif-mediated binding to mInsc and NuMA. This scaffolding activity requires the release of the autoinhibited conformation of LGN by binding of Gαi to its C-terminal GoLoco (GL) motifs. The interaction between the GL and TPR motifs of LGN represents a distinct GL/target binding mode with an unknown mechanism. Here, we show that two consecutive GL motifs of LGN form a minimal TPR-motif-binding unit. GL12 and GL34 bind to TPR0–3 and TPR4–7, respectively. The crystal structure of a truncated LGN reveals that GL34 forms a pair of parallel α helices and binds to the concave surface of TPR4–7, thereby preventing LGN from binding to other targets. Importantly, the GLs bind to TPR motifs with a mode distinct from that observed in the GL/Gαi·GDP complexes. Our results also indicate that multiple and orphan GL motif proteins likely respond to G proteins with distinct mechanisms.

ATP-binding cassette transporter BtuCD mediating vitamin B12 uptake in Escherichia coli couples the energy of ATP hydrolysis to the translocation of vitamin B12 across the membrane into the cell. Elastic normal mode analysis of BtuCD demonstrates that the simultaneous substrate trapping at periplasmic cavity and ATP binding at the ATP-binding cassette (BtuD) dimer proceeds readily along the lowest energy pathway. The transport power stroke is attributed to ATP-hydrolysis-induced opening of the nucleotide-binding domain dimer, which is coupled to conformational rearrangement of transmembrane domain (BtuC) helices leading to the closing at the periplasmic side and opening at the cytoplasmic gate. Simultaneous hydrolysis of two ATP is supported by the fact that antisymmetric movement of BtuD dimer implying alternating hydrolysis cannot induce effective conformational change of the translocation pathway. A plausible mechanism of translocation cycle is proposed in which the possible effect of the association of periplasmic binding protein BtuF to the transporter is also considered.

Netrin receptor DCC plays critical roles in many cellular processes, including axonal outgrowth and migration, angiogenesis, and apoptosis, but the molecular basis of DCC-mediated signaling is largely unclear. ERK2, a member of the MAPK family, is one of the few proteins known to be involved in DCC-mediated signaling. Here, we report that ERK2 directly interacts with DCC, and the ERK2-binding region was mapped to the conserved intracellular P1 domain of the receptor. The structure of ERK2 in complex with the P1 domain of DCC reveals that DCC contains a MAPK docking motif. The docking of the P1 domain onto ERK2 physically positions several phosphorylation sites of DCC in the vicinity of the kinase active site. We further show that the docking interaction between the P1 domain and ERK2 is essential for the ERK2-mediated phosphorylation of DCC. We conclude that DCC signaling is directly coupled with MAPK signaling cascades.Graphical AbstractDownload high-res image (320KB)Download full-size imageHighlights► The cytoplasmic domain of the netrin receptor DCC physically binds to ERK2 ► The P1 domain contains an ERK2-docking motif ► The structure of ERK2/DCC P1 complex reveals the detailed interaction mode of the interaction ► ERK2 directly phosphorylate DCC

Maltose transporter MalFGK2 is a type-I importer in the ATP-binding cassette (ABC) transporter superfamily. Upon the binding of its periplasmic binding protein, MalE, the ATPase activity of MalFGK2 can be greatly enhanced. Crystal structures of the MalFGK2–MalE–maltose complex in a so-called “pretranslocation” (“pre-T”) state with a partially closed conformation suggest that the formation of this MalE-stabilized intermediate state is a key step leading to the outward-facing catalytic state. On the contrary, crosslinking and fluorescence studies suggest that ATP binding alone is sufficient to promote the outward-facing catalytic state, thereby doubting the role of MalE binding. To clarify the role of MalE binding and to gain deeper understanding of the molecular mechanisms of MalFGK2, we calculated the free energy surfaces (FESs) related to the lateral motion in the presence and absence of MalE using atomistic metadynamics simulations. The results showed that, in the absence of MalE, laterally closing motion was energetically forbidden but, upon MalE binding, more closed conformations similar to the pre-T state become more stable. The significant effect of MalE binding on the free energy landscapes was in agreement with crystallographic studies and confirmed the important role of MalE in stabilizing the pre-T state. Our simulations also revealed that the allosteric effect of MalE stimulation originates from the MalE-binding-promoted vertical motion between MalF and MalG cores, which was further supported by MD simulation of the MalE-independent mutant MalF500.