Deciphering collaborative sidechain motions in proteins during molecular dynamics simulations

Molecular mechanisms for stabilizing biologics in the solid state

Protein drugs exhibit challenges of biophysical and biochemical instability due to their structural complexity and rich dynamics. Solid-state biologics aim to enhance stability by increasing molecular rigidity within the formulation matrix, representing a primary category of drug products alongside sterile liquid formulations. Understanding the molecular mechanisms behind the stabilization and destabilization of protein drugs, influenced by formulation composition and drying processes, provides scientific rationale for drug product design. This review aims to elaborate on the two primary models of water-to-sugar substitution and matrix vitrification, respectively, via thermodynamic and kinetic stabilization. It offers an up-to-date review of experimental investigations into these hypotheses, specifically elucidating protein structure and protein-excipient interactions at the molecular level, molecular dynamics across a broad range of motion regimes, and microscopic attributes such as protein-sugar and protein-salt miscibility and microenvironmental acidity, in relevant liquid, frozen, and solid states, using advanced biophysical techniques for solid-state analysis. Moreover, we discuss how these mechanistic understandings facilitate the investigation and prediction of critical stability behaviors and enables the design of solid biological drug products.

Mechanical stress and anionic lipids synergistically stabilize an atypical structure of the angiotensin II type 1 receptor (AT1)

Environmental factors, including mechanical stress and surrounding lipids, can influence the response of GPCRs, such as the mechanosensitive angiotensin II type 1 receptor (AT1). To investigate the impact of these factors on AT1 activation, we developed a steered molecular dynamics simulations protocol based on quaternion formalism. In this protocol, a pulling force was applied to the N-terminus of transmembrane helix 6 (TM6) to induce the TM6 opening characteristic of activation. Subsequently, the simulations were continued without constraints to allow the receptor to relax around the novel TM6 conformation under different conditions. We analyzed the responses of AT1 to membrane stretching, modeled by applying surface tension, in different bilayers. In phosphocholine bilayers without surface tension, we could observe a transient atypical structure of AT1, with an outward TM7 conformation, at the beginning of the activation process. This atypical structure then evolved toward a pre-active structure with outward TM6 and inward TM7. Strikingly, the presence of anionic phosphoglycerol lipids and application of surface tension synergistically favored the atypical structure, which led to an increase in the cross-section area of the receptor intracellular domain. Lipid internalization and H-bonds between lipid heads and the receptor C-terminus increased in phosphoglycerol vs phosphocholine bilayers, but did not depend on surface tension. The difference in the cross-section area of the atypical and pre-active conformations makes the conformational transition sensitive to lateral pressure, and favors the atypical conformation upon surface tension. Anionic lipids act as allosteric modulators of the conformational transition, by stabilizing the atypical conformation. These findings contribute to decipher the mechanisms underlying AT1 activation, highlighting the influence of environmental factors on GPCR responses. Moreover, our results reveal the existence of intermediary conformations that depend on receptor environment and could be targeted in drug design efforts.

Identifying and Assessing Putative Allosteric Sites and Modulators for CXCR4 Predicted through Network Modeling and Site Identification by Ligand Competitive Saturation

The chemokine receptor CXCR4 is a critical target for the treatment of several cancer types and HIV-1 infections. While orthosteric and allosteric modulators have been developed targeting its extracellular or transmembrane regions, the intramembrane region of CXCR4 may also include allosteric binding sites suitable for the development of allosteric drugs. To investigate this, we apply the Gaussian Network Model (GNM) to the monomeric and dimeric forms of CXCR4 to identify residues essential for its local and global motions located in the hinge regions of the protein. Residue interaction network (RIN) analysis suggests hub residues that participate in allosteric communication throughout the receptor. Mutual residues from the network models reside in regions with a high capacity to alter receptor dynamics upon ligand binding. We then investigate the druggability of these potential allosteric regions using the site identification by ligand competitive saturation (SILCS) approach, revealing two putative allosteric sites on the monomer and three on the homodimer. Two screening campaigns with Glide and SILCS-Monte Carlo docking using FDA-approved drugs suggest 20 putative hit compounds including antifungal drugs, anticancer agents, HIV protease inhibitors, and antimalarial drugs. In vitro assays considering mAB 12G5 and CXCL12 demonstrate both positive and negative allosteric activities of these compounds, supporting our computational approach. However, in vivo functional assays based on the recruitment of β-arrestin to CXCR4 do not show significant agonism and antagonism at a single compound concentration. The present computational pipeline brings a new perspective to computer-aided drug design by combining conformational dynamics based on network analysis and cosolvent analysis based on the SILCS technology to identify putative allosteric binding sites using CXCR4 as a showcase.

Cyclic peptides discriminate BCL-2 and its clinical mutants from BCL-XL by engaging a single-residue discrepancy

Overexpressed pro-survival B-cell lymphoma-2 (BCL-2) family proteins BCL-2 and BCL-XL can render tumor cells malignant. Leukemia drug venetoclax is currently the only approved selective BCL-2 inhibitor. However, its application has led to an emergence of resistant mutations, calling for drugs with an innovative mechanism of action. Herein we present cyclic peptides (CPs) with nanomolar-level binding affinities to BCL-2 or BCL-XL, and further reveal the structural and functional mechanisms of how these CPs target two proteins in a fashion that is remarkably different from traditional small-molecule inhibitors. In addition, these CPs can bind to the venetoclax-resistant clinical BCL-2 mutants with similar affinities as to the wild-type protein. Furthermore, we identify a single-residue discrepancy between BCL-2 D111 and BCL-XL A104 as a molecular "switch" that can differently engage CPs. Our study suggests that CPs may inhibit BCL-2 or BCL-XL by delicately modulating protein-protein interactions, potentially benefiting the development of next-generation therapeutics.

Probing the competitive inhibitor efficacy of frog-skin alpha helical AMPs identified against ACE2 binding to SARS-CoV-2 S1 spike protein as therapeutic scaffold to prevent COVID-19

In COVID-19 infection, the SARS-CoV-2 spike protein S1 interacts to the ACE2 receptor of human host, instigating the viral infection. To examine the competitive inhibitor efficacy of broad spectrum alpha helical AMPs extracted from frog skin, a comparative study of intermolecular interactions between viral S1 and AMPs was performed relative to S1-ACE2p interactions. The ACE2 binding region with S1 was extracted as ACE2p from the complex for ease of computation. Surprisingly, the Spike-Dermaseptin-S9 complex had more intermolecular interactions than the other peptide complexes and importantly, the S1-ACE2p complex. We observed how atomic displacements in docked complexes impacted structural integrity of a receptor-binding domain in S1 through conformational sampling analysis. Notably, this geometry-based sampling approach confers the robust interactions that endure in S1-Dermaseptin-S9 complex, demonstrating its conformational transition. Additionally, QM calculations revealed that the global hardness to resist chemical perturbations was found more in Dermaseptin-S9 compared to ACE2p. Moreover, the conventional MD through PCA and the torsional angle analyses indicated that Dermaseptin-S9 altered the conformations of S1 considerably. Our analysis further revealed the high structural stability of S1-Dermaseptin-S9 complex and particularly, the trajectory analysis of the secondary structural elements established the alpha helical conformations to be retained in S1-Dermaseptin-S9 complex, as substantiated by SMD results. In conclusion, the functional dynamics proved to be significant for viral Spike S1 and Dermaseptin-S9 peptide when compared to ACE2p complex. Hence, Dermaseptin-S9 peptide inhibitor could be a strong candidate for therapeutic scaffold to prevent infection of SARS-CoV-2.

Evolutionary information helps understand distinctive features of the angiotensin II receptors AT1 and AT2 in amniota

In vertebrates, the octopeptide angiotensin II (AngII) is an important in vivo regulator of the cardiovascular system. It acts mainly through two G protein-coupled receptors, AT1 and AT2. To better understand distinctive features of these receptors, we carried out a phylogenetic analysis that revealed a mirror evolution of AT1 and AT2, each one split into two clades, separating fish from terrestrial receptors. It also revealed that hallmark mutations occurred at, or near, the sodium binding site in both AT1 and AT2. Electrostatics computations and molecular dynamics simulations support maintained sodium binding to human AT1 with slow ingress from the extracellular side and an electrostatic component of the binding free energy around -3kT, to be compared to around -2kT for human AT2 and the δ opioid receptor. Comparison of the sodium binding modes in wild type and mutated AT1 and AT2 from humans and eels indicates that the allosteric control by sodium in both AT1 and AT2 evolved during the transition from fish to amniota. The unusual S7.46N mutation in AT1 is mirrored by a L3.36M mutation in AT2. In the presence of sodium, the N7.46 pattern in amniota AT1 stabilizes the inward orientation of N3.35 in the apo receptor, which should contribute to efficient N3.35 driven biased signaling. The M3.36 pattern in amniota AT2 favours the outward orientation of N3.35 and the receptor promiscuity. Both mutations have physiological consequences for the regulation of the renin-angiotensin system.

The analysis of protein sequences from different species can reveal interesting trends in the structural and functional evolution of a protein family. Here, we analyze the evolution of two G protein-coupled receptors, AT1 and AT2, which bind the angiotensin II peptide and are important regulators of the cardiovascular system. We show that these receptors underwent a mirror evolution. Specific mutations at, or near, the sodium binding pocket occurred in both AT1 and AT2 during the transition to terrestrial life. We carried out electrostatics computations and molecular dynamics simulations to decipher the details of the sodium binding mode in eel and human receptors, as prototypes of fish and amniota receptors. Our results indicate that sodium binding is kinetically slow but thermodynamically stable. Comparison of the sodium binding modes in eel and human receptors reveals that an unusual mutation in the sodium binding pocket of AT1 is critical for biased signaling of amniota AT1 whereas a mutation in AT2 promotes promiscuity of amniota AT2. In turn, these data indicate that a few mutations at a strategic position (here the sodium binding pocket) are an efficient way to gain functional evolution.

Experimental diagnostic of sequence-variant dynamic perturbations revealed by broadband dielectric spectroscopy

Genetic diversity leads to protein robustness, adaptability, and failure. Some sequence variants are structurally robust but functionally disturbed because mutations bring the protein onto unfolding/refolding routes resulting in misfolding diseases (e.g., Parkinson). We assume dynamic perturbations introduced by mutations foster the alternative unfolding routes and test this possibility by comparing the unfolding dynamics of the heat-labile enterotoxin B pentamers and the cholera toxin B pentamers, two pentamers structurally and functionally related and robust to 17 sequence variations. The B-subunit thermal unfolding dynamics are monitored by broadband dielectric spectroscopy in nanoconfined and weakly hydrated conditions. Distinct dielectric signals reveal the different B-subunits unfolding dynamics. Combined with network analyses, the experiments pinpoint the role of three mutations A1T, E7D, and E102A, in diverting LTB₅ to alternative unfolding routes that protect LTB₅ from dissociation. Altogether, the methodology diagnoses dynamics faults that may underlie functional disorder, drug resistance, or higher virulence of sequence variants.

Charge-Free, Stabilizing Amide-π Interactions Can Be Used to Control Collagen Triple-Helix Self-Assembly

There is a noted lack of understood, controllable interactions for directing the organization of collagen triple helices. While the field has had success using charge-pair interactions and cation-π interactions in helix design, these alone are not adequate for achieving the degree of specificity desirable for these supramolecular structures. Furthermore, because of the reliance on electrostatic interactions, designed heterotrimeric systems have been heavily charged, a property undesirable in some applications. Amide-π interactions are a comparatively understudied class of charge-free interactions, which could potentially be harnessed for triple-helix design. Herein, we propose, validate, and utilize pairwise amino acid amide-π interactions in collagen triple-helix design. Glutamine-phenylalanine pairs, when arranged in an axial geometry, are found to exhibit a moderately stabilizing effect, while in the lateral geometry, this pair is destabilizing. Together this allows glutamine-phenylalanine pairs to effectively set the register of triple helices. In contrast, interactions between asparagine and phenylalanine appear to have little effect on triple-helical stability. After deconvoluting the contributions of these amino acids to triple-helix stability, we demonstrate these new glutamine-phenylalanine interactions in the successful design of a heterotrimeric triple helix. The results of all of these analyses are used to update our collagen triple-helix thermal stability prediction algorithm, Scoring function for Collagen Emulating Peptides' Temperature of Transition (SCEPTTr).

Bios2cor: an R package integrating dynamic and evolutionary correlations to identify functionally important residues in proteins

SUMMARY

Both dynamic correlations in protein sidechain motions during molecular dynamics (MD) simulations and evolutionary correlations in multiple sequence alignments (MSA) of homologous proteins may reveal functionally important residues. We developed the R package Bios2cor that provides a unique framework to investigate and, possibly, integrate both analyses. Bios2cor starts with an MSA or a MD trajectory and computes correlation/covariation scores between positions in the MSA or between sidechain dihedral angles or rotamers in the MD trajectory. In addition, Bios2cor provides a variety of tools for the analysis, the visualization and the interpretation of the data.

AVAILABILITY

The R package Bios2cor is available from the Comprehensive R Archive Network, at http://cran.r-project.org/ web/packages/Bios2cor/index.html.