Independent Metrics for Protein Backbone and Side-Chain Flexibility: Time Scales and Effects of Ligand Binding

Atomistic simulations reveal impacts of missense mutations on the structure and function of SynGAP1

De novo mutations in the synaptic GTPase activating protein (SynGAP) are associated with neurological disorders like intellectual disability, epilepsy, and autism. SynGAP is also implicated in Alzheimer's disease and cancer. Although pathogenic variants are highly penetrant in neurodevelopmental conditions, a substantial number of them are caused by missense mutations that are difficult to diagnose. Hence, in silico mutagenesis was performed for probing the missense effects within the N-terminal region of SynGAP structure. Through extensive molecular dynamics simulations, encompassing three 150-ns replicates for 211 variants, the impact of missense mutations on the protein fold was assessed. The effect of the mutations on the folding stability was also quantitatively assessed using free energy calculations. The mutations were categorized as potentially pathogenic or benign based on their structural impacts. Finally, the study introduces wild-type-SynGAP in complex with RasGTPase at the inner membrane, while considering the potential effects of mutations on these key interactions. This study provides structural perspective to the clinical assessment of SynGAP missense variants and lays the foundation for future structure-based drug discovery.

Ion-pair interactions between voltage-sensing domain IV and pore domain I regulate CaV1.1 gating

The voltage-gated calcium channel CaV1.1 belongs to the family of pseudo-heterotetrameric cation channels, which are built of four structurally and functionally distinct voltage-sensing domains (VSDs) arranged around a common channel pore. Upon depolarization, positive gating charges in the S4 helices of each VSD are moved across the membrane electric field, thus generating the conformational change that prompts channel opening. This sliding helix mechanism is aided by the transient formation of ion-pair interactions with countercharges located in the S2 and S3 helices within the VSDs. Recently, we identified a domain-specific ion-pair partner of R1 and R2 in VSD IV of CaV1.1 that stabilizes the activated state of this VSD and regulates the voltage dependence of current activation in a splicing-dependent manner. Structure modeling of the entire CaV1.1 in a membrane environment now revealed the participation in this process of an additional putative ion-pair partner (E216) located outside VSD IV, in the pore domain of the first repeat (IS5). This interdomain interaction is specific for CaV1.1 and CaV1.2 L-type calcium channels. Moreover, in CaV1.1 it is sensitive to insertion of the 19 amino acid peptide encoded by exon 29. Whole-cell patch-clamp recordings in dysgenic myotubes reconstituted with wild-type or E216 mutants of GFP-CaV1.1e (lacking exon 29) showed that charge neutralization (E216Q) or removal of the side chain (E216A) significantly shifted the voltage dependence of activation (V1/2) to more positive potentials, suggesting that E216 stabilizes the activated state. Insertion of exon 29 in the GFP-CaV1.1a splice variant strongly reduced the ionic interactions with R1 and R2 and caused a substantial right shift of V1/2, whereas no further shift of V1/2 was observed on substitution of E216 with A or Q. Together with our previous findings, these results demonstrate that inter- and intradomain ion-pair interactions cooperate in the molecular mechanism regulating VSD function and channel gating in CaV1.1.

Systemic lupus erythematosus overlapping dermatomyositis owing to a heterozygous TREX1 Asp130Asn missense mutation

The 3' repair exonuclease 1 (TREX1) gene encodes a nuclear protein with 3' exonuclease activity, and the mutations have been associated with autoimmune diseases. Herein, we performed genetic analysis for the TREX1 gene in 55 patients with systemic lupus erythematosus (SLE). We identified one SLE patient with overlapping dermatomyositis having a heterozygous p.Asp130Asn mutation in the TREX1 gene. The patient had a high level of serum interferon (IFN)-α compared with that in healthy controls and other patients with SLE. In addition, the patient expressed elevated IFN signature genes compared with healthy controls. Our molecular dynamics simulation of the TREX1 protein in a complex with double-stranded DNA revealed that the D130N mutant causes significant changes in the active site's interaction network. One of our cases exhibited a heterozygous TREX1 p.Asp130Asn mutation that contributed to the type I IFN pathway, which may lead to the development of a severe SLE phenotype.

Decision Making in Structure-Based Drug Discovery: Visual Inspection of Docking Results

Molecular docking is a computational method widely used in drug discovery. Due to the inherent inaccuracies of molecular docking, visual inspection of binding modes is a crucial routine in the decision making process of computational medicinal chemists. Despite its apparent importance for medicinal chemistry projects, guidelines for the visual docking pose assessment have been hardly discussed in the literature. Here, we review the medicinal chemistry literature with the aim of identifying consistent principles for visual inspection, highlighting cases of its successful application, and discussing its limitations. In this context, we conducted a survey reaching experts in both academia and the pharmaceutical industry, which also included a challenge to distinguish native from incorrect poses. We were able to collect 93 expert opinions that offer valuable insights into visually supported decision-making processes. This perspective shall motivate discussions among experienced computational medicinal chemists and guide young scientists new to the field to stratify their compounds.

Antibodies exhibit multiple paratope states influencing VH-VL domain orientations

In the last decades, antibodies have emerged as one of the most important and successful classes of biopharmaceuticals. The highest variability and diversity of an antibody is concentrated on six hypervariable loops, also known as complementarity determining regions (CDRs) shaping the antigen-binding site, the paratope. Whereas it was assumed that certain sequences can only adopt a limited set of backbone conformations, in this study we present a kinetic classification of several paratope states in solution. Using molecular dynamics simulations in combination with experimental structural information we capture the involved conformational transitions between different canonical clusters and additional dominant solution structures occurring in the micro-to-millisecond timescale. Furthermore, we observe a strong correlation of CDR loop movements. Another important aspect when characterizing different paratope states is the relative VH/VL orientation and the influence of the distinct CDR loop states on the VH/VL interface. Conformational rearrangements of the CDR loops do not only have an effect on the relative VH/VL orientations, but also influence in some cases the elbow-angle dynamics and shift the respective distributions. Thus, our results show that antibodies exist as several interconverting paratope states, each contributing to the antibody's properties.

Idiopathic multicentric Castleman disease with novel heterozygous Ile729Met mutation in exon 10 of familial Mediterranean fever gene

Objective

While the aetiology of idiopathic multicentric Castleman disease (iMCD) remains unclear, the involvement of autoinflammatory mechanisms has been suggested. Herein we report a Japanese patient with iMCD with a novel heterozygous Ile729Met mutation in exon 10 of the Mediterranean fever (MEFV) gene.

Methods

We performed genetic analysis via targeted next-generation sequencing analysis and Sanger sequencing and conducted molecular dynamics simulations to investigate the hydrophobic interactions around the 729th amino acid in human pyrin.

Results

In February 2011, a 59-year-old man was diagnosed with IgG4-related disease at our department based on the findings of cervical and hilar lymphadenopathies, typical lung lesions and cervical lymph node biopsy. The patient was followed up without treatment, as he was asymptomatic. However, he had been experiencing prolonged fatigue and fever with high levels of CRP since June 2017. Axillary lymph node biopsy findings led to the diagnosis of iMCD. He was successfully treated with an IL-6 inhibitor and has been in remission for 12 months. Genetic analyses for hereditary autoinflammatory disease–related genes were performed, revealing a novel heterozygous Ile729Met mutation in exon 10 of the MEFV gene. We identified that this novel mutation significantly altered the local interaction of the human pyrin B30.2 domain by molecular dynamics simulation analysis and experimentally had the potential for inflammasome activation with increased inflammatory cytokines.

Conclusion

The abnormal function of pyrin due to a mutation in the MEFV gene in this patient may have contributed to the development of MCD by inducing IL-6 production via inflammasome signalling.

Large-Scale, Quantitative Protein Assays on a High-Throughput DNA Sequencing Chip

High-throughput DNA sequencing techniques have enabled diverse approaches for linking DNA sequence to biochemical function. In contrast, assays of protein function have substantial limitations in terms of throughput, automation, and widespread availability. We have adapted an Illumina high-throughput sequencing chip to display an immense diversity of ribosomally translated proteins and peptides and then carried out fluorescence-based functional assays directly on this flow cell, demonstrating that a single, widely available high-throughput platform can perform both sequencing-by-synthesis and protein assays. We quantified the binding of the M2 anti-FLAG antibody to a library of 1.3 × 10⁴ variant FLAG peptides, exploring non-additive effects of combinations of mutations and discovering a "superFLAG" epitope variant. We also measured the enzymatic activity of 1.56 × 10⁵ molecular variants of full-length human O⁶-alkylguanine-DNA alkyltransferase (SNAP-tag). This comprehensive corpus of catalytic rates revealed amino acid interaction networks and cooperativity, linked positive cooperativity to structural proximity, and revealed ubiquitous positively cooperative interactions with histidine residues.

An unexpected switch in peptide binding mode: from simulation to substrate specificity

A ten microsecond molecular dynamics simulation of a kallikrein-related peptidase 7 peptide complex revealed an unexpected change in binding mode. After more than two microseconds unrestrained sampling we observe a spontaneous transition of the binding pose including a 180° rotation around the P1 residue. Subsequently, the substrate peptide occupies the prime side region rather than the cognate non-prime side in a stable conformation. We characterize the unexpected binding mode in terms of contacts, solvent-accessible surface area, molecular interactions and energetic properties. We compare the new pose to inhibitor-bound structures of kallikreins with occupied prime side and find that a similar orientation is adopted. Finally, we apply in silico mutagenesis based on the alternative peptide binding position to explore the prime side specificity of kallikrein-related peptidase 7 and compare it to available experimental data. Our study provides the first microsecond time scale simulation data on a kallikrein protease and shows previously unexplored prime side interactions. Therefore, we expect our study to advance the rational design of inhibitors targeting kallikrein-related peptidase 7, an emerging drug target involved in several skin diseases as well as cancer.

Effect of PDGF-B aptamer on PDGFRβ/PDGF-B interaction: Molecular dynamics study

PDGFRβ/PDGF-B interaction plays a role in angiogenesis, and is mandatory in wound healing and cancer treatment. It has been reported that the PDGF-B aptamer was able to bind to PDGF-B, thus regulating the angiogenesis. However, the binding interaction between the aptamer and the growth factor, including the binding sites, has not been well investigated. This study applied a molecular dynamics (MD) simulation to investigate the aptamer-growth factor interaction in the presence or absence of a receptor (PDGFRβ). Characterization of the structure of an aptamer-growth factor complex revealed binding sites from each section in the complex. Upon the complex formation, PDGF-B and its aptamer exhibited less flexibility in their molecular movement, as indicated by the minimum values of RMSD, RMSF, loop-to-loop distance, and the summation of PCA eigenvalues. Our study of residue pairwise interaction demonstrated that the binding interaction was mainly contributed by electrostatic interaction between the positively-charged amino acid and the negatively-charged phosphate backbone. The role of the PDGF-B aptamer in PDGFRβ/PDGF-B interaction was also investigated. We demonstrated that the stability of the Apt-PDGF-B complex could prevent the presence of a competitor, of PDGFRβ, interrupting the binding process. Because the aptamer was capable of binding with PDGF-B, and blocking the growth factor from the PDGFRβ, it could down regulate the consequent signaling pathway. We provide evidence that the PDGF-BB aptamer is a promising molecule for regulation of angiogenesis. The MD study provides a molecular understanding to modification of the aptamer binding interaction, which could be used in a number of medical applications.

Structure-Based Target-Specific Screening Leads to Small-Molecule CaMKII Inhibitors

Target-specific scoring methods are more commonly used to identify small-molecule inhibitors among compounds docked to a target of interest. Top candidates that emerge from these methods have rarely been tested for activity and specificity across a family of proteins. In this study we docked a chemical library into CaMKIIδ, a member of the Ca2+ /calmodulin (CaM)-dependent protein kinase (CaMK) family, and re-scored the resulting protein-compound structures using Support Vector Machine SPecific (SVMSP), a target-specific method that we developed previously. Among the 35 selected candidates, three hits were identified, such as quinazoline compound 1 (KIN-1; N4-[7-chloro-2-[(E)-styryl]quinazolin-4-yl]-N1,N1-diethylpentane-1,4-diamine), which was found to inhibit CaMKIIδ kinase activity at single-digit micromolar IC50 . Activity across the kinome was assessed by profiling analogues of 1, namely 6 (KIN-236; N4-[7-chloro-2-[(E)-2-(2-chloro-4,5-dimethoxyphenyl)vinyl]quinazolin-4-yl]-N1,N1-diethylpentane-1,4-diamine), and an analogue of hit compound 2 (KIN-15; 2-[4-[(E)-[(5-bromobenzofuran-2-carbonyl)hydrazono]methyl]-2-chloro-6-methoxyphenoxy]acetic acid), namely 14 (KIN-332; N-[(E)-[4-(2-anilino-2-oxoethoxy)-3-chlorophenyl]methyleneamino]benzofuran-2-carboxamide), against 337 kinases. Interestingly, for compound 6, CaMKIIδ and homologue CaMKIIγ were among the top ten targets. Among the top 25 targets of 6, IC50 values ranged from 5 to 22 μm. Compound 14 was found to be not specific toward CaMKII kinases, but it does inhibit two kinases with sub-micromolar IC50 values among the top 25. Derivatives of 1 were tested against several kinases including several members of the CaMK family. These data afforded a limited structure-activity relationship study. Molecular dynamics simulations with explicit solvent followed by end-point MM-GBSA free-energy calculations revealed strong engagement of specific residues within the ATP binding pocket, and also changes in the dynamics as a result of binding. This work suggests that target-specific scoring approaches such as SVMSP may hold promise for the identification of small-molecule kinase inhibitors that exhibit some level of specificity toward the target of interest across a large number of proteins.

Determinants of Macromolecular Specificity from Proteomics-Derived Peptide Substrate Data

BACKGROUND

Recent advances in proteomics methodologies allow for high throughput profiling of proteolytic cleavage events. The resulting substrate peptide distributions provide deep insights in the underlying macromolecular recognition events, as determinants of biomolecular specificity identified by proteomics approaches may be compared to structure-based analysis of corresponding protein-protein interfaces.

METHOD

Here, we present an overview of experimental and computational methodologies and tools applied in the area and provide an outlook beyond the protein class of proteases.

RESULTS AND CONCLUSION

We discuss here future potential, synergies and needs of the emerging overlap disciplines of proteomics and structure-based modelling.

Analysis of Differential Efficacy and Affinity of GABAA (α1/α2) Selective Modulators

Selective modulators of the γ-amino butyric acid (GABA_A) family of receptors have the potential to treat a range of disease states related to cognition, pain, and anxiety. While the development of various α subunit-selective modulators is currently underway for the treatment of anxiety disorders, a mechanistic understanding of the correlation between their bioactivity and efficacy, based on ligand-target interactions, is currently still lacking. In order to alleviate this situation, in the current study we have analyzed, using ligand- and structure-based methods, a data set of 5440 GABA_A modulators. The Spearman correlation (ρ) between binding activity and efficacy of compounds was calculated to be 0.008 and 0.31 against the α1 and α2 subunits of GABA receptor, respectively; in other words, the compounds had little diversity in structure and bioactivity, but they differed significantly in efficacy. Two compounds were selected as a case study for detailed interaction analysis due to the small difference in their structures and affinities (ΔpK_{i(comp1_α1 - comp2_α1)} = 0.45 log units, ΔpK_{i(comp1_α2 - comp2_α2)} = 0 log units) as compared to larger relative efficacies (ΔRE_{(comp1_α1 - comp2_α1)} = 1.03, ΔRE_{(comp1_α2 - comp2_α2)} = 0.21). Docking analysis suggested that His-101 is involved in a characteristic interaction of the α1 receptor with both compounds 1 and 2. Residues such as Phe-77, Thr-142, Asn-60, and Arg-144 of the γ chain of the α1γ2 complex also showed interactions with heterocyclic rings of both compounds 1 and 2, but these interactions were disturbed in the case of α2γ2 complex docking results. Binding pocket stability analysis based on molecular dynamics identified three substitutions in the loop C region of the α2 subunit, namely, G200E, I201T, and V202I, causing a reduction in the flexibility of α2 compared to α1. These amino acids in α2, as compared to α1, were also observed to decrease the vibrational and dihedral entropy and to increase the hydrogen bond content in α2 in the apo state. However, freezing of both α1 and α2 was observed in the ligand-bound state, with an increased number of internal hydrogen bonds and increased entropy. Therefore, we hypothesize that the amino acid differences in the loop C region of α2 are responsible for conformational changes in the protein structure compared to α1, as well as for the binding modes of compounds and hence their functional signaling.

Oncogenic Mutations Differentially Affect Bax Monomer, Dimer, and Oligomeric Pore Formation in the Membrane

Dysfunction of Bax, a pro-apoptotic regulator of cellular metabolism is implicated in neurodegenerative diseases and cancer. We have constructed the first atomistic models of the Bax oligomeric pore consisting with experimental residue-residue distances. The models are stable, capturing well double electron-electron resonance (DEER) spectroscopy measurements and provide structural details in line with the DEER data. Comparison with the latest experimental results revealed that our models agree well with both Bax and Bak pores, pointed to a converged structural arrangement for Bax and Bak pore formation. Using multi-scale molecular dynamics simulations, we probed mutational effects on Bax transformation from monomer → dimer → membrane pore formation at atomic resolution. We observe that two cancer-related mutations, G40E and S118I, allosterically destabilize the monomer and stabilize an off-pathway swapped dimer, preventing productive pore formation. This observation suggests a mechanism whereby the mutations may work mainly by over-stabilizing the monomer → dimer transformation toward an unproductive off-pathway swapped-dimer state. Our observations point to misfolded Bax states, shedding light on the molecular mechanism of Bax mutation-elicited cancer. Most importantly, the structure of the Bax pore facilitates future study of releases cytochrome C in atomic detail.

Localization of Millisecond Dynamics: Dihedral Entropy from Accelerated MD

Here, we demonstrate a method to capture local dynamics on a time scale 3 orders of magnitude beyond state-of-the-art simulation approaches. We apply accelerated molecular dynamics simulations for conformational sampling and extract reweighted backbone dihedral distributions. Local dynamics are characterized by torsional probabilities, resulting in residue-wise dihedral entropies. Our approach is successfully validated for three different protein systems of increasing size: alanine dipeptide, bovine pancreatic trypsin inhibitor (BPTI), and the major birch pollen allergen Bet v 1a. We demonstrate excellent agreement of flexibility profiles with both large-scale computer simulations and NMR experiments. Thus, our method provides efficient access to local protein dynamics on extended time scales of high biological relevance.

Effect of BET Missense Mutations on Bromodomain Function, Inhibitor Binding and Stability

Lysine acetylation is an important epigenetic mark regulating gene transcription and chromatin structure. Acetylated lysine residues are specifically recognized by bromodomains, small protein interaction modules that read these modification in a sequence and acetylation dependent way regulating the recruitment of transcriptional regulators and chromatin remodelling enzymes to acetylated sites in chromatin. Recent studies revealed that bromodomains are highly druggable protein interaction domains resulting in the development of a large number of bromodomain inhibitors. BET bromodomain inhibitors received a lot of attention in the oncology field resulting in the rapid translation of early BET bromodomain inhibitors into clinical studies. Here we investigated the effects of mutations present as polymorphism or found in cancer on BET bromodomain function and stability and the influence of these mutants on inhibitor binding. We found that most BET missense mutations localize to peripheral residues in the two terminal helices. Crystal structures showed that the three dimensional structure is not compromised by these mutations but mutations located in close proximity to the acetyl-lysine binding site modulate acetyl-lysine and inhibitor binding. Most mutations affect significantly protein stability and tertiary structure in solution, suggesting new interactions and an alternative network of protein-protein interconnection as a consequence of single amino acid substitution. To our knowledge this is the first report studying the effect of mutations on bromodomain function and inhibitor binding.

Kinetic barriers in the isomerization of substituted ureas: implications for computer-aided drug design

Urea derivatives are ubiquitously found in many chemical disciplines. N,N'-substituted ureas may show different conformational preferences depending on their substitution pattern. The high energetic barrier for isomerization of the cis and trans state poses additional challenges on computational simulation techniques aiming at a reproduction of the biological properties of urea derivatives. Herein, we investigate energetics of urea conformations and their interconversion using a broad spectrum of methodologies ranging from data mining, via quantum chemistry to molecular dynamics simulation and free energy calculations. We find that the inversion of urea conformations is inherently slow and beyond the time scale of typical simulation protocols. Therefore, extra care needs to be taken by computational chemists to work with appropriate model systems. We find that both knowledge-driven approaches as well as physics-based methods may guide molecular modelers towards accurate starting structures for expensive calculations to ensure that conformations of urea derivatives are modeled as adequately as possible.

Sequence diversity of NanA manifests in distinct enzyme kinetics and inhibitor susceptibility

Streptococcus pneumoniae is the leading pathogen causing bacterial pneumonia and meningitis. Its surface-associated virulence factor neuraminidase A (NanA) promotes the bacterial colonization by removing the terminal sialyl residues from glycoconjugates on eukaryotic cell surface. The predominant role of NanA in the pathogenesis of pneumococci renders it an attractive target for therapeutic intervention. Despite the highly conserved activity of NanA, our alignment of the 11 NanAs revealed the evolutionary diversity of this enzyme. The amino acid substitutions we identified, particularly those in the lectin domain and in the insertion domain next to the catalytic centre triggered our special interest. We synthesised the representative NanAs and the mutagenized derivatives from E. coli for enzyme kinetics study and neuraminidase inhibitor susceptibility test. Via molecular docking we got a deeper insight into the differences between the two major variants of NanA and their influence on the ligand-target interactions. In addition, our molecular dynamics simulations revealed a prominent intrinsic flexibility of the linker between the active site and the insertion domain, which influences the inhibitor binding. Our findings for the first time associated the primary sequence diversity of NanA with the biochemical properties of the enzyme and with the inhibitory efficiency of neuraminidase inhibitors.

Dynamics of an Active-Site Flap Contributes to Catalysis in a JAMM Family Metallo Deubiquitinase

The endosome-associated deubiquitinase (DUB) AMSH is a member of the JAMM family of zinc-dependent metallo isopeptidases with high selectivity for Lys63-linked polyubiquitin chains, which play a key role in endosomal-lysosomal sorting of activated cell surface receptors. The catalytic domain of the enzyme features a flexible flap near the active site that opens and closes during its catalytic cycle. Structural analysis of its homologues, AMSH-LP (AMSH-like protein) and the fission yeast counterpart, Sst2, suggests that a conserved Phe residue in the flap may be critical for substrate binding and/or catalysis. To gain insight into the contribution of this flap in substrate recognition and catalysis, we generated mutants of Sst2 and characterized them using a combination of enzyme kinetics, X-ray crystallography, molecular dynamics simulations, and isothermal titration calorimetry (ITC). Our analysis shows that the Phe residue in the flap contributes key interactions during the rate-limiting step but not to substrate binding, since mutants of Phe403 exhibit a defect only in kcat but not in KM. Moreover, ITC studies show Phe403 mutants have similar KD for ubiquitin compared to the wild-type enzyme. The X-ray structures of both Phe403Ala and the Phe403Trp, in both the free and ubiquitin bound form, reveal no appreciable structural change that might impair substrate or alter product binding. We observed that the side chain of the Trp residue is oriented identically with respect to the isopeptide moiety of the substrate as the Phe residue in the wild-type enzyme, so the loss of activity seen in this mutant cannot be explained by the absence of a group with the ability to provide van der Waals interactions that facilitate the hyrdolysis of the Lys63-linked diubiquitin. Molecular dynamics simulations indicate that the flap in the Trp mutant is quite flexible, allowing almost free rotation of the indole side chain. Therefore, it is possible that these different dynamic properties of the flap in the Trp mutant, compared to the wild-type enzyme, manifest as a defect in interactions that facilitate the rate-limiting step. Consistent with this notion, the Trp mutant was able to cleave Lys48-linked and Lys11-linked diubiquitin better than the wild-type enzyme, indicating altered mobility and hence reduced selectivity.

Experimental and computational evidence on conformational fluctuations as a source of catalytic defects in genetic diseases

Catalytic mutants causing inherited type III galactosemia alter active site structural dynamics and shift the native conformational equilibrium towards inactive conformations.