Search

Malaria pathology is driven by the accumulation of Plasmodium falciparum-infected erythrocytes in microvessels 1 . This process is mediated by the polymorphic erythrocyte membrane protein 1 (PfEMP1) adhesion proteins of the parasite. A subset of PfEMP1 variants that bind to human endothelial protein C receptor (EPCR) through their CIDRα1 domains is responsible for severe malaria pathogenesis 2 . A longstanding question is whether individual antibodies can recognize the large repertoire of circulating PfEMP1 variants. Here we describe two broadly reactive and inhibitory human monoclonal antibodies to CIDRα1. The antibodies isolated from two different individuals exhibited similar and consistent EPCR-binding inhibition of diverse CIDRα1 domains, representing five of the six subclasses of CIDRα1. Both antibodies inhibited EPCR binding of both recombinant full-length and native PfEMP1 proteins, as well as parasite sequestration in bioengineered 3D human brain microvessels under physiologically relevant flow conditions. Structural analyses of the two antibodies in complex with three different CIDRα1 antigen variants reveal similar binding mechanisms that depend on interactions with three highly conserved amino acid residues of the EPCR-binding site in CIDRα1. These broadly reactive antibodies are likely to represent a common mechanism of acquired immunity to severe malaria and offer novel insights for the design of a vaccine or treatment targeting severe malaria., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Recent breakthroughs in protein structure prediction have enhanced the precision and speed at which protein configurations can be determined. Additionally, molecular dynamics (MD) simulations serve as a crucial tool for capturing the conformational space of proteins, providing valuable insights into their structural fluctuations. However, the scope of MD simulations is often limited by the accessible timescales and the computational resources available, posing challenges to comprehensively exploring protein behaviors. Recently emerging approaches have focused on expanding the capability of AlphaFold2 (AF2) to predict conformational substates of protein. Here, we benchmark the performance of various workflows that have adapted AF2 for ensemble prediction and compare the obtained structures with ensembles obtained from MD simulations and NMR. We provide an overview of the levels of performance and accessible timescales that can currently be achieved with machine learning (ML) based ensemble generation. Significant minima of the free energy surfaces remain undetected., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024. Published by Elsevier Inc.)

The humanization of camelid-derived variable domain heavy chain antibodies (VHHs) poses challenges including immunogenicity, stability, and potential reduction of affinity. Critical to this process are complementarity-determining regions (CDRs), Vernier and Hallmark residues, shaping the three-dimensional fold and influencing VHH structure and function. Additionally, the presence of non-canonical disulfide bonds further contributes to conformational stability and antigen binding. In this study, we systematically humanized two camelid-derived VHHs targeting the natural cytotoxicity receptor NKp30. Key structural positions in Vernier and Hallmark regions were exchanged with residues from the most similar human germline sequences. The resulting variants were characterized for binding affinities, yield, and purity. Structural binding modes were elucidated through crystal structure determination and AlphaFold2 predictions, providing insights into differences in binding affinity. Comparative structural and molecular dynamics characterizations of selected variants were performed to rationalize their functional properties and elucidate the role of specific sequence motifs in antigen binding. Furthermore, systematic analyses of next-generation sequencing (NGS) and Protein Data Bank (PDB) data was conducted, shedding light on the functional significance of Hallmark motifs and non-canonical disulfide bonds in VHHs in general. Overall, this study provides valuable insights into the structural determinants governing the functional properties of VHHs, offering a roadmap for their rational design, humanization, and optimization for therapeutic applications., (© 2024 The Author(s). Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

Background and Purpose: Pathogenic gain-of-function mutations in Ca v 1.3 L-type voltage-gated Ca 2+ -channels (CACNA1D) cause neurodevelopmental disorders with or without endocrine symptoms. We aimed to confirm a pathogenic gain-of function phenotype of CACNA1D de novo missense mutations A749T and L271H, and investigated the molecular mechanism causing their enhanced sensitivity for the Ca 2+ -channel blocker isradipine, a potential therapeutic for affected patients., Experimental Approach: Wildtype and mutant channels were expressed in tsA-201 cells and their gating analysed using whole-cell and single-channel patch-clamp recordings. The voltage-dependence of isradipine action was quantified using protocols inducing variable fractions of inactivated channels. The molecular basis for altered channel gating in the mutants was investigated using in silico modelling and molecular dynamics simulations., Key Results: Both mutations were confirmed pathogenic due to characteristic shifts of voltage-dependent activation and inactivation towards negative potentials (~20 mV). At negative holding potentials both mutations showed significantly higher isradipine sensitivity compared to wildtype. The affinity for wildtype and mutant channels increased with channel inactivation as predicted by the modulated receptor hypothesis (30- to 40-fold). The IC 50 was indistinguishable for wildtype and mutants when >50% of channels were inactivated., Conclusions and Implications: Mutations A749T and L271H induce pathogenic gating changes. Like wildtype, isradipine inhibition is strongly voltage-dependent. Our data explains their apparent higher drug sensitivity at a given negative voltage by the availability of more inactivated channels due to their more negative inactivation voltage range. Low nanomolar isradipine concentrations will only inhibit Ca v 1.3 channels in neurons during prolonged depolarized states without selectivity for mutant channels., (© 2024 The Author(s). British Journal of Pharmacology published by John Wiley & Sons Ltd on behalf of British Pharmacological Society.)

Immunoglobulin G (IgG) antibodies that bind their cognate antigen in a pH-dependent manner (acid-switched antibodies) can release their bound antigen for degradation in the acidic environment of endosomes, while the IgGs are rescued by the neonatal Fc receptor (FcRn). Thus, such IgGs can neutralize multiple antigens over time and therefore be used at lower doses than their non-pH-responsive counterparts. Here, we show that light-chain shuffling combined with phage display technology can be used to discover IgG1 antibodies with increased pH-dependent antigen binding properties, using the snake venom toxins, myotoxin II and α-cobratoxin, as examples. We reveal differences in how the selected IgG1s engage their antigens and human FcRn and show how these differences translate into distinct cellular handling properties related to their pH-dependent antigen binding phenotypes and Fc-engineering for improved FcRn binding. Our study showcases the complexity of engineering pH-dependent antigen binding IgG1s and demonstrates the effects on cellular antibody-antigen recycling., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Voltage-gated L-type Cav1.3 Ca 2+ channels support numerous physiological functions including neuronal excitability, sinoatrial node pacemaking, hearing, and hormone secretion. De novo missense mutations in the gene of their pore-forming α1-subunit (CACNA1D) induce severe gating defects which lead to autism spectrum disorder and a more severe neurological disorder with and without endocrine symptoms. The number of CACNA1D variants reported is constantly rising, but their pathogenic potential often remains unclear, which complicates clinical decision-making. Since functional tests are time-consuming and not always available, bioinformatic tools further improving pathogenicity potential prediction of novel variants are needed. Here we employed evolutionary analysis considering sequences of the Cav1.3 α1-subunit throughout the animal kingdom to predict the pathogenicity of human disease-associated CACNA1D missense variants. Co-variation analyses of evolutionary information revealed residue-residue couplings and allowed to generate a score, which correctly predicted previously identified pathogenic variants, supported pathogenicity in variants previously classified as likely pathogenic and even led to the re-classification or re-examination of 18 out of 80 variants previously assessed with clinical and electrophysiological data. Based on the prediction score, we electrophysiologically tested one variant (V584I) and found significant gating changes associated with pathogenic risks. Thus, our co-variation model represents a valuable addition to complement the assessment of the pathogenicity of CACNA1D variants completely independent of clinical diagnoses, electrophysiology, structural or biophysical considerations, and solely based on evolutionary analyses., (© 2024. The Author(s).)

Skeletal muscle contractions are initiated by action potentials, which are sensed by the voltage-gated calcium channel (Ca V 1.1) and are conformationally coupled to calcium release from intracellular stores. Notably, Ca V 1.1 contains four separate voltage-sensing domains (VSDs), which activate channel gating and excitation-contraction (EC-) coupling at different voltages and with distinct kinetics. Here we show that a single VSD of Ca V 1.1 controls skeletal muscle EC-coupling. Whereas mutations in VSDs I, II and IV affect the current properties but not EC-coupling, only mutations in VSD III alter the voltage-dependence of depolarization-induced calcium release. Molecular dynamics simulations reveal comprehensive, non-canonical state transitions of VSD III in response to membrane depolarization. Identifying the voltage sensor that activates EC-coupling and detecting its unique conformational changes opens the door to unraveling the downstream events linking VSD III motion to the opening of the calcium release channel, and thus resolving the signal transduction mechanism of skeletal muscle EC-coupling., (© 2024. The Author(s).)

Single-domain antibodies (sdAbs), such as V H Hs, are increasingly being developed for gastrointestinal (GI) applications against pathogens to strengthen gut health. However, what constitutes a suitable developability profile for applying these proteins in a gastrointestinal setting remains poorly explored. Here, we describe an in vitro methodology for the identification of sdAb derivatives, more specifically divalent V H H constructs, that display extraordinary developability properties for oral delivery and functionality in the GI environment. We showcase this by developing a heterodivalent V H H construct that cross-inhibits the toxic activity of the glycosyltransferase domains (GTDs) from three different toxinotypes of cytotoxin B (TcdB) from lineages of Clostridium difficile. We show that the V H H construct possesses high stability and binding activity under gastric conditions, in the presence of bile salts, and at high temperatures. We suggest that the incorporation of early developability assessment could significantly aid in the efficient discovery of V H Hs and related constructs fit for oral delivery and GI applications., (© 2024 The Author(s). Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

Protein folding is a fascinating, not fully understood phenomenon in biology. Molecular dynamics (MD) simulations are an invaluable tool to study conformational changes in atomistic detail, including folding and unfolding processes of proteins. However, the accuracy of the conformational ensembles derived from MD simulations inevitably relies on the quality of the underlying force field in combination with the respective water model. Here, we investigate protein folding, unfolding, and misfolding of fast-folding proteins by examining different force fields with their recommended water models, i.e., ff14SB with the TIP3P model and ff19SB with the OPC model. To this end, we generated long conventional MD simulations highlighting the perks and pitfalls of these setups. Using Markov state models, we defined kinetically independent conformational substates and emphasized their distinct characteristics, as well as their corresponding state probabilities. Surprisingly, we found substantial differences in thermodynamics and kinetics of protein folding, depending on the combination of the protein force field and water model, originating primarily from the different water models. These results emphasize the importance of carefully choosing the force field and the respective water model as they determine the accuracy of the observed dynamics of folding events. Thus, the findings support the hypothesis that the water model is at least equally important as the force field and hence needs to be considered in future studies investigating protein dynamics and folding in all areas of biophysics.

In peanut allergy, Arachis hypogaea 2 (Ara h 2) and Arachis hypogaea 6 (Ara h 6) are two clinically relevant peanut allergens with known structural and sequence homology and demonstrated cross-reactivity. We have previously utilized X-ray crystallography and epitope binning to define the epitopes on Ara h 2. We aimed to quantitatively characterize the cross-reactivity between Ara h 2 and Ara h 6 on a molecular level using human monoclonal antibodies (mAbs) and structural characterization of allergenic epitopes. We utilized mAbs cloned from Ara h 2 positive single B cells isolated from peanut-allergic, oral immunotherapy-treated patients to quantitatively analyze cross-reactivity between recombinant Ara h 2 (rAra h 2) and Ara h 6 (rAra h 6) proteins using biolayer interferometry and indirect inhibitory ELISA. Molecular dynamics simulations assessed time-dependent motions and interactions in the antibody-antigen complexes. Three epitopes-conformational epitopes 1.1 and 3, and the sequential epitope KRELRNL/KRELMNL-are conserved between Ara h 2 and Ara h 6, while two more conformational and three sequential epitopes are not. Overall, mAb affinity was significantly lower to rAra h 6 than it was to rAra h 2. This difference in affinity was primarily due to increased dissociation of the antibodies from rAra h 6, a phenomenon explained by the higher conformational flexibility of the Ara h 6-antibody complexes in comparison to Ara h 2-antibody complexes. Our results further elucidate the cross-reactivity of peanut 2S albumins on a molecular level and support the clinical immunodominance of Ara h 2., (© The Author(s) 2024. Published by Oxford University Press on behalf of the British Society for Immunology. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com. Elements of the work have been written by employees of the US Government.)

Recycling IgG antibodies bind to their target antigen at physiological pH in the blood stream and release them upon endocytosis when pH levels drop, allowing the IgG antibodies to be recycled into circulation via FcRn-mediated cellular pathways, while the antigens undergo lysosomal degradation. This enables recycling antibodies to achieve comparable therapeutic effect at lower doses than their non-recycling counterparts. The development of such antibodies is typically achieved by histidine doping of their variable regions or by performing in vitro antibody selection campaigns utilizing histidine doped libraries. Both are strategies that may introduce sequence liabilities. Here, we present a methodology that employs a naïve antibody phage display library, consisting of natural variable domains, to discover antibodies that bind α-cobratoxin from the venom of Naja kaouthia in a pH-dependent manner. As a result, an antibody was discovered that exhibits a 7-fold higher off-rate at pH 5.5 than pH 7.4 in bio-layer interferometry experiments. Interestingly, no histidine residues were found in its variable domains, and in addition, the antibody showed pH-dependent binding to a histidine-devoid antigen mutant. As such, the results demonstrate that pH-dependent antigen-antibody binding may not always be driven by histidine residues. By employing molecular dynamics simulations, different protonation states of titratable residues were found, which potentially could be responsible for the observed pH-dependent antigen binding properties of the antibody. Finally, given the typically high diversity of naïve antibody libraries, the methodology presented here can likely be applied to discover recycling antibodies against different targets ab initio without the need for histidine doping., (© 2023 The Authors. Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

Cav1.4 L-type calcium channels are predominantly expressed at the photoreceptor terminals and in bipolar cells, mediating neurotransmitter release. Mutations in its gene, CACNA1F , can cause congenital stationary night-blindness type 2 (CSNB2). Due to phenotypic variability in CSNB2, characterization of pathological variants is necessary to better determine pathological mechanism at the site of action. A set of known mutations affects conserved gating charges in the S4 voltage sensor, two of which have been found in male CSNB2 patients. Here, we describe two disease-causing Cav1.4 mutations with gating charge neutralization, exchanging an arginine 964 with glycine (RG) or arginine 1288 with leucine (RL). In both, charge neutralization was associated with a reduction channel expression also reflected in smaller ON gating currents. In RL channels, the strong decrease in whole-cell current densities might additionally be explained by a reduction of single-channel currents. We further identified alterations in their biophysical properties, such as a hyperpolarizing shift of the activation threshold and an increase in slope factor of activation and inactivation. Molecular dynamic simulations in RL substituted channels indicated water wires in both, resting and active, channel states, suggesting the development of omega ( ω )currents as a new pathological mechanism in CSNB2. This sum of the respective channel property alterations might add to the differential symptoms in patients beside other factors, such as genomic and environmental deviations.

We evaluated a workflow to reliably sample the conformational space of a set of 47 peptidic macrocycles. Starting from SMILES strings, we use accelerated molecular dynamics simulations to overcome high energy barriers, in particular, the cis-trans isomerization of peptide bonds. We find that our approach performs very well in polar solvents like water and dimethyl sulfoxide. Interestingly, the protonation state of a secondary amine in the ring only slightly influences the conformational ensembles of our test systems. For several of the macrocycles, determining the conformational distribution in chloroform turns out to be considerably more challenging. Especially, the choice of partial charges crucially influences the ensembles in chloroform. We address these challenges by modifying initial structures and the choice of partial charges. Our results suggest that special care has to be taken to understand the configurational distribution in apolar solvents, which is a key step toward a reliable prediction of membrane permeation of macrocycles and their chameleonic properties.

The electrostatic properties of proteins arise from the number and distribution of polar and charged residues. Electrostatic interactions in proteins play a critical role in numerous processes such as molecular recognition, protein solubility, viscosity, and antibody developability. Thus, characterizing and quantifying electrostatic properties of a protein are prerequisites for understanding these processes. Here, we present PEP-Patch, a tool to visualize and quantify the electrostatic potential on the protein surface in terms of surface patches, denoting separated areas of the surface with a common physical property. We highlight its applicability to elucidate protease substrate specificity and antibody-antigen recognition and predict heparin column retention times of antibodies as an indicator of pharmacokinetics.

Allergies related to kiwi consumption have become a growing health concern, with their prevalence on the rise. Many of these allergic reactions are attributed to cross-reactivity, particularly with the major allergen found in birch pollen. This cross-reactivity is associated with proteins belonging to the pathogenesis-related class 10 (PR-10) protein family. In our study, we determined the three-dimensional structures of the two PR-10 proteins in gold and green kiwi fruits, Act c 8 and Act d 8, using nuclear magnetic resonance (NMR) spectroscopy. The structures of both kiwi proteins closely resemble the major birch pollen allergen, Bet v 1, providing a molecular explanation for the observed immunological cross-reactivity between kiwi and birch pollen. Compared to Act d 11, however, a kiwi allergen that shares the same architecture as PR-10 proteins, structural differences are apparent. Moreover, despite both Act c 8 and Act d 8 containing multiple cysteine residues, no disulfide bridges are present within their structures. Instead, all the cysteines are accessible on the protein's surface and exposed to the surrounding solvent, where they are available for reactions with components of the natural food matrix. This structural characteristic sets Act c 8 and Act d 8 apart from other kiwi proteins with a high cysteine content. Furthermore, we demonstrate that pyrogallol, the most abundant phenolic compound found in kiwi, binds into the internal cavities of these two proteins, albeit with low affinity. Our research offers a foundation for further studies aimed at understanding allergic reactions associated with this fruit and exploring how interactions with the natural food matrix might be employed to enhance food safety.

Antibodies and other new antibody-like formats have emerged as one of the most rapidly growing classes of biotherapeutic proteins. Understanding the structural features that drive antibody function and, consequently, their molecular recognition is critical for engineering antibodies. Here, we present the structural architecture of conventional IgG antibodies alongside other formats. We emphasize the importance of considering antibodies as conformational ensembles in solution instead of focusing on single-static structures because their functions and properties are strongly governed by their dynamic nature. Thus, in this review, we provide an overview of the unique structural and dynamic characteristics of antibodies with respect to their antigen recognition, biophysical properties, and effector functions. We highlight the numerous technical advances in antibody structure prediction and design, enabled by the vast number of experimentally determined high-quality structures recorded with cryo-EM, NMR, and X-ray crystallography. Lastly, we assess antibody and vaccine design strategies in the context of structure and dynamics.

In the last years, antibodies have emerged as a promising new class of therapeutics, due to their combination of high specificity with long serum half-life and low risk of side-effects. Diabodies are a popular novel antibody format, consisting of two F v domains connected with short linkers. Like IgG antibodies, they simultaneously bind two target proteins. However, they offer altered properties, given their smaller size and higher rigidity. In this study, we conducted the-to our knowledge-first molecular dynamics (MD) simulations of diabodies and find a surprisingly high conformational flexibility in the relative orientation of the two F v domains. We observe rigidifying effects through the introduction of disulfide bonds in the F v -F v interface and characterize the effect of different disulfide bond locations on the conformation. Additionally, we compare V H -V L orientations and paratope dynamics between diabodies and an antigen binding fragment (Fab) of the same sequence. We find mostly consistent structures and dynamics, indicating similar antigen binding properties. The most significant differences can be found within the CDR-H2 loop dynamics. Of all CDR loops, the CDR-H2 is located closest to the artificial F v -F v interface. All examined diabodies show similar V H -V L orientations, F v -F v packing and CDR loop conformations. However, the variant with a P14C-K64C disulfide bond differs most from the Fab in our measures, including the CDR-H3 loop conformational ensemble. This suggests altered antigen binding properties and underlines the need for careful validation of the disulfide bond locations in diabodies., (© 2023 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

BCRs (Abs) and TCRs (or adaptive immune receptors [AIRs]) are the means by which the adaptive immune system recognizes foreign and self-antigens, playing an integral part in host defense, as well as the emergence of autoimmunity. Importantly, the interaction between AIRs and their cognate Ags defies a simple key-in-lock paradigm and is instead a complex many-to-many mapping between an individual's massively diverse AIR repertoire, and a similarly diverse antigenic space. Understanding how adaptive immunity balances specificity with epitopic coverage is a key challenge for the field, and terms such as broad specificity, cross-reactivity, and polyreactivity remain ill-defined and are used inconsistently. In this Immunology Notes and Resources article, a group of experimental, structural, and computational immunologists define commonly used terms associated with AIR binding, describe methodologies to study these binding modes, as well as highlight the implications of these different binding modes for therapeutic design., (Copyright © 2023 by The American Association of Immunologists, Inc.)

The use of volatile biomarkers in exhaled breath as predictors to individual drug response would advance the field of personalised medicine by providing direct information on enzyme activity. This would result in enormous benefits, both for patients and for the healthcare sector. Non-invasive breath tests would also gain a high acceptance by patients. Towards this goal, differences in metabolism resulting from extensive polymorphisms in a major group of drug-metabolizing enzymes, the cytochrome P450 (CYP) family, need to be determined and quantified. CYP2C9 is responsible for metabolising many crucial drugs (e.g., diclofenac) and food ingredients (e.g., limonene). In this paper, we provide a proof-of-concept study that illustrates the in vitro bioconversion of diclofenac in recombinant HEK293T cells overexpressing CYP2C9 to 4'-hydroxydiclofenac. This in vitro approach is a necessary and important first step in the development of breath tests to determine and monitor metabolic processes in the human body. By focusing on the metabolic conversion of diclofenac, we have been able to establish a workflow using a cell-based system for CYP2C9 activity. Furthermore, we illustrate how the bioconversion of diclofenac is limited in the presence of limonene, which is another CYP2C9 metabolising substrate. We show that increasing limonene levels continuously reduce the production of 4'-hydroxydiclofenac. Michaelis-Menten kinetics were performed for the diclofenac 4'-hydroxylation with and without limonene, giving a kinetic constant of the reaction, K M , of 103 µ M and 94.1 µ M, respectively, and a maximum reaction rate, V max , of 46.8 pmol min -1 10 6 cells -1 and 56.0 pmol min -1 10 6 cells -1 with and without the inhibitor, respectively, suggesting a non-competitive or mixed inhibition type. The half-maximal inhibitory concentration value (IC 50 ) for the inhibition of the formation of 4'-hydroxydiclofenace by limonene is determined to be 1413 µ M., (Creative Commons Attribution license.)

A primary objective in malaria vaccine design is the generation of high-quality antibody responses against the circumsporozoite protein of the malaria parasite, Plasmodium falciparum (PfCSP). To enable rational antigen design, we solved a cryo-EM structure of the highly potent anti-PfCSP antibody L9 in complex with recombinant PfCSP. We found that L9 Fab binds multivalently to the minor (NPNV) repeat domain, which is stabilized by a unique set of affinity-matured homotypic, antibody-antibody contacts. Molecular dynamics simulations revealed a critical role of the L9 light chain in integrity of the homotypic interface, which likely impacts PfCSP affinity and protective efficacy. These findings reveal the molecular mechanism of the unique NPNV selectivity of L9 and emphasize the importance of anti-homotypic affinity maturation in protective immunity against P. falciparum., (© 2023. The Author(s).)

The adaptive immune system arose around 500 million years ago in jawed fish, and, since then, it has mediated the immune defense against pathogens in all vertebrates. Antibodies play a central role in the immune reaction, recognizing and attacking external invaders. During the evolutionary process, several immunoglobulin isotypes emerged, each having a characteristic structural organization and dedicated function. In this work, we investigate the evolution of the immunoglobulin isotypes, in order to highlight the relevant features that were preserved over time and the parts that, instead, mutated. The residues that are coupled in the evolution process are often involved in intra- or interdomain interactions, meaning that they are fundamental to maintaining the immunoglobulin fold and to ensuring interactions with other domains. The explosive growth of available sequences allows us to point out the evolutionary conserved residues and compare the biophysical properties among different animal classes and isotypes. Our study offers a general overview of the evolution of immunoglobulin isotypes and advances the knowledge of their characteristic biophysical properties, as a first step in guiding protein design from evolution.

The proto-oncogene myc has been intensively studied primarily in vertebrate cell culture systems. Myc transcription factors control fundamental cellular processes such as cell proliferation, cell cycle control and stem cell maintenance. Myc interacts with the Max protein and Myc/Max heterodimers regulate thousands of target genes. The genome of the freshwater polyp Hydra encodes four myc genes ( myc1-4 ). Previous structural and biochemical characterization showed that the Hydra Myc1 and Myc2 proteins share high similarities with vertebrate c-Myc, and their expression patterns suggested a function in adult stem cell maintenance. In contrast, an additional Hydra Myc protein termed Myc3 is highly divergent, lacking the common N-terminal domain and all conserved Myc-boxes. Single cell transcriptome analysis revealed that the myc3 gene is expressed in a distinct population of interstitial precursor cells committed to nerve- and gland-cell differentiation, where the Myc3 protein may counteract the stemness actions of Myc1 and Myc2 and thereby allow the implementation of a differentiation program. In vitro DNA binding studies showed that Myc3 dimerizes with Hydra Max, and this dimer efficiently binds to DNA containing the canonical Myc consensus motif (E-box). In vivo cell transformation assays in avian fibroblast cultures further revealed an unexpected high potential for oncogenic transformation in the conserved Myc3 C-terminus, as compared to Hydra Myc2 or Myc1. Structure modeling of the Myc3 protein predicted conserved amino acid residues in its bHLH-LZ domain engaged in Myc3/Max dimerization. Mutating these amino acid residues in the human c-Myc (MYC) sequence resulted in a significant decrease in its cell transformation potential. We discuss our findings in the context of oncogenic transformation and cell differentiation, both relevant for human cancer, where Myc represents a major driver.

Cavβ subunits are essential for surface expression of voltage-gated calcium channel complexes and crucially modulate biophysical properties like voltage-dependent inactivation. Here, we describe the discovery and characterization of a novel Cavβ 2 variant with distinct features that predominates in the retina. We determined spliced exons in retinal transcripts of the Cacnb2 gene, coding for Cavβ 2 , by RNA-Seq data analysis and quantitative PCR. We cloned a novel Cavβ 2 splice variant from mouse retina, which we are calling β 2i , and investigated biophysical properties of calcium currents with this variant in a heterologous expression system as well as its intrinsic membrane interaction when expressed alone. Our data showed that β 2i predominated in the retina with expression in photoreceptors and bipolar cells. Furthermore, we observed that the β 2i N-terminus exhibited an extraordinary concentration of hydrophobic residues, a distinct feature not seen in canonical variants. The biophysical properties resembled known membrane-associated variants, and β 2i exhibited both a strong membrane association and a propensity for clustering, which depended on hydrophobic residues in its N-terminus. We considered available Cavβ structure data to elucidate potential mechanisms underlying the observed characteristics but resolved N-terminus structures were lacking and thus, precluded clear conclusions. With this description of a novel N-terminus variant of Cavβ 2 , we expand the scope of functional variation through N-terminal splicing with a distinct form of membrane attachment. Further investigation of the molecular mechanisms underlying the features of β 2i could provide new angles on the way Cavβ subunits modulate Ca 2+ channels at the plasma membrane., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Therapeutic antibodies should not only recognize antigens specifically, but also need to be free from developability issues, such as poor stability. Thus, the mechanistic understanding and characterization of stability are critical determinants for rational antibody design. In this study, we use molecular dynamics simulations to investigate the melting process of 16 antigen binding fragments (Fabs). We describe the Fab dissociation mechanisms, showing a separation in the V H -V L and in the C H 1-C L domains. We found that the depths of the minima in the free energy curve, corresponding to the bound states, correlate with the experimentally determined melting temperatures. Additionally, we provide a detailed structural description of the dissociation mechanism and identify key interactions in the CDR loops and in the C H 1-C L interface that contribute to stabilization. The dissociation of the V H -V L or C H 1-C L domains can be represented by conformational changes in the bend angles between the domains. Our findings elucidate the melting process of antigen binding fragments and highlight critical residues in both the variable and constant domains, which are also strongly germline dependent. Thus, our proposed mechanisms have broad implications in the development and design of new and more stable antigen binding fragments., (© 2023. The Author(s).)

Germline gain-of-function missense variants in the pore-forming Cav1.3 α1-subunit (CACNA1D gene) confer high risk for a severe neurodevelopmental disorder with or without endocrine symptoms. Here, we report a 4-week-old new-born with the novel de novo missense variant F747S with a so far not described prominent jittering phenotype in addition to symptoms previously reported for CACNA1D mutations including developmental delay, elevated aldosterone level and transient hypoglycemia. We confirmed the pathogenicity of this variant in whole-cell patch-clamp experiments with wild-type and F747S mutant channels heterologously expressed together with α2δ1 and cytosolic β3 or membrane-bound β2a subunits. Mutation F747S caused the quantitatively largest shift in the voltage dependence of activation (-28 mV) reported so far for CACNA1D germline mutations. It also shifted inactivation to more negative voltages, slowed the time course of current inactivation and slowed current deactivation upon repolarization with both co-expressed β-subunits. In silico modelling and molecular docking, simulations revealed that this gain-of-function phenotype can be explained by formation of a novel inter-domain hydrogen bond between mutant residues S747 (IIS6) with N1145 (IIIS6) stabilizing selectively the activated open channel state. F747S displayed 2-6-fold increased sensitivity for the L-type Ca2+ channel blocker isradipine compared to wild type. Our data confirm the pathogenicity of the F747S variant with very strong gain-of-function gating changes, which may contribute to the novel jittering phenotype. Increased sensitivity for isradipine suggests this drug for potential symptomatic off-label treatment for carriers of this mutation., (© The Author(s) 2022. Published by Oxford University Press.)

Camelid heavy-chain antibody variable domains (V H H), nanobodies, are the smallest-known functional antibody fragments with high therapeutic potential. In this study, we investigate a V H H binding to hen egg-white lysozyme (HEL). We structurally and dynamically characterized the conformational diversity of four V H H variants to elucidate the antigen-binding process. For two of these antibodies, not only are the dissociation constants known, but also the experimentally determined crystal structures of the V H H in complex with HEL are available. We performed well-tempered metadynamics simulations in combination with molecular dynamics simulations to capture a broad conformational space and to reconstruct the thermodynamics and kinetics of conformational transitions in the antigen-binding site, the paratope. By kinetically characterizing the loop movements of the paratope, we found that, with an increase in affinity, the state populations shift towards the binding competent conformation. The contacts contributing to antigen binding, and those who contribute to the overall stability, show a clear trend towards less variable but more intense contacts. Additionally, these investigated nanobodies clearly follow the conformational selection paradigm, as the binding competent conformation pre-exists within the structural ensembles without the presence of the antigen.

Advances in structural biology and the exponential increase in the amount of high-quality experimental structural data available in the Protein Data Bank has motivated numerous studies to tackle the grand challenge of predicting protein structures. In 2020 AlphaFold2 revolutionized the field using a combination of artificial intelligence and the evolutionary information contained in multiple sequence alignments. Antibodies are one of the most important classes of biotherapeutic proteins. Accurate structure models are a prerequisite to advance biophysical property predictions and consequently antibody design. Specialized tools used to predict antibody structures based on different principles have profited from current advances in protein structure prediction based on artificial intelligence. Here, we emphasize the importance of reliable protein structure models and highlight the enormous advances in the field, but we also aim to increase awareness that protein structure models, and in particular antibody models, may suffer from structural inaccuracies, namely incorrect cis-amide bonds, wrong stereochemistry or clashes. We show that these inaccuracies affect biophysical property predictions such as surface hydrophobicity. Thus, we stress the importance of carefully reviewing protein structure models before investing further computing power and setting up experiments. To facilitate the assessment of model quality, we provide a tool "TopModel" to validate structure models.

Beyond potency, a good developability profile is a key attribute of a biological drug. Selecting and screening for such attributes early in the drug development process can save resources and avoid costly late-stage failures. Here, we review some of the most important developability properties that can be assessed early on for biologics. These include the influence of the source of the biologic, its biophysical and pharmacokinetic properties, and how well it can be expressed recombinantly. We furthermore present in silico, in vitro , and in vivo methods and techniques that can be exploited at different stages of the discovery process to identify molecules with liabilities and thereby facilitate the selection of the most optimal drug leads. Finally, we reflect on the most relevant developability parameters for injectable versus orally delivered biologics and provide an outlook toward what general trends are expected to rise in the development of biologics.

The family of profilin allergens is a common class of proteins found in plants, viruses and various eukaryotes including mammals. Profilins are characterized by an evolutionary conserved structural fold, which is responsible for their cross-reactive nature of Immunoglobulin E (IgE) antibodies. Despite their high overall structural similarity, they exhibit substantial differences in their biophysical properties, such as thermal and pH stability. To understand the origin of these functional differences of Amb a 8, Art v 4 and Bet v 2, we performed constant pH molecular dynamics simulation in combination with Gaussian accelerated MD simulations. Depending on the respective protonation at different pH levels, we find distinct differences in conformational flexibility, which are consistent with experimentally determined melting temperatures. These variations in flexibility are accompanied by ensemble shifts in the conformational landscape and quantified and localized by residue-wise B-factors and dihedral entropies. These findings strengthen the link between flexibility of profilin allergens and their thermal stability. Thus, our results clearly show the importance of considering protonation dependent conformational ensembles in solution to elucidate biophysical differences between these structurally similar allergens., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (© 2022 Hofer, Fischer, Kamenik, Waibl, Fernández-Quintero and Liedl.)

The skeletal muscle voltage-gated calcium channel (CaV1.1) primarily functions as a voltage sensor for excitation-contraction coupling. Conversely, its ion-conducting function is modulated by multiple mechanisms within the pore-forming α1S subunit and the auxiliary α2δ-1 and γ1 subunits. In particular, developmentally regulated alternative splicing of exon 29, which inserts 19 amino acids in the extracellular IVS3-S4 loop of CaV1.1a, greatly reduces the current density and shifts the voltage dependence of activation to positive potentials outside the physiological range. We generated new HEK293 cell lines stably expressing α2δ-1, β3, and STAC3. When the adult (CaV1.1a) and embryonic (CaV1.1e) splice variants were expressed in these cells, the difference in the voltage dependence of activation observed in muscle cells was reproduced, but not the reduced current density of CaV1.1a. Only when we further coexpressed the γ1 subunit was the current density of CaV1.1a, but not that of CaV1.1e, reduced by >50%. In addition, γ1 caused a shift of the voltage dependence of inactivation to negative voltages in both variants. Thus, the current-reducing effect of γ1, unlike its effect on inactivation, is specifically dependent on the inclusion of exon 29 in CaV1.1a. Molecular structure modeling revealed several direct ionic interactions between residues in the IVS3-S4 loop and the γ1 subunit. However, substitution of these residues by alanine, individually or in combination, did not abolish the γ1-dependent reduction of current density, suggesting that structural rearrangements in CaV1.1a induced by inclusion of exon 29 may allosterically empower the γ1 subunit to exert its inhibitory action on CaV1.1 calcium currents., (© 2022 El Ghaleb et al.)

Sharks and other cartilaginous fish produce new antigen receptor (IgNAR) antibodies, as key part of their humoral immune response and are the phylogenetically oldest living organisms that possess an immunoglobulin (Ig)-based adaptive immune system. IgNAR antibodies are naturally occurring heavy-chain-only antibodies, that recognize antigens with their single domain variable regions (VNARs). In this study, we structurally and biophysically elucidate the effect of antibody humanization of a previously published spiny dogfish VNAR (parent E06), which binds with high affinity to the human serum albumin (HSA). We analyze different humanization variants together with the parental E06 VNAR and the human Vκ1 light chain germline DPK9 antibody to characterize the influence of point mutations in the framework and the antigen binding site on the specificity of VNARs as reported by Kovalenko et al. We find substantially higher flexibility in the humanized variants, reflected in a broader conformational space and a higher conformational entropy, as well as population shifts of the dominant binding site ensembles in solution. A further variant, in which some mutations are reverted, largely restores the conformational stability and the dominant binding minimum of the parent E06. We also identify differences in surface hydrophobicity between the human Vκ1 light chain germline DPK9 antibody, the parent VNAR E06 and the humanized variants. Additional simulations of VNAR-HSA complexes of the parent E06 VNAR and a humanized variant reveal that the parent VNAR features a substantially stronger network of stabilizing interactions. Thus, we conclude that a structural and dynamic understanding of the VNAR binding site upon humanization is a key aspect in antibody humanization., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Fernández-Quintero, Fischer, Kokot, Waibl, Seidler and Liedl.)

While antibody-based therapeutics have grown to be one of the major classes of novel medicines, some antibody development candidates face significant challenges regarding expression levels, solubility, as well as stability and aggregation, under physiological and storage conditions. A major determinant of those properties is surface hydrophobicity, which promotes unspecific interactions and has repeatedly proven problematic in the development of novel antibody-based drugs. Multiple computational methods have been devised for in-silico prediction of antibody hydrophobicity, often using hydrophobicity scales to assign values to each amino acid. Those approaches are usually validated by their ability to rank potential therapeutic antibodies in terms of their experimental hydrophobicity. However, there is significant diversity both in the hydrophobicity scales and in the experimental methods, and consequently in the performance of in-silico methods to predict experimental results. In this work, we investigate hydrophobicity of monoclonal antibodies using hydrophobicity scales. We implement several scoring schemes based on the solvent-accessibility and the assigned hydrophobicity values, and compare the different scores and scales based on their ability to predict retention times from hydrophobic interaction chromatography. We provide an overview of the strengths and weaknesses of several commonly employed hydrophobicity scales, thereby improving the understanding of hydrophobicity in antibody development. Furthermore, we test several datasets, both publicly available and proprietary, and find that the diversity of the dataset affects the performance of hydrophobicity scores. We expect that this work will provide valuable guidelines for the optimization of biophysical properties in future drug discovery campaigns., Competing Interests: HK and GG are employees of Roche Diagnostics GmbH. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Waibl, Fernández-Quintero, Wedl, Kettenberger, Georges and Liedl.)

PR-10 proteins constitute a major cause of food allergic reactions. Birch-pollen-related food allergies are triggered by the immunologic cross-reactivity of IgE antibodies with structurally homologous PR-10 proteins that are present in birch pollen and various food sources. While the three-dimensional structures of PR-10 food allergens have been characterized in detail, only a few experimental studies have addressed the structural flexibility of these proteins. In this study, we analyze the millisecond-timescale structural flexibility of thirteen PR-10 proteins from prevalent plant food sources by NMR relaxation-dispersion spectroscopy, in a comparative manner. We show that all the allergens in this study have inherently flexible protein backbones in solution, yet the extent of the structural flexibility appears to be strikingly protein-specific (but not food-source-specific). Above-average flexibility is present in the two short helices, α1 and α2, which form a V-shaped support for the long C-terminal helix α3, and shape the internal ligand-binding cavity, which is characteristic for PR-10 proteins. An in-depth analysis of the NMR relaxation-dispersion data for the PR-10 allergen from peanut reveals the presence of at least two subglobal conformational transitions on the millisecond timescale, which may be related to the release of bound low-molecular-weight ligands from the internal cavity.

Grid Inhomogeneous Solvation Theory (GIST) has proven useful to calculate localized thermodynamic properties of water around a solute. Numerous studies have leveraged this information to enhance structure-based binding predictions. We have recently extended GIST toward chloroform as a solvent to allow the prediction of passive membrane permeability. Here, we further generalize the GIST algorithm toward all solvents that can be modeled as rigid molecules. This restriction is inherent to the method and is already present in the inhomogeneous solvation theory. Here, we show that our approach can be applied to various solvent molecules by comparing the results of GIST simulations with thermodynamic integration (TI) calculations and experimental results. Additionally, we analyze and compare a matrix consisting of 100 entries of ten different solvent molecules solvated within each other. We find that the GIST results are highly correlated with TI calculations as well as experiments. For some solvents, we find Pearson correlations of up to 0.99 to the true entropy, while others are affected by the first-order approximation more strongly. The enthalpy-entropy splitting provided by GIST allows us to extend a recently published approach, which estimates higher order entropies by a linear scaling of the first-order entropy, to solvents other than water. Furthermore, we investigate the convergence of GIST in different solvents. We conclude that our extension to GIST reliably calculates localized thermodynamic properties for different solvents and thereby significantly extends the applicability of this widely used method.

Variable domains of camelid antibodies (so-called nanobodies or V H H) are the smallest antibody fragments that retain complete functionality and therapeutic potential. Understanding of the nanobody-binding interface has become a pre-requisite for rational antibody design and engineering. The nanobody-binding interface consists of up to three hypervariable loops, known as the CDR loops. Here, we structurally and dynamically characterize the conformational diversity of an anti-GFP-binding nanobody by using molecular dynamics simulations in combination with experimentally derived data from nuclear magnetic resonance (NMR) spectroscopy. The NMR data contain both structural and dynamic information resolved at various timescales, which allows an assessment of the quality of protein MD simulations. Thus, in this study, we compared the ensembles for the anti-GFP-binding nanobody obtained from MD simulations with results from NMR. We find excellent agreement of the NOE-derived distance maps obtained from NMR and MD simulations and observe similar conformational spaces for the simulations with and without NOE time-averaged restraints. We also compare the measured and calculated order parameters and find generally good agreement for the motions observed in the ps-ns timescale, in particular for the CDR3 loop. Understanding of the CDR3 loop dynamics is especially critical for nanobodies, as this loop is typically critical for antigen recognition.

Plasmalogens are an abundant class of glycerophospholipids in the mammalian body, with special occurrence in the brain and in immune cell membranes. Plasmanylethanolamine desaturase (PEDS1) is the final enzyme of plasmalogen biosynthesis, which introduces the characteristic 1-O-alk-1'-enyl double bond. The recent sequence identification of PEDS1 as transmembrane protein 189 showed that its protein sequence is related to a special class of plant desaturases (FAD4), with whom it shares a motif of 8 conserved histidines, which are essential for the enzymatic activity. In the present work, we wanted to gain more insight into the sequence-function relationship of this enzyme and mutated to alanine additional 28 amino acid residues of murine plasmanylethanolamine desaturase including those 20 residues, which are also totally conserved-in addition to the eight-histidine-motif-among the animal PEDS1 and plant FAD4 plant desaturases. We measured the enzymatic activity by transient transfection of tagged murine PEDS1 expression clones to a PEDS1-deficient human HAP1 cell line by monitoring of labeled plasmalogens formed from supplemented 1-O-pyrenedecyl-sn-glycerol in relation to recombinant protein expression. Surprisingly, only a single mutation, namely aspartate 100, led to a total loss of PEDS1 activity. The second strongest impact on enzymatic activity had mutation of phenylalanine 118, leaving only 6% residual activity. A structural model obtained by homology modelling to available structures of stearoyl-CoA reductase predicted that this aspartate 100 residue interacts with histidine 96, and phenylalanine 118 interacts with histidine 187, both being essential histidines assumed to be involved in the coordination of the di-metal center of the enzyme., (© 2022. The Author(s).)

Structure-based antibody design and accurate predictions of antibody-antigen interactions remain major challenges in computational biology. By using molecular dynamics simulations, we show that a single static X-ray structure is not sufficient to identify determinants of antibody-antigen recognition. Here, we investigate antibodies that undergo substantial conformational changes upon antigen binding and have been classified as difficult cases in an extensive benchmark for antibody-antigen docking. We present thermodynamics and transition kinetics of these conformational rearrangements and show that paratope states can be used to improve antibody-antigen docking. By using the unbound antibody X-ray structure as starting structure for molecular dynamics simulations, we retain a binding competent conformation substantially different to the unbound antibody X-ray structure. We also observe that the kinetically dominant antibody paratope conformations are chosen by the bound antigen conformation with the highest probability. Thus, we show that paratope states in solution can improve antibody-antigen docking and structure prediction., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

A new format of therapeutic proteins is bispecific antibodies, in which two different heavy chains heterodimerize to obtain two different binding sites. Therefore, it is crucial to understand and optimize the third constant domain (CH3-CH3) interface to favor heterodimerization over homodimerization, and to preserve the physicochemical properties, as thermal stability. Here, we use molecular dynamics simulations to investigate the dissociation process of 19 CH3-CH3 crystal structures that differ from each other in few point mutations. We describe the dissociation of the dimeric interface as a two-steps mechanism. As confirmed by a Markov state model, apart from the bound and the dissociated state, we observe an additional intermediate state, which corresponds to an encounter complex. The analysis of the interdomain contacts reveals key residues that stabilize the interface. We expect that our results will improve the understanding of the CH3-CH3 interface interactions and thus advance the developability and design of new antibodies formats., (© The Author(s) 2022. Published by Oxford University Press.)

Broadly neutralizing antibodies that target epitopes of haemagglutinin on the influenza virus have the potential to provide near universal protection against influenza virus infection 1 . However, viral mutants that escape broadly neutralizing antibodies have been reported 2,3 . The identification of broadly neutralizing antibody classes that can neutralize viral escape mutants is critical for universal influenza virus vaccine design. Here we report a distinct class of broadly neutralizing antibodies that target a discrete membrane-proximal anchor epitope of the haemagglutinin stalk domain. Anchor epitope-targeting antibodies are broadly neutralizing across H1 viruses and can cross-react with H2 and H5 viruses that are a pandemic threat. Antibodies that target this anchor epitope utilize a highly restricted repertoire, which encodes two public binding motifs that make extensive contacts with conserved residues in the fusion peptide. Moreover, anchor epitope-targeting B cells are common in the human memory B cell repertoire and were recalled in humans by an oil-in-water adjuvanted chimeric haemagglutinin vaccine 4,5 , which is a potential universal influenza virus vaccine. To maximize protection against seasonal and pandemic influenza viruses, vaccines should aim to boost this previously untapped source of broadly neutralizing antibodies that are widespread in the human memory B cell pool., (© 2021. The Author(s).)

Hydration thermodynamics play a fundamental role in fields ranging from the pharmaceutical industry to environmental research. Numerous methods exist to predict solvation thermodynamics of compounds ranging from small molecules to large biomolecules. Arguably the most precise methods are those based on molecular dynamics (MD) simulations in explicit solvent. One theory that has seen increased use is inhomogeneous solvation theory (IST). However, while many applications require accurate description of salt-water mixtures, no implementation of IST is currently able to estimate solvation properties involving more than one solvent species. Here, we present an extension to grid inhomogeneous solvation theory (GIST) that can take salt contributions into account. At the example of carbazole in 1 M NaCl solution, we compute the solvation energy as well as first and second order entropies. While the effect of the first order ion entropy is small, both the water-water and water-ion entropies contribute strongly. We show that the water-ion entropies are efficiently approximated using the Kirkwood superposition approximation. However, this approach cannot be applied to the water-water entropy. Furthermore, we test the quantitative validity of our method by computing salting-out coefficients and comparing them to experimental data. We find a good correlation to experimental salting-out constants, while the absolute values are overpredicted due to the approximate second order entropy. Since ions are frequently used in MD, either to neutralize the system or as a part of the investigated process, our method greatly extends the applicability of GIST. The use-cases range from biopharmaceuticals, where many assays require high salt concentrations, to environmental research, where solubility in sea water is important to model the fate of organic substances., (© 2022. The Author(s).)

As the current biotherapeutic market is dominated by antibodies, the design of different antibody formats, like bispecific antibodies and other new formats, represent a key component in advancing antibody therapy. When designing new formats, a targeted modulation of pairing preferences is key. Several existing approaches are successful, but expanding the repertoire of design possibilities would be desirable. Cognate immunoglobulin G antibodies depend on homodimerization of the fragment crystallizable regions of two identical heavy chains. By modifying the dimeric interface of the third constant domain (C H 3-C H 3), with different mutations on each domain, the engineered Fc fragments form rather heterodimers than homodimers. The first constant domain (C H 1-C L ) shares a very similar fold and interdomain orientation with the C H 3-C H 3 dimer. Thus, numerous well-established design efforts for C H 3-C H 3 interfaces, have also been applied to C H 1-C L dimers to reduce the number of mispairings in the Fabs. Given the high structural similarity of the C H 3-C H 3 and C H 1-C L domains we want to identify additional opportunities in comparing the differences and overlapping interaction profiles. Our vision is to facilitate a toolkit that allows for the interchangeable usage of different design tools from crosslinking the knowledge between these two interface types. As a starting point, here, we use classical molecular dynamics simulations to identify differences of the C H 3-C H 3 and C H 1-C L interfaces and already find unexpected features of these interfaces shedding new light on possible design variations. Apart from identifying clear differences between the similar C H 3-C H 3 and C H 1-C L dimers, we structurally characterize the effects of point-mutations in the C H 3-C H 3 interface on the respective dynamics and interface interaction patterns. Thus, this study has broad implications in the field of antibody engineering as it provides a structural and mechanistical understanding of antibody interfaces and thereby presents a crucial aspect for the design of bispecific antibodies., Competing Interests: Authors AB, GG and HK were employed by the company Roche. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Fernández-Quintero, Quoika, Wedl, Seidler, Kroell, Loeffler, Pomarici, Hoerschinger, Bujotzek, Georges, Kettenberger and Liedl.)

As the current biotherapeutic market is dominated by antibodies, the design of different antibody formats, like bispecific antibodies, is critical to the advancement of the field. In contrast to monovalent antibodies, which consist of two identical antigen-binding sites, bispecific antibodies can target two different epitopes by containing two different antigen-binding sites. Thus, the rise of new formats as successful therapeutics has reignited the interest in advancing and facilitating the efficient production of bispecific antibodies. Here, we investigate the influence of point mutations in the antigen-binding site, the paratope, on heavy and light chain pairing preferences by using molecular dynamics simulations. In agreement with experiments, we find that specific residues in the antibody variable domain (Fv), i.e., the complementarity-determining region (CDR) L3 and H3 loops, determine heavy and light chain pairing preferences. Excitingly, we observe substantial population shifts in CDR-H3 and CDR-L3 loop conformations in solution accompanied by a decrease in bispecific IgG yield. These conformational changes in the CDR3 loops induced by point mutations also influence all other CDR loop conformations and consequentially result in different CDR loop states in solution. However, besides their effect on the obtained CDR loop ensembles, point mutations also lead to distinct interaction patterns in the V H -V L interface. By comparing the interaction patterns among all investigated variants, we observe specific contacts in the interface that drive heavy and light chain pairing. Thus, these findings have broad implications in the field of antibody engineering and design because they provide a mechanistic understanding of antibody interfaces, by identifying critical factors driving the pairing preferences, and thus can help to advance the design of bispecific antibodies.

The voltage-gated calcium channel Ca V 1.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 Ca V 1.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 Ca V 1.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 Ca V 1.1 and Ca V 1.2 L-type calcium channels. Moreover, in Ca V 1.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-Ca V 1.1e (lacking exon 29) showed that charge neutralization (E216Q) or removal of the side chain (E216A) significantly shifted the voltage dependence of activation (V 1/2 ) to more positive potentials, suggesting that E216 stabilizes the activated state. Insertion of exon 29 in the GFP-Ca V 1.1a splice variant strongly reduced the ionic interactions with R1 and R2 and caused a substantial right shift of V 1/2 , whereas no further shift of V 1/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 Ca V 1.1., (Copyright © 2021 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

T-type calcium channels (Cav3.1 to Cav3.3) regulate low-threshold calcium spikes, burst firing and rhythmic oscillations of neurons and are involved in sensory processing, sleep, and hormone and neurotransmitter release. Here, we examined four heterozygous missense variants in CACNA1I, encoding the Cav3.3 channel, in patients with variable neurodevelopmental phenotypes. The p.(Ile860Met) variant, affecting a residue in the putative channel gate at the cytoplasmic end of the IIS6 segment, was identified in three family members with variable cognitive impairment. The de novo p.(Ile860Asn) variant, changing the same amino acid residue, was detected in a patient with severe developmental delay and seizures. In two additional individuals with global developmental delay, hypotonia, and epilepsy, the variants p.(Ile1306Thr) and p.(Met1425Ile), substituting residues at the cytoplasmic ends of IIIS5 and IIIS6, respectively, were found. Because structure modelling indicated that the amino acid substitutions differentially affect the mobility of the channel gate, we analysed possible effects on Cav3.3 channel function using patch-clamp analysis in HEK293T cells. The mutations resulted in slowed kinetics of current activation, inactivation, and deactivation, and in hyperpolarizing shifts of the voltage-dependence of activation and inactivation, with Cav3.3-I860N showing the strongest and Cav3.3-I860M the weakest effect. Structure modelling suggests that by introducing stabilizing hydrogen bonds the mutations slow the kinetics of the channel gate and cause the gain-of-function effect in Cav3.3 channels. The gating defects left-shifted and increased the window currents, resulting in increased calcium influx during repetitive action potentials and even at resting membrane potentials. Thus, calcium toxicity in neurons expressing the Cav3.3 variants is one likely cause of the neurodevelopmental phenotype. Computer modelling of thalamic reticular nuclei neurons indicated that the altered gating properties of the Cav3.3 disease variants lower the threshold and increase the duration and frequency of action potential firing. Expressing the Cav3.3-I860N/M mutants in mouse chromaffin cells shifted the mode of firing from low-threshold spikes and rebound burst firing with wild-type Cav3.3 to slow oscillations with Cav3.3-I860N and an intermediate firing mode with Cav3.3-I860M, respectively. Such neuronal hyper-excitability could explain seizures in the patient with the p.(Ile860Asn) mutation. Thus, our study implicates CACNA1I gain-of-function mutations in neurodevelopmental disorders, with a phenotypic spectrum ranging from borderline intellectual functioning to a severe neurodevelopmental disorder with epilepsy., (© The Author(s) (2021). Published by Oxford University Press on behalf of the Guarantors of Brain.)

Antibodies have emerged as one of the fastest growing classes of biotherapeutic proteins. To improve the rational design of antibodies, we investigate the conformational diversity of 16 different germline combinations, which are composed of 4 different kappa light chains paired with 4 different heavy chains. In this study, we systematically show that different heavy and light chain pairings strongly influence the paratope, interdomain interaction patterns and the relative V H -V L interface orientations. We observe changes in conformational diversity and substantial population shifts of the complementarity determining region (CDR) loops, resulting in distinct dominant solution structures and differently favored canonical structures. Additionally, we identify conformational changes in the structural diversity of the CDR-H3 loop upon different heavy and light chain pairings, as well as upon changes in sequence and structure of the neighboring CDR loops, despite having an identical CDR-H3 loop amino acid sequence. These results can also be transferred to all CDR loops and to the relative V H -V L orientation, as certain paratope states favor distinct interface angle distributions. Furthermore, we directly compare the timescales of sidechain rearrangements with the well-described transition kinetics of conformational changes in the backbone of the CDR loops. We show that sidechain flexibilities are strongly affected by distinct heavy and light chain pairings and decipher germline-specific structural features co-determining stability. These findings reveal that all CDR loops are strongly correlated and that distinct heavy and light chain pairings can result in different paratope states in solution, defined by a characteristic combination of CDR loop conformations and V H -V L interface orientations. Thus, these results have broad implications in the field of antibody engineering, as they clearly show the importance of considering paired heavy and light chains to understand the antibody binding site, which is one of the key aspects in the design of therapeutics., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Fernández-Quintero, Kroell, Bacher, Loeffler, Quoika, Georges, Bujotzek, Kettenberger and Liedl.)

Skeletal muscle excitation-contraction (EC) coupling roots in Ca 2+ -influx-independent inter-channel signaling between the sarcolemmal dihydropyridine receptor (DHPR) and the ryanodine receptor (RyR1) in the sarcoplasmic reticulum. Although DHPR Ca 2+ influx is irrelevant for EC coupling, its putative role in other muscle-physiological and developmental pathways was recently examined using two distinct genetically engineered mouse models carrying Ca 2+ non-conducting DHPRs: DHPR(N617D) (Dayal et al., 2017) and DHPR(E1014K) (Lee et al., 2015). Surprisingly, despite complete block of DHPR Ca 2+ -conductance, histological, biochemical, and physiological results obtained from these two models were contradictory. Here, we characterize the permeability and selectivity properties and henceforth the mechanism of Ca 2+ non-conductance of DHPR(N617). Our results reveal that only mutant DHPR(N617D) with atypical high-affinity Ca 2+ pore-binding is tight for physiologically relevant monovalent cations like Na + and K + . Consequently, we propose a molecular model of cooperativity between two ion selectivity rings formed by negatively charged residues in the DHPR pore region., Competing Interests: AD, MF, KL, MG No competing interests declared, (© 2021, Dayal et al.)

We revived and implemented a method developed by Kuhn in 1934, originally only published in German, that is, the so-called "freely jointed chain" model. This approach turned out to be surprisingly useful for analyzing state-of-the-art computer simulations of the thermosensitive coil-globule transition of N -Isopropylacrylamide 20-mer. Our atomistic computer simulations are orders of magnitude longer than those of previous studies and lead to a reliable description of thermodynamics and kinetics at many different temperatures. The freely jointed chain model provides a coordinate system, which allows us to construct a Markov state model of the conformational transitions. Furthermore, this guarantees a reliable reconstruction of the kinetics in back-and-forth directions. In addition, we obtain a description of the high diversity and variability of both conformational states. Thus, we gain a detailed understanding of the coil-globule transition. Surprisingly, conformational entropy turns out to play only a minor role in the thermodynamic balance of the process. Moreover, we show that the radius of gyration is an unexpectedly unsuitable coordinate to comprehend the transition kinetics because it does not capture the high conformational diversity within the different states. Consequently, the approach presented here allows for an exhaustive description and resolution of the conformational ensembles of arbitrary linear polymer chains.