Phosphorylation of Fibronectin Influences the Structural Stability of the Predicted Interchain Domain

Molecular interplay of ADAMTS13-MDTCS and von willebrand Factor-A2: deepened insights from extensive atomistic simulations

Thrombotic thrombocytopenic purpura (TTP) is a rare and life-threatening disease. One hallmark is severe ADAMTS13 deficiency, causing ultra-large von Willebrand factor (VWF) multimers to accumulate, leading to microthrombi and lastly to microangiopathic hemolytic anemia and severe thrombocytopenia. Despite great success in recent decades, the molecular picture of the interaction between VWF and ADAMTS13 remains vague. Here, we utilized modern replica-exchange molecular dynamics simulations with the TIGER2h method to sample a vast configurational space of the isolated ADAMTS13-MDTCS domains and the exposure to its substrate and activating cofactor - the unraveled VWF-A2 domain. The sampling of binding sites and conformations was guided and filtered in agreement with available experimental evidence. We provide comprehensive information on exosites for each domain and direct pairs of interacting amino acids, for the first time. The major binding cluster for the active site of the MP domain contrasts the previous mapping of VWF-A2 residues and reciprocal binding pockets. Two major binding modes are revealed and provide access to conformational changes of an extended gatekeeper tetrad upon overcoming local latency during substrate binding and to a dedicated recruitment mechanism. Our work adds the first molecular interaction model that places previous experimental results in perspective to better understand disease-related mutations towards improved therapies. Numerous empirical targets are proposed to verify the given binding modes, to refine the overall picture of MP binding pockets, the role of Dis binding in MP activation and the passage of the Cys-rich domain through VWF-A2, thus deepening the understanding of a highly dynamic interplay.Communicated by Ramaswamy H. Sarma.

Fundamental Redesign of the TIGER2hs Kernel to Address Severe Parameter Sensitivity

Replica exchange molecular dynamics simulations are one of the most popular approaches to enhance conformational sampling of molecular systems. Applications range from protein folding to protein-protein or other host-guest interactions, as well as binding free energy calculations. While these methods are computationally expensive, highly accurate results can be obtained. We recently developed TIGER2hs, an improved version of the temperature intervals with global exchange of replicas (TIGER2) algorithm. This method combines the replica-based enhanced sampling in an explicit solvent with a hybrid solvent energy evaluation. During the exchange attempts, bulk water is replaced by an implicit solvent model, allowing sampling with significantly less replicas than parallel tempering (REMD). This enables accurate enhanced sampling calculations with only a fraction of computational resources compared to REMD. Our latest results highlight several issues with sampling imbalance and parameter sensitivity within the original TIGER2 exchange algorithms that affect the overall state populations. A high sensitivity on replica number and maximum temperature is eliminated by changing to a pairwise exchange kernel (PE) without additional sorting. Simulations are controlled by adjusting the average temperature change per exchange ⟨ΔT/χ⟩ to below 30 K to mimic a controlled temperature mixing of replicas similar to REMD. Thus, this parameter provides an applicable property for selecting combinations of replica number and maximum temperature to adjust simulations for best accuracy, with flexible resource investment. This increases the robustness of the method and ensures results in excellent agreement with REMD, as demonstrated for three different peptides.

Structure-Property Relationships of Elastin-like Polypeptides: A Review of Experimental and Computational Studies

Elastin is a structural protein with outstanding mechanical properties (e.g., elasticity and resilience) and biologically relevant functions (e.g., triggering responses like cell adhesion or chemotaxis). It is formed from its precursor tropoelastin, a 60-72 kDa water-soluble and temperature-responsive protein that coacervates at physiological temperature, undergoing a phenomenon termed lower critical solution temperature (LCST). Inspired by this behavior, many scientists and engineers are developing recombinantly produced elastin-inspired biopolymers, usually termed elastin-like polypeptides (ELPs). These ELPs are generally comprised of repetitive motifs with the sequence VPGXG, which corresponds to repeats of a small part of the tropoelastin sequence, X being any amino acid except proline. ELPs display LCST and mechanical properties similar to tropoelastin, which renders them promising candidates for the development of elastic and stimuli-responsive protein-based materials. Unveiling the structure-property relationships of ELPs can aid in the development of these materials by establishing the connections between the ELP amino acid sequence and the macroscopic properties of the materials. Here we present a review of the structure-property relationships of ELPs and ELP-based materials, with a focus on LCST and mechanical properties and how experimental and computational studies have aided in their understanding.

Unveiling the N-Terminal Homodimerization of BCL11B by Hybrid Solvent Replica-Exchange Simulations

Transcription factors play a crucial role in regulating biological processes such as cell growth, differentiation, organ development and cellular signaling. Within this group, proteins equipped with zinc finger motifs (ZFs) represent the largest family of sequence-specific DNA-binding transcription regulators. Numerous studies have proven the fundamental role of BCL11B for a variety of tissues and organs such as central nervous system, T cells, skin, teeth, and mammary glands. In a previous work we identified a novel atypical zinc finger domain (CCHC-ZF) which serves as a dimerization interface of BCL11B. This domain and formation of the dimer were shown to be critically important for efficient regulation of the BCL11B target genes and could therefore represent a promising target for novel drug therapies. Here, we report the structural basis for BCL11B-BCL11B interaction mediated by the N-terminal ZF domain. By combining structure prediction algorithms, enhanced sampling molecular dynamics and fluorescence resonance energy transfer (FRET) approaches, we identified amino acid residues indispensable for the formation of the single ZF domain and directly involved in forming the dimer interface. These findings not only provide deep insight into how BCL11B acquires its active structure but also represent an important step towards rational design or selection of potential inhibitors.

Tethering of the IgG1 Antibody to Amorphous Silica for Immunosensor Development: A Molecular Dynamics Study

A key factor for improving the sensitivity and performance of immunosensors based on mechanical-plasmonic methods is the orientation of the antibody proteins immobilized on the inorganic surface. Although experimental techniques fail to determine surface phenomena at the molecular level, modern simulations open the possibility for improving our understanding of protein-surface interactions. In this work, replica exchange molecular dynamics (REMD) simulations have been used to model the IgG1 protein tethered onto the amorphous silica surface by considering a united-atom model and a relatively large system (2500 nm² surface). Additional molecular dynamics (MD) simulations have been conducted to derive an atomistic model for the amorphous silica surface using the cristobalite crystal structure as a starting point and to examine the structure of the free IgG1 antibody in the solution for comparison when immobilized. Analyses of the trajectories obtained for the tethered IgG1, which was sampled considering 32 different temperatures, have been used to define the geometry of the protein with respect to the inorganic surface. The tilt angle of the protein with respect to the surface plane increases with temperature, the most populated values being 24, 66, and 87° at the lowest (250 K), room (298 K), and the highest (380 K) temperatures. This variation indicates that the importance of protein-surface interactions decreases with increasing temperature. The influence of the surface on the structure of the antibody is very significant in the constant region, which is directly involved in the tethering process, while it is relatively unimportant for the antigen-binding fragments, which are farthest from the surface. These results are expected to contribute to the development of improved mechanical-plasmonic sensor microarrays in the near future.