Buried charged surface in proteins

Structures of the ribosome bound to EF-Tu-isoleucine tRNA elucidate the mechanism of AUG avoidance

The frequency of errors upon decoding of messenger RNA by the bacterial ribosome is low, with one misreading event per 1 × 104 codons. In the universal genetic code, the AUN codon box specifies two amino acids, isoleucine and methionine. In bacteria and archaea, decoding specificity of the AUA and AUG codons relies on the wobble avoidance strategy that requires modification of C34 in the anticodon loop of isoleucine transfer RNAIleCAU (tRNAIleCAU). Bacterial tRNAIleCAU with 2-lysylcytidine (lysidine) at the wobble position deciphers AUA while avoiding AUG. Here we report cryo-electron microscopy structures of the Escherichia coli 70S ribosome complexed with elongation factor thermo unstable (EF-Tu) and isoleucine-tRNAIleLAU in the process of decoding AUA and AUG. Lysidine in tRNAIleLAU excludes AUG by promoting the formation of an unusual Hoogsteen purine-pyrimidine nucleobase geometry at the third position of the codon, weakening the interactions with the mRNA and destabilizing the EF-Tu ternary complex. Our findings elucidate the molecular mechanism by which tRNAIleLAU specifically decodes AUA over AUG.

Water will Find Its Way: Transport through Narrow Tunnels in Hydrolases

An aqueous environment is vital for life as we know it, and water is essential for nearly all biochemical processes at the molecular level. Proteins utilize water molecules in various ways. Consequently, proteins must transport water molecules across their internal network of tunnels to reach the desired action sites, either within them or by functioning as molecular pipes to control cellular osmotic pressure. Despite water playing a crucial role in enzymatic activity and stability, its transport has been largely overlooked, with studies primarily focusing on water transport across membrane proteins. The transport of molecules through a protein's tunnel network is challenging to study experimentally, making molecular dynamics simulations the most popular approach for investigating such events. In this study, we focused on the transport of water molecules across three different α/β-hydrolases: haloalkane dehalogenase, epoxide hydrolase, and lipase. Using a 5 μs adaptive simulation per system, we observed that only a few tunnels were responsible for the majority of water transport in dehalogenase, in contrast to a higher diversity of tunnels in other enzymes. Interestingly, water molecules could traverse narrow tunnels with subangstrom bottlenecks, which is surprising given the commonly accepted water molecule radius of 1.4 Å. Our analysis of the transport events in such narrow tunnels revealed a markedly increased number of hydrogen bonds formed between the water molecules and protein, likely compensating for the steric penalty of the process. Overall, these commonly disregarded narrow tunnels accounted for ∼20% of the total water transport observed, emphasizing the need to surpass the standard geometrical limits on the functional tunnels to properly account for the relevant transport processes. Finally, we demonstrated how the obtained insights could be applied to explain the differences in a mutant of the human soluble epoxide hydrolase associated with a higher incidence of ischemic stroke.

Molecular investigation of antimicrobial peptides against Helicobacter pylori proteins using a peptide-protein docking approach

The impact of H. pylori resistance on patient's treatment failure is a major concern. Therefore, the development of novel or alternative therapies for H. pylori is urgently needed. The purpose of this study was to investigate the molecular interactions of various antimicrobial peptides (AMPs) to H. pylori proteins. We performed the peptide-protein molecular docking using HADDOCK 2.4 webserver. Fourteen AMPs were tested for their binding efficacy against four H. pylori proteins. Simulation of the peptide-protein complex was performed using molecular dynamic software package AMBER20. From molecular docking analysis, five peptides (LL-37, Tilapia piscidin 4, napin, snakin-1 and EcAMP1) showed strong binding interactions against H. pylori proteins. The strongest binding affinity was observed in the interactions between Snakin-1 and PBP2, TP4 and type I HopQ and EcAMP1 and type I HopQ with -11.1, -13.6 and -13.8 kcal/mol, respectively. The dynamic simulation was performed for two complexes (snakin1-PBP2 and EcAMP1-HopQ). Results of the dynamics simulation showed that EcAMP1 had stable interaction and binding to type I HopQ protein without significant structural changes. In conclusion, both results of docking and simulation showed that EcAMP1 might be useful as a potential therapeutic agent for H. pylori treatment. This molecular approach provides deep understanding of the interaction insights between AMPs and H. pylori proteins. It paves the way for the development of novel anti-H. pylori using antimicrobial peptides.

Living in trinity of extremes: Genomic and proteomic signatures of halophilic, thermophilic, and pH adaptation

Since nucleic acids and proteins of unicellular prokaryotes are directly exposed to extreme environmental conditions, it is possible to explore the genomic-proteomic compositional determinants of molecular mechanisms of adaptation developed by them in response to harsh environmental conditions. Using a wealth of currently available complete genomes/proteomes we were able to explore signatures of adaptation to three environmental factors, pH, salinity, and temperature, observing major trends in compositions of their nucleic acids and proteins. We derived predictors of thermostability, halophilic, and pH adaptations and complemented them by the principal components analysis. We observed a clear difference between thermophilic and salinity/pH adaptations, whereas latter invoke seemingly overlapping mechanisms. The genome-proteome compositional trade-off reveals an intricate balance between the work of base paring and base stacking in stabilization of coding DNA and r/tRNAs, and, at the same time, universal requirements for the stability and foldability of proteins regardless of the nucleotide biases. Nevertheless, we still found hidden fingerprints of ancient evolutionary connections between the nucleotide and amino acid compositions indicating their emergence, mutual evolution, and adjustment. The evolutionary perspective on the adaptation mechanisms is further studied here by means of the comparative analysis of genomic/proteomic traits of archaeal and bacterial species. The overall picture of genomic/proteomic signals of adaptation obtained here provides a foundation for future engineering and design of functional biomolecules resistant to harsh environments.

Low-sugar universal mRNA vaccine against coronavirus variants with deletion of glycosites in the S2 or stem of SARS-CoV-2 spike messenger RNA (mRNA)

Since the outbreak of Severe Acute Respiratory Syndrome Virus-2 (SARS-CoV-2) in 2019, more than 15 million spike protein sequences have been identified, raising a new challenge for the development of a broadly protective vaccine against the various emerging variants. We found that the virus, like most other human viruses, depends on host-made glycans to shield the conserved epitopes on spike protein from immune response and demonstrated that deletion of the glycan shields exposed highly conserved epitopes and elicited broadly protective immune responses. In this study, we identified 17 conserved epitopes from 14 million spike protein sequences and 11 of the conserved epitopes are in the S2 domain, including the six most conserved epitopes in the stem region. We also demonstrated that deletion of the glycosites in the spike messenger RNA (mRNA) S2 domain or the stem region exposed the highly conserved epitopes and elicited broadly protective immune responses, particularly CD-8+ T cell response against various SARS-CoV-2 variants, and other human coronaviruses including MERS, SARS viruses, and those causing common cold.

Conformational features and interaction mechanisms of VH H antibodies with β-hairpin CDR3: A case of Nb8-HigB2 interaction

The β-hairpin conformation is regarded as an important basic motif to form and regulate protein-protein interactions. Single-domain VH H antibodies are potential therapeutic and diagnostic tools, and the third complementarity-determining regions of the heavy chains (CDR3s) of these antibodies are critical for antigen recognition. Although the sequences and conformations of the CDR3s are diverse, CDR3s sometimes adopt β-hairpin conformations. However, characteristic features and interaction mechanisms of β-hairpin CDR3s remain to be fully elucidated. In this study, we investigated the molecular recognition of the anti-HigB2 VH H antibody Nb8, which has a CDR3 that forms a β-hairpin conformation. The interaction was analyzed by evaluation of alanine-scanning mutants, molecular dynamics simulations, and hydrogen/deuterium exchange mass spectrometry. These experiments demonstrated that positions 93 and 94 (Chothia numbering) in framework region 3, which is right outside CDR3 by definition, play pivotal roles in maintaining structural stability and binding properties of Nb8. These findings will facilitate the design and optimization of single-domain antibodies.

From in-silico screening to in-vitro evaluation: Enhancing the detection of Microcystins with engineered PP1 mutant variants

Cyanotoxins produced during harmful algal blooms (CyanoHABs) have become a worldwide issue of concern. Microcystins (MC) are the most ubiquitous group of cyanotoxins and have known carcinogenic and hepatotoxic effects. The protein phosphatase inhibition assays (PPIAs), based on the inhibition of Protein Phosphatase 1/2A (PP1/PP2A) by MC, are one of the most cost-effective options for detecting MC. In this work, we aimed to design in-silico and evaluate in-vitro mutant variants of the PP1 protein, in order to enhance their capabilities as a MC biosensor. To this end, we performed an in-silico active site-saturated mutagenesis screening, followed by stability and docking affinity calculation with the MCLR cyanotoxin. Candidates with improved both affinity and stability were further tested in a fully flexible active-site docking. The best-scored mutations (19) were individually analysed regarding their locations and interactions. Four of them (p.D197F; p.Q249Y; p.S129W; p.D220Q) were selected for in-vitro expression and evaluation. Mutant p.D197F, exhibited a significant increment in inhibition by MCLR with respect to the WT, while showing a non-significant difference in stability nor activity. This successful PP1 inhibition enhancement suggests the potential of the p.D197F variant for practical MC detection applications.

Alloreactivity and autoreactivity converge to support B cell epitope targeting in transplant rejection

Antibody (Ab) responses against human leukocyte antigen (HLA) proteins mismatched between donor and recipient are leading cause of allograft loss in kidney transplantation. However, therapies targeting alloreactive B cell and Ab-secreting cell (ASC) are lacking, motivating the need to understand how to prevent and abrogate these alloresponses. Using molecular, structural, and proteomic techniques, we profiled the B cell response in a kidney transplant recipient with antibody-mediated rejection and graft loss. We found that this response spanned the rejected organ and peripheral blood, stimulated the differentiation of multiple B cell subsets, and produced a high-affinity, donor-specific, anti-HLA response. We found epitopic immunodominance that relied on highly exposed, solvent-accessible mismatched HLA residues as well as structural and biomolecular evidence of autoreactivity against the recipient's self-HLA allele. These alloreactive and autoreactive signatures converged in the recipient's circulating donor-specific Ab repertoire, suggesting that rejection requires both the recognition of non-self and breaches of tolerance to lead to alloinjury and graft loss.

Structures of Plasmodium falciparum Chloroquine Resistance Transporter (PfCRT) Isoforms and Their Interactions with Chloroquine

Using a recently elucidated atomic-resolution cryogenic electron microscopy (cryo-EM) structure for the Plasmodium falciparum chloroquine resistance transporter (PfCRT) protein 7G8 isoform as template [Kim, J.; Nature 2019, 576, 315-320], we use Monte Carlo molecular dynamics (MC/MD) simulations of PfCRT embedded in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane to solve energy-minimized structures for 7G8 PfCRT and two additional PfCRT isoforms that harbor 5 or 7 amino acid substitutions relative to 7G8 PfCRT. Guided by drug binding previously defined using chloroquine (CQ) photoaffinity probe labeling, we also use MC/MD energy minimization to elucidate likely CQ binding geometries for the three membrane-embedded isoforms. We inventory salt bridges and hydrogen bonds in these structures and summarize how the limited changes in primary sequence subtly perturb local PfCRT isoform structure. In addition, we use the "AlphaFold" artificial intelligence AlphaFold2 (AF2) algorithm to solve for domain structure that was not resolved in the previously reported 7G8 PfCRT cryo-EM structure, and perform MC/MD energy minimization for the membrane-embedded AF2 structures of all three PfCRT isoforms. We compare energy-minimized structures generated using cryo-EM vs AF2 templates. The results suggest how amino acid substitutions in drug resistance-associated isoforms of PfCRT influence PfCRT structure and CQ transport.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Molecular interactions of hesperidin with DMPC/cholesterol bilayers

Since cell membranes are complex systems, the use of model lipid bilayers is quite important for the study of their interactions with bioactive molecules. Mammalian cell membranes require cholesterol (CHOL) for their structure and function. For this reason, the mixtures of phospholipid and cholesterol are necessary to use in model membrane studies to better simulate the real systems. In the present study, we investigated the effect of the incorporation of hesperidin in model membranes consisting of dimyristoylphosphatidylcholine (DMPC) and CHOL by using differential scanning calorimetry (DSC), attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy, and atomic force microscopy (AFM). ATR-FTIR results demonstrated that hesperidin increases the fluidity of the DMPC/CHOL binary system. DSC findings indicated that the presence of 5 mol% hesperidin induces a broadening of the main phase transition consisting of three overlapping components. AFM experiments showed that hesperidin increases the thickness of DMPC/CHOL lipid bilayer model membranes. In addition to experimental results, molecular docking studies were conducted with hesperidin and human lanosterol synthase (LS), which is an enzyme found in the final step of cholesterol synthesis, to characterize hesperidin's interactions with its surrounding via its hydroxyl and oxygen groups. Then, hesperidin's ADME/Tox (absorption, distribution, metabolism, excretion and toxicity) profile was computed to see the potential impact on living system. In conclusion, considering the data obtained from experimental studies, this work ensures molecular insights in the interaction between a flavonoid, as an antioxidant drug model, and lipids mimicking those found in mammalian membranes. Moreover, computational studies demonstrated that hesperidin may be a great potential for use as a therapeutic agent for hypercholesterolemia due to its antioxidant property.

A comprehensive in silico exploration of the impacts of missense variants on two different conformations of human pirin protein

BACKGROUND

Pirin, a member of the cupin superfamily, is an iron-binding non-heme protein. It acts as a coregulator of several transcription factors, especially the members of NFκB transcription factor family. Based on the redox state of its iron cofactor, it can assume two different conformations and thereby act as a redox sensor inside the nucleus. Previous studies suggested that pirin may be associated with cancer, inflammatory diseases as well as COVID-19 severities. Hence, it is important to explore the pathogenicity of its missense variants. In this study, we used a number of in silico tools to investigate the effects of missense variants of pirin on its structure, stability, metal cofactor binding affinity and interactions with partner proteins. In addition, we used protein dynamics simulation to elucidate the effects of selected variants on its dynamics. Furthermore, we calculated the frequencies of haplotypes containing pirin missense variants across five major super-populations (African, Admixed American, East Asian, European and South Asian).

RESULTS

Among a total of 153 missense variants of pirin, 45 were uniformly predicted to be pathogenic. Of these, seven variants can be considered for further experimental studies. Variants R59P and L116P were predicted to significantly destabilize and damage pirin structure, substantially reduce its affinity to its binding partners and alter pirin residue fluctuation profile via changing the flexibility of several key residues. Additionally, variants R59Q, F78V, G98D, V151D and L220P were found to impact pirin structure and function in multiple ways. As no haplotype was identified to be harboring more than one missense variant, further interrogation of the individual effects of these seven missense variants is highly recommended.

CONCLUSIONS

Pirin is involved in the transcriptional regulation of several genes and can play an important role in inflammatory responses. The variants predicted to be pathogenic in this study may thus contribute to a better understanding of the underlying molecular mechanisms of various inflammatory diseases. Future studies should be focused on clarifying if any of these variants can be used as disease biomarkers.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s42269-022-00917-7.

Large protein complex interfaces have evolved to promote cotranslational assembly

Assembly pathways of protein complexes should be precise and efficient to minimise misfolding and unwanted interactions with other proteins in the cell. One way to achieve this efficiency is by seeding assembly pathways during translation via the cotranslational assembly of subunits. While recent evidence suggests that such cotranslational assembly is widespread, little is known about the properties of protein complexes associated with the phenomenon. Here, using a combination of proteome-specific protein complex structures and publicly available ribosome profiling data, we show that cotranslational assembly is particularly common between subunits that form large intermolecular interfaces. To test whether large interfaces have evolved to promote cotranslational assembly, as opposed to cotranslational assembly being a non-adaptive consequence of large interfaces, we compared the sizes of first and last translated interfaces of heteromeric subunits in bacterial, yeast, and human complexes. When considering all together, we observe the N-terminal interface to be larger than the C-terminal interface 54% of the time, increasing to 64% when we exclude subunits with only small interfaces, which are unlikely to cotranslationally assemble. This strongly suggests that large interfaces have evolved as a means to maximise the chance of successful cotranslational subunit binding.

Finding alternatives to 5-fluorouracil: application of ensemble-based virtual screening for drug repositioning against human thymidylate synthase

5-fluorouracil and analogs are used in the treatment of many solid tumours. However, there are many cases of resistance and high toxicity associated with 5-fluorouracil chemotherapy. Repurposing FDA drugs against human thymidylate synthase revealed a number of FDA drugs that have a potential to be further developed for the treatment of various cancers for which 5-fluorouracil and analogs have been used for chemotherapy. Four FDA drugs prioritized for further validation included Erismodegib, Irinotecan, Conivaptan and Ergotamine. The role of water in mediating drug interactions and its contribution to the total binding energy was also shown. MM-PBSA calculations revealed that the binding affinity was the lowest for the hTS-Ergotamine complex (-66.702 ± 1.807 kJ/mol) suggesting moderate inhibition despite a large energetic contribution from van der Waal interactions (-190.889 ± 1.027 kJ/mol).Communicated by Ramaswamy H. Sarma.

Allosteric control of ACE2 peptidase domain dynamics

The Angiotensin Converting Enzyme 2 (ACE2) assists the regulation of blood pressure and is the main target of the coronaviruses responsible for SARS and COVID19. The catalytic function of ACE2 relies on the opening and closing motion of its peptidase domain (PD). In this study, we investigated the possibility of allosterically controlling the ACE2 PD functional dynamics. After confirming that ACE2 PD binding site opening-closing motion is dominant in characterizing its conformational landscape, we observed that few mutations in the viral receptor binding domain fragments were able to impart different effects on the binding site opening of ACE2 PD. This showed that binding to the solvent exposed area of ACE2 PD can effectively alter the conformational profile of the protein, and thus likely its catalytic function. Using a targeted machine learning model and relative entropy-based statistical analysis, we proposed the mechanism for the allosteric perturbation that regulates the ACE2 PD binding site dynamics at atomistic level. The key residues and the source of the allosteric regulation of ACE PD dynamics are also presented.

Vaccination with SARS-CoV-2 spike protein lacking glycan shields elicits enhanced protective responses in animal models

A major challenge to end the pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is to develop a broadly protective vaccine that elicits long-term immunity. As the key immunogen, the viral surface spike (S) protein is frequently mutated, and conserved epitopes are shielded by glycans. Here, we revealed that S protein glycosylation has site-differential effects on viral infectivity. We found that S protein generated by lung epithelial cells has glycoforms associated with increased infectivity. Compared to the fully glycosylated S protein, immunization of S protein with N-glycans trimmed to the mono-GlcNAc-decorated state (S_MG) elicited stronger immune responses and better protection for human angiotensin-converting enzyme 2 (hACE2) transgenic mice against variants of concern (VOCs). In addition, a broadly neutralizing monoclonal antibody was identified from S_MG-immunized mice that could neutralize wild-type SARS-CoV-2 and VOCs with subpicomolar potency. Together, these results demonstrate that removal of glycan shields to better expose the conserved sequences has the potential to be an effective and simple approach for developing a broadly protective SARS-CoV-2 vaccine.

Electronic Polarization Is Essential for the Stabilization and Dynamics of Buried Ion Pairs in Staphylococcal Nuclease Mutants

Buried charged residues play important roles in the modulation of protein stabilities and conformational dynamics and make crucial contributions to protein functions. Considering the generally nonpolar nature of protein interior, a key question concerns the contribution of electronic polarization to the stabilization and properties of buried charges. We answer this question by conducting free energy simulations using the latest polarizable CHARMM force field based on Drude oscillators for a series of Staphylococcal nuclease mutants that involve a buried Glu-Lys pair in different titration states and orientations. While a nonpolarizable model suggests that the ionized form of the buried Glu-Lys pair is more than 40 kcal/mol less stable than the charge-neutral form, the two titration states are comparable in stability when electronic polarization is included explicitly, a result better reconcilable with available experimental data. Analysis of free energy components suggests that additional stabilization of the ionized Glu-Lys pair has contributions from both the enhanced salt-bridge strength and stronger interaction between the ion-pair and surrounding protein residues and penetrated water. Despite the stronger direct interaction between Glu and Lys, the ion-pair exhibits considerably larger and faster structural fluctuations when polarization is included, due to compensation of interactions in the cavity. Collectively, observations from this work provide compelling evidence that electronic polarization is essential to the stability, hydration, dynamics, and therefore function of buried charges in proteins. Therefore, our study advocates for the explicit consideration of electronic polarization for mechanistic and engineering studies that implicate buried charged residues, such as enzymes and ion transporters.

Glycosite-deleted mRNA of SARS-CoV-2 spike protein as a broad-spectrum vaccine

Development of the messenger RNA (mRNA) vaccine has emerged as an effective and speedy strategy to control the spread of new pathogens. After vaccination, the mRNA is translated into the real protein vaccine, and there is no need to manufacture the protein in vitro. However, the fate of mRNA and its posttranslational modification inside the cell may affect immune response. Here, we showed that the mRNA vaccine of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein with deletion of glycosites in the receptor-binding domain (RBD) or especially the subunit 2 (S2) domain to expose more conserved epitopes elicited stronger antibody and CD8+ T cell responses with broader protection against the alpha, beta, gamma, delta, and omicron variants, compared to the unmodified mRNA. Immunization of such mRNA resulted in accumulation of misfolded spike protein in the endoplasmic reticulum, causing the up-regulation of BiP/GRP78, XBP1, and p-eIF2α to induce cell apoptosis and strong CD8+ T cell response. In addition, dendritic cells (DCs) incubated with S2-glysosite deleted mRNA vaccine increased class I major histocompatibility complex (MHC I) expression. This study provides a direction for the development of broad-spectrum mRNA vaccines which may not be achieved with the use of expressed proteins as antigens.

Computational Simulations Highlight the IL2Rα Binding Potential of Polyphenol Stilbenes from Fenugreek

Widely used in global households, fenugreek is well known for its culinary and medicinal uses. The various reported medicinal properties of fenugreek are by virtue of the different natural phytochemicals present in it. Regarded as a promising target, interleukin 2 receptor subunit alpha (IL2Rα) has been shown to influence immune responses. In the present research, using in silico techniques, we have demonstrated the potential IL2Rα binding properties of three polyphenol stilbenes (desoxyrhaponticin, rhaponticin, rhapontigenin) from fenugreek. As the first step, molecular docking was performed to assess the binding potential of the fenugreek phytochemicals with IL2Rα. All three phytochemicals demonstrated interactions with active site residues. To confirm the reliability of our molecular docking results, 100 ns molecular dynamics simulations studies were undertaken. As discerned by the RMSD and RMSF analyses, IL2Rα in complex with the desoxyrhaponticin, rhaponticin, and rhapontigenin indicated stability. The RMSD analysis of the phytochemicals alone also demonstrated no significant structural changes. Based on the stable molecular interactions and comparatively slightly better MM/PBSA binding free energy, rhaponticin seems promising. Additionally, ADMET analysis performed for the stilbenes indicated that all of them obey the ADMET rules. Our computational study thus supports further in vitro IL2Rα binding studies on these stilbenes, especially rhaponticin.

Computational prediction of the potential target of SARS-CoV-2 inhibitor plitidepsin via molecular docking, dynamic simulations and MM-PBSA calculations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication depends on the interaction between the viral proteins and the human translation machinery. The cytotoxic peptide plitidepsin was found to inhibit CoV-2 up to 90% at a concentration of 0.88 nM. In vitro studies suggest that this activity may be attributed to the inhibition of the eukaryotic translation elongation factor 1A (eEF1A). However, recent reports raised the potential for other cellular targets which plitidepsin may use to exert its potent antiviral activity. The lack of data about these potential targets represents a major limitation for its structural optimization. This work describes the use of a molecular modeling approach to rationalize the in vitro antiviral activity of plitidepsin and to identify potential cellular targets. The developed protocol involves an initial molecular docking step followed by molecular dynamics and binding free energy calculations. The results reveal the potential for plitidepsin to bind to the active site of the key enzyme SARS-CoV-2 RdRp. The results also highlight the importance of van der Waals interactions for proper binding with the enzyme. We believe that the results presented in this study could provide the grounds for the optimization of plitidepsin analogues as SARS-CoV-2 inhibitors.

Effect of monovalent salt concentration and peptide secondary structure in peptide-micelle binding

Recently, we reported a cationic 14 residue peptide LL-14 (LKWLKKLLKWLKKL) with salt-sensitive broad-spectrum antimicrobial potency. However, the mechanism of its salt (NaCl) sensitivity remained unclear. In this study, we have reported computational (∼14.2 μs of MD) and experimental (CD, fluorescence) investigations to examine the salt-sensitivity and the role of peptide secondary structure on LL-14 binding to simple membrane mimetic (SDS, DPC) systems. LL-14 was shown to adopt a random coil (P^c) conformation in water and α-helical conformation (P^h) in the peptide:SDS micelle complex, accompanied by tryptophan burial, using both simulations and experiments. Simulations successfully deconvoluted the LL-14:micelle binding event in terms of secondary structure (random coil P^cversus helix P^h) and gave atomic insight into the initial and final LL-14:SDS complexes. Electrostatics drove the N-terminus (L1 and K2) of LL-14 (P^c or P^h) to bind the SDS micellar surface, initiating complex formation. LL-14 in amphipathic P^h conformation bound faster and buried deeper into the SDS micelle relative to P^c. Increasing NaCl concentration incrementally delayed LL-14:micelle binding by shielding the overall charges of the interacting partners. LL-14 binding to the SDS micelle was significantly faster relative to that of the zwitterionic DPC micelle due to electrostatic reasons. Cationic α-helical amphipathic peptides (with positively charged N-terminus) with low salt-ion concentration seemed to be ideal for faster SDS binding.

We report computational (∼14.2 μs of MD) and experimental (CD, fluorescence) investigations to examine the salt-sensitivity and the role of the peptide secondary structure on LL-14 binding to simple membrane mimetic systems.

Interaction and Inhibition of a Ganoderma lucidum Proteoglycan on PTP1B Activity for Anti-diabetes

Protein tyrosine phosphatase 1B (PTP1B) is a key negative regulator of insulin and an effective target for the treatment of type 2 diabetes (T2D). A natural hyperbranched proteoglycan extracted from Ganoderma lucidum, namely, Fudan-Yueyang G. Lucidum (FYGL), was demonstrated capable of inhibiting the activity of PTP1B. Here, to identify the effective active components of FYGL, three different components, the polysaccharide FYGL-1, proteoglycans FYGL-2, and FYGL-3, were isolated from FYGL, and then, the protein moiety of FYGL-3 was further separated, namely, FYGL-3-P. Their abilities to enhance the glucose uptake in cells and inhibit the activity of PTP1B were compared. The inhibitory mechanisms were systematically explored by spectroscopic methods and MD simulations. The results showed that FYGL-3 and FYGL-3-P significantly enhanced the insulin-provoked glucose uptake in insulin-resistant HepG2 cells, detected by the glucose oxidase method. Also, the FYGL-3-P protein moiety in FYGL played an essential role in inhibiting the activity of PTP1B. A strong, enthalpy-driven, and multitargeted interaction by electrostatic forces between PTP1B and FYGL-3-P dramatically inhibited the catalytic activity of PTP1B. These results provided deep insights into the molecular mechanisms of FYGL inhibiting the activity of PTP1B and structurally helped researchers seek natural PTP1B inhibitors.

Broadening the Mutation Spectrum in GJA8 and CHMP4B: Novel Missense Variants and the Associated Phenotypes in Six Chinese Han Congenital Cataracts Families

Purpose: To broaden the mutation and phenotype spectrum of the GJA8 and CHMP4B genes and to reveal genotype-phenotype correlations in a cohort of Chinese patients with congenital cataracts (CCs).

Methods: Six Chinese Han families with CCs inherited in an autosomal dominant (AD) pattern were recruited for this study. All patients underwent full ocular examinations. Genomic DNA was extracted from the leukocytes of peripheral blood collected from all available patients and their unaffected family members. Whole-exome sequencing (WES) was performed on all probands and at least one of their parents. Candidate variants were further confirmed by Sanger sequencing. Bioinformatic analysis with several computational predictive programs was performed to assess the impacts of the candidate variants on the structure and function of the proteins.

Results: Four heterozygous candidate variants in three different genes (CRYBB2, GJA8, and CHMP4B) were identified in affected individuals from the six families, including two novel missense variants (GJA8: c.64G > C/p. G22R, and CHMP4B: c.587C > G/p. S196C), one missense mutation (CRYBB2: c.562C > T/p. R188C), and one small deletion (GJA8: c.426_440delGCTGGAGGGGACCCT/p.143_147delLEGTL). The three missense mutations were predicted as deleterious in all four computational prediction programs. In the homologous model, the GJA8: p.143_147delLEGTL mutation showed a sequence deletion of five amino acids at the cytoplasmic loop of the Cx50 protein, close to the third transmembrane domain. Patients carrying mutations in the same gene showed similar cataract phenotypes at a young age, including total cataracts, Y-sutural with fetal nuclear cataracts, and subcapsular cataracts.

Conclusion: This study further expands the mutation spectrum and genotype-phenotype correlation of CRYBB2, GJA8, and CHMP4B underlying CCs. This study sheds light on the importance of comparing congenital cataract phenotypes in patients at the same age stage. It offers clues for the pathogenesis of CCs and allows for an early prenatal diagnosis for families carrying these genetic variants.

Mechanistic insights into COVID-19 by global analysis of the SARS-CoV-2 3CLpro substrate degradome

The main viral protease (3CLpro) is indispensable for SARS-CoV-2 replication. We delineate the human protein substrate landscape of 3CLpro by TAILS substrate-targeted N-terminomics. We identify more than 100 substrates in human lung and kidney cells supported by analyses of SARS-CoV-2-infected cells. Enzyme kinetics and molecular docking simulations of 3CLpro engaging substrates reveal how noncanonical cleavage sites, which diverge from SARS-CoV, guide substrate specificity. Cleaving the interactors of essential effector proteins, effectively stranding them from their binding partners, amplifies the consequences of proteolysis. We show that 3CLpro targets the Hippo pathway, including inactivation of MAP4K5, and key effectors of transcription, mRNA processing, and translation. We demonstrate that Spike glycoprotein directly binds galectin-8, with galectin-8 cleavage disengaging CALCOCO2/NDP52 to decouple antiviral-autophagy. Indeed, in post-mortem COVID-19 lung samples, NDP52 rarely colocalizes with galectin-8, unlike in healthy lungs. The 3CLpro substrate degradome establishes a foundational substrate atlas to accelerate exploration of SARS-CoV-2 pathology and drug design.

The Molecular Basis for Escherichia coli O157:H7 Phage FAHEc1 Endolysin Function and Protein Engineering to Increase Thermal Stability

Bacteriophage-encoded endolysins have been identified as antibacterial candidates. However, the development of endolysins as mainstream antibacterial agents first requires a comprehensive biochemical understanding. This study defines the atomic structure and enzymatic function of Escherichia coli O157:H7 phage FAHEc1 endolysin, LysF1. Bioinformatic analysis suggests this endolysin belongs to the T4 Lysozyme (T4L)-like family of proteins and contains a highly conserved catalytic triad. We then solved the structure of LysF1 with x-ray crystallography to 1.71 Å. LysF1 was confirmed to exist as a monomer in solution by sedimentation velocity experiments. The protein architecture of LysF1 is conserved between T4L and related endolysins. Comparative analysis with related endolysins shows that the spatial orientation of the catalytic triad is conserved, suggesting the catalytic mechanism of peptidoglycan degradation is the same as that of T4L. Differences in the sequence illustrate the role coevolution may have in the evolution of this fold. We also demonstrate that by mutating a single residue within the hydrophobic core, the thermal stability of LysF1 can be increased by 9.4 °C without compromising enzymatic activity. Overall, the characterization of LysF1 provides further insight into the T4L-like class of endolysins. Our study will help advance the development of related endolysins as antibacterial agents, as rational engineering will rely on understanding mutable positions within this protein fold.

Deciphering Molecular Virulence Mechanism of Mycobacterium tuberculosis Dop isopeptidase Based on Its Sequence-Structure-Function Linkage

The pupylation pathway marks proteins for prokaryotic ubiquitin-like protein (Pup)-proteasomal degradation and survival strategy of mycobacteria inside of the host macrophages. Deamidase of Pup (Dop) plays a central role in the pupylation pathway. It is still a matter of investigation to know the function of Dop in virulence of mycobacterial lineage. Hence, the present study was intended to describe the sequence-structure-function-virulence link of Dop for understanding the molecular virulence mechanism of Mycobacterium tuberculosis H37Rv (Mtb). Phylogenetic analysis of this study indicated that Dop has extensively diverged across the proteasome-harboring bacteria. The functional part of Dop was converged across the pathogenic mycobacterial lineage. The genome-wide analysis pointed out that the pupylation gene locus was identical to each other, but its genome neighborhood differed from species to species. Molecular modeling and dynamic studies proved that the predicted structure of Mtb Dop was energetically stable and low conformational freedom. Moreover, evolutionary constraints in Mtb Dop were intensively analyzed for inferring its sequence-structure-function relationships for the full virulence of Mtb. It indicated that evolutionary optimization was extensively required to stabilize its local structural environment at the side chains of mutable residues. The sequence-structure-function-virulence link of Dop might have retained in Mtb by reordering hydrophobic and hydrogen bonding patterns in the local structural environment. Thus, the results of our study provide a quest to understand the molecular virulence and pathogenesis mechanisms of Mtb during the infection process.

LONP1 de novo dominant mutation causes mitochondrial encephalopathy with loss of LONP1 chaperone activity and excessive LONP1 proteolytic activity

LONP1 is an ATP-dependent protease and chaperone that plays multiple vital roles in mitochondria. LONP1 is essential for mitochondrial homeostasis due to its role in maintenance of the mitochondrial genome and its central role in regulating mitochondrial processes such as oxidative phosphorylation, mitophagy, and heme biosynthesis. Bi-allelic LONP1 mutations have been reported to cause a constellation of clinical presentations. We report a patient heterozygous for a de novo mutation in LONP1: c.901C>T,p.R301W presenting as a neonate with seizures, encephalopathy, pachygyria and microcephaly. Assays of respiratory chain activity in muscle showed complex II-III function at 8% of control. Functional studies in patient fibroblasts showed a signature of dysfunction that included significant decreases in known proteolytic targets of LONP1 (TFAM, PINK1, phospho-PDH E1α) as well as loss of mitochondrial ribosome subunits MRPL44 and MRPL11 with concomitant decreased activity and level of protein subunits of oxidative phosphorylation complexes I and IV. These results indicate excessive LONP1 proteolytic activity and a loss of LONP1 chaperone activity. Further, we demonstrate that the LONP1 N-terminal domain is involved in hexamer stability of LONP1 and that the ability to make conformational changes is necessary for LONP1 to regulate proper functioning of both its proteolytic and chaperone activities.

Investigating the stability of dengue virus envelope protein dimer using well-tempered metadynamics simulations

We explored the stability of the dengue virus envelope (E) protein dimer since it is widely assumed that the E protein dimer is stabilized by drug ligands or antibodies in an acidic environment, neutralizing the virus's ability to fuse with human cells. During this process, a large conformational change of the E protein dimer is required. We performed Molecular Dynamics simulations to mimic the conformational change and stability of the dimer in neutral and acidic conditions with the well-tempered metadynamics method. Furthermore, as a few neutralizing antibodies discovered from dengue patients were reported, we used the same simulation method to examine the influence of a selected antibody on the dimer stability in both neutral and acidic conditions. We also investigated the antibody's influence on a point-mutated E protein that had been reported to interrupt the protein-antibody interaction and result in more than 95% loss of the antibody's binding ability. Our simulation results are highly consistent with the experimental conclusion that binding of the antibody to the E protein dimer neutralizes the virus, especially in a low pH condition, while the mutation of W101A or N153A significantly reduces the antibody's ability in stabilizing the E protein dimer. We demonstrate that well-tempered metadynamics can be used to accurately explore the antibody's interaction on large protein complexes such as the E protein dimer, and the computational approach in this work is promising in future antibody development.

Substrate-imprinted docking of Agrobacterium tumefaciens uronate dehydrogenase for increased substrate selectivity

Agrobacterium tumefaciens uronate dehydrogenase (AtuUdh) belongs to the short-chain dehydrogenase superfamily, specifically those acting on the CH-OH group of donor with NAD⁺ or NADP⁺ as acceptor. It is apparently required for the production of D-glucaric acid. AtuUdh-catalyzed reaction is reversible with dual substrate-specific activity (D-galacturonic acid and D-glucuronic acid) in nature. In our study, 34 mutants were pre-screened from 155 mutants generated from AtuUdh (wild-type) and selected 10 structurally stable mutants with increased substrate selectivity. The specificity, efficiency, and selectivity of these mutants for different substrates and cofactors were predicted from 121 docked models using a substrate-imprinted docking approach. Q14F, S36L, and S75T mutants have shown a high binding affinity to D-glucuronic acid and its substrate intermediates such as D-glucaro-1,4-lactone and D-glucaro-1,5-lactone. These mutants exhibited a low binding affinity to the substrate and cofactor required for D-galactaric acid. D34S, N112E and S165E mutants found to show a high selectivity of D-galacturonic acid and its substrate intermediates for D-galactaric acid production. Ser75, Ser165, and Arg174 are active residues playing an imperative role in the substrate selectivity and also contributed in the conjecture the mechanism of transition state stabilization catalyzed by AtuUdh mutants. The present approach was used to reveal the substrate binding mechanism of AtuUdh mutants for a better understanding of the structural basis for selectivity and function.

The Effect of Electrostatic Interactions on the Folding Kinetics of a 3-α-Helical Bundle Protein Family

The trio of protein segment repeats called spectrins diverges by more than 2 orders of magnitude in their folding and unfolding rates, despite having very similar stabilities and almost coincidental topologies. Experimental studies revealed that the mutation of five particular residues dramatically alters the kinetic rates in the slow folders, making them similar to the rates of the fast folder. This is considered to be an exceptional behavior which seems in principle to challenge the current understanding of the protein folding process. In this work, we analyze this scenario, using a simplified computational model, combined with state-of-the-art kinetic analysis techniques. Our model faithfully separates the kinetics of the fast and slow folders and captures the effect of the five mutations. We show that the inclusion of electrostatics in the model is necessary to explain the experimental findings.

Substitutions of a buried glutamate residue hinder the conformational change in horse liver alcohol dehydrogenase and yield a surprising complex with endogenous 3'-Dephosphocoenzyme A

Glu-267 is highly conserved in alcohol dehydrogenases and buried as a negatively-charged residue in a loop of the NAD coenzyme binding domain. Glu-267 might have a structural role and contribute to a rate-promoting vibration that facilitates catalysis. Substitutions of Glu-267 with histidine or asparagine residues increase the dissociation constants for the coenzymes (NAD⁺ by ∼40-fold, NADH by ∼200-fold) and significantly decrease catalytic efficiencies by 16-1200-fold various substrates and substituted enzymes. The turnover numbers modestly change with the substitutions, but hydride transfer is at least partially rate-limiting for turnover for alcohol oxidation. X-ray structures of the E267H and E267 N enzymes are similar to the apoenzyme (open) conformation of the wild-type enzyme, and the substitutions are accommodated by local changes in the structure. Surprisingly, the E267H and E267 N enzymes have endogenous (from the expression in E. coli) 3'-dephosphocoenzyme A bound in the active site with the ADP moiety in the NAD binding site and the pantethiene sulfhydryl bound to the catalytic zinc. The kinetics and crystallography show that the substitutions of Glu-267 hinder the conformational change, which occurs when wild-type enzyme binds coenzymes, and affect productive binding of substrates.

Immunomodulating protein aggregates in soy and whey hydrolysates and their resistance to digestion in an in vitro infant gastrointestinal model: new insights in the mechanism of immunomodulatory hydrolysates

Hydrolysates, which are used in hypoallergenic infant formulas, have been found to possess immune modulating effects. For an optimal utilization of hydrolysates, the working mechanisms and responsible proteins underlying the effects should be elucidated. In this study, the immunomodulating activity of whey and soy hydrolysates was studied by quantifying TLR activation and assessing cytokine production in hydrolysate stimulated dendritic cells. The responsible protein fraction was identified and characterized by gel electrophoresis. The immune effects under gastrointestinal conditions were studied by digesting the hydrolysates in an in vitro infant digestion model, after which the digests were analyzed. In both soy and whey hydrolysates, TLR activation and cytokine production in dendritic cells were induced by a fraction containing protein aggregates larger than 1000 kDa, which were formed by electrostatic interactions and disulfide bonds. Only soy aggregates remained intact during duodenal digestion, and maintained the TLR activating capacity. Soy and whey protein aggregates larger than 1000 kDa possess immunomodulatory properties, but only soy aggregates remain under intestinal digestion conditions. This knowledge is important for a better understanding of the effects of hydrolysates.

Differential protein structural disturbances and suppression of assembly partners produced by nonsense GABRG2 epilepsy mutations: implications for disease phenotypic heterogeneity

Mutations in GABA_A receptor subunit genes are frequently associated with epilepsy, and nonsense mutations in GABRG2 are associated with several epilepsy syndromes including childhood absence epilepsy, generalized tonic clonic seizures and the epileptic encephalopathy, Dravet syndrome. The molecular basis for the phenotypic heterogeneity of mutations is unclear. Here we focused on three nonsense mutations in GABRG2 (GABRG2(R136*), GABRG2(Q390*) and GABRG2(W429*)) associated with epilepsies of different severities. Structural modeling and structure-based analysis indicated that the surface of the wild-type γ2 subunit was naturally hydrophobic, which is suitable to be buried in the cell membrane. Different mutant γ2 subunits had different stabilities and different interactions with their wild-type subunit binding partners because they adopted different conformations and had different surface hydrophobicities and different tendency to dimerize. We utilized flow cytometry and biochemical approaches in combination with lifted whole cell patch-clamp recordings. We demonstrated that the truncated subunits had no to minimal surface expression and unchanged or reduced surface expression of wild-type partnering subunits. The amplitudes of GABA-evoked currents from the mutant α1β2γ2(R136*), α1β2γ2(Q390*) and α1β2γ2(W429*) receptors were reduced compared to the currents from α1β2γ2 receptors but with differentially reduced levels. This thus suggests differential protein structure disturbances are correlated with disease severity.

Structural constraints-based evaluation of immunogenic avirulent toxins from Clostridium botulinum C2 and C3 toxins as subunit vaccines

Clostridium botulinum (group-III) is an anaerobic bacterium producing C2 and C3 toxins in addition to botulinum neurotoxins in avian and mammalian cells. C2 and C3 toxins are members of bacterial ADP-ribosyltransferase superfamily, which modify the eukaryotic cell surface proteins by ADP-ribosylation reaction. Herein, the mutant proteins with lack of catalytic and pore forming function derived from C2 (C2I and C2II) and C3 toxins were computationally evaluated to understand their structure-function integrity. We have chosen many structural constraints including local structural environment, folding process, backbone conformation, conformational dynamic sub-space, NAD-binding specificity and antigenic determinants for screening of suitable avirulent toxins. A total of 20 avirulent mutants were identified out of 23 mutants, which were experimentally produced by site-directed mutagenesis. No changes in secondary structural elements in particular to α-helices and β-sheets and also in fold rate of all-β classes. Structural stability was maintained by reordered hydrophobic and hydrogen bonding patterns. Molecular dynamic studies suggested that coupled mutations may restrain the binding affinity to NAD(+) or protein substrate upon structural destabilization. Avirulent toxins of this study have stable energetic backbone conformation with a common blue print of folding process. Molecular docking studies revealed that avirulent mutants formed more favorable hydrogen bonding with the side-chain of amino acids near to conserved NAD-binding core, despite of restraining NAD-binding specificity. Thus, structural constraints in the avirulent toxins would determine their immunogenic nature for the prioritization of protein-based subunit vaccine/immunogens to avian and veterinary animals infected with C. botulinum.

Multiple Surface Regions on the Niemann-Pick C2 Protein Facilitate Intracellular Cholesterol Transport

The cholesterol storage disorder Niemann-Pick type C (NPC) disease is caused by defects in either of two late endosomal/lysosomal proteins, NPC1 and NPC2. NPC2 is a 16-kDa soluble protein that binds cholesterol in a 1:1 stoichiometry and can transfer cholesterol between membranes by a mechanism that involves protein-membrane interactions. To examine the structural basis of NPC2 function in cholesterol trafficking, a series of point mutations were generated across the surface of the protein. Several NPC2 mutants exhibited deficient sterol transport properties in a set of fluorescence-based assays. Notably, these mutants were also unable to promote egress of accumulated intracellular cholesterol from npc2(-/-) fibroblasts. The mutations mapped to several regions on the protein surface, suggesting that NPC2 can bind to more than one membrane simultaneously. Indeed, we have previously demonstrated that WT NPC2 promotes vesicle-vesicle interactions. These interactions were abrogated, however, by mutations causing defective sterol transfer properties. Molecular modeling shows that NPC2 is highly plastic, with several intense positively charged regions across the surface that could interact favorably with negatively charged membrane phospholipids. The point mutations generated in this study caused changes in NPC2 surface charge distribution with minimal conformational changes. The plasticity, coupled with membrane flexibility, probably allows for multiple cholesterol transfer routes. Thus, we hypothesize that, in part, NPC2 rapidly traffics cholesterol between closely appositioned membranes within the multilamellar interior of late endosomal/lysosomal proteins, ultimately effecting cholesterol egress from this compartment.

Structure of the human MLH1 N-terminus: implications for predisposition to Lynch syndrome

Mismatch repair prevents the accumulation of erroneous insertions/deletions and non-Watson-Crick base pairs in the genome. Pathogenic mutations in the MLH1 gene are associated with a predisposition to Lynch and Turcot's syndromes. Although genetic testing for these mutations is available, robust classification of variants requires strong clinical and functional support. Here, the first structure of the N-terminus of human MLH1, determined by X-ray crystallography, is described. The structure shares a high degree of similarity with previously determined prokaryotic MLH1 homologs; however, this structure affords a more accurate platform for the classification of MLH1 variants.

Carboxyl-peptide plane stacking is important for stabilization of buried E305 of Trichoderma reesei Cel5A

Hydrogen bonds or salt bridges are usually formed to stabilize the buried ionizable residues. However, such interactions do not exist for two buried residues D271 and E305 of Trichoderma reesei Cel5A, an endoglucanase. Mutating D271 to alanine or leucine improves the enzyme thermostability quantified by the temperature T50 due to the elimination of the desolvation penalty of the aspartic acid. However, the same mutations for E305 decrease the enzyme thermostability. Free energy calculations based on the molecular dynamics simulation predict the thermostability of D271A, D271L, and E305A (compared to WT) in line with the experimental observation but overestimate the thermostability of E305L. Quantum mechanical calculations suggest that the carboxyl-peptide plane stacking interactions occurring to E305 but not D271 are important for the carboxyl group stabilization. For the protonated carboxyl group, the interaction energy can be as much as about -4 kcal/mol for parallel stacking and about -7 kcal/mol for T-shaped stacking. For the deprotonated carboxyl group, the largest interaction energies for parallel stacking and T-shaped stacking are comparable, about -7 kcal/mol. The solvation effect generally weakens the interaction, especially for the charged system. A search of the carboxyl-peptide plane stacking in the PDB databank indicates that parallel stacking but not T-shaped stacking is quite common, and the most probable distance between the two stacking fragments is close to the value predicted by the QM calculations. This work highlights the potential role of carboxyl amide π-π stacking in the stabilization of aspartic acid and glutamic acid in proteins.

Driving Forces in Pressure-Induced Protein Transitions

The molecular mechanisms underlying pressure-induced protein denaturation can be analyzed based on the pressure-dependent differences in the apparent volume occupied by amino acids inside the protein and when exposed to water in an unfolded conformation. This chapter presents a volumetric analysis of the peptide group and the 20 naturally occurring amino acid side chains in the interior of the native state, the micelle-like interior of the pressure-induced denatured state, and in the unfolded conformation modeled by low-molecular analogs of proteins. The transfer of a peptide group from the protein interior to water becomes increasingly favorable as pressure increases. This observation classifies solvation of peptide groups as a major driving force in pressure-induced protein denaturation. Polar side chains do not appear to exhibit significant pressure-dependent changes in their preference for the protein interior or solvent. The transfer of nonpolar side chains from the protein interior to water becomes more unfavorable as pressure increases. An inference can be drawn that a sizeable population of nonpolar side chains remains buried inside a solvent-inaccessible core of the pressure-induced denatured state. At elevated pressures this core, owing to the absence of structural constraints, may become packed almost as tightly as the interior of the native state. The presence and partial disappearance of large intraglobular voids is another driving force facilitating pressure-induced protein denaturation. Volumetric data presented here have implications for the kinetics of protein folding and shed light on the nature of the folding transition state ensembles.

The structure of vanin 1: a key enzyme linking metabolic disease and inflammation

Although part of the coenzyme A pathway, vanin 1 (also known as pantetheinase) sits on the cell surface of many cell types as an ectoenzyme, catalyzing the breakdown of pantetheine to pantothenic acid (vitamin B5) and cysteamine, a strong reducing agent. Vanin 1 was initially discovered as a protein involved in the homing of leukocytes to the thymus. Numerous studies have shown that vanin 1 is involved in inflammation, and more recent studies have shown a key role in metabolic disease. Here, the X-ray crystal structure of human vanin 1 at 2.25 Å resolution is presented, which is the first reported structure from the vanin family, as well as a crystal structure of vanin 1 bound to a specific inhibitor. These structures illuminate how vanin 1 can mediate its biological roles by way of both enzymatic activity and protein-protein interactions. Furthermore, it sheds light on how the enzymatic activity is regulated by a novel allosteric mechanism at a domain interface.

Forces stabilizing proteins

The goal of this article is to summarize what has been learned about the major forces stabilizing proteins since the late 1980s when site-directed mutagenesis became possible. The following conclusions are derived from experimental studies of hydrophobic and hydrogen bonding variants. (1) Based on studies of 138 hydrophobic interaction variants in 11 proteins, burying a -CH2- group on folding contributes 1.1±0.5 kcal/mol to protein stability. (2) The burial of non-polar side chains contributes to protein stability in two ways: first, a term that depends on the removal of the side chains from water and, more importantly, the enhanced London dispersion forces that result from the tight packing in the protein interior. (3) Based on studies of 151 hydrogen bonding variants in 15 proteins, forming a hydrogen bond on folding contributes 1.1±0.8 kcal/mol to protein stability. (4) The contribution of hydrogen bonds to protein stability is strongly context dependent. (5) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (6) Polar group burial can make a favorable contribution to protein stability even if the polar group is not hydrogen bonded. (7) Hydrophobic interactions and hydrogen bonds both make large contributions to protein stability.

A coarse-grained protein model in a water-like solvent

Simulations employing an explicit atom description of proteins in solvent can be computationally expensive. On the other hand, coarse-grained protein models in implicit solvent miss essential features of the hydrophobic effect, especially its temperature dependence, and have limited ability to capture the kinetics of protein folding. We propose a free space two-letter protein (“H-P”) model in a simple, but qualitatively accurate description for water, the Jagla model, which coarse-grains water into an isotropically interacting sphere. Using Monte Carlo simulations, we design protein-like sequences that can undergo a collapse, exposing the “Jagla-philic” monomers to the solvent, while maintaining a “hydrophobic” core. This protein-like model manifests heat and cold denaturation in a manner that is reminiscent of proteins. While this protein-like model lacks the details that would introduce secondary structure formation, we believe that these ideas represent a first step in developing a useful, but computationally expedient, means of modeling proteins.

Protons as second messenger regulators of G protein signaling

In response to environmental stress, cells often generate pH signals that serve to protect vital cellular components and reprogram gene expression for survival. A major barrier to our understanding of this process has been the identification of signaling proteins that detect changes in intracellular pH. To identify candidate pH sensors, we developed a computer algorithm that searches proteins for networks of proton-binding sidechains. This analysis indicates that Gα subunits, the principal transducers of G protein-coupled receptor (GPCR) signals, are pH sensors. Our structure-based calculations and biophysical investigations reveal that Gα subunits contain networks of pH-sensing sidechains buried between their Ras and helical domains. Further, we show that proton binding induces changes in conformation that promote Gα phosphorylation and suppress receptor-initiated signaling. Together, our computational, biophysical, and cellular analyses reveal an unexpected function for G proteins as mediators of stress-response signaling.

Peptide salt bridge stability: from gas phase via microhydration to bulk water simulations

The salt bridge formation and stability in the terminated lysine-glutamate dipeptide is investigated in water clusters of increasing size up to the limit of bulk water. Proton transfer dynamics between the acidic and basic side chains is described by DFT-based Born-Oppenheimer molecular dynamics simulations. While the desolvated peptide prefers to be in its neutral state, already the addition of a single water molecule can trigger proton transfer from the glutamate side chain to the lysine side chain, leading to a zwitterionic salt bridge state. Upon adding more water molecules we find that stabilization of the zwitterionic state critically depends on the number of hydrogen bonds between side chain termini, the water molecules, and the peptidic backbone. Employing classical molecular dynamics simulations for larger clusters, we observed that the salt bridge is weakened upon additional hydration. Consequently, long-lived solvent shared ion pairs are observed for about 30 water molecules while solvent separated ion pairs are found when at least 40 or more water molecules hydrate the dipeptide. These results have implications for the formation and stability of salt bridges at partially dehydrated surfaces of aqueous proteins.

Protein-protein interactions: general trends in the relationship between binding affinity and interfacial buried surface area

Protein-protein interactions play key roles in many cellular processes and their affinities and specificities are finely tuned to the functions they perform. Here, we present a study on the relationship between binding affinity and the size and chemical nature of protein-protein interfaces. Our analysis focuses on heterodimers and includes curated structural and thermodynamic data for 113 complexes. We observe a direct correlation between binding affinity and the amount of surface area buried at the interface. For a given amount of surface area buried, the binding affinity spans four orders of magnitude in terms of the dissociation constant (Kd ). Across the entire dataset, we observe no obvious relationship between binding affinity and the chemical composition of the interface. We also calculate the free energy per unit surface area buried, or "surface energy density," of each heterodimer. For interfacial surface areas between 500 and 2000 Å(2) , the surface energy density decreases as the buried surface area increases. As the buried surface area increases beyond about 2000 Å(2) , the surface energy density levels off to a constant value. We believe that these analyses and data will be useful for researchers with an interest in understanding, designing or inhibiting protein-protein interfaces.