An Efficient Approach to the Accurate Prediction of Mutational Effects in Antigen Binding to the MHC1

The major histocompatibility complex (MHC) can recognize and bind to external peptides to generate effective immune responses by presenting the peptides to T cells. Therefore, understanding the binding modes of peptide-MHC complexes (pMHC) and predicting the binding affinity of pMHCs play a crucial role in the rational design of peptide vaccines. In this study, we employed molecular dynamics (MD) simulations and free energy calculations with an Alanine Scanning with Generalized Born and Interaction Entropy (ASGBIE) method to investigate the protein-peptide interaction between HLA-A*02:01 and the G9209 peptide derived from the melanoma antigen gp100. The energy contribution of individual residue was calculated using alanine scanning, and hotspots on both the MHC and the peptides were identified. Our study shows that the pMHC binding is dominated by the van der Waals interactions. Furthermore, we optimized the ASGBIE method, achieving a Pearson correlation coefficient of 0.91 between predicted and experimental binding affinity for mutated antigens. This represents a significant improvement over the conventional MM/GBSA method, which yields a Pearson correlation coefficient of 0.22. The computational protocol developed in this study can be applied to the computational screening of antigens for the MHC1 as well as other protein-peptide binding systems.

Targeting hemoglobin receptors IsdH and IsdB of Staphylococcus aureus with a single VHH antibody inhibits bacterial growth

Methicillin-resistant Staphylococcus aureus, or MRSA, is one of the major causative agents of hospital-acquired infections worldwide. Novel antimicrobial strategies efficient against antibiotic-resistant strains are necessary and not only against S. aureus. Among those, strategies that aim at blocking or dismantling proteins involved in the acquisition of essential nutrients, helping the bacteria to colonize the host, are intensively studied. A major route for S. aureus to acquire iron from the host organism is the Isd (iron surface determinant) system. In particular, the hemoglobin receptors IsdH and IsdB located on the surface of the bacterium are necessary to acquire the heme moiety containing iron, making them a plausible antibacterial target. Herein, we obtained an antibody of camelid origin that blocked heme acquisition. We determined that the antibody recognized the heme-binding pocket of both IsdH and IsdB with nanomolar order affinity through its second and third complementary-determining regions. The mechanism explaining the inhibition of acquisition of heme in vitro could be described as a competitive process in which the complementary-determining region 3 from the antibody blocked the acquisition of heme by the bacterial receptor. Moreover, this antibody markedly reduced the growth of three different pathogenic strains of MRSA. Collectively, our results highlight a mechanism for inhibiting nutrient uptake as an antibacterial strategy against MRSA.

Categorization of hotspots into three types - weak, moderate and strong to distinguish protein-protein versus protein-peptide interactions

Protein-protein and protein-peptide interactions (PPI and PPepI) belong to a similar category of interactions, yet seemingly subtle differences exist among them. To characterize differences between protein-protein (PP) and protein-peptide (PPep) interactions, we have focussed on two important classes of residues-hotspot and anchor residues. Using implicit solvation-based free energy calculations, a very large-scale alanine scanning has been performed on benchmark datasets, consisting of over 5700 interface residues. The differences in the two categories are more pronounced, if the data were divided into three distinct types, namely - weak hotspots (having binding free energy loss upon Ala mutation, ΔΔG, ∼2-10 kcal/mol), moderate hotspots (ΔΔG, ∼10-20 kcal/mol) and strong hotspots (ΔΔG ≥ ∼20 kcal/mol). The analysis suggests that for PPI, weak hotspots are predominantly populated by polar and hydrophobic residues. The distribution shifts towards charged and polar residues for moderate hotspot and charged residues (principally Arg) are overwhelmingly present in the strong hotspot. On the other hand, in the PPepI dataset, the distribution shifts from predominantly hydrophobic and polar (in the weak type) to almost similar preference for polar, hydrophobic and charged residues (in moderate type) and finally the charged residue (Arg) and Trp are mostly occupied in the strong type. The preferred anchor residues in both categories are Arg, Tyr and Leu, possessing bulky side chain and which also strike a delicate balance between side chain flexibility and rigidity. The present knowledge should aid in effective design of biologics, by augmentation or disruption of PPIs with peptides or peptidomimetics.Communicated by Ramaswamy H. Sarma.

Unraveling the impact of ORF3a Q57H mutation on SARS-CoV-2: insights from molecular dynamics

ORF3a is a conserved accessory protein of SARS-CoV-2, linked to viral infection and pathogenesis, with acquired mutations at various locations. Previous studies have shown that the occurrence of the Q57H mutation is higher in comparison to other positions in ORF3a. This mutation is known to induce conformational changes, yet the extent of structural alteration and its role in the viral adaptation process remain unknown. Here we performed molecular dynamics (MD) simulations of wt-ORF3a, Q57H, and Q57A mutants to analyze structural changes caused by mutations compared to the native protein. The MD analysis revealed that Q57H and Q57A mutants show significant structural changes in the dimer conformation than the wt-ORF3a. This dimer conformer narrows down the ion channel cavity, which reduces Na + or K + permeability leading to decrease the antigenic response that can help the virus to escape the host immune system. Non-bonding interaction analysis shows the Q57H mutant has more interacting residues, resulting in more stability within dimer conformation than the wt-ORF3a and Q57A. Moreover, both mutant dimers (Q57H and Q57A) form a novel salt-bridge interaction at the same position between A:Asp142 and B:Lys61, whereas such an interaction is absent in the wt-ORF3a dimer. We have also noticed that the TM3 domain's flexibility in Q57H is increased because of strong inter-domain interactions of TM1 and TM2 within the dimer conformation. These unusual interactions and flexibility of Q57H mutant can have significant impacts on the SARS-CoV-2 adaptations, virulence, transmission, and immune system evasion. Our findings are consistent with the previous experimental data and provided details information on the structural perturbation in ORF3a caused by mutations, which can help better understand the structural change at the molecular level as well as the reason for the high virulence properties of this variant.Communicated by Ramaswamy H. Sarma.

Epitope Mapping of Anti-Mouse CCR3 Monoclonal Antibodies (C3Mab-6 and C3Mab-7)

One of G protein-coupled receptors, CC chemokine receptor 3 (CCR3), is expressed in eosinophils, basophils, a subset of Th2 lymphocytes, mast cells, and airway epithelial cells. CCR3 levels in the serum of colorectal cancer patients are significantly higher than in control groups. Moreover, CCR3 is essential for recruiting eosinophils into the lung. Therefore, CCR3 is considered both a therapeutic target for colorectal cancer and allergic diseases. Previously, we established anti-mouse CCR3 (mCCR3) monoclonal antibodies (mAbs), C₃Mab-6 (rat IgG₁, kappa) and C₃Mab-7 (rat IgG₁, kappa), by immunizing a rat with an N-terminal peptide of mCCR3. These mAbs can be used in flow cytometry and enzyme-linked immunosorbent assays. In this study, we performed the epitope mapping of C₃Mab-6 and C₃Mab-7 using alanine scanning. The reactivity between these mAbs and point mutants of mCCR3 were analyzed using flow cytometry. The results indicated that Phe3, Asn4, Thr5, Asp6, Glu7, Lys9, Thr10, and Glu13 of mCCR3 are essential for C₃Mab-6 binding, whereas Phe15 and Glu16 are essential for C₃Mab-7 binding.

Effect of mutation of NS2B cofactor residues on Dengue 2 NS2B-NS3 protease complex - an insight to viral replication

Dengue fever is an endemic virus-borne disease that causes many severe ailments, including dengue hemorrhagic fever and dengue shock syndrome. NS2B-NS3 protease is present in all four strains of the dengue virus. NS2B-NS3 is a non-structural protein that performs three distinct functions: protease activity, helicase activity, and nucleoside triphosphatase activity. NS2B-NS3 pro-complex plays a crucial role in viral replication, and NS2B interacts with NS3 protease at a flat active site with an amino acid of the N-terminal region. NS2B acts as a cofactor for NS3 protease. In the current study, the conserved residues of NS2B were identified. Dengue virus-2 NS2B was mutated at the identified conserved amino acid region to investigate the role of NS2B on activation of NS3 pro. Molecular dynamics simulations were performed to investigate the mutated complex's changes in stability, conformation, and free energy. The EAG mutant complex exhibited more unstable conformation, less hydrogen bond formation, and high binding energy than wild type. This result suggests a vital role of E63, A65, G69 mutation in NS2B for the interruption of activation of the NS3.Communicated by Ramaswamy H. Sarma.

Epitope Mapping of an Anti-EpCAM Monoclonal Antibody (EpMab-37) Using the Alanine Scanning Method

The epithelial cell adhesion molecule (EpCAM) is a type I transmembrane glycoprotein, and plays critical roles in cell adhesion, proliferation, and tumorigenesis. EpCAM has been considered as a promising target for tumor diagnosis and therapy. Anti-EpCAM monoclonal antibodies (mAbs) have been developed for EpCAM-overexpressed tumors, and several clinical trials have demonstrated promising outcomes. We previously established an anti-EpCAM mAb, EpMab-37 (mouse IgG₁, kappa), using the Cell-Based Immunization and Screening method. EpMab-37 was revealed to recognize the conformational epitope of EpCAM. In this study, we determined the critical epitope of EpMab-37 by flow cytometry using the 1 × alanine scanning (1 × Ala-scan) and the 2 × alanine scanning (2 × Ala-scan) method. We first performed flow cytometry by 1 × Ala-scan using one alanine (or glycine)-substituted EpCAM mutants, which were expressed on Chinese hamster ovary-K1 cells, and found that the EpMab-37 did not recognize the R163A mutant of EpCAM. We next performed flow cytometry by 2 × Ala-scan using two alanine (or glycine) residues-substituted EpCAM mutants, and confirmed that EpMab-37 did not recognize R163A-including mutants of EpCAM. The results indicated that the critical binding epitope of EpMab-37 includes Arg163 of EpCAM.

Determination of the Binding Epitope of an Anti-Mouse CCR9 Monoclonal Antibody (C9Mab-24) Using the 1× Alanine and 2× Alanine-Substitution Method

C-C chemokine receptor 9 (CCR9) is a receptor for C-C-chemokine ligand 25 (CCL25). CCR9 is crucial in the chemotaxis of immune cells and inflammatory responses. Moreover, CCR9 is highly expressed in tumors, including several solid tumors and T-cell acute lymphoblastic leukemia. Several preclinical studies have shown that anti-CCR9 monoclonal antibodies (mAbs) exert antitumor activity. Therefore, CCR9 is an attractive target for tumor therapy. In this study, we conducted the epitope mapping of an anti-mouse CCR9 (mCCR9) mAb, C₉Mab-24 (rat IgG_2a, kappa), using the 1× alanine (1× Ala)- and 2× alanine (2× Ala)-substitution methods via enzyme-linked immunosorbent assay. We first performed the 1× Ala-substitution method using one alanine-substituted peptides of the mCCR9 N-terminus (amino acids 1-19). C₉Mab-24 did not recognize two peptides (F14A and F17A), indicating that Phe14 and Phe17 are critical for C₉Mab-24-binding to mCCR9. Furthermore, we conducted the 2× Ala-substitution method using two consecutive alanine-substituted peptides of the mCCR9 N-terminus, and showed that C₉Mab-24 did not react with four peptides (M13A-F14A, F14A-D15A, D16A-F17A, and F17A-S18A), indicating that _13-MFDDFS_-18 is involved in C₉Mab-24-binding to mCCR9. Overall, combining, the 1× Ala- or 2× Ala-scanning methods could be useful for understanding for target-antibody interaction.

Integrating deep mutational scanning and low-throughput mutagenesis data to predict the impact of amino acid variants

BACKGROUND

Evaluating the impact of amino acid variants has been a critical challenge for studying protein function and interpreting genomic data. High-throughput experimental methods like deep mutational scanning (DMS) can measure the effect of large numbers of variants in a target protein, but because DMS studies have not been performed on all proteins, researchers also model DMS data computationally to estimate variant impacts by predictors.

RESULTS

In this study, we extended a linear regression-based predictor to explore whether incorporating data from alanine scanning (AS), a widely used low-throughput mutagenesis method, would improve prediction results. To evaluate our model, we collected 146 AS datasets, mapping to 54 DMS datasets across 22 distinct proteins.

CONCLUSIONS

We show that improved model performance depends on the compatibility of the DMS and AS assays, and the scale of improvement is closely related to the correlation between DMS and AS results.

Insights into small molecule inhibitor bindings to PD-L1 with residue-specific binding free energy calculation

Targeting the immunological checkpoint PD-1/PD-L1 with antibodies has shown opportunities to improve cancer treatment in recent years. However, antibody therapy is a double-edged sword with high cost, low patient tolerance, lack of oral bioavailability, and a reaction to most solid tumors that prevents the adoption of antibodies. Advancement of small-molecule PD-1/PD-L1 inhibitors that could overwhelm these drawbacks is sluggish because of the poor pharmacodynamic properties and shallow pocket of the PD-1/PD-L1 binding interface. Recently, a number of compounds have been discovered to bind the PD-L1/PD-L1 dimer interface, providing an excellent alternative to inhibit the interaction between PD-1/PD-L1 and small molecules. Quantitative characterization of PD-L1 interactions with these inhibitors will advance the design of novel and efficient inhibitors in the future. Here, the binding free energies of 35 PD-L1 dimer inhibitors have been calculated using the alanine-scanning-interaction-entropy (AS-IE) method. Hotspot residues on PD-L1 and potential modification groups on the inhibitors were identified. The experimental results for the AS-IE method were better correlated than the classical MM/GBSA method. These results may set the stage for the design the more powerful PD-L1 inhibitors.Communicated by Ramaswamy H. Sarma.

Epitope Mapping of Anti-Mouse CCR3 Monoclonal Antibodies Using Flow Cytometry

The CC chemokine receptor 3 (CCR3) is a receptor for CC chemokines, including CCL5/RANTES, CCL7/MCP-3, and CCL11/eotaxin. CCR3 is expressed on the surface of eosinophils, basophils, a subset of Th2 lymphocytes, mast cells, and airway epithelial cells. CCR3 and its ligands are involved in airway hyperresponsiveness in allergic asthma, ocular allergies, and cancers. Therefore, CCR3 is an attractive target for those therapies. Previously, anti-mouse CCR3 (mCCR3) monoclonal antibodies (mAbs), C₃Mab-3 (rat IgG_2a, kappa), and C₃Mab-4 (rat IgG_2a, kappa) were developed using the Cell-Based Immunization and Screening (CBIS) method. In this study, the binding epitope of these mAbs was investigated using flow cytometry. A CCR3 extracellular domain-substituted mutant analysis showed that C₃Mab-3, C₃Mab-4, and a commercially available mAb (J073E5) recognized the N-terminal region (amino acids 1-38) of mCCR3. Next, alanine scanning was conducted in the N-terminal region. The results revealed that the Ala2, Phe3, Asn4, and Thr5 of mCCR3 are involved in C₃Mab-3 binding, whereas Ala2, Phe3, and Thr5 are essential to C₃Mab-4 binding, and Ala2 and Phe3 are crucial to J073E5 binding. These results reveal the involvement of the N-terminus of mCCR3 in the recognition of C₃Mab-3, C₃Mab-4, and J073E5.

[N-terminal truncation of prenyltransferase enhances the biosynthesis of prenylnaringenin]

8-prenylnaringenin (8-PN) is a potent estrogen with high medicinal values. It also serves as an important precursor for many prenylated flavonoids. Microbial synthesis of 8-PN is mainly hindered by the low catalytic activity of prenyltransferases (PTS) and insufficient supply of precursors. In this work, a SfN8DT-1 from Sophora flavescens was used to improve the efficiency of (2S)-naringenin prenylation. The predicted structure of SfN8DT-1 showed that its main body is comprised of 9 α-helices and 8 loops, along with a long side chain formed by nearly 120 amino acids. SfN8DT-1 mutants with different side-chain truncated were tested in Saccharomyces cerevisiae. A mutant expressing the truncated enzyme at K62 site, designated as SfND8T-1-t62, produced the highest 8-PN titer. Molecular docking of SfN8DT-1-t62 with (2S)-naringenin and dimethylallyl diphosphate (DMAPP) showed that K185 was a potentially crucial residue. Alanine scanning within a range of 0.5 nm around these two substrates showed that the mutant K185A may decrease its affinity to substrates, which also indicated K185 was a potentially critical residue. Besides, the mutant K185W enhanced the affinity to ligands implied by the simulated saturation mutation, while the saturated mutation of K185 showed a great decrease in 8-PN production, indicating K185 is vital for the activity of SfN8DT-1. Subsequently, overexpressing the key genes of Mevalonate (MVA) pathway further improved the titer of 8-PN to 31.31 mg/L, which indicated that DMAPP supply is also a limiting factor for 8-PN synthesis. Finally, 44.92 mg/L of 8-PN was produced in a 5 L bioreactor after 120 h, which is the highest 8-PN titer reported to date.

Cellular ATP Levels Determine the Stability of a Nucleotide Kinase

The energy currency of the cell ATP, is used by kinases to drive key cellular processes. However, the connection of cellular ATP abundance and protein stability is still under investigation. Using Fast Relaxation Imaging paired with alanine scanning and ATP depletion experiments, we study the nucleotide kinase (APSK) domain of 3'-phosphoadenosine-5'-phosphosulfate (PAPS) synthase, a marginally stable protein. Here, we show that the in-cell stability of the APSK is determined by ligand binding and directly connected to cellular ATP levels. The observed protein stability change for different ligand-bound states or under ATP-depleted conditions ranges from ΔGf 0 = -10.7 to +13.8 kJ/mol, which is remarkable since it exceeds changes measured previously, for example upon osmotic pressure, cellular stress or differentiation. The results have implications for protein stability during the catalytic cycle of APS kinase and suggest that the cellular ATP level functions as a global regulator of kinase activity.

Generation of a Novel High-Affinity Antibody Binding to PCSK9 Catalytic Domain with Slow Dissociation Rate by CDR-Grafting, Alanine Scanning and Saturated Site-Directed Mutagenesis for Favorably Treating Hypercholesterolemia

Inhibition of proprotein convertase subtilisin/kexin type 9 (PCSK9) has become an attractive therapeutic strategy for lowering low-density lipoprotein cholesterol (LDL-C). In this study, a novel high affinity humanized IgG1 mAb (named h5E12-L230G) targeting the catalytic domain of human PCSK9 (hPCSK9) was generated by using CDR-grafting, alanine-scanning mutagenesis, and saturated site-directed mutagenesis. The heavy-chain constant region of h5E12-L230G was modified to eliminate the cytotoxic effector functions and mitigate the heterogeneity. The biolayer interferometry (BLI) binding assay and molecular docking study revealed that h5E12-L230G binds to the catalytic domain of hPCSK9 with nanomolar affinity (K_D = 1.72 nM) and an extremely slow dissociation rate (k_off, 4.84 × 10⁻⁵ s⁻¹), which interprets its quite low binding energy (−54.97 kcal/mol) with hPCSK9. Additionally, h5E12-L230G elevated the levels of LDLR and enhanced the LDL-C uptake in HepG2 cells, as well as reducing the serum LDL-C and total cholesterol (TC) levels in hyperlipidemic mouse model with high potency comparable to the positive control alirocumab. Our data indicate that h5E12-L230G is a high-affinity anti-PCSK9 antibody candidate with an extremely slow dissociation rate for favorably treating hypercholesterolemia and relevant cardiovascular diseases.

HISNAPI: a bioinformatic tool for dynamic hot spot analysis in nucleic acid-protein interface with a case study

Protein-nucleic acid interactions play essential roles in many biological processes, such as transcription, replication and translation. In protein-nucleic acid interfaces, hotspot residues contribute the majority of binding affinity toward molecular recognition. Hotspot residues are commonly regarded as potential binding sites for compound molecules in drug design projects. The dynamic property is a considerable factor that affects the binding of ligands. Computational approaches have been developed to expedite the prediction of hotspot residues on protein-nucleic acid interfaces. However, existing approaches overlook hotspot dynamics, despite their essential role in protein function. Here, we report a web server named Hotspots In silico Scanning on Nucleic Acid and Protein Interface (HISNAPI) to analyze hotspot residue dynamics by integrating molecular dynamics simulation and one-step free energy perturbation. HISNAPI is capable of not only predicting the hotspot residues in protein-nucleic acid interfaces but also providing insights into their intensity and correlation of dynamic motion. Protein dynamics have been recognized as a vital factor that has an effect on the interaction specificity and affinity of the binding partners. We applied HISNAPI to the case of SARS-CoV-2 RNA-dependent RNA polymerase, a vital target of the antiviral drug for the treatment of coronavirus disease 2019. We identified the hotspot residues and characterized their dynamic behaviors, which might provide insight into the target site for antiviral drug design. The web server is freely available via a user-friendly web interface at http://chemyang.ccnu.edu.cn/ccb/server/HISNAPI/ and http://agroda.gzu.edu.cn:9999/ccb/server/HISNAPI/.

A MIF-Derived Cyclopeptide that Inhibits MIF Binding and Atherogenic Signaling via the Chemokine Receptor CXCR2

Macrophage migration inhibitory factor (MIF) is an inflammatory cytokine and atypical chemokine with a key role in inflammatory diseases including atherosclerosis. Key atherogenic functions of MIF are mediated by noncognate interaction with the chemokine receptor CXCR2. The MIF N-like loop comprising the sequence 47-56 is an important structural determinant of the MIF/CXCR2 interface and MIF(47-56) blocks atherogenic MIF activities. However, the mechanism and critical structure-activity information within this sequence have remained elusive. Here, we show that MIF(47-56) directly binds to CXCR2 to compete with MIF receptor activation. By using alanine scanning, essential and dispensable residues were identified. Moreover, MIF(cyclo10), a designed cyclized variant of MIF(47-56), inhibited key inflammatory and atherogenic MIF activities in vitro and in vivo/ex vivo, and exhibited strongly improved resistance to proteolytic degradation in human plasma in vitro, thus suggesting that it could serve as a promising basis for MIF-derived anti-atherosclerotic peptides.

Computational Analysis of Residue-Specific Binding Free Energies of Androgen Receptor to Ligands

Androgen receptor (AR) is an important therapeutic target for the treatment of diseases such as prostate cancer, hypogonadism, muscle wasting, etc. In this study, the complex structures of the AR ligand-binding domain (LBD) with fifteen ligands were analyzed by molecular dynamics simulations combined with the alanine-scanning-interaction-entropy method (ASIE). The quantitative free energy contributions of the pocket residues were obtained and hotspot residues are quantitatively identified. Our calculation shows that that these hotspot residues are predominantly hydrophobic and their interactions with binding ligands are mainly van der Waals interactions. The total binding free energies obtained by summing over binding contributions by individual residues are in good correlation with the experimental binding data. The current quantitative analysis of binding mechanism of AR to ligands provides important insight on the design of future inhibitors.

An in-silico approach to identify the potential hot spots in SARS-CoV-2 spike RBD to block the interaction with ACE2 receptor

A novel acute viral pneumonia induced by SARS-CoV-2 exploded at the end of 2019, causing a severe medical and economic crisis. For developing specific pharmacotherapy against SARS-CoV-2, an in silico virtual screening was developed for the available in-house molecules. The conserved domain analysis was performed to identify the highly conserved and exposed amino acid regions in the SARS-CoV-2-S RBD sites. The Protein-Protein interaction analyses demonstrated the higher affinity between the SARS-CoV-2-S and ACE2 due to varieties of significant interactions between them. The computational alanine scanning mutation study has recognized the highly stabilized amino acids in the SARS-CoV-2-S RBD/ACE2 complex. The cumulative sequence investigations have inferred that Lys417, Phe486, Asn487, Tyr489, and Gln493 are perhaps the iconic target amino acids to develop a drug molecule or vaccine against SARS-CoV-2 infection. Most of the selected compounds include luteolin, zhebeirine, 3-dehydroverticine, embelin, andrographolide, ophiopogonin D, crocin-1, sprengerinin A, B, C, peimine, etc. were exhibited distinguish drug actions through the strong hydrogen bonding with the hot spots of the RBD. Besides, the 100 ns molecular dynamics simulation and free energy binding analysis showed the significant efficacy of luteolin to inhibit the infection of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Modulation of co-stimulatory signal from CD2-CD58 proteins by a grafted peptide

Peptides were designed to inhibit the protein-protein interaction of CD2 and CD58 to modulate the immune response. This work involved the design and synthesis of eight different peptides by replacing each amino acid residue in peptide 6 with alanine as well as grafting the peptide to the sunflower trypsin-inhibitor framework. From the alanine scanning studies, mutation at position 2 of the peptide was shown to result in increased potency to inhibit cell adhesion interactions. The most potent peptide from the alanine scanning was further studied for its detailed three-dimensional structure and binding to CD58 protein using surface plasmon resonance and flow cytometry. This peptide was used to graft to the sunflower trypsin inhibitor to improve the stability of the peptide. The grafted peptide, SFTI-a1, was further studied for its potency as well as its thermal, chemical, and enzymatic stability. The grafted peptide exhibited improved activity compared to our previously grafted peptide and was stable against thermal and enzymatic degradation.

Coevolution, Dynamics and Allostery Conspire in Shaping Cooperative Binding and Signal Transmission of the SARS-CoV-2 Spike Protein with Human Angiotensin-Converting Enzyme 2

Binding to the host receptor is a critical initial step for the coronavirus SARS-CoV-2 spike protein to enter into target cells and trigger virus transmission. A detailed dynamic and energetic view of the binding mechanisms underlying virus entry is not fully understood and the consensus around the molecular origins behind binding preferences of SARS-CoV-2 for binding with the angiotensin-converting enzyme 2 (ACE2) host receptor is yet to be established. In this work, we performed a comprehensive computational investigation in which sequence analysis and modeling of coevolutionary networks are combined with atomistic molecular simulations and comparative binding free energy analysis of the SARS-CoV and SARS-CoV-2 spike protein receptor binding domains with the ACE2 host receptor. Different from other computational studies, we systematically examine the molecular and energetic determinants of the binding mechanisms between SARS-CoV-2 and ACE2 proteins through the lens of coevolution, conformational dynamics, and allosteric interactions that conspire to drive binding interactions and signal transmission. Conformational dynamics analysis revealed the important differences in mobility of the binding interfaces for the SARS-CoV-2 spike protein that are not confined to several binding hotspots, but instead are broadly distributed across many interface residues. Through coevolutionary network analysis and dynamics-based alanine scanning, we established linkages between the binding energy hotspots and potential regulators and carriers of signal communication in the virus-host receptor complexes. The results of this study detailed a binding mechanism in which the energetics of the SARS-CoV-2 association with ACE2 may be determined by cumulative changes of a number of residues distributed across the entire binding interface. The central findings of this study are consistent with structural and biochemical data and highlight drug discovery challenges of inhibiting large and adaptive protein-protein interfaces responsible for virus entry and infection transmission.

Mutations at the Serine Hydroxymethyltransferase Impact its Interaction with a Soluble NSF Attachment Protein and a Pathogenesis-Related Protein in Soybean

Resistance to soybean cyst nematodes (SCN) in “Peking-type” resistance is bigenic, requiring Rhg4-a and rhg1-a. Rhg4-a encodes a serine hydroxymethyltransferase (GmSHMT08) and rhg1-a encodes a soluble NSF attachment protein (GmSNAP18). Recently, it has been shown that a pathogenesis-related protein, GmPR08-Bet VI, potentiates the interaction between GmSHMT08 and GmSNAP18. Mutational analysis using spontaneously occurring and ethyl methanesulfonate (EMS)-induced mutations was carried out to increase our knowledge of the interacting GmSHMT08/GmSNAP18/GmPR08-Bet VI multi-protein complex. Mutations affecting the GmSHMT08 protein structure (dimerization and tetramerization) and interaction sites with GmSNAP18 and GmPR08-Bet VI proteins were found to impact the multi-protein complex. Interestingly, mutations affecting the PLP/THF substrate binding and catalysis did not affect the multi-protein complex, although they resulted in increased susceptibility to SCN. Most importantly, GmSHMT08 and GmSNAP18 from PI88788 were shown to interact within the cell, being potentiated in the presence of GmPR08-Bet VI. In addition, we have shown the presence of incompatibility between the GmSNAP18 (rhg1-b) of PI88788 and GmSHMT08 (Rhg4-a) from Peking. Components of the reactive oxygen species (ROS) pathway were shown to be induced in the SCN incompatible reaction and were mapped to QTLs for resistance to SCN using different mapping populations.

Molecular Mechanism of Selective Binding of NMS-P118 to PARP-1 and PARP-2: A Computational Perspective

The nuclear protein poly (ADP-ribose) polymerase-1 (PARP-1) inhibitors have been proven effective to potentiate both chemotherapeutic agents and radiotherapy. However, a major problem of most current PARP inhibitors is their lack of selectivity for PARP-1 and its closest isoform PARP-2. NMS-P118 is a highly selective PARP inhibitor that binds PARP-1 stronger than PARP-2 and has many advantages such as excellent pharmacokinetic profiles. In this study, molecular dynamics (MD) simulations of NMS-P118 in complex with PARP-1 and PARP-2 were performed to understand the molecular mechanism of its selectivity. Alanine scanning together with free energy calculation using MM/GBSA and interaction entropy reveal key residues that are responsible for the selectivity. Although the conformation of the binding pockets and NMS-P118 are very similar in PARP-1 and PARP-2, most of the hot-spot residues in PARP-1 have stronger binding free energy than the corresponding residues in PARP-2. Detailed analysis of the binding energy shows that the 4′4-difluorocyclohexyl ring on NMS-P118 form favorable hydrophobic interaction with Y889 in PARP-1. In addition, the H862 residue in PARP-1 has stronger binding free energy than H428 in PARP-2, which is due to shorter distance and stronger hydrogen bonds. Moreover, the negatively charged E763 residue in PARP-1 forms stronger electrostatic interaction energy with the positively charged NMS-P118 than the Q332 residue in PARP-2. These results rationalize the selectivity of NMS-P118 and may be useful for designing novel selective PARP inhibitors.

Understanding microcystin-LR antibody binding interactions using in silico docking and in vitro mutagenesis

Microcystins (MCs) are a group of highly potent cyanotoxins that are becoming more widely distributed due to increased global temperatures and climate change. Microcystin-leucine-arginine (MC-LR) is the most potent and most common variant, with a guideline limit of 1 μg/l in drinking water. We previously developed a novel avian single-chain fragment variable (scFv), designated 2G1, for use in an optical-planar waveguide detection system for microcystin determination. This current work investigates interactions between 2G1 and MC-LR at the molecular level through modelling with an avian antibody template and molecular docking by AutoDock Vina to identify key amino acid (AA) residues involved. These potential AA interactions were investigated in vitro by targeted mutagenesis, specifically, by alanine scanning mutations. Glutamic acid (E) was found to play a critical role in the 2G1-MC-LR binding interaction, with the heavy chain glutamic acid (E) 102 (H-E102) forming direct bonds with the arginine (R) residue of MC-LR. In addition, alanine mutation of light chain residue aspartic acid 57 (L-D57) led to an improvement in antigen-binding observed using enzyme-linked immunosorbent assay (ELISA), and was confirmed by surface plasmon resonance (SPR). This work will contribute to improving the binding of recombinant anti-MC-LR to its antigen and aid in the development of a higher sensitivity harmful algal toxin diagnostic.

Computer-guided library generation applied to the optimization of single-domain antibodies

Computer-guided library generation is a plausible strategy to optimize antibodies. Herein, we report the improvement of the affinity of a single-domain camelid antibody for its antigen using such approach. We first conducted experimental and computational alanine scanning to describe the precise energetic profile of the antibody-antigen interaction surface. Based on this characterization, we hypothesized that in-silico mutagenesis could be employed to guide the development of a small library for phage display with the goal of improving the affinity of an antibody for its antigen. Optimized antibody mutants were identified after three rounds of selection, in which an alanine residue at the core of the antibody-antigen interface was substituted by residues with large side-chains, generating diverse kinetic responses, and resulting in greater affinity (>10-fold) for the antigen.

Preprotein signature for full susceptibility to the co-translational translocation inhibitor cyclotriazadisulfonamide

Cyclotriazadisulfonamide (CADA) inhibits the co-translational translocation of human CD4 (huCD4) into the endoplasmic reticulum lumen in a signal peptide (SP)-dependent way. We propose that CADA binds the nascent huCD4 SP in a folded conformation within the translocon resembling a normally transitory state during translocation. Here, we used alanine scanning on the huCD4 SP to identify the signature for full susceptibility to CADA. In accordance with our previous work, we demonstrate that residues in the vicinity of the hydrophobic h-region are critical for sensitivity to CADA. In particular, exchanging Gln-15, Val-17 or Pro-20 in the huCD4 SP for Ala resulted in a resistant phenotype. Together with positively charged residues at the N-terminal portion of the mature protein, these residues mediate full susceptibility to the co-translational translocation inhibitory activity of CADA towards huCD4. In addition, sensitivity to CADA is inversely related to hydrophobicity in the huCD4 SP. In vitro translocation experiments confirmed that the general hydrophobicity of the h-domain and positive charges in the mature protein are key elements that affect both the translocation efficiency of huCD4 and the sensitivity towards CADA. Besides these two general SP parameters that determine the functionality of the signal sequence, unique amino acid pairs (L14/Q15 and L19/P20) in the SP hydrophobic core add specificity to the sensitivity signature for a co-translational translocation inhibitor.

Protein-RNA interactions: structural characteristics and hotspot amino acids

Structural information about protein-RNA complexes supports the understanding of crucial recognition processes in the cell, and it can allow the development of high affinity ligands to interfere with these processes. In this respect, the identification of amino acid hotspots is particularly important. In contrast to protein-protein interactions, in silico approaches for protein-RNA interactions lag behind in their development. Herein, we report an analysis of available protein-RNA structures. We assembled a data set of 322 crystal and NMR structures and analyzed them regarding interface properties. In addition, we describe a computational alanine-scanning approach which provides interaction scores for interface amino acids, allowing the identification of potential hotspots in protein-RNA interfaces. We have made the computational approach available as an online tool, which allows interaction scores to be calculated for any structure of a protein-RNA complex by uploading atomic coordinates to the PRI HotScore web server (https://pri-hotscore.labs.vu.nl).

Theoretical Study on Zearalenol Compounds Binding with Wild Type Zearalenone Hydrolase and V153H Mutant

Zearalenone hydrolase (ZHD) is the only reported α/β-hydrolase that can detoxify zearalenone (ZEN). ZHD has demonstrated its potential as a treatment for ZEN contamination that will not result in damage to cereal crops. Recent researches have shown that the V153H mutant ZHD increased the specific activity against α-ZOL, but decreased its specific activity to β-ZOL. To understand whyV153H mutation showed catalytic specificity for α-ZOL, four molecular dynamics simulations combining with protein network analysis for wild type ZHD α-ZOL, ZHD β-ZOL, V153H α-ZOL, and V153H β-ZOL complexes were performed using Gromacs software. Our theoretical results indicated that the V153H mutant could cause a conformational switch at the cap domain (residues Gly161–Thr190) to affect the relative position catalytic residue (H242). Protein network analysis illustrated that the V153H mutation enhanced the communication with the whole protein and residues with high betweenness in the four complexes, which were primarily assembled in the cap domain and residues Met241 to Tyr245 regions. In addition, the existence of α-ZOL binding to V153H mutation enlarged the distance from the O_AE atom in α-ZOL to the NE2 atom in His242, which prompted the side chain of H242 to the position with catalytic activity, thereby increasing the activity of V153H on the α-ZOL. Furthermore, α-ZOL could easily form a right attack angle and attack distance in the ZHD and α-ZOL complex to guarantee catalytic reaction. The alanine scanning results indicated that modifications of the residues in the cap domain produced significant changes in the binding affinity for α-ZOL and β-ZOL. Our results may provide useful theoretical evidence for the mechanism underlying the catalytic specificity of ZHD.

Identification of two distinct peptide-binding pockets in the SH3 domain of human mixed-lineage kinase 3

Mixed-lineage kinase 3 (MLK3; also known as MAP3K11) is a Ser/Thr protein kinase widely expressed in normal and cancerous tissues, including brain, lung, liver, heart, and skeletal muscle tissues. Its Src homology 3 (SH3) domain has been implicated in MLK3 autoinhibition and interactions with other proteins, including those from viruses. The MLK3 SH3 domain contains a six-amino-acid insert corresponding to the n-Src insert, suggesting that MLK3 may bind additional peptides. Here, affinity selection of a phage-displayed combinatorial peptide library for MLK3's SH3 domain yielded a 13-mer peptide, designated "MLK3 SH3-interacting peptide" (MIP). Unlike most SH3 domain peptide ligands, MIP contained a single proline. The 1.2-Å crystal structure of the MIP-bound SH3 domain revealed that the peptide adopts a β-hairpin shape, and comparison with a 1.5-Å apo SH3 domain structure disclosed that the n-Src loop in SH3 undergoes an MIP-induced conformational change. A 1.5-Å structure of the MLK3 SH3 domain bound to a canonical proline-rich peptide from hepatitis C virus nonstructural 5A (NS5A) protein revealed that it and MIP bind the SH3 domain at two distinct sites, but biophysical analyses suggested that the two peptides compete with each other for SH3 binding. Moreover, SH3 domains of MLK1 and MLK4, but not MLK2, also bound MIP, suggesting that the MLK1-4 family may be differentially regulated through their SH3 domains. In summary, we have identified two distinct peptide-binding sites in the SH3 domain of MLK3, providing critical insights into mechanisms of ligand binding by the MLK family of kinases.

Understanding Insulin Endocrinology in Decapod Crustacea: Molecular Modelling Characterization of an Insulin-Binding Protein and Insulin-Like Peptides in the Eastern Spiny Lobster, Sagmariasus verreauxi

The insulin signalling system is one of the most conserved endocrine systems of Animalia from mollusc to man. In decapod Crustacea, such as the Eastern spiny lobster, Sagmariasus verreauxi (Sv) and the red-claw crayfish, Cherax quadricarinatus (Cq), insulin endocrinology governs male sexual differentiation through the action of a male-specific, insulin-like androgenic gland peptide (IAG). To understand the bioactivity of IAG it is necessary to consider its bio-regulators such as the insulin-like growth factor binding protein (IGFBP). This work has employed various molecular modelling approaches to represent S. verreauxi IGFBP and IAG, along with additional Sv-ILP ligands, in order to characterise their binding interactions. Firstly, we present Sv- and Cq-ILP2: neuroendocrine factors that share closest homology with Drosophila ILP8 (Dilp8). We then describe the binding interaction of the N-terminal domain of Sv-IGFBP and each ILP through a synergy of computational analyses. In-depth interaction mapping and computational alanine scanning of IGFBP_N' highlight the conserved involvement of the hotspot residues Q₆₇, G₇₀, D₇₁, S₇₂, G₉₁, G₉₂, T₉₃ and D₉₄. The significance of the negatively charged residues D₇₁ and D₉₄ was then further exemplified by structural electrostatics. The functional importance of the negative surface charge of IGFBP is exemplified in the complementary electropositive charge on the reciprocal binding interface of all three ILP ligands. When examined, this electrostatic complementarity is the inverse of vertebrate homologues; such physicochemical divergences elucidate towards ligand-binding specificity between Phyla.

Functional Segregation of Overlapping Genes in HIV

Overlapping genes pose an evolutionary dilemma as one DNA sequence evolves under the selection pressures of multiple proteins. Here, we perform systematic statistical and mutational analyses of the overlapping HIV-1 genes tat and rev and engineer exhaustive libraries of non-overlapped viruses to perform deep mutational scanning of each gene independently. We find a "segregated" organization in which overlapped sites encode functional residues of one gene or the other, but never both. Furthermore, this organization eliminates unfit genotypes, providing a fitness advantage to the population. Our comprehensive analysis reveals the extraordinary manner in which HIV minimizes the constraint of overlapping genes and repurposes that constraint to its own advantage. Thus, overlaps are not just consequences of evolutionary constraints, but rather can provide population fitness advantages.

Alanine-scanning mutagenesis of human signal transducer and activator of transcription 1 to estimate loss- or gain-of-function variants

BACKGROUND

Germline heterozygous mutations in human signal transducer and activator of transcription 1 (STAT1) can cause loss of function (LOF), as in patients with Mendelian susceptibility to mycobacterial diseases, or gain of function (GOF), as in patients with chronic mucocutaneous candidiasis. LOF and GOF mutations are equally rare and can affect the same domains of STAT1, especially the coiled-coil domain (CCD) and DNA-binding domain (DBD). Moreover, 6% of patients with chronic mucocutaneous candidiasis with a GOF STAT1 mutation have mycobacterial disease, obscuring the functional significance of the identified STAT1 mutations. Current computational approaches, such as combined annotation-dependent depletion, do not distinguish LOF and GOF variants.

OBJECTIVE

We estimated variations in the CCD/DBD of STAT1.

METHODS

We mutagenized 342 individual wild-type amino acids in the CCD/DBD (45.6% of full-length STAT1) to alanine and tested the mutants for STAT1 transcriptional activity.

RESULTS

Of these 342 mutants, 201 were neutral, 30 were LOF, and 111 were GOF mutations in a luciferase assay. This assay system correctly estimated all previously reported LOF mutations (100%) and slightly fewer GOF mutations (78.1%) in the CCD/DBD of STAT1. We found that GOF alanine mutants occurred at the interface of the antiparallel STAT1 dimer, suggesting that they destabilize this dimer. This assay also precisely predicted the effect of 2 hypomorphic and dominant negative mutations, E157K and G250E, in the CCD of STAT1 that we found in 2 unrelated patients with Mendelian susceptibility to mycobacterial diseases.

CONCLUSION

The systematic alanine-scanning assay is a useful tool to estimate the GOF or LOF status and the effect of heterozygous missense mutations in STAT1 identified in patients with severe infectious diseases, including mycobacterial and fungal diseases.

Basigin Interacts with Plasmodium vivax Tryptophan-rich Antigen PvTRAg38 as a Second Erythrocyte Receptor to Promote Parasite Growth

Elucidating the molecular mechanisms of the host-parasite interaction during red cell invasion by Plasmodium is important for developing newer antimalarial therapeutics. Recently, we have characterized a Plasmodium vivax tryptophan-rich antigen PvTRAg38, which is expressed by its merozoites, binds to host erythrocytes, and interferes with parasite growth. Interaction of this parasite ligand with the host erythrocyte occurs through its two regions present at amino acid positions 167-178 (P₂) and 197-208 (P₄). Each region recognizes its own erythrocyte receptor. Previously, we identified band 3 as the chymotrypsin-sensitive erythrocyte receptor for the P₄ region, but the other receptor, binding to P₂ region, remained unknown. Here, we have identified basigin as the second erythrocyte receptor for PvTRAg38, which is resistant to chymotrypsin. The specificity of interaction between PvTRAg38 and basigin was confirmed by direct interaction where basigin was specifically recognized by P₂ and not by the P₄ region of this parasite ligand. Interaction between P₂ and basigin is stabilized through multiple amino acid residues, but Gly-171 and Leu-175 of P₂ were more critical. These two amino acids were also critical for parasite growth. Synthetic peptides P₂ and P₄ of PvTRAg38 interfered with the parasite growth independently but had an additive effect if combined together indicating involvement of both the receptors during red cell invasion. In conclusion, PvTRAg38 binds to two erythrocyte receptors basigin and band 3 through P₂ and P₄ regions, respectively, to facilitate parasite growth. This advancement in our knowledge on molecular mechanisms of host-parasite interaction can be exploited to develop therapeutics against P. vivax malaria.

Global mapping of antibody recognition of the hepatitis C virus E2 glycoprotein: Implications for vaccine design

The E2 envelope glycoprotein is the primary target of human neutralizing antibody response against hepatitis C virus (HCV), and is thus a major focus of vaccine and immunotherapeutics efforts. There is emerging evidence that E2 is a highly complex, dynamic protein with residues across the protein that are modulating antibody recognition, local and global E2 stability, and viral escape. To comprehensively map these determinants, we performed global E2 alanine scanning with a panel of 16 human monoclonal antibodies (hmAbs), resulting in an unprecedented dataset of the effects of individual alanine substitutions across the E2 protein (355 positions) on antibody recognition. Analysis of shared energetic effects across the antibody panel identified networks of E2 residues involved in antibody recognition and local and global E2 stability, as well as predicted contacts between residues across the entire E2 protein. Further analysis of antibody binding hotspot residues defined groups of residues essential for E2 conformation and recognition for all 14 conformationally dependent E2 antibodies and subsets thereof, as well as residues that enhance antibody recognition when mutated to alanine, providing a potential route to engineer E2 vaccine immunogens. By incorporating E2 sequence variability, we found a number of E2 polymorphic sites that are responsible for loss of neutralizing antibody binding. These data and analyses provide fundamental insights into antibody recognition of E2, highlighting the dynamic and complex nature of this viral envelope glycoprotein, and can serve as a reference for development and rational design of E2-targeting vaccines and immunotherapeutics.

Enabling the Discovery and Virtual Screening of Potent and Safe Antimicrobial Peptides. Simultaneous Prediction of Antibacterial Activity and Cytotoxicity

Antimicrobial peptides (AMPs) represent promising alternatives to fight against bacterial pathogens. However, cellular toxicity remains one of the main concerns in the early development of peptide-based drugs. This work introduces the first multitasking (mtk) computational model focused on performing simultaneous predictions of antibacterial activities, and cytotoxicities of peptides. The model was created from a data set containing 3592 cases, and it displayed accuracy higher than 96% for classifying/predicting peptides in both training and prediction (test) sets. The technique known as alanine scanning was computationally applied to illustrate the calculation of the quantitative contributions of the amino acids (in their respective positions of the sequence) to the biological effects of a defined peptide. A small library formed by 10 peptides was generated, where peptides were designed by considering the interpretations of the different descriptors in the mtk-computational model. All the peptides were predicted to exhibit high antibacterial activities against multiple bacterial strains, and low cytotoxicity against various cell types. The present mtk-computational model can be considered a very useful tool to support high throughput research for the discovery of potent and safe AMPs.

Structure-activity relationship of the peptide binding-motif mediating the BRCA2:RAD51 protein-protein interaction

RAD51 is a recombinase involved in the homologous recombination of double‐strand breaks in DNA. RAD51 forms oligomers by binding to another molecule of RAD51 via an ‘FxxA’ motif, and the same recognition sequence is similarly utilised to bind BRCA2. We have tabulated the effects of mutation of this sequence, across a variety of experimental methods and from relevant mutations observed in the clinic. We use mutants of a tetrapeptide sequence to probe the binding interaction, using both isothermal titration calorimetry and X‐ray crystallography. Where possible, comparison between our tetrapeptide mutational study and the previously reported mutations is made, discrepancies are discussed and the importance of secondary structure in interpreting alanine scanning and mutational data of this nature is considered.

Host-parasite interaction: multiple sites in the Plasmodium vivax tryptophan-rich antigen PvTRAg38 interact with the erythrocyte receptor band 3

Tryptophan‐rich antigens of malarial parasites interact with host molecules and play an important role in parasite survival. Merozoite expressed Plasmodium vivax tryptophan‐rich antigen PvTRAg38 binds to human erythrocytes and facilitates parasite growth in a heterlologous Plasmodium falciparum culture system. Recently, we identified band 3 in human erythrocytes as one of its receptors, although the receptor‐ligand binding mechanisms remain unknown. In the present study, using synthetic mutated peptides of PvTRAg38, we show that multiple amino acid residues of its 12 amino acid domain (KWVQWKNDKIRS) at position 197–208 interact with three different ectodomains of band 3 receptor on human erythrocytes. Our findings may help in the design of new therapeutic approaches for malaria.

Molecular dynamics study of biodegradation of azo dyes via their interactions with AzrC azoreductase

Azo dyes are one of the most important class of dyes, which have been widely used in industries. Because of the environmental pollution of azo dyes, many studies have been performed to study their biodegradation using bacterial systems. In present work, the AzrC of mesophilic gram-positive Bacillus sp. B29 has been considered to study its interaction with five common azo dyes (orange G, acid red 88, Sudan I, orange I, and methyl red). The molecular dynamics simulations have been employed to study the interaction between AzrC and azo dyes. The trajectory was confirmed using root mean square deviation and the root mean square fluctuation analyses. Then, the hydrogen bond and alanine scanning analyses were performed to reveal active site residues. Phe105 (A), Phe125 (B), Phe172 (B), and Pro132 (B) have been found as the most important hydrophobic residues whereas Asn104 (A), Tyr127 (B), and Asn187 (A) have key role in making hydrogen bond. The results of molecular mechanics Poisson-Boltzmann surface area and molecular mechanics generalized Born surface area calculations proved that the hydrophobic azo dyes like Acid red 88 binds more tightly to the AzrC protein. The calculated data suggested MR A 121 (B) I as a potential candidate for improving the AzrC-MR interactions.

Stabilization of G protein-coupled receptors by point mutations

G protein-coupled receptors (GPCRs) are flexible integral membrane proteins involved in transmembrane signaling. Their involvement in many physiological processes makes them interesting targets for drug development. Determination of the structure of these receptors will help to design more specific drugs, however, their structural characterization has so far been hampered by the low expression and their inherent instability in detergents which made protein engineering indispensable for structural and biophysical characterization. Several approaches to stabilize the receptors in a particular conformation have led to breakthroughs in GPCR structure determination. These include truncations of the flexible regions, stabilization by antibodies and nanobodies, fusion partners, high affinity and covalently bound ligands as well as conformational stabilization by mutagenesis. In this review we focus on stabilization of GPCRs by insertion of point mutations, which lead to increased conformational and thermal stability as well as improved expression levels. We summarize existing mutagenesis strategies with different coverage of GPCR sequence space and depth of information, design and transferability of mutations and the molecular basis for stabilization. We also discuss whether mutations alter the structure and pharmacological properties of GPCRs.

Mixture-based combinatorial libraries from small individual peptide libraries: a case study on α1-antitrypsin deficiency

The design, synthesis and screening of diversity-oriented peptide libraries using a "libraries from libraries" strategy for the development of inhibitors of α1-antitrypsin deficiency are described. The major buttress of the biochemical approach presented here is the use of well-established solid-phase split-and-mix method for the generation of mixture-based libraries. The combinatorial technique iterative deconvolution was employed for library screening. While molecular diversity is the general consideration of combinatorial libraries, exquisite design through systematic screening of small individual libraries is a prerequisite for effective library screening and can avoid potential problems in some cases. This review will also illustrate how large peptide libraries were designed, as well as how a conformation-sensitive assay was developed based on the mechanism of the conformational disease. Finally, the combinatorially selected peptide inhibitor capable of blocking abnormal protein aggregation will be characterized by biophysical, cellular and computational methods.

IFITMs restrict the replication of multiple pathogenic viruses

The interferon-inducible transmembrane protein (IFITM) family inhibits a growing number of pathogenic viruses, among them influenza A virus, dengue virus, hepatitis C virus, and Ebola virus. This review covers recent developments in our understanding of the IFITM's molecular determinants, potential mechanisms of action, and impact on pathogenesis.

News from the protein mutability landscape

Some mutations of protein residues matter more than others, and these are often conserved evolutionarily. The explosion of deep sequencing and genotyping increasingly requires the distinction between effect and neutral variants. The simplest approach predicts all mutations of conserved residues to have an effect; however, this works poorly, at best. Many computational tools that are optimized to predict the impact of point mutations provide more detail. Here, we expand the perspective from the view of single variants to the level of sketching the entire mutability landscape. This landscape is defined by the impact of substituting every residue at each position in a protein by each of the 19 non-native amino acids. We review some of the powerful conclusions about protein function, stability and their robustness to mutation that can be drawn from such an analysis. Large-scale experimental and computational mutagenesis experiments are increasingly furthering our understanding of protein function and of the genotype-phenotype associations. We also discuss how these can be used to improve predictions of protein function and pathogenicity of missense variants.

Computational studies of difference in binding modes of peptide and non-peptide inhibitors to MDM2/MDMX based on molecular dynamics simulations

Inhibition of p53-MDM2/MDMX interaction is considered to be a promising strategy for anticancer drug design to activate wild-type p53 in tumors. We carry out molecular dynamics (MD) simulations to study the binding mechanisms of peptide and non-peptide inhibitors to MDM2/MDMX. The rank of binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) method agrees with one of the experimental values. The results suggest that van der Waals energy drives two kinds of inhibitors to MDM2/MDMX. We also find that the peptide inhibitors can produce more interaction contacts with MDM2/MDMX than the non-peptide inhibitors. Binding mode predictions based on the inhibitor-residue interactions show that the π–π, CH–π and CH–CH interactions dominated by shape complimentarity, govern the binding of the inhibitors in the hydrophobic cleft of MDM2/MDMX. Our studies confirm the residue Tyr99 in MDMX can generate a steric clash with the inhibitors due to energy and structure. This finding may theoretically provide help to develop potent dual-specific or MDMX inhibitors.

Site2 binding energetics of the regulatory step of growth hormone-induced receptor homodimerization

Receptor signaling in the growth hormone (GH)-growth hormone receptor (GHR) system is controlled through a sequential two-step hormone-induced dimerization of two copies of the extracellular domain (ECD) of the receptor. The regulatory step of this process is the binding of the second ECD (ECD2) to the stable preassociated 1 : 1 GH/ECD1 complex on the cell surface. To determine the energetics that governs this step, the binding kinetics of 38 single- and double-alanine mutants in the hGH Site2 contact with ECD2 were measured by using trimolecular surface plasmon resonance (TM-SPR). We find that the Site2 interface of hGH does not have a distinct binding hot-spot region, and the most important residues are not spatially clustered, but rather are distributed over the whole binding surface. In addition, it was determined through analysis of a set of pairwise double alanine mutations that there is a significant degree of negative cooperativity among Site2 residues. Residues that show little effect or even improved binding on substitution with alanine, when paired with D116A-hGH, display significant negative cooperativity. Because most of these pairwise mutated residues are spatially separated by >or=10 A, this indicates that the Site2 binding interface of the hGH-hGHR ternary complex displays both structural and energetic malleability.