Nonpolar Residues Research Articles

Impact of F80M and F83M mutations on the functionality of fluoride ion channel elucidated in microsecond level molecular dynamic simulation

Fluoride ion channels of the Fluc family plays a critically important role in combating environmental fluoride toxicity. As per the crystal structure of these fluoride ion channels, the pore region is densely packed with a series of hydrogen bond donating residues arranged in a ladder fashion creating an ion conducting pathway. In earlier studies, it was revealed that although the ion conducting pathway polarity is highly conserved, however the functionality of the channel protein depends on several residues at particular positions. While, a threonine at end of the pore is critically important in forming initial interactions, two phenylalanines at the central region coordinate F- transportation through the channel. It was also revealed that these two phenylalanines cannot be substituted by any other aromatic, polar or non-polar residues without hindering the functionality with exception of methionine. In another study, it was revealed that these two phenylalanines F80 and F83 when mutated with methionine; F80M lead to active state, while the F83M has lead to inactivity of F- anion conductivity. However, the exact atomic level detailing on how exactly these mutations have impacted the conductivity remained elusive. In this scenario, in this present study, we have modeled these two mutations and performed a microsecond level simulation on each mutation compared with wild type towards understanding the atomic level detailing revealing several insights on what exactly happening at these residues responsible for the selective conductivity of F- ions. Communicated by Ramaswamy H. Sarma

Quantifying the Long‐Range Coupling of Electronic Properties in Proteins with ab initio Molecular Dynamics**

AbstractThe delicate interplay of covalent and non‐covalent interactions in proteins is inherently quantum mechanical and highly dynamic in nature. To directly interrogate the evolving nature of the electronic structure of proteins, we carry out 100‐ps‐scale ab initio molecular dynamics simulations of three representative small proteins with range‐separated hybrid density functional theory. We quantify the nature and length‐scale of the coupling of residue‐specific charge probability distributions in these proteins. While some nonpolar residues exhibit expectedly narrow charge distributions, most polar and charged residues exhibit broad, multimodal distributions. Even for nonpolar residues, we observe sequence‐specific deviations corresponding to charge accumulation or depletion that would be challenging to capture in a fixed charge force field. We quantify the effect of residue‐residue interactions on charge distributions first with linear cross‐correlations. We then show how additional insight can be gained from evaluating the mutual information of charge distributions. We show that a significant number of residues couple most strongly with residues that are distant in both sequence and space over a range of secondary structures including α‐helical, β‐sheet, disulfide bridging, and lasso motifs. The mutual information analysis is necessary to capture coupling between some polar and charged residues that would be otherwise missed.

Computational Simulation of Holin S105 in Membrane Bilayer and Its Dimerization Through a Helix-Turn-Helix Motif.

During the final step of the bacteriophage infection cycle, the cytoplasmic membrane of host cells is disrupted by small membrane proteins called holins. The function of holins in cell lysis is carried out by forming a highly ordered structure called lethal lesion, in which the accumulation of holins in the cytoplasmic membrane leads to the sudden opening of a hole in the middle of this oligomer. Previous studies showed that dimerization of holins is a necessary step to induce their higher order assembly. However, the molecular mechanism underlying the holin-mediated lesion formation is not well understood. In order to elucidate the functions of holin, we first computationally constructed a structural model for our testing system: the holin S105 from bacteriophage lambda. All atom molecular dynamic simulations were further applied to refine its structure and study its dynamics as well as interaction in lipid bilayer. Additional simulations on association between two holins provide supportive evidence to the argument that the C-terminal region of holin plays a critical role in regulating the dimerization. In detail, we found that the adhesion of specific nonpolar residues in transmembrane domain 3 (TMD3) in a polar environment serves as the driven force of dimerization. Our study therefore brings insights to the design of binding interfaces between holins, which can be potentially used to modulate the dynamics of lesion formation.

N-Terminally Added Tag Selectively Enhances Heterologous Expression of VacA Cytotoxin Variants from Helicobacter pylori.

Gastric pathogen Helicobacter pylori secretes VacA cytotoxin displaying a high degree of polymorphic variations of which the highest VacA pathogenicity correlates with m1-type variant followed by VacA-m2. To comparatively evaluate expression in Escherichia coli of the mature VacA variants (m1- and m2-types) and their 33- and 55/59-kDa domains fused with His(6) tag at N- or C-terminus. All VacA clones expressed in E. coli TOP10™ were analyzed by SDS-PAGE and Western blotting. VacA inclusions were solubilized under native conditions (~150-rpm shaking at 37°C for 2 h in 20 mM HEPES (pH7.4) and 150 mM NaCl). Membrane-perturbing and cytotoxic activities of solubilized VacA proteins were assessed via liposome-entrapped dye leakage and resazurin- based cell viability assays, respectively. VacA binding to human gastric adenocarcinoma cells was assessed by immunofluorescence microscopy. Side-chain hydrophobicity of VacA was analyzed through modeled structures constructed by homology- and ab initio-based modeling. Both full-length VacA-m1 and 33-kDa domain were efficiently expressed only in the presence of N-terminal extension while its 55-kDa domain was capably expressed with either N- or Cterminal extension. Selectively enhanced expression was also observed for VacA-m2. Protein expression profiles revealed a critical period in IPTG-induced production of the 55-kDa domain with N-terminal extension unlike its C-terminal extension showing relatively stable expression. Both VacA- m1 isolated domains were able to independently bind to cultured gastric cells similar to the full- length toxin, albeit the 33-kDa domain exhibited significantly higher activity of membrane perturbation than others. Membrane-perturbing and cytotoxic activities observed for VacA-m1 appeared to be higher than those of VacA-m2. Homology-based modeling and sequence analysis suggested a potential structural impact of non-polar residues located at the N-terminus of the mature VacA toxin and its 33-kDa domain. Our data provide molecular insights into selective influence of the N-terminally added tag on efficient expression of recombinant VacA variants, signifying biochemical and biological implications of the hydrophobic stretch within the N-terminal domain.

Inhibition of jack bean urease by amphiphilic peptides

In the current study, amphiphilic peptides were designed and screened against Jack bean urease by using computer aided drug discovery approach. The result showed that out of thirty-eight amphiphilic peptides 1, 3, 12, 18, 30, and 33 exhibit stronger binding affinity with the active site of the enzyme through chelation of charged amino acids with the nickel ions i.e., Ni+2 841 and Ni+2 842 as well as hydrophobic contacts of the nonpolar tail with the nonpolar residues in the active site. The selected amphiphilic peptides were synthesized by solid-phase peptide synthesis strategy, characterized by fast atomic bombardment mass spectroscopy (FAB-MS) and nuclear magnetic resonance spectroscopy (1H and 13C-NMR) and in vitro urease inhibitory activity of amphiphilic peptides was studied. Amphiphilic peptides 12 and 33 showed excellent urease inhibitory activity, (p < 0.001) with IC50 values 20.5 ± 0.01, and 28.1 ± 0.03 µM respectively, which was considerably better than thiourea used as positive control.

Interaction fingerprint of transmembrane segments in voltage sensor domains

Voltage sensor domain (VSD) in channel and non-channel membrane proteins shares a common function in the detection of changes in the transmembrane electric potential. The VSD is made of four helical transmembrane segments (S1–S4) that form a structurally conserved scaffold through inter-transmembrane residue-residue interactions. Details about these interactions are yet to be fully understood in the context of the unique structural and physical characteristics of the voltage sensor unit. In this study, molecular dynamics simulations were carried out to investigate transmembrane helix-helix interactions via residue-based nonbonding energies using the activated and resting state conformations of VSD from Hv1, CiVSP, KvAP and NavAb. Inter-transmembrane interaction energies within the VSD were determined. Analysis of electrostatic and van der Waals components revealed the strengths and weaknesses of the interactions between each pair of transmembrane segments. In all cases the S4 helix had the highest electrostatic contribution to favor the key role as the voltage sensitive segment. Electrostatic interactions for the S1–S2 pair as well as the S1–S3 pair were relatively weak. Van der Waal interaction energies between adjacent segments were on average greater than that between diagonally opposite segments. Salt bridge interactions between S4-arginines and the negatively charged residues in other segments appear to contribute more to stabilizing the energy than the van der Waals interactions between nonpolar residues. The overall behavior of residue-residue contacts is similar among the transmembrane domains, reflecting the common inter- transmembrane interaction pattern in the VSD. In addition, analysis of the residue positions suggested that subtle differences in the orientation of the salt-bridges can be attributed to the difference in the inter-transmembrane interaction strengths inside the VSDs.

Fluorescence Anisotropy Decays and Microscale-Volume Viscometry Reveal the Compaction of Ribosome-Bound Nascent Proteins.

This work introduces a technology that combines fluorescence anisotropy decay with microscale-volume viscometry to investigate the compaction and dynamics of ribosome-bound nascent proteins. Protein folding in the cell, especially when nascent chains emerge from the ribosomal tunnel, is poorly understood. Previous investigations based on fluorescence anisotropy decay determined that a portion of the ribosome-bound nascent protein apomyoglobin (apoMb) forms a compact structure. This work, however, could not assess the size of the compact region. The combination of fluorescence anisotropy with microscale-volume viscometry, presented here, enables identifying the size of compact nascent-chain subdomains using a single fluorophore label. Our results demonstrate that the compact region of nascent apoMb contains 57-83 amino acids and lacks residues corresponding to the two native C-terminal helices. These amino acids are necessary for fully burying the nonpolar residues in the native structure, yet they are not available for folding before ribosome release. Therefore, apoMb requires a significant degree of post-translational folding for the generation of its native structure. In summary, the combination of fluorescence anisotropy decay and microscale-volume viscometry is a powerful approach to determine the size of independently tumbling compact regions of biomolecules. This technology is of general applicability to compact macromolecules linked to larger frameworks.

Native to designed: microbial -amylases for industrial applications.

Background-amylases catalyze the endo-hydrolysis of -1,4-D-glycosidic bonds in starch into smaller moieties. While industrial processes are usually performed at harsh conditions, -amylases from mainly the bacteria, fungi and yeasts are preferred for their stabilities (thermal, pH and oxidative) and specificities (substrate and product). Microbial -amylases can be purified and characterized for industrial applications. While exploring novel enzymes with these properties in the nature is time-costly, the advancements in protein engineering techniques including rational design, directed evolution and others have privileged their modifications to exhibit industrially ideal traits. However, the commentary on the strategies and preferably mutated residues are lacking, hindering the design of new mutants especially for enhanced substrate specificity and oxidative stability. Thus, our review ensures wider accessibility of the previously reported experimental findings to facilitate the future engineering work.Survey methodology and objectivesA traditional review approach was taken to focus on the engineering of microbial -amylases to enhance industrially favoured characteristics. The action mechanisms of - and -amylases were compared to avoid any bias in the research background. This review aimed to discuss the advances in modifying microbial -amylases via protein engineering to achieve longer half-life in high temperature, improved resistance (acidic, alkaline and oxidative) and enhanced specificities (substrate and product). Captivating results were discussed in depth, including the extended half-life at 100C, pH 3.5 and 10, 1.8 M hydrogen peroxide as well as enhanced substrate (65.3%) and product (42.4%) specificities. These shed light to the future microbial -amylase engineering in achieving paramount biochemical traits ameliorations to apt in the industries.ConclusionsMicrobial -amylases can be tailored for specific industrial applications through protein engineering (rational design and directed evolution). While the critical mutation points are dependent on respective enzymes, formation of disulfide bridge between cysteine residues after mutations is crucial for elevated thermostability. Amino acids conversion to basic residues was reported for enhanced acidic resistance while hydrophobic interaction resulted from mutated hydrophobic residues in carbohydrate-binding module or surface-binding sites is pivotal for improved substrate specificity. Substitution of oxidation-prone methionine residues with non-polar residues increases the enzyme oxidative stability. Hence, this review provides conceptual advances for the future microbial -amylases designs to exhibit industrially significant characteristics. However, more attention is needed to enhance substrate specificity and oxidative stability since they are least reported.

Identification of Pyruvate Dehydrogenase E1 as a Potential Target against Magnaporthe oryzae through Experimental and Theoretical Investigation.

Magnaporthe oryzae (M. oryzae) is a typical cause of rice blast in agricultural production. Isobavachalcone (IBC), an active ingredient of Psoralea corylifolia L. extract, is an effective fungicide against rice blast. To determine the mechanism of IBC against M. oryzae, the effect of IBC on the metabolic pathway of M. oryzae was explored by transcriptome profiling. In M. oryzae, the expression of pyruvate dehydrogenase E1 (PDHE1), part of the tricarboxylic acid (TCA cycle), was significantly decreased in response to treatment with IBC, which was verified by qPCR and testing of enzyme activity. To further elucidate the interactions between IBC and PDHE1, the 3D structure model of the PDHE1 from M. oryzae was established based on homology modeling. The model was utilized to analyze the molecular interactions through molecular docking and molecular dynamics simulation, revealing that IBC has π-π stacking interactions with residue TYR139 and undergoes hydrogen bonding with residue ASP217 of PDHE1. Additionally, the nonpolar residues PHE111, MET174, ILE 187, VAL188, and MET250 form strong hydrophobic interactions with IBC. The above results reveal that PDHE1 is a potential target for antifungal agents, which will be of great significance for guiding the design of new fungicides. This research clarified the mechanism of IBC against M. oryzae at the molecular level, which will underpin further studies of the inhibitory mechanism of flavonoids and the discovery of new targets. It also provides theoretical guidance for the field application of IBC.

Computational optimization of the SARS-CoV-2 receptor-binding-motif affinity for human ACE2

The coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is responsible for the coronavirus disease 2019 pandemic, and the closely related SARS-CoV coronavirus enter cells by binding at the human angiotensin converting enzyme 2 (hACE2). The stronger hACE2 affinity of SARS-CoV-2 has been connected with its higher infectivity. In this work, we study hACE2 complexes with the receptor-binding domains (RBDs) of the human SARS-CoV-2 and human SARS-CoV viruses, using all-atom molecular dynamics simulations and computational protein design with a physics-based energy function. The molecular dynamics simulations identify charge-modifying substitutions between the CoV-2 and CoV RBDs, which either increase or decrease the hACE2 affinity of the SARS-CoV-2 RBD. The combined effect of these mutations is small, and the relative affinity is mainly determined by substitutions at residues in contact with hACE2. Many of these findings are in line and interpret recent experiments. Our computational protein design calculations redesign positions 455, 493, 494, and 501 of the SARS-CoV-2 receptor binding motif, which contact hACE2 in the complex and are important for ACE2 recognition. Sampling is enhanced by an adaptive importance sampling Monte Carlo method. Sequences with increased affinity replace CoV-2 glutamine by a negative residue at position 493; serine by a nonpolar or aromatic residue or an asparagine at position 494; and asparagine by valine or threonine at position 501. Substitutions at positions 455 and 501 have a smaller effect on affinity. Substitutions suggested by our design are seen in viral sequences encountered in other species, including bat and pangolin. Our results might be used to identify potential virus strains with higher human infectivity and assist in the design of peptide-based or peptidomimetic compounds with the potential to inhibit SARS-CoV-2 binding at hACE2.

Binding Mechanisms of Amyloid-like Peptides to Lipid Bilayers and Effects of Divalent Cations.

In several neurodegenerative diseases, cell toxicity can emerge from damage produced by amyloid aggregates to lipid membranes. The details accounting for this damage are poorly understood including how individual amyloid peptides interact with phospholipid membranes before aggregation. Here, we use all-atom molecular dynamics simulations to investigate the molecular mechanisms accounting for amyloid-membrane interactions and the role played by calcium ions in this interaction. Model peptides known to self-assemble into amyloid fibrils and bilayer made from zwitterionic and anionic lipids are used in this study. We find that both electrostatic and hydrophobic interactions contribute to peptide-bilayer binding. In particular, the attraction of peptides to lipid bilayers is dominated by electrostatic interactions between positive residues and negative phosphate moieties of lipid head groups. This attraction is stronger for anionic bilayers than for zwitterionic ones. Hydrophobicity drives the burial of nonpolar residues into the interior of the bilayer producing strong binding in our simulations. Moreover, we observe that the attraction of peptides to the bilayer is significantly reduced in the presence of calcium ions. This is due to the binding of calcium ions to negative phosphate moieties of lipid head groups, which leaves phospholipid bilayers with a net positive charge. Strong binding of the peptide to the membrane occurs less frequently in the presence of calcium ions and involves the formation of a "Ca2+ bridge".

Analysis of Binding Modes of Antigen-Antibody Complexes by Molecular Mechanics Calculation.

Antibodies are one of the most important protein molecules in biopharmaceutics. Due to the recent advance in technology for producing monoclonal antibodies, many structural data are available on the antigen-antibody complexes. To characterize the molecular interaction in antigen-antibody recognition, we computationally analyzed 500 complex structures by molecular mechanics calculations. The presence of Ser and Tyr is markedly large in the complementarity-determining regions (CDRs). Although Ser is abundant in CDRs, its contribution to the binding score is not large. Instead, Tyr, Asp, Glu, and Arg significantly contribute to the molecular interaction from the viewpoint of the binding score. The decomposition of the binding score suggests that the hydrophilic interaction is predominant in all CDRs compared with the hydrophobic one. The contribution of the heavy chain is larger than that of the light chain. In particular, H2 and H3 largely contribute to the binding interaction. Tyr is a main contributing residue both in H2 and H3. The positively charged residue Arg also significantly contributes to the binding score in H3, while the contribution of Lys is small. The appearance of Ser is remarkable in H2, and Asp is abundant in H3. The non-charged polar residues, Thr, Asn, and Gln, appear much in H2, compared to appearing in H3. The negatively charged residues Asp and Glu significantly contribute to the binding score in H3. The contributions of Phe and Trp are not large in spite that the aromatic residues are capable of making the π-π or CH-π interaction. Gly is commonly abundant both in H2 and H3. The average distance of the shortest direct hydrogen bond between the antigen and antibody is longer than that of the hydrogen bonds observed in the complexes between compounds and their target proteins. Therefore, the antigen-antibody interface is not so tight as the compound-target protein interface. The calculation of shape complementarity is consistent with the result of the hydrogen bonds in that the fitness of the antigen-antibody contact is not so high as that of the compound-target protein contact. There exist many water molecules at the antigen-antibody interface. These findings suggest that Tyr, Asp, Glu, and Arg are rich in H3 and work as major contributors for the interaction with the antigen. Ser, Thr, Asn, and Gln are rich in H2 and support the interaction with enhancing molecular fitness. Gly is helpful in increasing flexibility and geometrical diversity. Because the antigen-antibody binding is fundamentally hydrophilic-driven, the non-polar residues are unfavorable for mediating the contact even for the aromatic residues such as Phe and Trp.

Synthesis of SARS‐CoV‐2 Transmembrane Domains Using Solid Phase Peptide Chemistry and Analysis of Their Oligomerization for Novel Drug Design

The SARS‐CoV‐2 spike protein transmembrane (TM) domain is an important component of the protein that contributes to its oligomerization for its active configuration. An understanding of this function can yield valuable data to future research on how to inhibit the virus’ infection of a host. The TM region of the protein is not as well studied due to the complication that nonpolar peptides can have in large groups; therefore, an understanding of this region could lead to novel discoveries. A peptide of the TM domain must be synthesized by using Solid Phase Peptide Synthesis, in which a resin is utilized to build a peptide off the c‐terminus using Fmoc‐protected amino acids and a strong base. A fluorescence assay can be performed to determine its oligomerization state by looking at the change of fluorescence that occurs when a quencher is added to the solution. Once the oligomerization state of the spike proteins TM domain is determined, this information will be beneficial the development of drugs against the virus. The set up to run this synthesis was developed in the lab and is causing premature termination of the amino acid chain. Current work involves addressing a potential issue of hydrophobicity due to the volume non‐polar residues, which is currently being address by trying a TGR resin that has a polystyrene core which is suited for long peptides. It is anticipated that addressing this issue will lead to a peptide that can be purified. The future implications of this research is that it has the potential of being useful for drugs that would inhibit the active configuration of SARS‐CoV‐2.

The effect of acidic pH on the adsorption and lytic activity of the peptides Polybia-MP1 and its histidine-containing analog in anionic lipid membrane: a biophysical study by molecular dynamics and spectroscopy.

Antimicrobial peptides (AMPs) are part of the innate immune system of many species. AMPs are short sequences rich in charged and non-polar residues. They act on the lipid phase of the plasma membrane without requiring membrane receptors. Polybia-MP1 (MP1), extracted from a native wasp, is a broad-spectrum bactericide, an inhibitor of cancer cell proliferation being non-hemolytic and non-cytotoxic. MP1 mechanism of action and its adsorption mode is not yet completely known. Its adsorption to lipid bilayer and lytic activity is most likely dependent on the ionization state of its two acidic and three basic residues and consequently on the bulk pH. Here we investigated the effect of bulk acidic (pH 5.5) and neutral pH (7.4) solution on the adsorption, insertion, and lytic activity of MP1 and its analog H-MP1 to anionic (7POPC:3POPG) model membrane. H-MP1 is a synthetic analog of MP1 with lysines replaced by histidines. Bulk pH changes could modulate this peptide efficiency. The combination of different experimental techniques and molecular dynamics (MD) simulations showed that the adsorption, insertion, and lytic activity of H-MP1 are highly sensitive to bulk pH in opposition to MP1. The atomistic details, provided by MD simulations, showed peptides contact their N-termini to the bilayer before the insertion and then lay parallel to the bilayer. Their hydrophobic faces inserted into the acyl chain phase disturb the lipid-packing.

Current Understanding of the Structure, Stability and Dynamic Properties of Amyloid Fibrils.

Amyloid fibrils are supramolecular protein assemblies represented by a cross-β structure and fibrous morphology, whose structural architecture has been previously investigated. While amyloid fibrils are basically a main-chain-dominated structure consisting of a backbone of hydrogen bonds, side-chain interactions also play an important role in determining their detailed structures and physicochemical properties. In amyloid fibrils comprising short peptide segments, a steric zipper where a pair of β-sheets with side chains interdigitate tightly is found as a fundamental motif. In amyloid fibrils comprising longer polypeptides, each polypeptide chain folds into a planar structure composed of several β-strands linked by turns or loops, and the steric zippers are formed locally to stabilize the structure. Multiple segments capable of forming steric zippers are contained within a single protein molecule in many cases, and polymorphism appears as a result of the diverse regions and counterparts of the steric zippers. Furthermore, the β-solenoid structure, where the polypeptide chain folds in a solenoid shape with side chains packed inside, is recognized as another important amyloid motif. While side-chain interactions are primarily achieved by non-polar residues in disease-related amyloid fibrils, the participation of hydrophilic and charged residues is prominent in functional amyloids, which often leads to spatiotemporally controlled fibrillation, high reversibility, and the formation of labile amyloids with kinked backbone topology. Achieving precise control of the side-chain interactions within amyloid structures will open up a new horizon for designing useful amyloid-based nanomaterials.

Functional Selectivity of a Biased Cannabinoid-1 Receptor (CB1R) Antagonist.

Seven-transmembrane receptors signal via G-protein- and β-arrestin-dependent pathways. We describe a peripheral CB1R antagonist (MRI-1891) highly biased toward inhibiting CB1R-induced β-arrestin-2 (βArr2) recruitment over G-protein activation. In obese wild-type and βArr2-knockout (KO) mice, MRI-1891 treatment reduces food intake and body weight without eliciting anxiety even at a high dose causing partial brain CB1R occupancy. By contrast, the unbiased global CB1R antagonist rimonabant elicits anxiety in both strains, indicating no βArr2 involvement. Interestingly, obesity-induced muscle insulin resistance is improved by MRI-1891 in wild-type but not in βArr2-KO mice. In C2C12 myoblasts, CB1R activation suppresses insulin-induced akt-2 phosphorylation, preventable by MRI-1891, βArr2 knockdown or overexpression of CB1R-interacting protein. MRI-1891, but not rimonabant, interacts with nonpolar residues on the N-terminal loop, including F108, and on transmembrane helix-1, including S123, a combination that facilitates βArr2 bias. Thus, CB1R promotes muscle insulin resistance via βArr2 signaling, selectively mitigated by a biased CB1R antagonist at reduced risk of central nervous system (CNS) side effects.

C-terminal troponin-I residues trap tropomyosin in the muscle thin filament blocked-state

Tropomyosin and troponin regulate muscle contraction by participating in a macromolecular scale steric-mechanism to control myosin-crossbridge – actin interactions and consequently contraction. At low-Ca2+, the C-terminal 30% of troponin subunit-I (TnI) is proposed to trap tropomyosin in a position on thin filaments that sterically interferes with myosin-binding, thus causing muscle relaxation. In contrast, at high-Ca2+, inhibition is released after the C-terminal domains dissociate from F-actin-tropomyosin as its component switch-peptide domain binds to the N-lobe of troponin-C (TnC). Recent, paradigm-shifting, cryo-EM reconstructions by the Namba group have revealed density attributed to TnI along cardiac muscle thin filaments at both low- and high-Ca2+ concentration. Modeling the reconstructions showed expected high-Ca2+ hydrophobic interactions of the TnI switch-peptide and TnC. However, under low-Ca2+ conditions, sparse interactions of TnI and tropomyosin, and in particular juxtaposition of non-polar switch-peptide residues and charged tropomyosin amino acids in the published model seem difficult to reconcile with an expected steric-blocking conformation. This anomaly is likely due to inaccurate fitting of tropomyosin into the cryo-EM volume. In the current study, the low-Ca2+ cryo-EM volume was fitted with a more accurate tropomyosin model and representation of cardiac TnI. Our results show that at low-Ca2+ a cluster of hydrophobic residues at the TnI switch-peptide and adjacent H4 helix (Ala149, Ala151, Met 154, Leu159, Gly160, Ala161, Ala163, Leu167, Leu169, Ala171, Leu173) draw-in tropomyosin surface residues (Ile143, Ile146, Ala151, Ile154), presumably attracting the entire tropomyosin cable to its myosin-blocking position on actin. The modeling confirms that neighboring TnI “inhibitory domain” residues (Arg145, Arg148) bind to thin filaments at actin residue Asp25, as previously suggested. ClusPro docking of TnI residues 137–184 to actin-tropomyosin, including the TnI inhibitory-domain, switch-peptide and Helix H4, verified the modeled configuration. Our residue-to-residue contact-mapping of the TnI-tropomyosin association lends itself to experimental validation and functional localization of disease-bearing mutations.

His180 in the pore-lining α4 of the Bacillus thuringiensis Cry4Aa δ-endotoxin is crucial for structural arrangements of the α4-α5 transmembrane hairpin and hence biotoxicity

One proposed toxic mechanism of Bacillus thuringiensis Cry δ-endotoxins involves pore formation in target membranes by the α4-α5 transmembrane hairpin constituting their pore-forming domain. Here, nine selected charged and uncharged polar residues in the pore-lining α4 of the Cry4Aa mosquito-active toxin were substituted with Ala. All mutant toxins, i.e., D169A, R171A, Q173A, H178A, Y179A, H180A, Q182A, N183A and E187A, were over-expressed in Escherichia coli as 130-kDa protoxin inclusions at levels comparable to the wild-type toxin. Bioassays against Aedes aegypti larvae revealed that only H178A and H180A mutants displayed a drastic reduction in biotoxicity, albeit almost complete insolubility observed for H178A, but not for H180A inclusions. Further mutagenic analysis showed that replacements of His180 with charged (Arg, Lys, Asp, Glu), small uncharged polar (Ser, Cys) or small non-polar (Gly, Val) residues severely impaired the biotoxicity, unlike substitutions with relatively large uncharged (Asn, Gln, Leu) or aromatic (Phe, Tyr, Trp) residues. Similar to the trypsin-activated wild-type toxin, both bio-active and -inactive H180 mutants were still capable of releasing entrapped calcein from lipid vesicles and producing cation-selective channels with ~130-pS maximum conductance. Analysis of the Cry4Aa structure revealed the existence of a hydrophobic cavity near the critical His180 side-chain. Analysis of simulated structures revealed that His180-to-smaller residue conversions create a gap disrupting such cavity's hydrophobicity and hence structural arrangements of the α4-α5 hairpin. Altogether, our data disclose a critical involvement in Cry4Aa-biotoxicity of His180 exclusively present in the lumen-facing α4 for providing proper environment for the α4-α5 hairpin prior to membrane-inserted pore formation.

Profiling thimet oligopeptidase-mediated proteolysis in Arabidopsis thaliana.

Protein homeostasis (proteostasis) is crucial for proper cellular function, including the production of peptides with biological functions through controlled proteolysis. Proteostasis has roles in maintenance of cellular functions and plant interactions with the environment under physiological conditions. Plant stress continues to reduce agricultural yields causing substantial economic losses; thus, it is critical to understand how plants perceive stress signals to elicit responses for survival. As previously shown in Arabidopsis thaliana, thimet oligopeptidases (TOPs) TOP1 (also referred to as organellar oligopeptidase) and TOP2 (also referred to as cytosolic oligopeptidase) are essential components in plant response to pathogens, but further characterization of TOPs and their peptide substrates is required to understand their contributions to stress perception and defense signaling. Herein, label-free peptidomics via liquid chromatography-tandem mass spectrometry was used to differentially quantify 1111 peptides, originating from 369 proteins, between the Arabidopsis Col-0 wild type and top1top2 knock-out mutant. This revealed 350 peptides as significantly more abundant in the mutant, representing accumulation of these potential TOP substrates. Ten direct substrates were validated using in vitro enzyme assays with recombinant TOPs and synthetic candidate peptides. These TOP substrates are derived from proteins involved in photosynthesis, glycolysis, protein folding, biogenesis, and antioxidant defense, implicating TOP involvement in processes aside from defense signaling. Sequence motif analysis revealed TOP cleavage preference for non-polar residues in the positions surrounding the cleavage site. Identification of these substrates provides a framework for TOP signaling networks, through which the interplay between proteolytic pathways and defense signaling can be further characterized.

Temperature Dependency of CLC-ec1 Dimerization in Lipid Bilayers

Why do greasy membrane proteins form stable complexes via greasy binding interfaces in the similarly greasy lipid solvent? To answer this question, we must quantify the thermodynamic factors that drive membrane protein association in membranes. Previously, we found that the CLC-ec1 Cl-/H+ antiporter participates in an equilibrium dimerization reaction in 2:1 POPE:POPG lipid bilayers. The dimer assembles via a large membrane embedded interface lined by non-polar residues, and the complex is stable, exhibiting a free energy of dimerization of −11 kcal/mole relative to the subunit/lipid standard state. Thus, CLC-ec1 provides an ideal model system for studying this question. In the current study, we carried out a van't Hoff analysis of the CLC-ec1 (WT) dimerization equilibrium in lipid bilayers. To verify that we are studying the equilibrated system, we examined subunit-exchange kinetics between WT-Cy3, and WT-Cy5 EPL (E. coli polar lipid) proteoliposomes, fused to form multi-lamellar membranes. Mixing of subunits led to formation of heterodimers (WT-Cy3/WT-Cy5) and the resulting increase in bulk Förster Resonance Energy Transfer (FRET) was monitored as a function of time, and temperature (22-56°C). Additionally, temperature jump type experiments conducted for the co-labelled WT showed changes in FRET, with both fused and co-labelled samples converging to the same FRET- plateaus. We found that rate of subunit-exchange increases with temperature, and that FRET- plateau decreases above 44oC indicating a decrease in oligomerization, with the chloride transport function remaining unchanged across the same temperature range. To measure the full dimerization binding isotherm, we equilibrated WT-Cy5 in EPL membranes at various temperatures and carried out subunit capture single-molecule photobleaching analysis for measurement of Keq(T). Finally, a van't Hoff plot of ln(Keq) vs. 1/T is presented, allowing a thermodynamic dissection of enthalpy, entropy, and heat capacity changes associated with CLC dimerization in membranes.