Deciphering S100B Allosteric Signaling: The Role of a Peptide Target, TRTK-12, as an Ensemble Modulator.

Allostery is an essential biological phenomenon in which perturbation at one site in a biomolecule elicits a functional response at a distal location(s). It is integral to biological processes, such as cellular signaling, metabolism, and transcription regulation. Understanding allostery is also crucial for rational drug discovery. In this work, we focus on an allosteric S100B protein that belongs to the S100 class of EF-hand Ca2+-binding proteins. The Ca2+-binding affinity of S100B is modulated allosterically by TRTK-12 peptide binding 25 Å away from the Ca2+-binding site. We investigated S100B allostery by carrying out nuclear magnetic resonance (NMR) measurements along with microsecond-long molecular dynamics (MD) simulations on S100B/Ca2+ with/without TRTK-12 at different NaCl salt concentrations. NMR HSQC results show that TRTK-12 reorganizes how S100B/Ca2+ responds to different salt concentrations at both orthosteric and allosteric sites. The MD data suggest that TRTK-12 breaks the dynamic aromatic and hydrogen-bond interactions (not observed in X-ray crystallographic structures) between the hinge/helix and Ca2+-binding EF-hand loop of the two subunits in the homodimeric protein. This triggers rearrangement in the protein network architectures and leads to allosteric communication. Finally, computational studies of S100B at distinct ionic strengths suggest that ligand-bound species are more robust to the changing environment relative to the S100B/Ca2+ complex.

Opposing roles of organic salts on mini-protein structure.

We investigated the effects of 1-ethyl-3-methylimidazolium chloride ([EMIM][Cl]) and choline chloride ([Chol][Cl]) on the local environment and conformational landscapes of Trp-cage and Trpzip4 mini-proteins using experimental and computational approaches. Fluorescence experiments and computational simulations revealed distinct behaviors of the mini-proteins in the presence of these organic salts. [EMIM][Cl] showed a strong interaction with Trp-cage, leading to fluorescence quenching and destabilization of its native structural interactions. Conversely, [Chol][Cl] had a negligible impact on Trp-cage fluorescence at low concentrations but increased it at high concentrations, indicating a stabilizing role. Computational simulations elucidated that [EMIM][Cl] disrupted the hydrophobic core packing and decreased proline-aromatic residue contacts in Trp-cage, resulting in a more exposed environment for Trp residues. In contrast, [Chol][Cl] subtly influenced the hydrophobic core packing, creating a hydrophobic environment near the tryptophan residues. Circular dichroism experiments revealed that [Chol][Cl] stabilized the secondary structure of both mini-proteins, although computational simulations did not show significant changes in secondary content at the explored concentrations. The simulations also demonstrated a more rugged free energy landscape for Trp-cage and Trpzip4 in [EMIM][Cl], suggesting destabilization of the tertiary structure for Trp-cage and secondary structure for Trpzip4. Similar fluorescence trends were observed for Trpzip4, with [EMIM][Cl] quenching fluorescence and exhibiting stronger interaction, while [Chol][Cl] increased the fluorescence at high concentrations. These findings highlight the interplay between [EMIM][Cl] and [Chol][Cl] with the mini-proteins and provide a detailed molecular-level understanding of how these organic salts impact the nearby surroundings and structural variations. Understanding such interactions is valuable for diverse applications, from biochemistry to materials science.

Modulation of Aß 16-22 aggregation by glucose.

The self-assembly of amyloid-beta (Aß) peptides into fibrillar structures in the brain is a signature of Alzheimer's disease. Recent studies have reported correlations between Alzheimer's disease and type-2 diabetes. Structurally, hyperglycemia induces covalent protein crosslinkings by advanced glycation end products (AGE), which can affect the stability of Aß oligomers. In this work, we leverage physics-based coarse-grained molecular simulations to probe alternate thermodynamic pathways that affect peptide aggregation propensities at varying concentrations of glucose molecules. Similar to previous experimental reports, our simulations show a glucose concentration-dependent increase in Aß aggregation rates, without changes in the overall secondary structure content. We discovered that glucose molecules prefer partitioning onto the aggregate-water interface at a specific orientation, resulting in a loss of molecular rotational entropy. This effectively hastens the aggregation rates, as peptide self-assembly can reduce the available surface area for peptide-glucose interactions. This work introduces a new thermodynamic-driven pathway, beyond chemical cross-linking, that can modulate Aß aggregation.

Folding and modulation of the helical conformation of Glycophorin A by point mutations.

Transmembrane helix folding and self-association play important roles in biological signaling and transportation pathways across biomembranes. With molecular simulations, studies to explore the structural biochemistry of this process have been limited to focusing on individual fragments of this process - either helix formation or dimerization. While at an atomistic resolution, it can be prohibitive to access long spatio-temporal scales, at the coarse grained (CG) level, current methods either employ additional constraints to prevent spontaneous unfolding or have a low resolution on sidechain beads that restricts the study of dimer disruption caused by mutations. To address these research gaps, in this work, we apply our recent, in-house developed CG model (ProMPT) to study the folding and dimerization of Glycophorin A (GpA) and its mutants in the presence of Dodecyl-phosphocholine (DPC) micelles. Our results first validate the two-stage model that folding and dimerization are independent events for transmembrane helices and found a positive correlation between helix folding and DPC-peptide contacts. The wild type (WT) GpA is observed to be a right-handed dimer with specific GxxxG contacts, which agrees with experimental findings. Specific point mutations reveal several features responsible for the structural stability of GpA. While the T87L mutant forms anti-parallel dimers due to an absence of T87 interhelical hydrogen bonds, a slight loss in helicity and a hinge-like feature at the GxxxG region develops for the G79L mutant. We note that the local changes in the hydrophobic environment, affected by the point mutation, contribute to the development of this helical bend. This work presents a holistic overview of the structural stability of GpA in a micellar environment, while taking secondary structural fluctuations into account. Moreover, it presents opportunities for applications of computationally efficient CG models to study conformational alterations of transmembrane proteins that have physiological relevance.

Emergence of allostery through reorganization of protein residue network architecture.

Despite more than a century of study, consensus on the molecular basis of allostery remains elusive. A comparison of allosteric and non-allosteric members of a protein family can shed light on this important regulatory mechanism, and the bacterial biotin protein ligases, which catalyze post-translational biotin addition, provide an ideal system for such comparison. While the Class I bacterial ligases only function as enzymes, the bifunctional Class II ligases use the same structural architecture for an additional transcription repression function. This additional function depends on allosterically activated homodimerization followed by DNA binding. In this work, we used experimental, computational network, and bioinformatics analyses to uncover distinguishing features that enable allostery in the Class II biotin protein ligases. Experimental studies of the Class II Escherichia coli protein indicate that catalytic site residues are critical for both catalysis and allostery. However, allostery also depends on amino acids that are more broadly distributed throughout the protein structure. Energy-based community network analysis of representative Class I and Class II proteins reveals distinct residue community architectures, interactions among the communities, and responses of the network to allosteric effector binding. Bioinformatics mutual information analyses of multiple sequence alignments indicate distinct networks of coevolving residues in the two protein families. The results support the role of divergent local residue community network structures both inside and outside of the conserved enzyme active site combined with distinct inter-community interactions as keys to the emergence of allostery in the Class II biotin protein ligases.

Mechanistic insights into the inhibition of amyloid-ß aggregation by chitosan.

Neurodegeneration related to Alzheimer's disease has long been linked to the accumulation of abnormal aggregates of amyloid-ß (Aß) peptides. Pre-fibrillar oligomeric intermediates of Aß aggregation are considered the primary drivers of neurotoxicity, however, their targetting remains an unresolved challenge. In response, the effects of macromolecular components of the blood-brain barrier, artificial extracellular matrix mimics, and polymeric drug delivery particles, on the aggregation of Aß peptides are gaining interest. Multiple experimental studies have demonstrated the potential of one such macromolecule, chitosan (CHT) - a polysaccharide with acid induced cationicity (pKa 6.5) - to inhibit the aggregation of Aß, and reduce the associated neurotoxic effects. However, the mechanistic details of this inhibitory action, and the structural details of the emergent Aß complexes are not understood. In this work, we probed how CHT modulated the aggregation of Aß's central hydrophobic core fragment, K16LVFFAE22, using coarse-grained molecular dynamics simulations. CHT was found to bind and sequester Aß peptides, thus limiting their ultimate aggregation numbers. The intensity of this inhibitory action was enhanced by CHT concentration, as well as CHT's pH-dependent degree of cationicity, corroborating experimental observations. Furthermore, CHT was found to reshape the conformational landscapes of Aß peptides, enriching collapsed peptides at near-physiological conditions of pH 7.5, and extended peptides at slightly acidic conditions of pH 6.5, where the charge profile of K16LVFFAE22 peptides remained unchanged. These conformational changes were limited to peptides in direct contact in CHT, thus emphasizing the influence of local environments on Aß conformations. These findings add to basic knowledge of the aggregation behaviour of Aß peptides, and could potentially guide the development of advanced CHT-based materials for the treatment of Alzheimer's disease.

Discovery and characterization of a functional scFv for CCR2 inhibition via an extracellular loop.

The chemokine receptor CCR2 plays a key role in cellular migration and inflammatory processes. While tremendous progress has been made in elucidating CCR2 function and inhibition, the majority of approaches target its N-terminal domain and less is known about the function of the remaining extracellular loops and their potential as targets. Here, we used phage display to identify an antibody-derived scFv (single chain variable fragment) clone that specifically targets the second extracellular epitope of CCR2 (ECL2) for inhibition. Using in silico molecular docking, we identified six potential primary binding conformations of the novel scFv to the specified CCR2 epitope. In silico molecular dynamic analysis was used to determine conformational stability and identify protein-protein interactions. Umbrella sampling of a range of configurations with incrementally increasing separation of scFv and target generated by force pulling simulations was used to calculate binding energies. Downstream characterization by ELISA showed high binding affinity of the ECL2-scFv to CCR2. Furthermore, we showed that blocking the second extracellular loop inhibits macrophage migration and polarized macrophages towards M1 inflammatory cytokine production as potently as lipopolysaccharide (LPS). These studies highlight the applicability of epitope-specific targeting, emphasize the importance of in silico predictive modeling, and warrant further investigation into the role of the remaining epitopes of CCR2.

Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery.

EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 µM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 µM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.

Transferable and Polarizable Coarse Grained Model for Proteins─ProMPT.

The application of classical molecular dynamics (MD) simulations at atomic resolution (fine-grained level, FG), to most biomolecular processes, remains limited because of the associated computational complexity of representing all the atoms. This problem is magnified in the presence of protein-based biomolecular systems that have a very large conformational space, and MD simulations with fine-grained resolution have slow dynamics to explore this space. Current transferable coarse grained (CG) force fields in literature are either limited to only peptides with the environment encoded in an implicit form or cannot capture transitions into secondary/tertiary peptide structures from a primary sequence of amino acids. In this work, we present a transferable CG force field with an explicit representation of the environment for accurate simulations with proteins. The force field consists of a set of pseudoatoms representing different chemical groups that can be joined/associated together to create different biomolecular systems. This preserves the transferability of the force field to multiple environments and simulation conditions. We have added electronic polarization that can respond to environmental heterogeneity/fluctuations and couple it to protein's structural transitions. The nonbonded interactions are parametrized with physics-based features such as solvation and partitioning free energies determined by thermodynamic calculations and matched with experiments and/or atomistic simulations. The bonded potentials are inferred from corresponding distributions in nonredundant protein structure databases. We present validations of the CG model with simulations of well-studied aqueous protein systems with specific protein fold types─Trp-cage, Trpzip4, villin, WW-domain, and ß-α-ß. We also explore the applications of the force field to study aqueous aggregation of Aß 16-22 peptides.

Recovery of enzyme structure and activity following rehydration from ionic liquid.

Long-term preservation of proteins at room temperature continues to be a major challenge. Towards using ionic liquids (ILs) to address this challenge, here we present a combination of experiments and simulations to investigate changes in lysozyme upon rehydration from IL mixtures using two imidazolium-based ILs (1-ethyl-3-methylimidazolium ethylsulfate, [EMIM][EtSO4] and 1-ethyl-3-methylimidazolium diethylphosphate, [EMIM][Et2PO4]). Various spectroscopic experiments and molecular dynamics simulations are performed to ascertain the structure and activity of lysozyme. Circular dichroism spectroscopy confirms that lysozyme maintains its secondary structure upon rehydration, even after 295 days. Increasing the IL concentration decreases the activity of lysozyme and is ultimately quenched at sufficiently high IL concentrations, but the rehydration of lysozyme from high IL concentrations completely restores its activity. Such rehydration occurs in the most common lysozyme activity assay, but without careful attention, this effect on the IL concentration can be overlooked. From simulations we observe occupation of [EMIM+] ions near the vicinity of the active site and the ligand-lysozyme complex is less stable in the presence of ILs, which results in the reduction of lysozyme activity. Upon rehydration, fast leaving of [EMIM+] is observed and the availability of active site is restored. In addition, suppression of structural fluctuations is also observed when in high IL concentrations, which also explains the decrease of activity. This structure suppression is recovered after undergoing rehydration. The return of native protein structure and activity indicates that after rehydration lysozyme returns to its original state. Our results also suggest a simple route to protein recovery following extended storage.

Effects of applied surface-tension on membrane-assisted Aß aggregation.

Accumulation of protein-based (Aß) aggregates on cellular membranes with varying structural properties is commonly recognized as the key step in Alzheimer's pathogenesis. But experimental and computational challenges have made this biophysical characterization difficult. In particular, studies connecting biological membrane organization and Aß aggregation are limited. While experiments have suggested that an increased membrane curvature results in faster Aß peptide aggregation in the context of Alzheimer's disease, a mechanistic explanation for this relation is missing. In this work, we are leveraging molecular simulations with a physics-based coarse grained model to address and understand the relationships between curved cellular membranes and aggregation of a model template peptide Aß 16-22. In agreement with experimental results, our simulations also suggest a positive correlation between increased peptide aggregation and membrane curvature. More curved membranes have higher lipid packing defects that engage peptide hydrophobic groups and promote faster diffusion leading to peptide fibrillar structures. In addition, we curated the effects of peptide aggregation on the membrane's structure and organization. Interfacial peptide aggregation results in heterogeneous headgroup-peptide interactions and an induced crowding effect at the lipid headgroup region, leading to a more ordered headgroup region and disordered lipid-tails at the membrane core. This work presents a mechanistic and morphological overview of the relationships between the biomembrane local structure and organization, and Aß peptide aggregation.

Dual mechanism of ionic liquid-induced protein unfolding.

Ionic liquids (ILs) are gaining attention as protein stabilizers and refolding additives. However, varying degrees of success with this approach motivates the need to better understand fundamental IL-protein interactions. A combination of experiment and simulation is used to investigate the thermal unfolding of lysozyme in the presence of two imidazolium-based ILs (1-ethyl-3-methylimidazolium ethylsulfate, [EMIM][EtSO4] and 1-ethyl-3-methylimidazolium diethylphosphate, [EMIM][Et2PO4]). Both ILs reduce lysozyme melting temperature Tm, but more gradually than strong denaturants. [EMIM][Et2PO4] lowers lysozyme Tm more readily than [EMIM][EtSO4], as well as requiring less energy to unfold the protein, as determined by the calorimetric enthalpy ΔH. Intrinsic fluorescence measurements indicate that both ILs bind to tryptophan residues in a dynamic mode, and furthermore, molecular dynamics simulations show a high density of [EMIM]+ near lysozyme's Trp62 residue. For both ILs approximately half of the [EMIM]+ cations near Trp62 show perfect alignment of their respective rings. The [EMIM]+ cations, having a "local" effect in binding to tryptophan, likely perturb a critically important Arg-Trp-Arg bridge through favorable π-π and cation-π interactions. Simulations show that the anions, [EtSO4]- and [Et2PO4]-, interact in a "global" manner with lysozyme, due to this protein's strong net positive charge. The anions also determine the local distribution of ions surrounding the protein. [Et2PO4]- is found to have a closer first coordination shell around the protein and stronger Coulomb interactions with lysozyme than [EtSO4]-, which could explain why the former anion is more destabilizing. Patching of ILs to the protein surface is also observed, suggesting there is no universal IL solvent for proteins, and highlighting the complexity of the IL-protein environment.

Microscopic Picture of Calcium-Assisted Lipid Demixing and Membrane Remodeling Using Multiscale Simulations.

The specificity of anionic phospholipids-calcium ion interaction and lipid demixing has been established as a key regulatory mechanism in several cellular signaling processes. The mechanism and implications of this calcium-assisted demixing have not been elucidated from a microscopic point of view. Here, we present an overview of atomic interactions between calcium and phospholipids that can drive nonideal mixing of lipid molecules in a model lipid bilayer composed of zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS)) lipids with computer simulations at multiple resolutions. Lipid nanodomain formation and growth were driven by calcium-enabled lipid bridging of the charged phosphatidylserine (PS) headgroups, which were favored against inter-POPS dipole interactions. Consistent with several experimental studies of calcium-associated membrane sculpting, our analyses also suggest modifications in local membrane curvature and cross-leaflet couplings as a response to such induced lateral heterogeneity. In addition, reverse mapping to a complementary atomistic description revealed structural insights in the presence of anionic nanodomains, at timescales not accessed by previous computational studies. This work bridges information across multiple scales to reveal a mechanistic picture of calcium ion's impact on membrane biophysics.

Tuning Allostery through Integration of Disorder to Order with a Residue Network.

In allostery, a signal from one site in a protein is transmitted to a second site to alter its function. Due to its ubiquity in biology and the potential for its exploitation in drug and protein design, the molecular basis of allosteric communication continues to be the subject of intense research. Although allosterically coupled sites are frequently characterized by disorder, how communication between disordered segments occurs remains obscure. Allosteric activation of Escherichia coli BirA dimerization occurs via coupled distant disorder-to-order transitions. In this work, combined structural and computational studies reveal an extensive residue network in BirA. Substitution of several network residues yields large perturbations to allostery. Force distribution analysis reveals that disruptions to the disorder-to-order transitions through amino acid substitution are manifested in shifts in the energy experienced by network residues as well as alterations in packing of an α-helix that plays a critical role in allostery. The combined results reveal a highly distributed allosteric mechanism that is robust to sequence change.

Computational insights into lipid assisted peptide misfolding and aggregation in neurodegeneration.

Peptide misfolding and aberrant assembly in membranous micro-environments have been associated with numerous neurodegenerative diseases. The biomolecular mechanisms and biophysical implications of these amyloid membrane interactions have been under extensive research and can assist in understanding disease pathogenesis and potential development of rational therapeutics. But, the complex nature and diversity of biomolecular interactions, structural transitions, and dependence on local environmental conditions have made accurate microscopic characterization challenging. In this review, using cases of Alzheimer's disease (amyloid-beta peptide), Parkinson's disease (alpha-synuclein peptide) and Huntington's disease (huntingtin protein), we illustrate existing challenges in experimental investigations and summarize recent relevant numerical simulation studies into amyloidogenic peptide-membrane interactions. In addition we project directions for future in silico studies and discuss shortcomings of current computational approaches.

Pathways of amyloid-beta absorption and aggregation in a membranous environment.

Aggregation of misfolded oligomeric amyloid-beta (Aß) peptides on lipid membranes has been identified as a primary event in Alzheimer's pathogenesis. However, the structural and dynamical features of this membrane assisted Aß aggregation have not been well characterized. The microscopic characterization of dynamic molecular-level interactions in peptide aggregation pathways has been challenging both computationally and experimentally. In this work, we explore differential patterns of membrane-induced Aß 16-22 (K-L-V-F-F-A-E) aggregation from the microscopic perspective of molecular interactions. Physics-based coarse-grained molecular dynamics (CG-MD) simulations were employed to investigate the effect of lipid headgroup charge - zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine: POPC) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine: POPS) - on Aß 16-22 peptide aggregation. Our analyses present an extensive overview of multiple pathways for peptide absorption and biomechanical forces governing peptide folding and aggregation. In agreement with experimental observations, anionic POPS molecules promote extended configurations in Aß peptides that contribute towards faster emergence of ordered ß-sheet-rich peptide assemblies compared to POPC, suggesting faster fibrillation. In addition, lower cumulative rates of peptide aggregation in POPS due to higher peptide-lipid interactions and slower lipid diffusion result in multiple distinct ordered peptide aggregates that can serve as nucleation seeds for subsequent Aß aggregation. This study provides an in-silico assessment of experimentally observed aggregation patterns, presents new morphological insights and highlights the importance of lipid headgroup chemistry in modulating the peptide absorption and aggregation process.

Mesoporous Encapsulated Chiral Nanogold for Use in Enantioselective Reactions.

Although various nanomaterials have been designed for biocatalysis, few of them can accelerate chemical reactions with high selectivity and stereocontrol, which remains them from being perfect alternatives to nature enzymes. Herein, inspired by the natural enzymes, an enantioselective nanomaterial has been constructed, with gold nanoparticles (AuNPs) as active centers, chiral cysteine (Cys) as selectors for chiral recognition, and expanded mesoporous silica (EMSN) as a skeleton of the artificial enzyme. In the oxidation of chiral 3,4-dihydroxy-phenylalanine (DOPA), the nanozyme with d-Cys shows preference to l-DOPA while the artificial enzyme with l-Cys shows preference to d-DOPA. Subsequent calculation of apparent steady-state kinetic parameters and activation energies together with molecular dynamics (MD) simulations showed that the different affinity precipitated by hydrogen bonding formation between chiral Cys and DOPA is the origin of chiral selectivity.

How Hydrophobic Hydration Destabilizes Surfactant Micelles at Low Temperature: A Coarse-Grained Simulation Study.

Micelles are self-assembled aggregates of amphiphilic surfactant molecules that are important in a variety of applications, including drug delivery, detergency, and catalysis. It is known that the micellization process is driven by the same physiochemical forces that drive protein folding, aggregation, and biological membrane self-assembly. Nevertheless, the molecular details of how micelle stability changes in water at low temperature are not fully clear. We develop and use a coarse-grained model to investigate how the interplay between nonionic surfactants and the surrounding water at the nanoscale affects the stability of micelles at high and low temperatures. Simulations of preformed C12E5 micelles in explicit water at a range of temperatures reveal the existence of two distinct surfactant conformations within the micelle, a bent structure and an extended structure, the latter being more prevalent at low temperature. Favorable interactions of the surfactant with more ordered solvation water stabilizes the extended configuration, allowing nanoscale wetting of the dry, hydrophobic core of the micelle, leading to the micelle breaking. Taken together, our coarse-grained simulations unravel how energetic and structural changes of the surfactant and the surrounding water destabilize micelles at low temperature, which is a direct consequence of the weakened hydrophobicity. Our approach thus provides an effective mean for extracting the molecular-level changes during hydrophobicity-driven destabilization of molecular self-assembly, which is important in a wide range of fields, including biology, polymer science, and nanotechnology.