The N-terminal domain of RfaH plays an active role in protein fold-switching

The bacterial elongation factor RfaH promotes the expression of virulence factors by specifically binding to RNA polymerases (RNAP) paused at a DNA signal. This behavior is unlike that of its paralog NusG, the major representative of the protein family to which RfaH belongs. Both proteins have an N-terminal domain (NTD) bearing an RNAP binding site, yet NusG C-terminal domain (CTD) is folded as a β-barrel while RfaH CTD is forming an α-hairpin blocking such site. Upon recognition of the specific DNA exposed by RNAP, RfaH is activated via interdomain dissociation and complete CTD structural rearrangement into a β-barrel structurally identical to NusG CTD. Although RfaH transformation has been extensively characterized computationally, little attention has been given to the role of the NTD in the fold-switching process, as its structure remains unchanged. Here, we used Associative Water-mediated Structure and Energy Model (AWSEM) molecular dynamics to characterize the transformation of RfaH, spotlighting the sequence-dependent effects of NTD on CTD fold stabilization. Umbrella sampling simulations guided by native contacts recapitulate the thermodynamic equilibrium experimentally observed for RfaH and its isolated CTD. Temperature refolding simulations of full-length RfaH show a high success towards α-folded CTD, whereas the NTD interferes with βCTD folding, becoming trapped in a β-barrel intermediate. Meanwhile, NusG CTD refolding is unaffected by the presence of RfaH NTD, showing that these NTD-CTD interactions are encoded in RfaH sequence. Altogether, these results suggest that the NTD of RfaH favors the α-folded RfaH by specifically orienting the αCTD upon interdomain binding and by favoring β-barrel rupture into an intermediate from which fold-switching proceeds.     

A specific interdomain interaction preserves the structural and binding properties of the ModA protein from the phytopathogen Xanthomonas citri domain interaction and transport in ModA

The periplasmic-binding proteins in ATP-binding cassette systems (ABC Transporters) are responsible for the capture and delivery of ligands to their specific transporters, triggering a series of ATP-driven conformational changes that leads to the transport of the ligand. Structurally consisting of two lobes, the proteins change conformation after interaction with the ligand. The structure of the molybdate-binding protein (ModA) from Xanthomonas citri, bound to molybdate, was previously solved by our group and an interdomain interaction, mediated by a salt bridge between K127 and D59, apparently supports the binding properties and keeps the domains closed. To determinate the importance of this interaction, we built two ModA mutants, K127S and D59A, and analysed their functional and structural properties. Based on a set of spectroscopic experiments, crystallisation trials, structure determination and molecular dynamics (MD) simulations, we showed that the salt bridge is essential to maintain the structure and binding properties. Additionally, the MD simulations revealed that this mutant adopted a more compact structure that packed down the ligand-binding pocket. From the closed bound to open structure, the positioning of the helices forming the dipole and the salt bridge are essential to induce an intermediate state.     

β‐strand interactions at the domain interface critical for the stability of human lens γD‐crystallin

Human age‐onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens γD‐crystallin protein and its isolated domains. A partially unfolded folding intermediate of γD‐crystallin is detected in simulations with its C‐terminal domain (C‐td) folded and N‐terminal domain (N‐td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N‐td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C‐td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non‐native surface salt‐bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non‐native Glu135‐Arg142 salt‐bridge causes significant destabilization to the folding core of the isolated C‐td, which, in turn, induces unfolding of the N‐td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of γD‐crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the γD‐crystallin surface might lead to unfolding and subsequent aggregation.     

Coupled intra- and interdomain dynamics support domain cross-talk in Pin1

The functional mechanisms of multidomain proteins often exploit interdomain interactions, or “cross-talk.” An example is human Pin1, an essential mitotic regulator consisting of a Trp–Trp (WW) domain flexibly tethered to a peptidyl-prolyl isomerase (PPIase) domain, resulting in interdomain interactions important for Pin1 function. Substrate binding to the WW domain alters its transient contacts with the PPIase domain via means that are only partially understood. Accordingly, we have investigated Pin1 interdomain interactions using NMR paramagnetic relaxation enhancement (PRE) and molecular dynamics (MD) simulations. The PREs show that apo-Pin1 samples interdomain contacts beyond the range suggested by previous structural studies. They further show that substrate binding to the WW domain simultaneously alters interdomain separation and the internal conformation of the WW domain. A 4.5-μs all-atom MD simulation of apo-Pin1 suggests that the fluctuations of interdomain distances are correlated with fluctuations of WW domain interresidue contacts involved in substrate binding. Thus, the interdomain/WW domain conformations sampled by apo-Pin1 may already include a range of conformations appropriate for binding Pin1''s numerous substrates. The proposed coupling between intra-/interdomain conformational fluctuations is a consequence of the dynamic modular architecture of Pin1. Such modular architecture is common among cell-cycle proteins; thus, the WW–PPIase domain cross-talk mechanisms of Pin1 may be relevant for their mechanisms as well.     

Akt phosphorylation regulates the tumour-suppressor merlin through ubiquitination and degradation

The neurofibromatosis-2 (NF2) tumour-suppressor gene encodes an intracellular membrane-associated protein, called merlin, whose growth-suppressive function is dependent on its ability to form interactions through its intramolecular amino-terminal domain (NTD) and carboxy-terminal domain (CTD). Merlin phosphorylation plays a critical part in dictating merlin NTD/CTD interactions as well as in controlling binding to its effector proteins. Merlin is partially regulated by phosphorylation of Ser 518, such that hyperphosphorylated merlin is inactive and fails to form productive intramolecular and intermolecular interactions. Here, we show that the protein kinase Akt directly binds to and phosphorylates merlin on residues Thr 230 and Ser 315, which abolishes merlin NTD/CTD interactions and binding to merlin's effector protein PIKE-L and other binding partners. Furthermore, Akt-mediated phosphorylation leads to merlin degradation by ubiquitination. These studies demonstrate that Akt-mediated merlin phosphorylation regulates the function of merlin in the absence of an inactivating mutation.     

Targeting disordered-structured domain interactions in Galectin-3 based on NMR and enhanced MD

Intrinsically disordered regions (IDRs) are common and important functional domains in many proteins. However, IDRs are difficult to target for drug development due to the lack of defined structures that would facilitate the identification of possible drug-binding pockets. Galectin-3 is a carbohydrate-binding protein of which overexpression has been implicated in a wide variety of disorders, including cancer and inflammation. Apart from its carbohydrate-recognition/binding domain (CRD), Galectin-3 also contains a functionally important disordered N-terminal domain (NTD) that contacts the C-terminal domain (CTD) and could be a target for drug development. To overcome challenges involved in inhibitor design due to lack of structure and the highly dynamic nature of the NTD, we used a protocol combining nuclear magnetic resonance data from recombinant Galectin-3 with accelerated molecular dynamics (MD) simulations. This approach identified a pocket in the CTD with which the NTD makes frequent contact. In accordance with this model, mutation of residues L131 and L203 in this pocket caused loss of Galectin-3 agglutination ability, signifying the functional relevance of the cavity. In silico screening was used to design candidate inhibitory peptides targeting the newly discovered cavity, and experimental testing of only three of these yielded one peptide that inhibits the agglutination promoted by wild-type Galectin-3. NMR experiments further confirmed that this peptide indeed binds to a cavity in the CTD, not within the actual CRD. Our results show that it is possible to apply a combination of MD simulations and NMR experiments to precisely predict the binding interface of a disordered domain with a structured domain, and furthermore use this predicted interface for designing inhibitors. This procedure can potentially be extended to many other targets in which similar IDR interactions play a vital functional role.     

The conserved Arg241‐Glu439 salt bridge determines flexibility of the inositol 1,4,5‐trisphosphate receptor binding core in the ligand‐free state

Inositol 1,4,5‐trisphosphate receptor (InsP₃R) is an intracellular Ca²⁺‐release channel activated by binding of inositol 1,4,5‐trisphosphate (InsP₃) to the InsP₃ binding core (IBC). Structural change in the IBC upon InsP₃ binding is the key process in channel pore opening. In this study, we performed molecular dynamics (MD) simulations of the InsP₃‐free form of the IBC, starting with removal of InsP₃ from the InsP₃‐bound crystal structure, and obtained the structural ensemble of the InsP₃‐free form of the IBC. The simulation revealed that the two domains of the IBC largely fluctuate around the average structure with the hinge angle opened 17° more than in the InsP₃‐bound form, and the twist angle rotated by 45°, forming interdomain contacts that are different from those in the bound form. The InsP₃ binding loop was disordered. The InsP₃‐free form thus obtained was reproduced four times in simulations started from a fully extended configuration of the two domains. Simulations beginning with the fully extended form indicated that formation of a salt bridge between Arg241 and Glu439 is crucial for stabilizing the closed form of the two domains. Mutation of Arg241 to Gln prevented formation of the compact structure by the two domains, but the fully flexible domain arrangement was maintained. Thus, the Arg241‐Glu439 salt bridge determines the flexibility of the InsP₃‐free form of the IBC.Proteins 2013; 81:1699–1708. © 2013 Wiley Periodicals, Inc.     

The weak interdomain coupling observed in the 70 kDa subunit of human replication protein A is unaffected by ssDNA binding

Replication protein A (RPA) is a heterotrimeric, multi-functional protein that binds single-stranded DNA (ssDNA) and is essential for eukaryotic DNA metabolism. Using heteronuclear NMR methods we have investigated the domain interactions and ssDNA binding of a fragment from the 70 kDa subunit of human RPA (hRPA70). This fragment contains an N-terminal domain (NTD), which is important for hRPA70–protein interactions, connected to a ssDNA-binding domain (SSB1) by a flexible linker (hRPA70_1–326). Correlation analysis of the amide ¹H and ¹⁵N chemical shifts was used to compare the structure of the NTD and SSB1 in hRPA70_1–326 with two smaller fragments that corresponded to the individual domains. High correlation coefficients verified that the NTD and SSB1 maintained their structures in hRPA70_1–326, indicating weak interdomain coupling. Weak interdomain coupling was also suggested by a comparison of the transverse relaxation rates for hRPA70_1–326 and one of the smaller hRPA70 fragments containing the NTD and the flexible linker (hRPA70_1–168). We also examined the structure of hRPA70_1–326 after addition of three different ssDNA substrates. Each of these substrates induced specific amide ¹H and/or ¹⁵N chemical shift changes in both the NTD and SSB1. The NTD and SSB1 have similar topologies, leading to the possibility that ssDNA binding induced the chemical shift changes observed for the NTD. To test this hypothesis we monitored the amide ¹H and ¹⁵N chemical shift changes of hRPA70_1–168 after addition of ssDNA. The same amide ¹H and ¹⁵N chemical shift changes were observed for the NTD in hRPA70_1–168 and hRPA70_1–326. The NTD residues with the largest amide ¹H and/or ¹⁵N chemical shift changes were localized to a basic cleft that is important for hRPA70–protein interactions. Based on this relationship, and other available data, we propose a model where binding between the NTD and ssDNA interferes with hRPA70–protein interactions.     

Immunoglobulin kappa light chain and its amyloidogenic mutants: a molecular dynamics study

AL amyloidosis and LCDD are pathological conditions caused by extracellural deposition of monoclonal Ig light chain variable domains. In the former case, deposits have a form of amyloid fibrils, in the latter, amorphous aggregates. 1REI kappa light chain variable domain and its two point mutants, R61N and D82I, were chosen for the analysis in this work. Wild 1REI does not create deposits in vitro, while R61N aggregates as amyloid fibrils and D82I creates amorphous aggregates. Both mutated residues create a conserved salt bridge; thus, substitution of any of them should decrease V(L) domain stability. For these three proteins, 5 ns MD simulations were conducted in temperatures of 300 K and 400 K, with protonated and unprotonated acidic residues, mimicking acidic and neutral experimental pH conditions (3 sets: N300, N400, and A400). The analysis of trajectories focused on characterization of changes in conformational behavior and stability of Ig kappa light chain variable domain caused by single aminoacid substitutions that were experimentally proved to enhance aggregation propensity, both in the form of amyloid and amorphous aggregates. Residue D82 turns out to be involved not only in R61-D82 but also in K45-D82 interaction, which was not observed in the X-ray structure, but frequently populated simulations of 1REI. The substitution D82I excludes both interactions, resulting in substantial destabilization (i.e., easier aggregation). Examination of behavior of edge regions of V(L) beta-sandwich reveals significant alterations in D82I mutant compared to wild 1REI, while relatively small changes occur in R61N. This suggests that mild and slow destabilization is the reason of the conversion of V(L) to partially folded amyloidogenic intermediate structure.     

Coarse‐grained simulations of proton‐dependent conformational changes in lactose permease

During lactose/H⁺ symport, the Escherichia coli lactose permease (LacY) undergoes a series of global conformational transitions between inward‐facing (open to cytoplasmic side) and outward‐facing (open to periplasmic side) states. However, the exact local interactions and molecular mechanisms dictating those large‐scale structural changes are not well understood. All‐atom molecular dynamics simulations have been performed to investigate the molecular interactions involved in conformational transitions of LacY, but the simulations can only explore early or partial global structural changes because of the computational limits (< 100 ns). In this work, we implement a hybrid force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid, to investigate the proton‐dependent dynamics and conformational changes of LacY. The effects of the protonation states on two key glutamate residues (Glu325 and Glu269) have been studied. Our results on the salt‐bridge dynamics agreed with all‐atom simulations at early short time period, validating our simulations. From our microsecond simulations, we were able to observe the complete transition from inward‐facing to outward‐facing conformations of LacY. Our results showed that all helices have participated during the global conformational transitions and helical movements of LacY. The inter‐helical distances measured in our simulations were consistent with the double electron‐electron resonance experiments at both cytoplasmic and periplasmic sides. Our simulations indicated that the deprotonation of Glu325 induced the opening of the periplasmics side and partial closure of the cytoplasmic side of LacY, while protonation of the Glu269 caused a stable cross‐domain salt‐bridge (Glu130‐Arg344) and completely closed the cytoplasmic side. Proteins 2016; 84:1067–1074. © 2016 Wiley Periodicals, Inc.     

Investigating the inter‐subunit/subdomain interactions and motions relevant to disease mutations in the N‐terminal domain of ryanodine receptors by molecular dynamics simulation

The ryanodine receptors (RyR) are essential to calcium signaling in striated muscles, and numerous disease mutations have been identified in two RyR isoforms, RyR1 in skeletal muscle and RyR2 in cardiac muscle. A deep understanding of the activation/regulation mechanisms of RyRs has been hampered by the shortage of high‐resolution structures and dynamic information for this giant tetrameric complex in different functional states. Toward elucidating the molecular mechanisms of disease mutations in RyRs, we performed molecular dynamics simulation of the N‐terminal domain (NTD) which is not only the best‐resolved structural component of RyRs, but also a hotspot of disease mutations. First, we simulated the tetrameric NTD of wild‐type RyR1 and three disease mutants (K155E, R157Q, and R164Q) that perturb the inter‐subunit interfaces. Our simulations identified a dynamic network of salt bridges involving charged residues at the inter‐subunit/subdomain interfaces and disease‐mutation sites. By perturbing this key network, the above three mutations result in greater flexibility with the highest inter‐subunit opening probability for R157Q. Next, we simulated the monomeric NTD of RyR2 in the presence or absence of a central Cl^– anion which is known to stabilize the interfaces between the three NTD subdomains (A, B, and C). We found that the loss of Cl^– restructures the salt‐bridge network near the Cl^–‐binding site, leading to rotations of subdomain A/B relative to subdomain C and enhanced mobility between the subdomains. This finding supports a mechanism for disease mutations in the NTD of RyR2 via perturbation of the Cl^– binding. The rich structural and dynamic information gained from this study will guide future mutational and functional studies of the NTD of RyRs. Proteins 2017; 85:1633–1644. © 2017 Wiley Periodicals, Inc.     

Decoding the structural events in substrate‐gating mechanism of eukaryotic prolyl oligopeptidase using normal mode analysis and molecular dynamics simulations

Prolyl oligopeptidase (POP) is a serine protease, unique for its ability to cleave various small oligopeptides shorter than 30 amino acids. POP is an important drug target since it is implicated in various neurological disorders. Although there is structural evidence that bacterial POPs undergo huge interdomain movements and acquire an “open” state in the substrate‐unbound form, hitherto, no crystal structure is available in the substrate‐unbound domain‐open form of eukaryotic POPs. Indeed, there is no difference between the substrate‐unbound/bound states of eukaryotic POPs. This raises unanswered questions about whether difference in the substrate access pathway exists between bacterial and eukaryotic POPs. Here, we have used normal mode analysis and molecular dynamics to unravel the mechanism of substrate entry in mammalian POPs, which has been debated until now. Motions observed using normal modes of porcine and bacterial POPs were analyzed and compared, augmented by molecular dynamics of these proteins. Identical to bacterial POPs, interdomain opening was found to be the possible pathway for the substrate‐gating in mammals as well. On the basis of our analyses and evidences, a mechanistic model of substrate entry in POPs has been proposed. Up‐down movement of N‐terminal hydrolase domain resulted in twisting motion of two domains, followed by the conformational changes of interdomain loop regions, which facilitate interdomain opening. Similar to bacterial POPs, an open form of porcine POP is also proposed with domain‐closing motion. This work has direct implications for the development of novel inhibitors of mammalian POPs to understand the etiology of various neurological diseases. Proteins 2014; 82:1428–1443. © 2014 Wiley Periodicals, Inc.     

Exploration micromechanism of VP35 IID interaction and recognition dsRNA: A molecular dynamics simulation

Multifunctional viral protein (VP35) encoded by the highly pathogenic Ebola viruses (EBOVs) can antagonize host double‐stranded RNA (dsRNA) sensors and immune response because of the simultaneous recognition of dsRNA backbone and blunt ends. Mutation of select hydrophobic conserved basic residues within the VP35 inhibitory domain (IID) abrogates its dsRNA‐binding activity, and impairs VP35‐mediated interferon (IFN) antagonism. Herein the detailed binding mechanism between dsRNA and WT, single mutant, and double mutant were investigated by all‐atom molecular dynamics (MD) simulation and binding energy calculation. R312A/R322A double mutations results in a completely different binding site and orientation upon the structure analyses. The calculated binding free energy results reveal that R312A, R322A, and K339A single mutations decrease the binding free energies by 17.82, 13.18, and 13.68 kcal mol^?1, respectively. The binding energy decomposition indicates that the strong binding affinity of the key residues is mainly due to the contributions of electrostatic interactions in the gas phase, where come from the positively charged side chain and the negatively charged dsRNA backbone. R312A, R322A, and K339A single mutations have no significant effect on VP35 IID conformation, but the mutations influence the contributions of electrostatic interactions in the gas phase. The calculated results reveal that end‐cap residues which mainly contribute VDW interactions can recognize and capture dsRNA blunt ends, and the central basic residues (R312, R322, and K339) which mainly contribute favorable electrostatic interactions with dsRNA backbone can fix dsRNA binding site and orientation. Proteins 2017; 85:1008–1023. © 2017 Wiley Periodicals, Inc.     

KCTD5 is endowed with large,functionally relevant,interdomain motions

The KCTD family is an emerging class of proteins that are involved in important biological processes whose biochemical and structural properties are rather poorly characterized or even completely undefined. We here used KCTD5, the only member of the family with a known three-dimensional structure, to gain insights into the intrinsic structural stability of the C-terminal domain (CTD) and into the mutual dynamic interplay between the two domains of the protein. Molecular dynamics (MD) simulations indicate that in the simulation timescale (120 ns), the pentameric assembly of the CTD is endowed with a significant intrinsic stability. Moreover, MD analyses also led to the identification of exposed β-strand residues. Being these regions intrinsically sticky, they could be involved in the substrate recognition. More importantly, simulations conducted on the full-length protein provide interesting information of the relative motions between the BTB domain and the CTD of the protein. Indeed, the dissection of the overall motion of the protein is indicative of a large interdomain twisting associated with limited bending movements. Notably, MD data indicate that the entire interdomain motion is pivoted by a single residue (Ser150) of the hinge region that connects the domains. The functional relevance of these motions was evaluated in the context of the functional macromolecular machinery in which KCTD5 is involved. This analysis indicates that the interdomain twisting motion here characterized may be important for the correct positioning of the substrate to be ubiquitinated with respect to the other factors of the ubiquitination machinery.     

Crystal structure of the phospholipase A and acyltransferase 4 (PLAAT4) catalytic domain

Phospholipase A and Acyltransferase 4 (PLAAT4) is a class II tumor suppressor, that also plays a role as a restrictor of intracellular Toxoplasma gondii infection through restriction of parasitic vacuole size. The catalytic N-terminal domain (NTD) interacts with the C-terminal domain (CTD), which is important for sub-cellular targeting and enzymatic function. The dynamics of the NTD main (L1) loop and the L2(B6) loop adjacent to the active site, have been shown to be important regulators of enzymatic activity. Here, we present the crystal structure of PLAAT4 NTD, determined from severely intergrown crystals using automated, laser-based crystal harvesting and data reduction technologies. The structure showed the L1 loop in two distinct conformations, highlighting a complex network of interactions likely influencing its conformational flexibility. Ensemble refinement of the crystal structure recapitulates the major correlated motions observed in solution by NMR. Our analysis offers useful insights on millisecond dynamics based on the crystal structure, complementing NMR studies which preclude structural information at this time scale.     

Specificity and cooperativity at β-lactamase position 104 in TEM-1/BLIP and SHV-1/BLIP interactions

Establishing a quantitative understanding of the determinants of affinity in protein–protein interactions remains challenging. For example, TEM‐1/β‐lactamase inhibitor protein (BLIP) and SHV‐1/BLIP are homologous β‐lactamase/β‐lactamase inhibitor protein complexes with disparate K_d values (3 nM and 2 μM, respectively), and a single substitution, D104E in SHV‐1, results in a 1000‐fold enhancement in binding affinity. In TEM‐1, E104 participates in a salt bridge with BLIP K74, whereas the corresponding SHV‐1 D104 does not in the wild type SHV‐1/BLIP co‐structure. Here, we present a 1.6 Å crystal structure of the SHV‐1 D104E/BLIP complex that demonstrates that this point mutation restores this salt bridge. Additionally, mutation of a neighboring residue, BLIP E73M, results in salt bridge formation between SHV‐1 D104 and BLIP K74 and a 400‐fold increase in binding affinity. To understand how this salt bridge contributes to complex affinity, the cooperativity between the E/K or D/K salt bridge pair and a neighboring hot spot residue (BLIP F142) was investigated using double mutant cycle analyses in the background of the E73M mutation. We find that BLIP F142 cooperatively stabilizes both interactions, illustrating how a single mutation at a hot spot position can drive large perturbations in interface stability and specificity through a cooperative interaction network. Proteins 2011. © 2011 Wiley‐Liss, Inc.