Inactivation induced by pathogenic Cav1.3 L-type Ca2+-channel variants enhances sensitivity for dihydropyridine Ca2+ channel blockers.

BACKGROUND AND PURPOSE: Pathogenic gain-of-function mutations in Cav1.3 L-type voltage-gated Ca2+-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 Ca2+-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 IC50 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 Cav1.3 channels in neurons during prolonged depolarized states without selectivity for mutant channels.

On the humanization of VHHs: Prospective case studies, experimental and computational characterization of structural determinants for functionality.

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.

CosolvKit: a Versatile Tool for Cosolvent MD Preparation and Analysis.

Cosolvent molecular dynamics (MDs) are an increasingly popular form of simulations where small molecule cosolvents are added to water-solvated protein systems. These simulations can perform diverse target characterization tasks, including cryptic and allosteric pocket identification and pharmacophore profiling and supplement suites of enhanced sampling methods to explore protein conformational landscapes. The behavior of these systems is tied to the cosolvents used, so the ability to define diverse and complex mixtures is critical in dictating the outcome of the simulations. However, existing methods for preparing cosolvent simulations only support a limited number of predefined cosolvents and concentrations. Here, we present CosolvKit, a tool for the preparation and analysis of systems composed of user-defined cosolvents and concentrations. This tool is modular, supporting the creation of files for multiple MD engines, as well as direct access to OpenMM simulations, and offering access to a variety of generalizable small-molecule force fields. To the best of our knowledge, CosolvKit represents the first generalized approach for the construction of these simulations.

Assessing AF2's ability to predict structural ensembles of proteins.

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.

Engineering of pH-dependent antigen binding properties for toxin-targeting IgG1 antibodies using light-chain shuffling.

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.

CaV1.1 voltage-sensing domain III exclusively controls skeletal muscle excitation-contraction coupling.

Skeletal muscle contractions are initiated by action potentials, which are sensed by the voltage-gated calcium channel (CaV1.1) and are conformationally coupled to calcium release from intracellular stores. Notably, CaV1.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 CaV1.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.

Fit-for-purpose heterodivalent single-domain antibody for gastrointestinal targeting of toxin B from Clostridium difficile.

Single-domain antibodies (sdAbs), such as VHHs, 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 VHH constructs, that display extraordinary developability properties for oral delivery and functionality in the GI environment. We showcase this by developing a heterodivalent VHH 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 VHH 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 VHHs and related constructs fit for oral delivery and GI applications.

Pathogenicity of de novo CACNA1D Ca2+ channel variants predicted from sequence co-variation.

Voltage-gated L-type Cav1.3 Ca2+ 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.

The Role of Force Fields and Water Models in Protein Folding and Unfolding Dynamics.

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.

Defining the cross-reactivity between peanut allergens Ara h 2 and Ara h 6 using monoclonal antibodies.

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.

Asymmetric contribution of a selectivity filter gate in triggering inactivation of CaV1.3 channels.

Voltage-dependent and Ca2+-dependent inactivation (VDI and CDI, respectively) of CaV channels are two biologically consequential feedback mechanisms that fine-tune Ca2+ entry into neurons and cardiomyocytes. Although known to be initiated by distinct molecular events, how these processes obstruct conduction through the channel pore remains poorly defined. Here, focusing on ultrahighly conserved tryptophan residues in the interdomain interfaces near the selectivity filter of CaV1.3, we demonstrate a critical role for asymmetric conformational changes in mediating VDI and CDI. Specifically, mutagenesis of the domain III-IV interface, but not others, enhanced VDI. Molecular dynamics simulations demonstrate that mutations in distinct selectivity filter interfaces differentially impact conformational flexibility. Furthermore, mutations in distinct domains preferentially disrupt CDI mediated by the N- versus C-lobes of CaM, thus uncovering a scheme of structural bifurcation of CaM signaling. These findings highlight the fundamental importance of the asymmetric arrangement of the pseudotetrameric CaV pore domain for feedback inhibition.

Addressing Challenges of Macrocyclic Conformational Sampling in Polar and Apolar Solvents: Lessons for Chameleonicity.

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.

PEP-Patch: Electrostatics in Protein-Protein Recognition, Specificity, and Antibody Developability.

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.

Structural Basis of the Immunological Cross-Reactivity between Kiwi and Birch Pollen.

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.

Structure and Dynamics Guiding Design of Antibody Therapeutics and Vaccines.

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.

Asymmetric Contribution of a Selectivity Filter Gate in Triggering Inactivation of CaV1.3 Channels.

Voltage-dependent and Ca2+-dependent inactivation (VDI and CDI, respectively) of CaV channels are two biologically consequential feedback mechanisms that fine-tune Ca2+ entry into neurons and cardiomyocytes. Although known to be initiated by distinct molecular events, how these processes obstruct conduction through the channel pore remains poorly defined. Here, focusing on ultra-highly conserved tryptophan residues in the inter-domain interfaces near the selectivity filter of CaV1.3, we demonstrate a critical role for asymmetric conformational changes in mediating VDI and CDI. Specifically, mutagenesis of the domain III-IV interface, but not others, enhanced VDI. Molecular dynamics simulations demonstrate that mutations in distinct selectivity filter interfaces differentially impact conformational flexibility. Furthermore, mutations in distinct domains preferentially disrupt CDI mediated by the N- versus C-lobes of CaM, thus uncovering a scheme of structural bifurcation of CaM signaling. These findings highlight the fundamental importance of the asymmetric arrangement of the pseudo-tetrameric CaV pore domain for feedback inhibition.

Phage display assisted discovery of a pH-dependent anti-α-cobratoxin antibody from a natural variable domain library.

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.

Defining and Studying B Cell Receptor and TCR Interactions.

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.

Establishing a cell-based screening workflow for determining the efficiency of CYP2C9 metabolism: moving towards the use of breath volatiles in personalised medicine.

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 thein vitrobioconversion of diclofenac in recombinant HEK293T cells overexpressing CYP2C9 to 4'-hydroxydiclofenac. Thisin vitroapproach 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,KM, of 103µM and 94.1µM, respectively, and a maximum reaction rate,Vmax, of 46.8 pmol min-1106cells-1and 56.0 pmol min-1106cells-1with and without the inhibitor, respectively, suggesting a non-competitive or mixed inhibition type. The half-maximal inhibitory concentration value (IC50) for the inhibition of the formation of 4'-hydroxydiclofenace by limonene is determined to be 1413µM.

Cross-linking disulfide bonds govern solution structures of diabodies.

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 Fv 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 Fv domains. We observe rigidifying effects through the introduction of disulfide bonds in the Fv -Fv interface and characterize the effect of different disulfide bond locations on the conformation. Additionally, we compare VH -VL 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 Fv -Fv interface. All examined diabodies show similar VH -VL orientations, Fv -Fv 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.