Skeletal Dysplasia Mutations Effect on Human Filamins' Structure and Mechanosensing

Atomistic Insights into gp82 Binding: A Microsecond, Million-Atom Exploration of Trypanosoma cruzi Host-Cell Invasion

Chagas disease, caused by the protozoan Trypanosoma cruzi, affects millions globally, leading to severe cardiac and gastrointestinal complications in its chronic phase. The invasion of host cells by T. cruzi is mediated by the interaction between the parasite's glycoprotein gp82 and the human receptor lysosome-associated membrane protein 2 (LAMP2). While experimental studies have identified a few residues involved in this interaction, a comprehensive molecular-level understanding has been lacking. In this study, we present a 1.44-million-atom computational model of the gp82 complex, including over 3300 lipids, glycosylation sites, and full molecular representations of gp82 and LAMP2, making it the most complete model of a parasite-host interaction to date. Using microsecond-long molecular dynamics simulations and dynamic network analysis, we identified critical residue interactions, including novel regions of contact that were previously uncharacterized. Our findings also highlight the significance of the transmembrane domain of LAMP2 in stabilizing the complex. These insights extend beyond traditional hydrogen bond interactions, revealing a complex network of cooperative motions that facilitate T. cruzi invasion. This study not only confirms key experimental observations but also uncovers new molecular targets for therapeutic intervention, offering a potential pathway to disrupt T. cruzi infection and combat Chagas disease.

Integrating Dynamic Network Analysis with AI for Enhanced Epitope Prediction in PD-L1:Affibody Interactions

Understanding binding epitopes involved in protein-protein interactions and accurately determining their structure are long-standing goals with broad applicability in industry and biomedicine. Although various experimental methods for binding epitope determination exist, these approaches are typically low throughput and cost-intensive. Computational methods have potential to accelerate epitope predictions; however, recently developed artificial intelligence (AI)-based methods frequently fail to predict epitopes of synthetic binding domains with few natural homologues. Here we have developed an integrated method employing generalized-correlation-based dynamic network analysis on multiple molecular dynamics (MD) trajectories, initiated from AlphaFold2Multimer structures, to unravel the structure and binding epitope of the therapeutic PD-L1:Affibody complex. Both AlphaFold2 and conventional molecular dynamics trajectory analysis were ineffective in distinguishing between two proposed binding models, parallel and perpendicular. However, our integrated approach, utilizing dynamic network analysis, demonstrated that the perpendicular mode was significantly more stable. These predictions were validated using a suite of experimental epitope mapping protocols, including cross-linking mass spectrometry and next-generation sequencing-based deep mutational scanning. Conversely, AlphaFold3 failed to predict a structure bound in the perpendicular pose, highlighting the necessity for exploratory research in the search for binding epitopes and challenging the notion that AI-generated protein structures can be accepted without scrutiny. Our research underscores the potential of employing dynamic network analysis to enhance AI-based structure predictions for more accurate identification of protein-protein interaction interfaces.

Uncovering and engineering the mechanical properties of the adhesion GPCR ADGRG1 GAIN domain

Key cellular functions depend on the transduction of extracellular mechanical signals by specialized membrane receptors including adhesion G-protein coupled receptors (aGPCRs). While recently solved structures support aGPCR activation through shedding of the extracellular GAIN domain, the molecular mechanisms underpinning receptor mechanosensing remain poorly understood. When probed using single-molecule atomic force spectroscopy and molecular simulations, ADGRG1 GAIN dissociated from its tethered agonist at forces significantly higher than other reported signaling mechanoreceptors. Strong mechanical resistance was achieved through specific structural deformations and force propagation pathways under mechanical load. ADGRG1 GAIN variants computationally designed to lock the alpha and beta subdomains and rewire mechanically-induced structural deformations were found to modulate the GPS-Stachel rupture forces. Our study provides unprecedented insights into the molecular underpinnings of GAIN mechanical stability and paves the way for engineering mechanosensors, better understanding aGPCR function, and informing drug-discovery efforts targeting this important receptor class.

May the force be with you: The role of hyper-mechanostability of the bone sialoprotein binding protein during early stages of Staphylococci infections

The bone sialoprotein-binding protein (Bbp) is a mechanoactive MSCRAMM protein expressed on the surface of Staphylococcus aureus that mediates adherence of the bacterium to fibrinogen-α (Fgα), a component of the bone and dentine extracellular matrix of the host cell. Mechanoactive proteins like Bbp have key roles in several physiological and pathological processes. Particularly, the Bbp: Fgα interaction is important in the formation of biofilms, an important virulence factor of pathogenic bacteria. Here, we investigated the mechanostability of the Bbp: Fgα complex using in silico single-molecule force spectroscopy (SMFS), in an approach that combines results from all-atom and coarse-grained steered molecular dynamics (SMD) simulations. Our results show that Bbp is the most mechanostable MSCRAMM investigated thus far, reaching rupture forces beyond the 2 nN range in typical experimental SMFS pulling rates. Our results show that high force-loads, which are common during initial stages of bacterial infection, stabilize the interconnection between the protein's amino acids, making the protein more "rigid". Our data offer new insights that are crucial on the development of novel anti-adhesion strategies.

Intervertebral disc degeneration is rescued by TGFβ/BMP signaling modulation in an ex vivo filamin B mouse model

Spondylocarpotarsal syndrome (SCT) is a rare musculoskeletal disorder characterized by short stature and vertebral, carpal, and tarsal fusions resulting from biallelic nonsense mutations in the gene encoding filamin B (FLNB). Utilizing a FLNB knockout mouse, we showed that the vertebral fusions in SCT evolved from intervertebral disc (IVD) degeneration and ossification of the annulus fibrosus (AF), eventually leading to full trabecular bone formation. This resulted from alterations in the TGFβ/BMP signaling pathway that included increased canonical TGFβ and noncanonical BMP signaling. In this study, the role of FLNB in the TGFβ/BMP pathway was elucidated using in vitro, in vivo, and ex vivo treatment methodologies. The data demonstrated that FLNB interacts with inhibitory Smads 6 and 7 (i-Smads) to regulate TGFβ/BMP signaling and that loss of FLNB produces increased TGFβ receptor activity and decreased Smad 1 ubiquitination. Through the use of small molecule inhibitors in an ex vivo spine model, TGFβ/BMP signaling was modulated to design a targeted treatment for SCT and disc degeneration. Inhibition of canonical and noncanonical TGFβ/BMP pathway activity restored Flnb-/- IVD morphology. These most effective improvements resulted from specific inhibition of TGFβ and p38 signaling activation. FLNB acts as a bridge for TGFβ/BMP signaling crosstalk through i-Smads and is key for the critical balance in TGFβ/BMP signaling that maintains the IVD. These findings further our understanding of IVD biology and reveal new molecular targets for disc degeneration as well as congenital vertebral fusion disorders.

Mechanisms of Nanonewton Mechanostability in a Protein Complex Revealed by Molecular Dynamics Simulations and Single-Molecule Force Spectroscopy

Can molecular dynamics simulations predict the mechanical behavior of protein complexes? Can simulations decipher the role of protein domains of unknown function in large macromolecular complexes? Here, we employ a wide-sampling computational approach to demonstrate that molecular dynamics simulations, when carefully performed and combined with single-molecule atomic force spectroscopy experiments, can predict and explain the behavior of highly mechanostable protein complexes. As a test case, we studied a previously unreported homologue from Ruminococcus flavefaciens called X-module-Dockerin (XDoc) bound to its partner Cohesin (Coh). By performing dozens of short simulation replicas near the rupture event, and analyzing dynamic network fluctuations, we were able to generate large simulation statistics and directly compare them with experiments to uncover the mechanisms involved in mechanical stabilization. Our single-molecule force spectroscopy experiments show that the XDoc-Coh homologue complex withstands forces up to 1 nN at loading rates of 10⁵ pN/s. Our simulation results reveal that this remarkable mechanical stability is achieved by a protein architecture that directs molecular deformation along paths that run perpendicular to the pulling axis. The X-module was found to play a crucial role in shielding the adjacent protein complex from mechanical rupture. These mechanisms of protein mechanical stabilization have potential applications in biotechnology for the development of systems exhibiting shear enhanced adhesion or tunable mechanics.

Critical Structural Defects Explain Filamin A Mutations Causing Mitral Valve Dysplasia

Mitral valve diseases affect ∼3% of the population and are the most common reasons for valvular surgery because no drug-based treatments exist. Inheritable genetic mutations have now been established as the cause of mitral valve insufficiency, and four different missense mutations in the filamin A gene (FLNA) have been found in patients suffering from nonsyndromic mitral valve dysplasia (MVD). The filamin A (FLNA) protein is expressed, in particular, in endocardial endothelia during fetal valve morphogenesis and is key in cardiac development. The FLNA-MVD-causing mutations are clustered in the N-terminal region of FLNA. How the mutations in FLNA modify its structure and function has mostly remained elusive. In this study, using NMR spectroscopy and interaction assays, we investigated FLNA-MVD-causing V711D and H743P mutations. Our results clearly indicated that both mutations almost completely destroyed the folding of the FLNA5 domain, where the mutation is located, and also affect the folding of the neighboring FLNA4 domain. The structure of the neighboring FLNA6 domain was not affected by the mutations. These mutations also completely abolish FLNA's interactions with protein tyrosine phosphatase nonreceptor type 12, which has been suggested to contribute to the pathogenesis of FLNA-MVD. Taken together, our results provide an essential structural and molecular framework for understanding the molecular bases of FLNA-MVD, which is crucial for the development of new therapies to replace surgery.

Non-syndromic Mitral Valve Dysplasia Mutation Changes the Force Resilience and Interaction of Human Filamin A

Filamin A (FLNa), expressed in endocardial endothelia during fetal valve morphogenesis, is key in cardiac development. Missense mutations in FLNa cause non-syndromic mitral valve dysplasia (FLNA-MVD). Here, we aimed to reveal the currently unknown underlying molecular mechanism behind FLNA-MVD caused by the FLNa P637Q mutation. The solved crystal structure of the FLNa3-5 P637Q revealed that this mutation causes only minor structural changes close to mutation site. These changes were observed to significantly affect FLNa's ability to transmit cellular force and to interact with its binding partner. The performed steered molecular dynamics simulations showed that significantly lower forces are needed to split domains 4 and 5 in FLNA-MVD than with wild-type FLNa. The P637Q mutation was also observed to interfere with FLNa's interactions with the protein tyrosine phosphatase PTPN12. Our results provide a crucial step toward understanding the molecular bases behind FLNA-MVD, which is critical for the development of drug-based therapeutics.

A comprehensive analysis of the allergenicity and IgE epitopes of myosinogen allergens in Scylla paramamosain

BACKGROUND

Scylla paramamosain is one of the most common and serious food allergens in Asia. Therefore, research on its prevalence, accurate diagnosis, and IgE-binding pattern of the allergens is crucial.

OBJECTIVE

To identify the IgE epitopes of the myosinogen allergens in S. paramamosain using phage peptide library.

METHODS

The prevalence of allergy to crabs (AC) and of sensitization was analysed using a questionnaire and a serological assay. BAT was performed by flow cytometry, and its diagnostic performance was evaluated in relation to allergens purified from crab myosinogen. IgE-binding epitopes were identified by phage display using the IgE from patients with AC. Sequence- and structure-based bioinformatics analyses were performed to identify allergenic epitopes.

RESULTS

Crab was the most common cause of food allergies in this study. Subjects with AC (n = 30) with clear clinical symptoms were identified by immunoblotting and BAT. All of the myosinogen allergens triggered basophil activation; surface expression of CD63 and CD203c was higher in patients allergic to AK and FLN c than in patients allergic to SCP and TIM. In addition to six conformational epitopes of SCP, six linear epitopes and eight conformational epitopes of AK were identified. Five linear epitopes and three conformational epitopes of TIM, nine linear and ten conformational epitopes of FLN c were also identified, and the sequence VH(I/T) L was appeared in epitopes of both TIM and FLN c. The number of epitopes showed consistency with the value of BAT.

CONCLUSIONS AND CLINICAL RELEVANCE

BAT can be used for accurate diagnosis of AC. Identification of particular allergenic motifs could be a valuable tool for prevention, diagnosis, and treatment of food allergies.

The filamin-B–refilin axis – spatiotemporal regulators of the actin-cytoskeleton in development and disease

During development, cycles of spatiotemporal remodeling of higher-order networks of actin filaments contribute to control cell fate specification and differentiation. Programs for controlling these dynamics are hard-wired into actin-regulatory proteins. The filamin family of actin-binding proteins exert crucial mechanotransduction and signaling functions in tissue morphogenesis. Filamin-B (FLNB) is a key player in chondrocyte progenitor differentiation for endochondral ossification. Biallelic loss-of-function mutations or gain-of-function mutations in FLNB cause two groups of skeletal disorders that can be attributed to either the loss of repressive function on TGF-β signaling or a disruption in mechanosensory properties, respectively. In this Review, we highlight a unique family of vertebrate-specific short-lived filamin-binding proteins, the refilins (refilin-A and refilin-B), that modulate filamin-dependent actin crosslinking properties. Refilins are downstream TGF-β effectors in epithelial cells. Double knockout of both refilin-A and refilin-B in mice results in precocious ossification of some axial skeletal elements, leading to malformations that are similar to those seen in FLNB-deficient mice. Based on these findings, we present a model summarizing the role of refilins in regulating the mechanosensory functions of FLNB during skeletal development. We also discuss the possible contribution of refilins to FLNB-related skeletal pathologies that are associated with gain-of-function mutations.

PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations

Molecular dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth has led to the production of a large amount of simulation data that need to be filtered for relevant information to address specific biomedical and biochemical questions. One of the most relevant molecular properties that can be investigated by all-atom MD simulations is the time-dependent evolution of the complex noncovalent interaction networks governing such fundamental aspects as molecular recognition, binding strength, and mechanical and structural stability. Extracting, evaluating, and visualizing noncovalent interactions is a key task in the daily work of structural biologists. We have developed PyContact, an easy-to-use, highly flexible, and intuitive graphical user interface-based application, designed to provide a toolkit to investigate biomolecular interactions in MD trajectories. PyContact is designed to facilitate this task by enabling identification of relevant noncovalent interactions in a comprehensible manner. The implementation of PyContact as a standalone application enables rapid analysis and data visualization without any additional programming requirements, and also preserves full in-program customization and extension capabilities for advanced users. The statistical analysis representation is interactively combined with full mapping of the results on the molecular system through the synergistic connection between PyContact and VMD. We showcase the capabilities and scientific significance of PyContact by analyzing and visualizing in great detail the noncovalent interactions underlying the ion permeation pathway of the human P2X3 receptor. As a second application, we examine the protein-protein interaction network of the mechanically ultrastable cohesin-dockering complex.

Combining in Vitro and in Silico Single-Molecule Force Spectroscopy to Characterize and Tune Cellulosomal Scaffoldin Mechanics

Cellulosomes are polyprotein machineries that efficiently degrade cellulosic material. Crucial to their function are scaffolds consisting of highly homologous cohesin domains, which serve a dual role by coordinating a multiplicity of enzymes as well as anchoring the microbe to its substrate. Here we combined two approaches to elucidate the mechanical properties of the main scaffold ScaA of Acetivibrio cellulolyticus. A newly developed parallelized one-pot in vitro transcription-translation and protein pull-down protocol enabled high-throughput atomic force microscopy (AFM)-based single-molecule force spectroscopy (SMFS) measurements of all cohesins from ScaA with a single cantilever, thus promising improved relative force comparability. Albeit very similar in sequence, the hanging cohesins showed considerably lower unfolding forces than the bridging cohesins, which are subjected to force when the microbe is anchored to its substrate. Additionally, all-atom steered molecular dynamics (SMD) simulations on homology models offered insight into the process of cohesin unfolding under force. Based on the differences among the individual force propagation pathways and their associated correlation communities, we designed mutants to tune the mechanical stability of the weakest hanging cohesin. The proposed mutants were tested in a second high-throughput AFM SMFS experiment revealing that in one case a single alanine to glycine point mutation suffices to more than double the mechanical stability. In summary, we have successfully characterized the force induced unfolding behavior of all cohesins from the scaffoldin ScaA, as well as revealed how small changes in sequence can have large effects on force resilience in cohesin domains. Our strategy provides an efficient way to test and improve the mechanical integrity of protein domains in general.