RESUMEN

G protein-coupled receptors (GPCRs) are membrane-bound signaling proteins that play an essential role in cellular signaling processes. Due to their intrinsic function of transmitting internal signals in response to external cues, these receptors are adapted to be highly dynamic in nature. The ß2 -adrenergic receptor (ß2 AR) is a representative member of the family that has been extensively analyzed in terms of its structure and activation. Although the structure of the transmembrane domain has been characterized in the different functional states of the receptor, the conformational dynamics of the extramembrane domains, especially the intrinsically disordered regions are still emerging. In this study, we analyze the state-dependent dynamics of extramembrane domains of ß2 AR using atomistic molecular dynamics simulations. We introduce a parameter, the residue excess dynamics that allows us to better quantify receptor dynamics. Using this measure, we show that the dynamics of the extramembrane domains are sensitive to the receptor state. Interestingly, the ligand-bound intermediate R ' state shows the maximal dynamics compared to either the active R*G or inactive R states. Ligand binding appears to be correlated with high residue excess dynamics that are dampened upon G protein coupling. The intracellular loop-3 (ICL3) domain has a tendency to flip towards the membrane upon ligand binding, which could contribute to receptor "priming." We highlight an important ICL1-helix-8 interplay that is broken in the ligand-bound state but is retained in the active state. Overall, our study highlights the importance of characterizing the functional dynamics of the GPCR loop domains.

Asunto(s)

RESUMEN

Human P-glycoprotein (P-gp) is a kind of ATP-binding cassette (ABC) transporters. Once human P-gp is overexpressed in tumor cells, which can lead to tumor multidrug resistance (MDR). However, the present experimental methods are difficult to obtain the large-scale conformational transition process of human P-gp. In this work, we explored the allosteric pathway of human P-gp from the inward-facing (IF) to the outward-facing (OF) state in the substrate transport process with the two-state anisotropic network model (tANM). These results suggest that the allosteric transitions proceed in a coupled way. The conformational changes of nucleotide-binding domains (NBDs) finally make the transmembrane domains (TMDs) to the OF state via the role of the allosteric propagation of the intracellular helices IH1 and IH2. Additionally, this allosteric pathway is advantageous in energy compared with other methods. This study reveals the conformational transition of P-gp, which contributes to an understanding of the allosteric mechanism of ABC exporters.

RESUMEN

The Arabidopsis Serine/Threonine Kinase 1 (SIK1) is a Sterile 20 (STE20)/Hippo orthologue that is also categorized as a Mitogen-Activated Protein Kinase Kinase Kinase Kinase (MAP4K). Like its animal and fungi orthologues, SIK1 is required for cell cycle exit, cell expansion, polarity establishment, as well as pathogenic response. The catalytic activity of SIK1, like other MAPKs, is presumably regulated by its phosphorylation states. Since no crystal structure for SIK1 has been reported yet, we built structural models for SIK1 kinase domain in different phosphorylation states with different pocket conformation to see how this kinase may be regulated. Using computational structural biology methods, we outlined a conduction path in which a phosphorylation site on the A-loop regulates the catalytic activity of SIK1 by controlling the closing or opening of the catalytic pocket at the G-loop. Furthermore, with analyses on the dynamic motions and in vitro kinase assay, we confirmed that three key residues in this conduction path, Lys278, Glu295, and Arg370, are indeed important for the kinase activity of SIK1. Since these residues are conserved in all STE20 kinases examined, the regulatory mechanism that we discovered may be common in STE20 kinases.

RESUMEN

Protein evolution proceeds by a complex response of organismal fitness to mutations that can simultaneously affect protein stability, structure, and enzymatic activity. To probe the relationship between genotype and phenotype, we chose a fundamental paradigm for protein evolution, folding, and design, the (ßα)8 TIM barrel fold. Here, we demonstrate the role of long-range allosteric interactions in the adaptation of an essential hyperthermophilic TIM barrel enzyme to mesophilic conditions in a yeast host. Beneficial fitness effects observed with single and double mutations of the canonical ßα-hairpin clamps and the α-helical shell distal to the active site revealed an underlying energy network between opposite faces of the cylindrical ß-barrel. We experimentally determined the fitness of multiple mutants in the energetic phase plane, contrasting the energy barrier of the chemical reaction and the folding free energy of the protein. For the system studied, the reaction energy barrier was the primary determinant of organism fitness. Our observations of long-range epistatic interactions uncovered an allosteric pathway in an ancient and ubiquitous enzyme that may provide a novel way of designing proteins with a desired activity and stability profile.

Asunto(s)

RESUMEN

With the increasing difficulty to develop new drugs and the emergence of resistance to traditional orthosteric-site inhibitors, the search for alternatives is finally approaching the focus on allosteric sites. Allosteric sites offer opportunities to regulate many pharmacologically targeted pathways by inhibition or activation. In addition, allosteric sites tend to be less conserved than the functional site, which may facilitate the design of specific effectors in the protein families for which specific orthosteric inhibitors have proved difficult to design. Furthermore, recent evidence suggests that all proteins might be susceptible of allosteric regulation, increasing the space of druggable targets. Computational identification of allosteric sites has therefore become an active field of research. The problem can be approached from two sides: (1) the identification of allosteric-communication pathways between the functional site and potential allosteric sites and (2) the functional-site-independent identification of allosteric sites. While the first approach tends to be more laborious and thus restricted to a single protein, the second tends to be more amenable to larger-scale analysis, thus providing tools for the two drug discovery scenarios: the analysis of known targets and the screening for new potential targets. Here, I show some basic concepts and methods useful to the identification of allosteric sites and pathways, in line with these two approaches. I describe them in some detail to build a clear framework, at the risk of losing the interest of experts. Examples of recent studies involving these methods are also illustrated, focusing on the techniques rather than on their findings on allosterism.

Asunto(s)

RESUMEN

Glycogen synthase kinase 3ß (GSK3ß), a serine/threonine protein kinase, is involved in several human diseases, including type II diabetes, mood disorders, prostate cancer, and Alzheimer's disease, representing a potential therapeutic target. GSK3ß has a unique specificity, with its primed substrates binding to the primed phosphate binding site, which is critical for the catalytic activity of GSK3ß. An L343R mutation located at the C-lobe of GSK3ß, remote from the catalytic site, causes kinase inactivation. However, the detailed mechanism of this remains unclear. Here, microsecond molecular dynamics (MD) simulations and network analysis were performed to elucidate the allosteric inactivation of GSK3ß triggered by the L343R mutation. Large-scale MD simulations of wild-type and the L343R mutant revealed that the L343R mutation caused disruption of the chemical environment near the mutation site, which propagated remotely to affect the conformational dynamics of the activation loop (A-loop). The resulting conformational rearrangement of the A-loop in the L343R mutant disrupted the primed phosphate binding site, thereby abrogating the catalytic activity of GSK3ß. Furthermore, network analysis identified the allosteric pathway from R343 to the primed phosphate binding site in the L343R mutant. Collectively, the results of this study provide a mechanistic explanation of how the L343R mutation allosterically affects the functional activity of GSK3ß, which contributes to our understanding of GSK3ß biology.

Asunto(s)

RESUMEN

Ras proteins, as small GTPases, mediate cell proliferation, survival and differentiation. Ras mutations have been associated with a broad spectrum of human cancers and thus targeting Ras represents a potential way forward for cancer therapy. A recently reported monobody NS1 allosterically disrupts the Ras-mediated signaling pathway, but its efficacy is reduced by R135K mutation in H-Ras. However, the detailed mechanism is unresolved. Here, using molecular dynamics (MD) simulations and dynamic network analysis, we explored the molecular mechanism for the unbinding of NS1 to H-Ras and shed light on the underlying allosteric network in H-Ras. MD simulations revealed that the overall structures of the two complexes did not change significantly, but the H-Ras-NS1 interface underwent significant conformational alteration in the mutant Binding free energy analysis showed that NS1 binding was unfavored after R135K mutation, which resulted in the unfavorable binding of NS1. Furthermore, the critical residues on H-Ras responsible for the loss of binding of NS1 were identified. Importantly, the allosteric networks for these important residues were revealed, which yielded a novel insight into the allosteric regulatory mechanism of H-Ras.

Asunto(s)

RESUMEN

This study identifies dynamical properties of maltose-binding protein (MBP) useful in unveiling active site residues susceptible to ligand binding. The described methodology has been previously used in support of novel topological techniques of persistent homology and statistical inference in complex, multi-scale, high-dimensional data often encountered in computational biophysics. Here we outline a computational protocol that is based on the anisotropic elastic network models of 14 all-atom three-dimensional protein structures. We introduce the notion of dynamical distance matrices as a measure of correlated interactions among 370 amino acid residues that constitute a single protein. The dynamical distance matrices serve as an input for a persistent homology suite of codes to further distinguish a small subset of residues with high affinity for ligand binding and allosteric activity. In addition, we show that ligand-free closed MBP structures require lower deformation energies than open MBP structures, which may be used in categorization of time-evolving molecular dynamics structures. Analysis of the most probable allosteric coupling pathways between active site residues and the protein exterior is also presented.

Asunto(s)

RESUMEN

Methylglyoxal synthase (MGS) is a homohexameric enzyme responsible for converting dihydroxyacetone phosphate (DHAP) to methylglyoxal and phosphate in the methylglyoxal bypass of glycolysis. Phosphate acts as an allosteric inhibitor and strong regulator for this enzyme. Previous studies on MGS from Thermus sp. GH5 (TMGS) had indicated a pathway for transmitting the signal through Pro82, Arg97 and Val101 to the active site. The necessity of these residues for heterotropic negative cooperativity between subunits of TMGS were also proposed. In this study, it has been shown that a path via a salt bridge between Arg80 and Asp100 in the narrow dimer interface provides an alternative pathway for transmission of the allosteric inhibitory signal through subunit interfaces.

Asunto(s)

RESUMEN

The transporter MsbA is a kind of multidrug resistance ATP-binding cassette transporter that can transport lipid A, lipopolysaccharides, and some amphipathic drugs from the cytoplasmic to the periplasmic side of the inner membrane. In this work, we explored the allosteric pathway of MsbA from the inward- to outward-facing states during the substrate transport process with the adaptive anisotropic network model. The results suggest that the allosteric transitions proceed in a coupled way. The large-scale closing motions of the nucleotide-binding domains occur first, accompanied with a twisting motion at the same time, which becomes more obvious in middle and later stages, especially for the later. This twisting motion plays an important role for the rearrangement of transmembrane helices and the opening of transmembrane domains on the periplasmic side that mainly take place in middle and later stages respectively. The topological structure plays an important role in the motion correlations above. The conformational changes of nucleotide-binding domains are propagated to the transmembrane domains via the intracellular helices IH1 and IH2. Additionally, the movement of the transmembrane domains proceeds in a nonrigid body, and the two monomers move in a symmetrical way, which is consistent with the symmetrical structure of MsbA. These results are helpful for understanding the transport mechanism of the ATP-binding cassette exporters.

Asunto(s)

RESUMEN

The amino acid sequences of apolipoprotein E (apoE) from 63 different mammalian species have been downloaded from the protein database. The sequences were compared to human apoE4 to determine conserved and non-conserved sequences of amino acids. ApoE4 is the major risk factor for the development of late onset Alzheimer's disease while apoE3, which differs from apoE4 by a single amino acid change at position 112, poses little or no risk for the development of this disease. Thus, the two proteins appear to be structurally and functionally different. Seven highly conserved regions, representing approximately 47 amino acids (of 299) have been found. These regions are distributed throughout the protein and reflect ligand binding sites as well as regions proposed to be involved in the propagation of the cysteine-arginine change at position 112 to distant regions of the protein in the N- and C-terminal domains. Highly non-conserved regions are at the N- and C-terminal ends of the apoE protein.

Asunto(s)

RESUMEN

The homohexameric enzyme methylglyoxal synthase (MGS) converts dihydroxyacetone phosphate (DHAP) to methylglyoxal and phosphate. This enzyme is allosterically inhibited by phosphate. The allosteric signal induced by phosphate in MGS from Thermus sp. GH5 (TMGS) has been tracked by site-directed mutagenesis, from the binding site of phosphate to the pathways that transmit the signal, and finally to the active site which is the receiver of the signal. In TMGS, Ser-55 distinguishes the inhibitory phosphate from the phosphoryl group of the substrate, DHAP, and transmits the allosteric signal through Pro-82, Arg-97 and Val-101 to the active site. Furthermore, the addition of a C-terminal tail to TMGS reinforces the allosteric signal by introducing a new salt bridge between Asp-10 and an Arg in this tail. Lastly, the active site amino acid, Gly-56, is shown to be involved in both allostery and phosphate elimination step from DHAP by TMGS. Interestingly, some of the mutations also trigger homotropic allostery, supporting the hypothesis that allostery is an intrinsic property of all dynamic proteins. The details of the TMGS allosteric network discussed in this study can serve as a model system for understanding the enigmatic allosteric mechanism of other proteins.

Asunto(s)