Detection of Side Chain Rearrangements Mediating the Motions of Transmembrane Helices in Molecular Dynamics Simulations of G Protein-Coupled Receptors

Conformational dynamics of complement protease C1r inhibitor proteins from Lyme disease- and relapsing fever-causing spirochetes

Borrelial pathogens are vector-borne etiological agents known to cause Lyme disease, relapsing fever, and Borrelia miyamotoi disease. These spirochetes each encode several surface-localized lipoproteins that bind components of the human complement system to evade host immunity. One borrelial lipoprotein, BBK32, protects the Lyme disease spirochete from complement-mediated attack via an alpha helical C-terminal domain that interacts directly with the initiating protease of the classical complement pathway, C1r. In addition, the B. miyamotoi BBK32 orthologs FbpA and FbpB also inhibit C1r, albeit via distinct recognition mechanisms. The C1r-inhibitory activities of a third ortholog termed FbpC, which is found exclusively in relapsing fever-causing spirochetes, remains unknown. Here, we report the crystal structure of the C-terminal domain of Borrelia hermsii FbpC to a limiting resolution of 1.5 Å. We used surface plasmon resonance and assays of complement function to demonstrate that FbpC retains potent BBK32-like anticomplement activities. Based on the structure of FbpC, we hypothesized that conformational dynamics of the complement inhibitory domains of borrelial C1r inhibitors may differ. To test this, we utilized the crystal structures of the C-terminal domains of BBK32, FbpA, FbpB, and FbpC to carry out molecular dynamics simulations, which revealed borrelial C1r inhibitors adopt energetically favored open and closed states defined by two functionally critical regions. Taken together, these results advance our understanding of how protein dynamics contribute to the function of bacterial immune evasion proteins and reveal a surprising plasticity in the structures of borrelial C1r inhibitors.

"Conformational dynamics of C1r inhibitor proteins from Lyme disease and relapsing fever spirochetes"

Borrelial pathogens are vector-borne etiological agents of Lyme disease, relapsing fever, and Borrelia miyamotoi disease. These spirochetes each encode several surface-localized lipoproteins that bind to components of the human complement system. BBK32 is an example of a borrelial lipoprotein that protects the Lyme disease spirochete from complement-mediated attack. The complement inhibitory activity of BBK32 arises from an alpha helical C-terminal domain that interacts directly with the initiating protease of the classical pathway, C1r. Borrelia miyamotoi spirochetes encode BBK32 orthologs termed FbpA and FbpB, and these proteins also inhibit C1r, albeit via distinct recognition mechanisms. The C1r-inhibitory activities of a third ortholog termed FbpC, which is found exclusively in relapsing fever spirochetes, remains unknown. Here we report the crystal structure of the C-terminal domain of B. hermsii FbpC to a limiting resolution of 1.5 Å. Surface plasmon resonance studies and assays of complement function demonstrate that FbpC retains potent BBK32-like anti-complement activities. Based on the structure of FbpC, we hypothesized that conformational dynamics of the complement inhibitory domains of borrelial C1r inhibitors may differ. To test this, we utilized the crystal structures of the C-terminal domains of BBK32, FbpA, FbpB, and FbpC to carry out 1 µs molecular dynamics simulations, which revealed borrelial C1r inhibitors adopt energetically favored open and closed states defined by two functionally critical regions. This study advances our understanding of how protein dynamics contribute to the function of bacterial immune evasion proteins and reveals a surprising plasticity in the structures of borrelial C1r inhibitors.

DHFR Mutants Modulate Their Synchronized Dynamics with the Substrate by Shifting Hydrogen Bond Occupancies

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize C_α-C_β bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-μs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

In silico analysis and in planta production of recombinant ccl21/IL1β protein and characterization of its in vitro anti-tumor and immunogenic activity

CCL21 has an essential role in anti-tumor immune activity. Epitopes of IL1β have adjuvant activity without causing inflammatory responses. CCR7 and its ligands play a vital role in the immune balance; specifically, in transport of T lymphocytes and antigen-presenting cells such as dendritic cells to the lymph nodes. This study aimed to produce epitopes of CCL21 and IL1β as a recombinant protein and characterize its in vitro anti-tumor and immunogenic activity. A codon-optimized ccl21/IL1β gene was designed and synthesized from human genes. Stability and binding affinity of CCL21/IL1β protein and CCR7 receptor were examined through in silico analyses. The construct was introduced into N. tabacum to produce this recombinant protein and the structure and function of CCL21/IL1β were examined. Purified protein from transgenic leaves generated a strong signal in SDS PAGE and western blotting assays. FTIR measurement and MALDI-TOF/TOF mass spectrography showed that ccl21/IL-1β was correctly expressed in tobacco plants. Potential activity of purified CCL21/IL1β in stimulating the proliferation and migration of MCF7 cancer cell line was investigated using the wound healing method. The results demonstrated a decrease in survival rate and metastasization of cancer cells in the presence of CCL21/IL1β, and IC₅₀ of CCL21 on MCF7 cells was less than that of non-recombinant protein. Agarose assay on PBMCsCCR7+ showed that CCL21/IL1β has biological activity and there is a distinguishable difference between chemokinetic (CCL21) and chemotactic (FBS) movements. Overall, the results suggest that CCL21/IL1β could be considered an effective adjuvant in future in vivo and clinical tests.

Mutual information analysis of the dynamic correlation between side chains in proteins

Protein dynamics play an essential role in function regulation. In recent years, many experimental and theoretical studies have shown that changes in protein fluctuations in the backbone and side chains fulfill a pivotal role associated with amino acid mutations, chemical modifications, and ligand binding. The dynamic correlations between protein side chains have not been sufficiently studied, and no reliable analysis method has been available so far. Therefore, we developed a method to evaluate the dynamic correlation between protein side chains using mutual information and molecular dynamics simulations. To eliminate the structural superposition errors dealing with conventional analysis methods, and to accurately extract the intrinsic fluctuation properties of the side chains, we employed distance principal component analysis (distPCA). The motion of the side chain was then projected onto the eigenvector space obtained by distPCA, and the mutual information between the projected motions was calculated. The proposed method was then applied to a small protein "eglin c" and the mutants. The results show that even a single mutation significantly changed the dynamic correlations and also suggest that the dynamic change is deeply related to the stability. Those results indicate that our developed method could be useful for analyzing the molecular mechanism of allosteric communication in proteins.

Conserved Conformational Hierarchy across Functionally Divergent Glycosyltransferases of the GT-B Structural Superfamily as Determined from Microsecond Molecular Dynamics

It has long been understood that some proteins undergo conformational transitions en route to the Michaelis Complex to allow chemistry. Examination of crystal structures of glycosyltransferase enzymes in the GT-B structural class reveals that the presence of ligand in the active site triggers an open-to-closed conformation transition, necessary for their catalytic functions. Herein, we describe microsecond molecular dynamics simulations of two distantly related glycosyltransferases that are part of the GT-B structural superfamily, HepI and GtfA. Simulations were performed using the open and closed conformations of these unbound proteins, respectively, and we sought to identify the major dynamical modes and communication networks that interconnect the open and closed structures. We provide the first reported evidence within the scope of our simulation parameters that the interconversion between open and closed conformations is a hierarchical multistep process which can be a conserved feature of enzymes of the same structural superfamily. Each of these motions involves of a collection of smaller molecular reorientations distributed across both domains, highlighting the complexities of protein dynamic involved in the interconversion process. Additionally, dynamic cross-correlation analysis was employed to explore the potential effect of distal residues on the catalytic efficiency of HepI. Multiple distal nonionizable residues of the C-terminal domain exhibit motions anticorrelated to positively charged residues in the active site in the N-terminal domain involved in substrate binding. Mutations of these residues resulted in a reduction in negatively correlated motions and an altered enzymatic efficiency that is dominated by lower K_m values with k_cat effectively unchanged. The findings suggest that residues with opposing conformational motions involved in the opening and closing of the bidomain HepI protein can allosterically alter the population and conformation of the “closed” state, essential to the formation of the Michaelis complex. The stabilization effects of these mutations likely equally influence the energetics of both the ground state and the transition state of the catalytic reaction, leading to the unaltered k_cat. Our study provides new insights into the role of conformational dynamics in glycosyltransferase’s function and new modality to modulate enzymatic efficiency.

Structural and Functional Analysis of Gly212 Mutants Reveals the Importance of Intersubunit Interactions in ASIC1a Channel Function

Acid-sensing ion channels (ASICs) act as pH sensors in neurons. ASICs contribute to pain sensation, learning, fear behavior and to neuronal death after ischemic stroke. Extracellular acidification induces a transient activation and subsequent desensitization of these Na⁺-selective channels. ASICs are trimeric channels made of identical or homologous subunits. We have previously shown that mutation of the highly conserved Gly212 residue of human ASIC1a to Asp affects the channel function. Gly212 is located in the proximity of a predicted Cl^– binding site at a subunit interface. Here, we have measured the function of a series of Gly212 mutants. We show that substitution of Gly212 affects the ASIC1a pH dependence and current decay kinetics. Intriguingly, the mutations to the acidic residues Asp and Glu have opposing effects on the pH dependence and the current decay kinetics. Analysis of molecular dynamics simulation trajectories started with the coordinates of the closed conformation indicates that the immediate environment of residue 212 in G212E, which shifts the pH dependence to more alkaline values, adopts a conformation closer to the open state. The G212D and G212E mutants have a different pattern of intersubunit salt bridges, that, in the case of G212E, leads to an approaching of neighboring subunits. Based on the comparison of crystal structures, the conformational changes in this zone appear to be smaller during the open-desensitized transition. Nevertheless, MD simulations highlight differences between mutants, suggesting that the changed function upon substitution of residue 212 is due to differences in intra- and intersubunit interactions in its proximity.

Cross-linking, DEER-spectroscopy and molecular dynamics confirm the inward facing state of P-glycoprotein in a lipid membrane

The drug efflux pump P-glycoprotein (P-gp) displays a complex transport mechanism involving multiple drug binding sites and two centres for nucleotide hydrolysis. Elucidating the molecular mechanism of transport remains elusive and the availability of P-gp structures in distinct natural and ligand trapped conformations will accelerate our understanding. The present investigation sought to provide biochemical data to validate specific features of these structures; with particular focus on the transmembrane domain that provides the transport conduit. Hence our focus was on transmembrane helices six and twelve (TM6/TM12), which are believed to participate in drug binding, as they line the central transport conduit and provide a direct link to the catalytic centres. A series of P-gp mutants were generated with a single cysteine in both TM6 and TM12 to facilitate measurement of inter-helical distances using cross-linking and DEER strategies. Experimental results were compared to published structures per se and those refined by MD simulations. This analysis revealed that the refined inward-facing murine structure (4M1M) of P-gp provides a good representation of the proximity, topography and relative motions of TM6 and TM12 in reconstituted human P-gp.

Mechanistic insights into the deleterious roles of Nasu-Hakola disease associated TREM2 variants

Recently, the critical roles played by genetic variants of TREM2 (Triggering Receptor Expressed on Myeloid cells 2) in Alzheimer's disease have been aggressively highlighted. However, few studies have focused on the deleterious roles of Nasu-Hakola disease (NHD) associated TREM2 variants. In order to get insights into the contributions made by these variants to neurodegeneration, we investigated the influences of four NHD associated TREM2 mutations (Y38C, W50C, T66M, and V126G) on loss-of-function, and followed this with in silico prediction and conventional molecular dynamics simulation. NHD mutations were predicted to be highly deleterious by eight different in silico bioinformatics tools and found to induce conformational changes by molecular dynamics simulation. As compared with the wild-type, the four variants produced substantial differences in the collective motions of loop regions, which not only promoted structural remodeling in the CDR2 (complementarity-determining region 2) loop but also in the CDR1 loop, by changing inter- and intra-loop hydrogen bonding networks. In addition, structural studies in a free energy landscape analysis showed that Y38, T66, and V126 are crucial for maintaining the structural features of CDR1 and CDR2 loops, and that mutations in these positions produced steric clashes and loss of ligand binding. These results showed the presence of mutations in the TREM2 ectodomain induced flexibility and caused structural alterations. Dynamical scenarios, as provided by the present study, may be critical to our understanding of the roles of these TREM2 mutations in neurodegenerative diseases.