Zinc-mediated conformational preselection mechanism in the allosteric control of DNA binding to the zinc transcriptional regulator (ZitR)

Modeling Zinc Complexes Using Neural Networks

Understanding the energetic landscapes of large molecules is necessary for the study of chemical and biological systems. Recently, deep learning has greatly accelerated the development of models based on quantum chemistry, making it possible to build potential energy surfaces and explore chemical space. However, most of this work has focused on organic molecules due to the simplicity of their electronic structures as well as the availability of data sets. In this work, we build a deep learning architecture to model the energetics of zinc organometallic complexes. To achieve this, we have compiled a configurationally and conformationally diverse data set of zinc complexes using metadynamics to overcome the limitations of traditional sampling methods. In terms of the neural network potentials, our results indicate that for zinc complexes, partial charges play an important role in modeling the long-range interactions with a neural network. Our developed model outperforms semiempirical methods in predicting the relative energy of zinc conformers, yielding a mean absolute error (MAE) of 1.32 kcal/mol with reference to the double-hybrid PWPB95 method.

Ligand-Dependent Modulation of the Dynamics of Intracellular Loops Dictates Functional Selectivity of 5-HT2AR

The serotonin 2A receptor (5-HT_2AR) subtype of the G protein-coupled receptor (GPCR) family is involved in a plethora of neuromodulatory functions (e.g., neurogenesis, sleep, and cognitive processes). 5-HT_2AR is the target of pharmacologically distinct classes of ligands, binding of which either activate or inactivate the receptor. Although high-resolution structures of 5-HT_2AR as well as several other 5-HT GPCRs provided snapshots of both active and inactive conformational states, these structures, representing a truncated form of the receptor, cannot fully explain the mechanism of conformational transitions during their function. Importantly, biochemical studies have suggested the importance of intracellular loops in receptor functions. In our previous study, a model of the ligand-free form of 5-HT_2AR with the third intracellular loop (ICL3) has been meticulously built. Here, we have investigated the functional regulation of 5-HT_2AR with intact intracellular loops in ligand-free and five distinct ligand-bound configurations using unbiased atomistic molecular dynamics (MD) simulations. The selected ligands belong to either of the full, partial, or inverse agonist classes, which exert distinct pharmacological responses. We have observed significant structural, dynamic, and thermodynamic differences within ligand-bound complexes. Our results revealed, for the first time, that either activation or inactivation of the receptor upon specific ligand binding is primarily achieved through conformational transitions of its second and third intracellular loops (ICL2 and ICL3). A remarkable allosteric cross-talk between the ligand-binding site and the distal intracellular parts of the receptor, where binding of a specific ligand thermodynamically controls (either stabilizes or destabilizes) the intracellular region, consisting of crucial dynamic elements ICL2 and ICL3, and differential conformational transitions of these loops determine ligand-dependent functional selectivity.

Transcription factor allosteric regulation through substrate coordination to zinc

The development of new synthetic biology circuits for biotechnology and medicine requires deeper mechanistic insight into allosteric transcription factors (aTFs). Here we studied the aTF UxuR, a homodimer of two domains connected by a highly flexible linker region. To explore how ligand binding to UxuR affects protein dynamics we performed molecular dynamics simulations in the free protein, the aTF bound to the inducer D-fructuronate or the structural isomer D-glucuronate. We then validated our results by constructing a sensor plasmid for D-fructuronate in Escherichia coli and performed site-directed mutagenesis. Our results show that zinc coordination is necessary for UxuR function since mutation to alanines prevents expression de-repression by D-fructuronate. Analyzing the different complexes, we found that the disordered linker regions allow the N-terminal domains to display fast and large movements. When the inducer is bound, UxuR can sample an open conformation with a more pronounced negative charge at the surface of the N-terminal DNA binding domains. In opposition, in the free and D-glucuronate bond forms the protein samples closed conformations, with a more positive character at the surface of the DNA binding regions. These molecular insights provide a new basis to harness these systems for biological systems engineering.

Deciphering the selective binding mechanisms of anaplastic lymphoma kinase-derived neuroblastoma tumor neoepitopes to human leukocyte antigen

Neuroblastoma (NB), as a metastatic form of solid tumor, has a high fatality rate found in early childhood. The two anaplastic lymphoma kinase (ALK) neoepitopes nonamer and decamer used in cancer immunotherapy against NB cancer can selectively bind to the human leukocyte antigen (HLA-B*15:01) groove with high affinities, whereas the native self-peptide is unable to interact with the HLA-B*15:01. Here, we performed molecular dynamics (MD) simulations and subsequent molecular mechanics-generalized born surface area (MM-GBSA) binding free energy calculations to explore the selective binding mechanisms of nonamer and decamer to the HLA-B*15:01 against the self-peptide. MD simulations revealed the significant conformational dynamics of the self-peptide in the HLA-B*15:01 groove against the nonamer and decamer. Binding free energy calculations showed that the binding affinities of HLA-B*15:01-neoepitope complexes were followed in the order decamer > nonamer > self-peptide. Detailed analysis of HLA-B*15:01-neoepitope structural complexes showed that compared to the nonamer, the self-peptide tended to move outward to the solvent, whereas the decamer moved deep to the HLA-B*15:01 groove. These different dynamic observations of the ALK neoepitopes can explain the distinct binding affinities of self-peptide, nonamer, and decamer to the HLA-B*15:01. The results may be useful for the design of more selective ALK neoepitopes.

How Parkinson's disease-related mutations disrupt the dimerization of WD40 domain in LRRK2: a comparative molecular dynamics simulation study

The multidomain kinase enzyme leucine-rich-repeat kinase 2 (LRRK2), activated through a homodimerization manner, has been identified as an important pathogenic factor in Parkinson's disease (PD), the second most common neurodegenerative disease wordwide. The Trp-Asp-40 (WD40) domain, located in the C-terminal LRRK2, harbours one of the most frequent PD-related variants, G2385R. However, the detailed dynamics of WD40 during LRRK2 dimerization and the underlying mechanism through which the pathogenic mutations disrupt the formation of the WD40 dimer have remained elusive. Here, microsecond-scale molecular dynamics simulations were employed to provide a mechanistic view underlying the WD40 dimerization and unveil the structural basis by which the interface-based mutations G2385R, H2391D and R2394E compromise the corresponding process. The simulation results identified important residues, D2351, R2394, E2395, R2413, and R2443, involved in establishing the complex binding network along the dimerization interface, which was significantly weakened in the presence of interfacial mutations. A "sag-bulge" model was proposed to explain the unfavorable dimer formation in the mutant systems. In addition, mutations altered the community configuration in the wild-type system in which inter-monomeric interplay is prominent, leading to the destabilization of the WD40 dimer under mutation.