RESUMEN

Sec-pathway is the main protein secretion pathway in prokaryotes and is essential for their survival. The motor protein SecA is the main coordinator of the pathway in bacteria as it is has evolved to perform multiple tasks, acting like a "swiss army knife", from binding pre-proteins to altering its oligomeric and conformational states. This study focuses on the role of its Preprotein Binding Domain (PBD), which is a key protein module that identified in three conformational states (WO, O and C). A thorough analysis was conducted to identify PBD's inter- and intra-protomeric interactions, highlighting the most significant and conserved ones. Both crystallographic and biophysical data indicate that the WO state is the main during dimerization, while the monomeric structure can adopt all three states. C-tail, StemPBD and 3ß-tipPBD are important elements for the stabilization of different oligomeric and conformational states, as they offer specific interactions. Alterations in the lipophilicity of the StemPBD causes increased proteins dynamics or/and Prl phenotype. In the C state, 3ß-tipPBD interacts and opens the ATPase motor. We hypothesize that this partial opening of the motor with the increased dynamics describes the Prl phenotype.

RESUMEN

The Fragile X messenger ribonucleoprotein (FMRP) is a multi-domain protein involved in interactions with various macromolecules, including proteins and coding/non-coding RNAs. The three KH domains (KH0, KH1 and KH2) within FMRP are recognized for their roles in mRNA binding. In the context of Fragile X syndrome (FXS), over-and-above CGG triplet repeats expansion, three specific point mutations have been identified, each affecting one of the three KH domains (R138QKH0, G266EKH1, and I304NKH2) resulting in the expression of non-functional FMRP. This study aims to elucidate the molecular mechanism underlying the loss of function associated with the G266EKH1 pathological variant. We investigate the conformational and dynamic properties of the isolated KH1 domain and the two KH1 site-directed mutants G266EKH1 and G266AKH1. Employing a combined in vitro and in silico approach, we reveal that the G266EKH1 variant lacks the characteristic features of a folded domain. This observation provides an explanation for functional impairment observed in FMRP carrying the G266E mutation within the KH1 domain, as it renders the domain unable to fold properly. Molecular Dynamics simulations suggest a pivotal role for residue 266 in regulating the structural stability of the KH domains, primarily through stabilizing the α-helices of the domain. Overall, these findings enhance our comprehension of the molecular basis for the dysfunction associated with the G266EKH1 variant in FMRP.

Asunto(s)

Proteína de la Discapacidad Intelectual del Síndrome del Cromosoma X Frágil , Síndrome del Cromosoma X Frágil , Proteína de la Discapacidad Intelectual del Síndrome del Cromosoma X Frágil/genética , Proteína de la Discapacidad Intelectual del Síndrome del Cromosoma X Frágil/metabolismo , Proteína de la Discapacidad Intelectual del Síndrome del Cromosoma X Frágil/química , Síndrome del Cromosoma X Frágil/genética , Síndrome del Cromosoma X Frágil/metabolismo , Humanos , Dominios Proteicos , Simulación de Dinámica Molecular , Conformación Proteica , Mutagénesis Sitio-Dirigida

RESUMEN

Peptide-lipid interactions play an important role in maintaining the integrity and function of the cell membrane. Even slight changes in these interactions can induce the development of various diseases. Specifically, peptide misfolding and aggregation in the membrane is considered to be one of the triggers of Alzheimer's disease (AD), however its exact mechanism is still unclear. To this end, an increase of amyloid-beta (Aß) peptide concentration in the human brain is widely accepted to gradually produce cytotoxic Aß aggregates (plaques). These plaques initiate a sequence of pathogenic events ending up in observable symptoms of dementia. Understanding the mechanism of the Aß interaction with cells is crucial for early detection and prevention of Alzheimer's disease. Hence, in this work, a comprehensive Raman analysis of the Aß42 conformational dynamics in water and in liposomes and lipodiscs that mimic the membrane system is presented. The obtained results show that the secondary structure of Aß42 in liposomes is dominated by the α-helix conformation, which remains stable over time. However, it comes as a surprise to reveal that the lipodisc environment induces the transformation of the Aß42 secondary structure to a ß-turn/random coil. Our Raman spectroscopy findings are supported with molecular dynamics (MD) and density functional theory (DFT) simulations, showing their good agreement.Communicated by Ramaswamy H. Sarma.

RESUMEN

Oncogenic Ras proteins are known to present multiple conformational states, as reported by the great variety of crystallographic structures. The GTP-bound states are grouped into two main states: the "inactive" state 1 and the "active" state 2. Recent reports on H-Ras have shown that state 2 exhibits two substates, directly related to the orientation of Tyr32: toward the GTP-bound pocket and outwards. In this paper, we show that N-Ras exhibits another substate of state 2, related to a third orientation of Tyr32, toward Ala18 and parallel to the GTP-bound pocket. We also show that this substate is highly sampled in the G12V mutation of N-Ras and barely present in its wild-type form, and that the G12V mutation prohibits the sampling of the GTPase-activating protein (GAP) binding substate, rendering this mutation oncogenic. Furthermore, using molecular dynamics simulations, we explore the importance of the membrane on N-Ras' conformational state dynamics and its strong influence on Ras protein stability. Moreover, the membrane has a significant influence on the conformational (sub)states sampling of Ras. This, in turn, is of crucial importance in the activation/deactivation cycle of Ras, due to the binding of guanine nucleotide exchange factor proteins (GEFs)/GTPase-activating proteins (GAPs).

Asunto(s)

Factores de Intercambio de Guanina Nucleótido , Mutación Puntual , Proteínas Proto-Oncogénicas p21(ras) , Factores de Intercambio de Guanina Nucleótido/genética , Guanosina Trifosfato/metabolismo , Mutación , Proteínas ras/metabolismo , Proteínas Proto-Oncogénicas p21(ras)/genética , Proteínas Proto-Oncogénicas p21(ras)/fisiología , Simulación de Dinámica Molecular

RESUMEN

Prolidase (EC 3.4.13.9) is an enzyme that specifically hydrolyzes Xaa-Pro dipeptides into free amino acids. We previously studied kinetic behaviours and solved the crystal structure of wild-type (WT) Lactococcus lactis prolidase (Llprol), showing that this homodimeric enzyme has unique characteristics: allosteric behaviour and substrate inhibition. In this study, we focused on solving the crystal structures of three Llprol mutants (D36S, H38S, and R293S) which behave differently in v-S plots. The D36S and R293S Llprol mutants do not show allosteric behaviour, and the Llprol mutant H38S has allosteric behaviour comparable to the WT enzyme (Hill constant 1.52 and 1.58, respectively). The crystal structures of Llprol variants suggest that the active site of Llprol formed with amino acid residues from both monomers, i.e., located in an interfacial area of dimer. The comparison between the structure models of Llprol indicated that the two monomers in the dimers of Llprol variants have different relative positions among Llprol variants. They showed different interatomic distances between the amino acid residues bridging the two monomers and varied sizes of the solvent-accessible interface areas in each Llprol variant. These observations indicated that Llprol could adapt to different conformational states with distinctive substrate affinities. It is strongly speculated that the domain movements required for productive substrate binding are restrained in allosteric Llprol (WT and H38S). At low substrate concentrations, only one out of the two active sites at the dimer interface could accept substrate; as a result, the asymmetrical activated dimer leads to allosteric behaviour.

Asunto(s)

Dipeptidasas , Lactococcus lactis , Regulación Alostérica , Lactococcus lactis/genética , Lactococcus lactis/metabolismo , Especificidad por Sustrato , Modelos Moleculares , Aminoácidos/metabolismo

RESUMEN

Full-length Bruton's tyrosine kinase (BTK) has been refractory to structural analysis. The nearest full-length structure of BTK to date consists of the autoinhibited SH3-SH2-kinase core. Precisely how the BTK N-terminal domains (the Pleckstrin homology/Tec homology [PHTH] domain and proline-rich regions [PRR] contain linker) contribute to BTK regulation remains unclear. We have produced crystals of full-length BTK for the first time but despite efforts to stabilize the autoinhibited state, the diffraction data still reveal only the SH3-SH2-kinase core with no electron density visible for the PHTH-PRR segment. Cryo-electron microscopy (cryoEM) data of full-length BTK, on the other hand, provide the first view of the PHTH domain within full-length BTK. CryoEM reconstructions support conformational heterogeneity in the PHTH-PRR region wherein the globular PHTH domain adopts a range of states arrayed around the autoinhibited SH3-SH2-kinase core. On the way to activation, disassembly of the SH3-SH2-kinase core opens a new autoinhibitory site on the kinase domain for PHTH domain binding that is ultimately released upon interaction of PHTH with phosphatidylinositol (3,4,5)-trisphosphate. Membrane-induced dimerization activates BTK and we present here a crystal structure of an activation loop swapped BTK kinase domain dimer that likely represents the conformational state leading to trans-autophosphorylation. Together, these data provide the first structural elucidation of full-length BTK and allow a deeper understanding of allosteric control over the BTK kinase domain during distinct stages of activation.

Asunto(s)

Agammaglobulinemia Tirosina Quinasa , Microscopía por Crioelectrón , Dominios Proteicos , Fosforilación , Dimerización

RESUMEN

The breakthrough of AlphaFold2 and the publication of AlphaFold DB represent a significant advance in the field of predicting static protein structures. However, AlphaFold2 models tend to represent a single static structure, and multiple-conformation prediction remains a challenge. In this work, we proposed a method named MultiSFold, which uses a distance-based multi-objective evolutionary algorithm to predict multiple conformations. To begin, multiple energy landscapes are constructed using different competing constraints generated by deep learning. Subsequently, an iterative modal exploration and exploitation strategy is designed to sample conformations, incorporating multi-objective optimization, geometric optimization and structural similarity clustering. Finally, the final population is generated using a loop-specific sampling strategy to adjust the spatial orientations. MultiSFold was evaluated against state-of-the-art methods using a benchmark set containing 80 protein targets, each characterized by two representative conformational states. Based on the proposed metric, MultiSFold achieves a remarkable success ratio of 56.25% in predicting multiple conformations, while AlphaFold2 only achieves 10.00%, which may indicate that conformational sampling combined with knowledge gained through deep learning has the potential to generate conformations spanning the range between different conformational states. In addition, MultiSFold was tested on 244 human proteins with low structural accuracy in AlphaFold DB to test whether it could further improve the accuracy of static structures. The experimental results demonstrate the performance of MultiSFold, with a TM-score better than that of AlphaFold2 by 2.97% and RoseTTAFold by 7.72%. The online server is at http://zhanglab-bioinf.com/MultiSFold .

Asunto(s)

Algoritmos , Conformación Proteica , Proteínas , Proteínas/química , Humanos , Biología Computacional/métodos , Modelos Moleculares , Aprendizaje Profundo

RESUMEN

Artificial intelligence has made significant advances in the field of protein structure prediction in recent years. In particular, DeepMind's end-to-end model, AlphaFold2, has demonstrated the capability to predict three-dimensional structures of numerous unknown proteins with accuracy levels comparable to those of experimental methods. This breakthrough has opened up new possibilities for understanding protein structure and function as well as accelerating drug discovery and other applications in the field of biology and medicine. Despite the remarkable achievements of artificial intelligence in the field, there are still some challenges and limitations. In this Review, we discuss the recent progress and some of the challenges in protein structure prediction. These challenges include predicting multidomain protein structures, protein complex structures, multiple conformational states of proteins, and protein folding pathways. Furthermore, we highlight directions in which further improvements can be conducted.

Asunto(s)

Inteligencia Artificial , Descubrimiento de Drogas , Pliegue de Proteína , Proyectos de Investigación

RESUMEN

By assuming that a stepwise outward movement of the four S4 segments of the hERG potassium channel determines a concomitant progressive increase in the flow of the permeant potassium ions, the inward and outward potassium currents can be simulated by using only one or two adjustable (i.e., free) parameters. This deterministic kinetic model differs from the stochastic models of hERG available in the literature, which usually require more than 10 free parameters. The K+ outward current of hERG contributes to the repolarization of the cardiac action potential. On the other hand, the K+ inward current increases with a positive shift in the transmembrane potential, in apparent contrast to both the electric and osmotic forces, which would concur in moving K+ ions outwards. This peculiar behavior can be explained by the appreciable constriction of the central pore midway along its length, with a radius < 1 Å and hydrophobic sacks surrounding it, as reported in an open conformation of the hERG potassium channel. This narrowing raises a barrier to the outward movement of K+ ions, inducing them to move increasingly inwards under a gradually more positive transmembrane potential.

RESUMEN

The spike protein of SARS-CoV-2 exists in two major conformational states, namely the 'open' and 'closed' states which are also known as the 'up' and 'down' states, respectively. In its open state, the receptor binding domain (RBD) of the protein is exposed for binding with ACE2, whereas the spike RBD is inaccessible to ACE2 in the closed state of the protein. In the current work, we have performed all-atom microsecond simulations of the full-length trimeric spike protein solvated in explicit aqueous medium with an average system size of ~0.7 million atoms to understand the molecular nature of intra- and inter-chain interactions, water-bridged interactions between different residues that contribute to the stability of the open and closed states of the protein, and also the free energy landscape for transition between the open and closed states of the protein. We have also examined the changes of such interactions that are associated with switching from one state to the other through both unbiased and biased simulations at all-atom level with total run length of 4 µs. Interestingly, after about 0.8 µs of unbiased molecular dynamics run of the spike system in the open state, we observed a gradual transition of the monomeric chain (B) from open to its partially closed or down state. Initially the residues at the interface of chain B RBD in the open state spike protein were at non-hydrogen-bonding distances from the residues of chain C RBD. However, the two RBDs gradually came closer and finally the residue S459 of the RBD of chain B made a hydrogen bond with F374 of chain C in the last 200 ns of the simulation along with formation of a few more hydrogen bonds involving other residues. Since no transition from closed to the open state of the protein is observed in the present 1 µs unbiased simulation of the closed state protein, the current study seems to suggest that the closed conformational state is preferred for the spike protein of SARS-CoV-2 in aqueous medium. Furthermore, calculations of the free energy surface of the conformational transition from open (up) to the closed (down) state using a biased simulation method reveal a free energy barrier of ~3.20 kcal/mol for the transition of RBD from open to the closed state, whereas the barrier for the reverse process is found to be significantly higher.

Asunto(s)

COVID-19 , SARS-CoV-2 , Humanos , Enzima Convertidora de Angiotensina 2 , Glicoproteína de la Espiga del Coronavirus , Unión Proteica , Simulación de Dinámica Molecular

RESUMEN

We unveil the intimate relationship between protein dynamics and allostery by following the trajectories of model proteins in their conformational and sequence spaces. Starting from a nonallosteric hyperthermophilic malate dehydrogenase, we have tracked the role of protein dynamics in the evolution of the allosteric capacity. Based on a large phylogenetic analysis of the malate (MalDH) and lactate dehydrogenase (LDH) superfamily, we identified two amino acid positions that could have had a major role for the emergence of allostery in LDHs, which we targeted for investigation by site-directed mutagenesis. Wild-type MalDH and the single and double mutants were tested with respect to their substrate recognition profiles. The double mutant displayed a sigmoid-shaped profile typical of homotropic activation in LDH. By using molecular dynamics simulations, we showed that the mutations induce a drastic change in the protein sampling of its conformational landscape, making transiently T-like (inactive) conformers, typical of allosteric LDHs, accessible. Our data fit well with the seminal key concept linking protein dynamics and evolvability. We showed that the selection of a new phenotype can be achieved by a few key dynamics-enhancing mutations causing the enrichment of low-populated conformational substates.

Asunto(s)

Malato Deshidrogenasa , Malatos , Regulación Alostérica , Aminoácidos/genética , L-Lactato Deshidrogenasa/química , L-Lactato Deshidrogenasa/genética , L-Lactato Deshidrogenasa/metabolismo , Malato Deshidrogenasa/genética , Mutación , Filogenia

RESUMEN

This work analyses the thermostability of a membrane protein, the gastric H,K-ATPase, by means of a detailed kinetic characterization of its inactivation process, which showed to exhibit first-order kinetics. We observed parallel time courses for the decrease of ATPase activity, the decrease of the autophosphorylation capacity and the loss of tertiary structure at 49 °C. Higher temperatures were required to induce a significant change in secondary structure. The correspondence between the kinetics of Trp fluorescence measured at 49 °C and the decrease of the residual activity after heating at that temperature, proves the irreversibility of the inactivation process. Inactivation proceeds at different rates in E1 or E2 conformations. The K+-induced E2 state exhibits a lower inactivation rate; the specific effect is exerted with a K0.5 similar to that found at 25 °C, providing a further inkling that K+ occlusion by the H,K-ATPase is not really favoured. Increasing [H+] from pH 8 to pH 7, which possibly shifts the protein to E1, produces a subtle destabilizing effect on the H,K-ATPase. We performed a prediction of potential intramolecular interactions and found that the differential stability between E1 and E2 may be mainly explained by the higher number of hydrophobic interactions in the α- and ß-subunits of E2 conformation.

Asunto(s)

ATPasa Intercambiadora de Sodio-Potasio , Sodio , Cationes/metabolismo , ATPasa Intercambiadora de Hidrógeno-Potásio/metabolismo , Cinética , Sodio/metabolismo , ATPasa Intercambiadora de Sodio-Potasio/metabolismo

RESUMEN

A significant proportion of proteins comprise multiple domains. Domain-domain docking is a tool that predicts multi-domain protein structures when individual domain structures can be accurately predicted but when domain orientations cannot be predicted accurately. GalaxyDomDock predicts an ensemble of domain orientations from given domain structures by docking. Such information would also be beneficial in elucidating the functions of proteins that have multiple states with different domain orientations. GalaxyDomDock is an ab initio domain-domain docking method based on GalaxyTongDock, a previously developed protein-protein docking method. Infeasible domain orientations for the given linker are effectively screened out from the docked conformations by a geometric filter, using the Dijkstra algorithm. In addition, domain linker conformations are predicted by adopting a loop sampling method FALC. The proposed GalaxyDomDock outperformed existing ab initio domain-domain docking methods, such as AIDA and Rosetta, in performance tests on the Rosetta benchmark set of two-domain proteins. GalaxyDomDock also performed better than or comparable to AIDA on the AIDA benchmark set of two-domain proteins and two-domain proteins containing discontinuous domains, including the benchmark set in which each domain of the set was modeled by the recent version of AlphaFold. The GalaxyDomDock web server is freely available as a part of GalaxyWEB at http://galaxy.seoklab.org/domdock.

Asunto(s)

Uso de Internet , Dominios Proteicos , Proteínas , Programas Informáticos , Algoritmos , Simulación del Acoplamiento Molecular , Proteínas/química

RESUMEN

Human α2-macroglobulin (hα2M) is a multidomain protein with a plethora of essential functions, including transport of signaling molecules and endopeptidase inhibition in innate immunity. Here, we dissected the molecular mechanism of the inhibitory function of the ∼720-kDa hα2M tetramer through eight cryo­electron microscopy (cryo-EM) structures of complexes from human plasma. In the native complex, the hα2M subunits are organized in two flexible modules in expanded conformation, which enclose a highly porous cavity in which the proteolytic activity of circulating plasma proteins is tested. Cleavage of bait regions exposed inside the cavity triggers rearrangement to a compact conformation, which closes openings and entraps the prey proteinase. After the expanded-to-compact transition, which occurs independently in the four subunits, the reactive thioester bond triggers covalent linking of the proteinase, and the receptor-binding domain is exposed on the tetramer surface for receptor-mediated clearance from circulation. These results depict the molecular mechanism of a unique suicidal inhibitory trap.

Asunto(s)

Péptido Hidrolasas , alfa-Macroglobulinas , Microscopía por Crioelectrón , Endopeptidasas/metabolismo , Humanos , Péptido Hidrolasas/metabolismo , Conformación Proteica , Factores de Transcripción , alfa-Macroglobulinas/química , alfa-Macroglobulinas/metabolismo

RESUMEN

The cellular functions are regulated by a complex interplay of diffuse and local signals. Studying the latter is challenging, but experimental work in cell physiology has led to recognize that understanding a cell's dynamics requires a deep comprehension of local fluctuations of cytosolic regulators. Macromolecular complexes are major determinants of local signaling. Multienzyme assemblies limit the diffusion restriction to reaction kinetics by direct exchange of metabolites. Likewise, close coupling of ion channels and transporters modulates the ion concentration around a channel mouth or transporter binding site. Extreme signal locality is brought about by conformational coupling between membrane proteins, as is typical of mechanotransduction. A paradigmatic case is integrin-mediated cell adhesion. Sensing the extracellular microenvironment and providing an appropriate response are essential in growth and development and have innumerable pathological implications. The process involves bidirectional signal transduction by complex supramolecular structures that link integrin receptors to ion channels and transporters, growth factor receptors, cytoskeletal elements, and other regulatory elements. The dynamics of such complexes are only beginning to be understood. A thoroughly studied example is the association between integrin receptors and the voltage-gated K+ channels Kv11.1. These channels are widely expressed in early embryos, where their physiological roles are poorly understood and apparently different from the shaping of action potential firing in the adult. Hints about these roles come from studies in cancer cells, where Kv11.1 is often overexpressed and appears to reassume functions it presumably exerts during embryogenesis, such as controlling cell proliferation/differentiation, apoptosis, and migration. Kv11.1 is implicated in these processes through its linking to integrin subunits, which in turn regulates channel expression. Specific cellular functions, such as proliferation and migration, appear to be modulated by distinct conformational states of the channel (e.g., open and closed), whose balance is affected by the link with integrin subunits.

Asunto(s)

Canales Iónicos , Mecanotransducción Celular , Integrinas/metabolismo , Canales Iónicos/metabolismo , Cinética , Transducción de Señal/fisiología

RESUMEN

Dopamine transporter mediates the neurotransmitter dopamine homeostasis in a sodium-dependent manner. The transport process involves an alternating access of a substrate to the extracellular and intracellular spaces, which is associated with different conformational states of the transporter. However, the underlying mechanism of modulation of the state transition remains elusive. Here we present a computational simulation study of human dopamine transporter to explore its two end states (outward-facing open and inward-facing open) that have not been determined experimentally. We show that the full-length transporter may tend to adopt the inward-facing open state in its free state. The binding of an amphetamine may not trap the transporter in the outward-facing open state with increasing length of the N-terminal. Furthermore, we identify distinct patterns in the interaction networks between the N-terminal and the intracellular region that could stabilize the state of the transporter, independent of substrate binding and phosphorylation. Our results reveal the essential role of the N-terminal dynamics in modulating the functional states of the dopamine transporter, providing molecular insights into the coupling of conformational transition and substrate passage in neurotransmitter transporters.

Asunto(s)

Proteínas de Transporte de Dopamina a través de la Membrana Plasmática , Simulación de Dinámica Molecular , Dopamina/metabolismo , Proteínas de Transporte de Dopamina a través de la Membrana Plasmática/química , Proteínas de Transporte de Dopamina a través de la Membrana Plasmática/metabolismo , Humanos , Conformación Proteica , Sodio

RESUMEN

Nucleic acids do not fold into a single conformation, and dynamic ensembles are needed to describe their propensities to cycle between different conformations when performing cellular functions. We review recent advances in solution-state nuclear magnetic resonance (NMR) methods and their integration with computational techniques that are improving the ability to probe the dynamic ensembles of DNA and RNA. These include computational approaches for predicting chemical shifts from structure and generating conformational libraries from sequence, measurements of exact nuclear Overhauser effects, development of new probes to study chemical exchange using relaxation dispersion, faster and more sensitive real-time NMR techniques, and new NMR approaches to tackle large nucleic acid assemblies. We discuss how these advances are leading to new mechanistic insights into gene expression and regulation.

Asunto(s)

Ácidos Nucleicos , ADN , Espectroscopía de Resonancia Magnética , Resonancia Magnética Nuclear Biomolecular , ARN

RESUMEN

INTRODUCTION: E-selectin is a member of the selectin family of cell adhesion molecules expressed on the plasma membrane of inflamed endothelium and facilitates initial leukocyte tethering and subsequent cell rolling during the early stages of the inflammatory response via binding to glycoproteins expressing sialyl LewisX and sialyl LewisA (sLeX/A). Existing crystal structures of the extracellular lectin/EGF-like domain of E-selectin complexed with sLeX have revealed that E-selectin can exist in two conformation states, a low affinity (bent) conformation, and a high affinity (extended) conformation. The differentiating characteristic of the two conformations is the interdomain angle between the lectin and the EGF-like domain. METHODS: Using molecular dynamics (MD) simulations we observed that in the absence of tensile force E-selectin undergoes spontaneous switching between the two conformational states at equilibrium. A single amino acid substitution at residue 2 (serine to tyrosine) on the lectin domain favors the extended conformation. RESULTS: Steered molecular dynamics (SMD) simulations of E-selectin and PSGL-1 in conjunction with experimental cell adhesion assays show a longer binding lifetime of E-selectin (S2Y) to PSGL-1 compared to wildtype protein. CONCLUSIONS: The findings in this study advance our understanding into how the structural makeup of E-selectin allosterically influences its adhesive dynamics.

RESUMEN

Bromodomains (BDs) are small protein modules that interact with acetylated marks in histones. These posttranslational modifications are pivotal to regulate gene expression, making BDs promising targets to treat several diseases. While the general structure of BDs is well known, their dynamical features and their interplay with other macromolecules are poorly understood, hampering the rational design of potent and selective inhibitors. Here, we combine extensive molecular dynamics simulations, Markov state modeling, and available structural data to reveal a transiently formed state that is conserved across all BD families. It involves the breaking of two backbone hydrogen bonds that anchor the ZA-loop with the αA helix, opening a cryptic pocket that partially occludes the one associated to histone binding. By analyzing more than 1,900 experimental structures, we unveil just two adopting the hidden state, explaining why it has been previously unnoticed and providing direct structural evidence for its existence. Our results suggest that this state is an allosteric regulatory switch for BDs, potentially related to a recently unveiled BD-DNA-binding mode.

Asunto(s)

Proteínas de Ciclo Celular/química , Proteínas Co-Represoras/química , Proteínas de Unión al ADN/química , Histona Acetiltransferasas/química , Péptidos y Proteínas de Señalización Intracelular/química , Factores Generales de Transcripción/química , Factores de Transcripción/química , Proteína 28 que Contiene Motivos Tripartito/química , Secuencia de Aminoácidos , Sitios de Unión , Proteínas de Ciclo Celular/genética , Proteínas de Ciclo Celular/metabolismo , Proteínas Co-Represoras/genética , Proteínas Co-Represoras/metabolismo , Cristalografía por Rayos X , ADN/química , ADN/genética , ADN/metabolismo , Proteínas de Unión al ADN/genética , Proteínas de Unión al ADN/metabolismo , Regulación de la Expresión Génica , Histona Acetiltransferasas/genética , Histona Acetiltransferasas/metabolismo , Humanos , Péptidos y Proteínas de Señalización Intracelular/genética , Péptidos y Proteínas de Señalización Intracelular/metabolismo , Cadenas de Markov , Simulación de Dinámica Molecular , Unión Proteica , Conformación Proteica en Hélice alfa , Conformación Proteica en Lámina beta , Dominios y Motivos de Interacción de Proteínas , Alineación de Secuencia , Homología de Secuencia de Aminoácido , Termodinámica , Factores de Transcripción/genética , Factores de Transcripción/metabolismo , Factores Generales de Transcripción/genética , Factores Generales de Transcripción/metabolismo , Proteína 28 que Contiene Motivos Tripartito/genética , Proteína 28 que Contiene Motivos Tripartito/metabolismo

RESUMEN

The COVID-19 pandemic has spread worldwide. However, as soon as the first vaccines-the only scientifically verified and efficient therapeutic option thus far-were released, mutations combined into variants of SARS-CoV-2 that are more transmissible and virulent emerged, raising doubts about their efficiency. This study aims to explain possible molecular mechanisms responsible for the increased transmissibility and the increased rate of hospitalizations related to the new variants. A combination of theoretical methods was employed. Constant-pH Monte Carlo simulations were carried out to quantify the stability of several spike trimeric structures at different conformational states and the free energy of interactions between the receptor-binding domain (RBD) and angiotensin-converting enzyme II (ACE2) for the most worrying variants. Electrostatic epitopes were mapped using the PROCEEDpKa method. These analyses showed that the increased virulence is more likely to be due to the improved stability to the S trimer in the opened state, in which the virus can interact with the cellular receptor, ACE2, rather than due to alterations in the complexation RBD-ACE2, since the difference observed in the free energy values was small (although more attractive in general). Conversely, the South African/Beta variant (B.1.351), compared with the SARS-CoV-2 wild type (wt), is much more stable in the opened state with one or two RBDs in the up position than in the closed state with three RBDs in the down position favoring the infection. Such results contribute to understanding the natural history of disease and indicate possible strategies for developing new therapeutic molecules and adjusting the vaccine doses for higher B-cell antibody production.