Dynamics of the trimeric AcrB transporter protein inferred from a B-factor analysis of the crystal structure

An Ion-Pair Induced Intermediate Complex Captured in Class D Carbapenemase Reveals Chloride Ion as a Janus Effector Modulating Activity

Antibiotic-resistant Enterobacterales that produce oxacillinase (OXA)-48-like Class D β-lactamases are often linked to increased clinical mortality. Though the catalytic mechanism of OXA-48 is known, the molecular origin of its biphasic kinetics has been elusive. We here identify selective chloride binding rather than decarbamylation of the carbamylated lysine as the source of biphasic kinetics, utilizing isothermal titration calorimetry (ITC) to monitor the complete reaction course with the OXA-48 variant having a chemically stable N-acetyl lysine. Further structural investigation enables us to capture an unprecedented inactive acyl intermediate wedged in place by a halide ion paired with a conserved active site arginine. Supported by mutagenesis and mathematical simulation, we identify chloride as a "Janus effector" that operates by allosteric activation of the burst phase and by inhibition of the steady state in kinetic assays of β-lactams. We show that chloride-induced biphasic kinetics directly affects antibiotic efficacy and facilitates the differentiation of clinical isolates encoding Class D from Class A and B carbapenemases. As chloride is present in laboratory and clinical procedures, our discovery greatly expands the roles of chloride in modulating enzyme catalysis and highlights its potential impact on the pharmacokinetics and efficacy of antibiotics during in vivo treatment.

Competition-cooperation mechanism between Escherichia coli and Staphylococcus aureus based on systems mapping

Introduction

Interspecies interactions are a crucial driving force of species evolution. The genes of each coexisting species play a pivotal role in shaping the structure and function within the community, but how to identify them at the genome-wide level has always been challenging.

Methods

In this study, we embed the Lotka-Volterra ordinary differential equations in the theory of community ecology into the systems mapping model, so that this model can not only describe how the quantitative trait loci (QTL) of a species directly affects its own phenotype, but also describe the QTL of the species how to indirectly affect the phenotype of its interacting species, and how QTL from different species affects community behavior through epistatic interactions.

Results

By designing and implementing a co-culture experiment for 100 pairs of Escherichia coli (E. coli) and Staphylococcus aureus (S. aureus), we mapped 244 significant QTL combinations in the interaction process of the two bacteria using this model, including 69 QTLs from E. coli and 59 QTLs from S. aureus, respectively. Through gene annotation, we obtained 57 genes in E. coli, among which the genes with higher frequency were ypdC, nrfC, yphH, acrE, dcuS, rpnE, and ptsA, while we obtained 43 genes in S. aureus, among which the genes with higher frequency were ebh, SAOUHSC_00172, capF, gdpP, orfX, bsaA, and phnE1.

Discussion

By dividing the overall growth into independent growth and interactive growth, we could estimate how QTLs modulate interspecific competition and cooperation. Based on the quantitative genetic model, we can obtain the direct genetic effect, indirect genetic effect, and genome-genome epistatic effect related to interspecific interaction genes, and then further mine the hub genes in the QTL networks, which will be particularly useful for inferring and predicting the genetic mechanisms of community dynamics and evolution. Systems mapping can provide a tool for studying the mechanism of competition and cooperation among bacteria in co-culture, and this framework can lay the foundation for a more comprehensive and systematic study of species interactions.

A Combined DFT, Energy Decomposition, and Data Analysis Approach to Investigate the Relationship Between Noncovalent Interactions and Selectivity in a Flexible DABCOnium/Chiral Anion Catalyst System

Developing strategies to study reactivity and selectivity in flexible catalyst systems has become an important topic of research. Herein, we report a combined experimental and computational study aimed at understanding the mechanistic role of an achiral DABCOnium cofactor in a regio- and enantiodivergent bromocyclization reaction. It was found that electron-deficient aryl substituents enable rigidified transition states via an anion-π interaction with the catalyst, which drives the selectivity of the reaction. In contrast, electron-rich aryl groups on the DABCOnium result in significantly more flexible transition states, where interactions between the catalyst and substrate are more important. An analysis of not only the lowest-energy transition state structures but also an ensemble of low-energy transition state conformers via energy decomposition analysis and machine learning was crucial to revealing the dominant noncovalent interactions responsible for observed changes in selectivity in this flexible system.

Reviving B-Factors: Activating ALK Mutations Increase Protein Dynamics of the Unphosphorylated Kinase

Anaplastic lymphoma kinase (ALK) is a receptor tyrosine kinase that can become oncogenic by activating mutations or overexpression. Full kinetic characterization of both phosphorylated and nonphosphorylated wildtype and mutant ALK kinase domain was done. Our structure-based drug design programs directed at ALK allowed us to interrogate whether X-ray crystallography data could be used to support the hypothesis that activation of ALK by mutation occurs due to increased protein dynamics. Crystallographic B-factors were converted to normalized B-factors, which allowed analysis of wildtype ALK, ALK-C1156Y, and ALK-L1196M. This data suggests that mobility of the P-loop, αC-helix, and activation loop (A-loop) may be important in catalytic activity increases, with or without phosphorylation. Both molecular dynamics simulations and hydrogen-deuterium exchange experimental data corroborated the normalized B-factors data.

X-ray refinement significantly underestimates the level of microscopic heterogeneity in biomolecular crystals

Biomolecular X-ray structures typically provide a static, time- and ensemble-averaged view of molecular ensembles in crystals. In the absence of rigid-body motions and lattice defects, B-factors are thought to accurately reflect the structural heterogeneity of such ensembles. In order to study the effects of averaging on B-factors, we employ molecular dynamics simulations to controllably manipulate microscopic heterogeneity of a crystal containing 216 copies of villin headpiece. Using average structure factors derived from simulation, we analyse how well this heterogeneity is captured by high-resolution molecular-replacement-based model refinement. We find that both isotropic and anisotropic refined B-factors often significantly deviate from their actual values known from simulation: even at high 1.0 Å resolution and Rfree of 5.9%, B-factors of some well-resolved atoms underestimate their actual values even sixfold. Our results suggest that conformational averaging and inadequate treatment of correlated motion considerably influence estimation of microscopic heterogeneity via B-factors, and invite caution in their interpretation.

An Analysis of MIF Structural Features that Control Functional Activation of CD74

For more than 15 years, the tautomerase active site of macrophage migration inhibitory factor (MIF) and its catalytic residue Pro1 have been being targeted for the development of therapeutics that block activation of its cell surface receptor, CD74. Neither the biological role of the MIF catalytic site nor the mechanistic details of CD74 activation are well understood. The inherently unstable structure of CD74 remains the biggest obstacle in structural studies with MIF for understanding the basis of CD74 activation. Using a novel approach, we elucidate the mechanistic details that control activation of CD74 by MIF surface residues and identify structural parameters of inhibitors that reduce CD74 biological activation. We also find that N-terminal mutants located deep in the catalytic site affect surface residues immediately outside the catalytic site, which are responsible for reduction of CD74 activation.

Functional rotation induced by alternating protonation states in the multidrug transporter AcrB: all-atom molecular dynamics simulations

The multidrug transporter AcrB actively exports a wide variety of noxious compounds using proton-motive force as an energy source in Gram-negative bacteria. AcrB adopts an asymmetric structure comprising three protomers with different conformations that are sequentially converted during drug export; these cyclic conformational changes during drug export are referred to as functional rotation. To investigate functional rotation driven by proton-motive force, all-atom molecular dynamics simulations were performed. Using different protonation states for the titratable residues in the middle of the transmembrane domain, our simulations revealed the correlation between the specific protonation states and the side-chain configurations. Changing the protonation state for Asp408 induced a spontaneous structural transition, which suggests that the proton translocation stoichiometry may be one proton per functional rotation cycle. Furthermore, our simulations demonstrate that alternating the protonation states in the transmembrane domain induces functional rotation in the porter domain, which is primarily responsible for drug transport.

Resolution of structural heterogeneity in dynamic crystallography

Dynamic behavior of proteins is critical to their function. X-ray crystallography, a powerful yet mostly static technique, faces inherent challenges in acquiring dynamic information despite decades of effort. Dynamic `structural changes' are often indirectly inferred from `structural differences' by comparing related static structures. In contrast, the direct observation of dynamic structural changes requires the initiation of a biochemical reaction or process in a crystal. Both the direct and the indirect approaches share a common challenge in analysis: how to interpret the structural heterogeneity intrinsic to all dynamic processes. This paper presents a real-space approach to this challenge, in which a suite of analytical methods and tools to identify and refine the mixed structural species present in multiple crystallographic data sets have been developed. These methods have been applied to representative scenarios in dynamic crystallography, and reveal structural information that is otherwise difficult to interpret or inaccessible using conventional methods.

Molecular Dynamics Computer Simulations of Multidrug RND Efflux Pumps

Over-expression of multidrug efflux pumps of the Resistance Nodulation Division (RND) protein super family counts among the main causes for microbial resistance against pharmaceuticals. Understanding the molecular basis of this process is one of the major challenges of modern biomedical research, involving a broad range of experimental and computational techniques. Here we review the current state of RND transporter investigation employing molecular dynamics simulations providing conformational samples of transporter components to obtain insights into the functional mechanism underlying efflux pump-mediated antibiotics resistance in Escherichia coli and Pseudomonas aeruginosa.

Porter domain opening and closing motions in the multi-drug efflux transporter AcrB

Acriflavine resistance protein B acts as the active transporter in the multi-drug efflux pump Acriflavine resistance proteins A / B - Tolerance to colicins protein in Escherichia coli. Within the same reaction cycle intermediate all Acriflavine resistance protein B X-ray structures display highly similar conformations of the substrate-recruiting and transporting porter domain. To assess if this structural homogeneity is an intrinsic feature of Acriflavine resistance protein B or stems from other causes we performed a series of six independent, unbiased 100 ns molecular dynamics simulations of membrane-embedded, asymmetric, substrate-free wild type Acriflavine resistance protein B in a 150 mM NaCl solution. We find the porter domain more flexible than previously assumed displaying clear opening and closing motions of the proximal binding pocket (L and T-state) and the exit of the drug transport channels (O-intermediate). Concurrently the hydrophobic binding pocket favors a closed conformation in all three protomers. Our findings suggest that the conformational homogeneity seen in the crystal structures is likely an effect of bound but structurally unresolved substrate. Our simulations further imply that each of the known three reaction cycle intermediates occurs in at least two variants, the Thr676 loop independently regulates porter domain access and likely plays a key role in substrate transport. On a 100 ns time scale we find no evidence supporting the proposed LLL resting state in the absence of substrate. If the proximal binding pocket dynamics have an inhibiting effect on Acriflavine resistance protein B pump activity lowering the life time of substrate-accessible conformations, the observed dynamics could provide a structural explanation for the Acriflavine resistance protein B activity-enhancing effect of the adaptor protein Acriflavine resistance protein A stabilizing PC1 and PC2 subdomain orientations.

Three ways in, one way out: water dynamics in the trans-membrane domains of the inner membrane translocase AcrB

Powered by proton-motive force, the inner membrane translocase AcrB is the engine of the AcrAB-TolC efflux pump in Escherichia coli. As proton conduction in proteins occurs along hydrogen-bonded networks of polar residues and water molecules, knowledge of the protein-internal water distribution and water-interacting residues allows drawing conclusions to possible pathways of proton conduction. Here, we report a series of 6× 50 ns independent molecular dynamics simulations of asymmetric AcrB embedded in a phospholipid/water environment. Simulating each monomer in its proposed protonation state, we calculated for each trans-membrane domain the average water distribution, identified residues interacting with these waters and quantified each residue's frequency of water hydrogen bond contact. Combining this information we find three possible routes of proton transfer connecting a continuously hydrated region of known key residues in the TMD interior to bulk water by one cytoplasmic and up to three periplasm water channels in monomer B and A. We find that water access of the trans-membrane domains is regulated by four groups of residues in a combination of side chain re-orientations and shifts of trans-membrane helices. Our findings support a proton release event via Arg971 during the C intermediate or in the transition to A, and proton uptake occurring in the A or B state or during a so far unknown intermediate in between B and C where cytoplasmic water access is still possible. Our simulations suggest experimentally testable hypotheses, which have not been investigated so far.

Membrane protein dynamics from femtoseconds to seconds

Membrane proteins play a key role in energy conversion, transport, signal recognition, transduction, and other fundamental biological processes. Despite considerable progress in experimental techniques, the determination of structure and dynamics of membrane proteins still represents a great challenge. Computer simulation methods are becoming an increasingly important tool not only in the interpretation of experiments but also in the prediction of membrane protein dynamics. In the present review, we give a brief introduction to molecular modeling techniques currently used to explore protein dynamics on time scales ranging from femtoseconds to microseconds. We then describe a few recent example applications of these techniques to membrane proteins. In conclusion, we also discuss some of the newest developments in simulation methodology that have the potential to further extend the time scale accessible to explore (membrane) protein dynamics.

Determination of ensemble-average pairwise root mean-square deviation from experimental B-factors

Root mean-square deviation (RMSD) after roto-translational least-squares fitting is a measure of global structural similarity of macromolecules used commonly. On the other hand, experimental x-ray B-factors are used frequently to study local structural heterogeneity and dynamics in macromolecules by providing direct information about root mean-square fluctuations (RMSF) that can also be calculated from molecular dynamics simulations. We provide a mathematical derivation showing that, given a set of conservative assumptions, a root mean-square ensemble-average of an all-against-all distribution of pairwise RMSD for a single molecular species, <RMSD(2)>(1/2), is directly related to average B-factors (<B>) and <RMSF(2)>(1/2). We show this relationship and explore its limits of validity on a heterogeneous ensemble of structures taken from molecular dynamics simulations of villin headpiece generated using distributed-computing techniques and the Folding@Home cluster. Our results provide a basis for quantifying global structural diversity of macromolecules in crystals directly from x-ray experiments, and we show this on a large set of structures taken from the Protein Data Bank. In particular, we show that the ensemble-average pairwise backbone RMSD for a microscopic ensemble underlying a typical protein x-ray structure is approximately 1.1 A, under the assumption that the principal contribution to experimental B-factors is conformational variability.

Bend-twist-stretch model for coarse elastic network simulation of biomolecular motion

The empirical harmonic potential function of elastic network models (ENMs) is augmented by three- and four-body interactions as well as by a parameter-free connection rule. In the new bend-twist-stretch (BTS) model the complexity of the parametrization is shifted from the spatial level of detail to the potential function, enabling an arbitrary coarse graining of the network. Compared to distance cutoff-based Hookean springs, the approach yields a more stable parametrization of coarse-grained ENMs for biomolecular dynamics. Traditional ENMs give rise to unbounded zero-frequency vibrations when (pseudo)atoms are connected to fewer than three neighbors. A large cutoff is therefore chosen in an ENM (about twice the average nearest-neighbor distance), resulting in many false-positive connections that reduce the spatial detail that can be resolved. More importantly, the required three-neighbor connectedness also limits the coarse graining, i.e., the network must be dense, even in the case of low-resolution structures that exhibit few spatial features. The new BTS model achieves such coarse graining by extending the ENM potential to include three-and four-atom interactions (bending and twisting, respectively) in addition to the traditional two-atom stretching. Thus, the BTS model enables reliable modeling of any three-dimensional graph irrespective of the atom connectedness. The additional potential terms were parametrized using continuum elastic theory of elastic rods, and the distance cutoff was replaced by a competitive Hebb connection rule, setting all free parameters in the model. We validate the approach on a carbon-alpha representation of adenylate kinase and illustrate its use with electron microscopy maps of E. coli RNA polymerase, E. coli ribosome, and eukaryotic chaperonin containing T-complex polypeptide 1, which were difficult to model with traditional ENMs. For adenylate kinase, we find excellent reproduction (>90% overlap) of the ENM modes and B factors when BTS is applied to the carbon-alpha representation as well as to coarser descriptions. For the volumetric maps, coarse BTS yields similar motions (70%-90% overlap) to those obtained from significantly denser representations with ENM. Our Python-based algorithms of ENM and BTS implementations are freely available.

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.

Nucleotide-dependent allostery within the ABC transporter ATP-binding cassette: a computational study of the MJ0796 dimer

ATP-binding cassette transporters perform energy-dependent transmembrane solute trafficking in all organisms. These proteins often mediate cellular resistance to therapeutic drugs and are involved in a range of human genetic diseases. Enzymological studies have implicated a helical subdomain within the ATP-binding cassette nucleotide-binding domain in coupling ATP hydrolysis to solute transport in the transmembrane domains. Consistent with this, structural and computational analyses have indicated that the helical subdomain undergoes nucleotide-dependent movement relative to the core of the nucleotide-binding domain fold. Here we use theoretical methods to examine the allosteric nucleotide dependence of helical subdomain transitions to further elucidate its role in interactions between the transmembrane and nucleotide-binding domains. Unrestrained 30-ns molecular dynamics simulations of the ATP-bound, ADP-bound, and apo states of the MJ0796 monomer support the idea that interaction of a conserved glutamine residue with the catalytic metal mediates the rotation of the helical subdomain in response to nucleotide binding and hydrolysis. Simulations of the nucleotide-binding domain dimer revealed that ATP hydrolysis induces a large transition of one helical subdomain, resulting in an asymmetric conformation of the dimer not observed previously. A coarse-grained elastic network analysis supports this finding, revealing the existence of corresponding dynamic modes intrinsic to the contact topology of the protein. The implications of these findings for the coupling of ATP hydrolysis to conformational changes in the transmembrane domains required for solute transport are discussed in light of recent whole transporter structures.

Conformation of the AcrB multidrug efflux pump in mutants of the putative proton relay pathway

We previously reported the X-ray structures of wild-type Escherichia coli AcrB, a proton motive force-dependent multidrug efflux pump, and its N109A mutant. These structures presumably reflect the resting state of AcrB, which can bind drugs. After ligand binding, a proton may bind to an acidic residue(s) in the transmembrane domain, i.e., Asp407 or Asp408, within the putative network of electrostatically interacting residues, which also include Lys940 and Thr978, and this may initiate a series of conformational changes that result in drug expulsion. Herein we report the X-ray structures of four AcrB mutants, the D407A, D408A, K940A, and T978A mutants, in which the structure of this tight electrostatic network is expected to become disrupted. These mutant proteins revealed remarkably similar conformations, which show striking differences from the previously known conformations of the wild-type protein. For example, the loop containing Phe386 and Phe388, which play a major role in the initial binding of substrates in the central cavity, becomes prominently extended into the center of the cavity, such that binding of large substrate molecules may become difficult. We believe that this new conformation may mimic, at least partially, one of the transient conformations of the transporter during the transport cycle.