Helical Regions Research Articles

Hydrophobic residues in disordered and helical regions mediate the oligomerization and phase separation of TDP-43 C-terminal domain.

TAR DNA binding protein 43 (TDP-43) is a nuclear, multidomain protein implicated in RNA metabolism. The cytoplasmic accumulation of TDP-43 can result in the formation of inclusion bodies which are a hallmark of neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The prion-like, disordered C-terminal domain (CTD) of TDP-43 (spanning residues 267-414) is aggregation-prone and harbors the majority (∼90%) of all ALS-related mutations. Several studies report that CTD can undergo liquid-liquid phase separation (LLPS). The conserved, hydrophobic region (CR, spanning residues 319-341) in CTD populates an α-helical structure which is essential for phase separation. The requirement of CR for phase separation is attributed to its tendency to undergo helix-helix self-association and promote oligomerization. In addition to CR, aromatic residues in the flanking disordered regions (IDR1/2) also play a critical role in CTD self-assembly. However, the relative contributions of critical residue types in CR and flanking regions on the overall phase separation propensity of CTD are unclear. Here, we utilized multiscale simulations coupled with biophysical experiments to determine the role of hydrophobic residues in the CR and IDR1/2 regions on CTD phase separation. Single chain and condensed phase simulations collectively identified the contribution of aromatic residues while also suggesting significant contributions to LLPS by other residue types, which were previously not recognized as essential for phase separation within and outside the CR region. The computational observations were successfully corroborated by in vitro microscopy and saturation concentration measurements. Finally, NMR experiments determined that CTD oligomerization was critically dependent on these residues within CR and IDR1/2 regions. Overall, our work uncovers the collective importance of hydrophobic residues across both helical and disordered regions in modulating the oligomerization and phase separation of CTD.

Revisiting the Ramachandran plot based on statistical analysis of static and dynamic characteristics of protein structures

Ramachandran plots, which describe protein structures by plotting the dihedral angle pairs of the backbone on a two-dimensional plane, have played an important role in structural biology over the past few decades. However, despite continued discovery of new protein structures to date, the Ramachandran plot is still constructed by only a small number of data points, and further it cannot reflect the steric information of proteins. Here, we investigated the secondary structure of proteins in terms of static and dynamic characteristics. As for static feature, the Ramachandran plot was revisited for the dataset consisting of 9,148 non-redundant high-resolution protein structures released in the protein data bank until April 1, 2022. By calculating amino acid propensities, it was found that the proportion of secondary structures with respect to residue depth is directly related to their hydrophobicity. As for dynamic feature, normal mode analysis (NMA) based on an elastic network model (ENM) was carried out for the dataset using our KOSMOS web server (http://bioengineering.skku.ac.kr/kosmos/). All ENM-based NMA results were stored in the KOSMOS database, allowing researchers to use them in various ways. In this process, it was commonly found that high B-factors appeared at the edge of the alpha helix region, which was elucidated by introducing residue depth. In addition, by investigating the change in dihedral angle, it was possible to quantitatively survey the contribution of structural change of protein on the Ramachandran plot. In conclusion, our statistical analysis of protein characteristics will provide insight into a range of protein structural studies.

Coenzyme Engineering of Glucose-6-phosphate Dehydrogenase on a Nicotinamide-Based Biomimic and Its Application as a Glucose Biosensor

Nicotinamide adenine dinucleotide (NAD(P)+)-dependent oxidoreductases have been widely employed as biocatalysts for numerous applications, such as in vitro biomanufacturing and biosensors. The application of biomimetic nicotinamide coenzymes (BNCs) in an enzymatic redox cascade constitutes a promising alternative that can eliminate the need for expensive natural NAD(P)+ coenzymes. Herein, we demonstrated that the coenzyme engineering of glucose-6-phosphate dehydrogenase from Zymomonas mobilis (ZmG6PDH) enhanced its catalytic efficiency (kcat/Km) on oxidized nicotinamide mononucleotide (NMN+). Compared with the wild-type (WT) enzyme, the optimal mutant R4 exhibited a 112-fold enhancement in catalytic efficiency on NMN+, with 4.7 × 103 and 2.6 × 103-fold change in the two coenzyme specificity values ([kcat/Km]NMN+/[kcat/Km]NAD+ and [kcat/Km]NMN+/[kcat/Km]NADP+), respectively. To obtain insights into the structure–function relationship of the enzyme and its mutant, we determined the crystal structures of WT G6PDH (apo-WT) and the complex of mutant R4 and NMN+ (R4:NMN+) and performed molecular dynamics simulations of the ternary complexes of the G6PDH/substrate (glucose-6-phosphate, G6P)/coenzyme. The results revealed that the helix region of the coenzyme-binding Rossmann-like domain in mutant R4 was nearer to the active site than the wild-type enzyme, thus creating a more compact substrate/coenzyme-binding site to favor the binding of the substrate G6P and NMN+. Furthermore, a glucose biosensor based on the glucokinase (GK)/G6PDH-NMN+ biosystem without cryopreservation was constructed using NMN+ as the coenzyme. These results indicate that the combination of BNCs and oxidoreductases can be widely employed in biosensors, biocatalysis, and in vitro biomanufacturing.

Mapping glycation and glycoxidation sites in collagen I of human cortical bone

Accumulation of advanced glycation end products (AGEs), particularly in long-lived extracellular matrix proteins, has been implicated in pathogenesis of diabetic complications and in aging. Knowledge about specific locations of AGEs and their precursors within protein primary structure is critical for understanding their physiological and pathophysiological impact. However, the information on specific AGE sites is lacking. Here, we identified sequence positions of four major AGEs, carboxymethyllysine, carboxyethyllysine, 5-hydro-5-methyl imidazolone, and 5-hydro-imidazolone, and an AGE precursor fructosyllysine within the triple helical region of collagen I from cortical bone of human femurs. The presented map provides a basis for site-specific quantitation of AGEs and other non-enzymatic post-translational modifications and identification of those sites affected by aging, diabetes, and other diseases such as osteoporosis; it can also help in guiding future studies of AGE impact on structure and function of collagen I in bone.

On a mechanistic impact of transmembrane tetramerization in the pathological activation of RTKs

Constitutive activation of receptor tyrosine kinases (RTKs) via different mutations has a strong impact on the development of severe human disorders, including cancer. Here we propose a putative activation scenario of RTKs, whereby transmembrane (TM) mutations can also promote higher-order oligomerization of the receptors that leads to the subsequent ligand-free activation. We illustrate this scenario using a computational modelling framework comprising sequence-based structure prediction and all-atom 1 µs molecular dynamics (MD) simulations in a lipid membrane for a previously characterised oncogenic TM mutation V536E in platelet-derived growth factor receptor alpha (PDGFRA). We show that in the course of MD simulations the mutant TM tetramer retains stable and compact configuration strengthened by tight protein-protein interactions, while the wild type TM tetramer demonstrates looser packing and a tendency to dissociate. Moreover, the mutation affects the characteristic motions of mutated TM helical segments by introducing additional non-covalent crosslinks in the middle of the TM tetramer, which operate as mechanical hinges. This leads to dynamic decoupling of the C-termini from the rigidified N-terminal parts and facilitates more pronounced possible displacement between the C-termini of the mutant TM helical regions that can provide more freedom for mutual rearrangement of the kinase domains located downstream. Our results for the V536E mutation in the context of PDGFRA TM tetramer allow for the possibility that the effect of oncogenic TM mutations can go beyond alternating the structure and dynamics of TM dimeric states and might also promote the formation of higher-order oligomers directly contributing to ligand-independent signalling effectuated by PDGFRA and other RTKs.

Structural studies of bacteriorhodopsin in BC era.

It marked half a century since the discovery of bacteriorhodopsin two years ago. On this occasion, I have revisited historically important diffraction studies of this membrane protein, based on my recollections. X-ray diffraction and electron diffraction, and electron microscopy, described the low-resolution structure of bacteriorhodopsin within the purple membrane. Neutron diffraction was effective to assign the helical regions in the primary structure with 7 rods revealed by low-resolution structure as well as to describe the retinal position. Substantial conformational changes upon light illumination were clarified by the structures of various photointermediates. Early trials of time-resolved studies were also introduced. Models for the mechanism of light-driven proton pump based on the low-resolution structural studies are also described. Significantly, they are not far from the today's understanding. I believe that the spirit of the early research scientists in this field and the essence of their studies, which constitute the foundations of the field, still actively fertilizes current membrane protein research.

Deep learning methods of predicting RNA torsion angle

Modeling of RNA tertiary structure is one of the basic problems in molecular biophysics, and it is very important in understanding the biological function of RNA and designing new structures. RNA tertiary structure is mainly determined by seven torsions of main-chain and side-chain backbone, the accurate prediction of these torsion angles is the basis of modeling RNA tertiary structure. At present, there are only a few methods of using deep learning to predict RNA torsion angles, and the prediction accuracy needs further improving if it is used to model RNA tertiary structure. In this study, we also develop a deep learning method, 1dRNA, to predict RNA backbone torsions and pseudotorsion angles, including two different deep learning models, the convolution model (DRCNN) that considers the features of adjacent nucleotides and the Hyper-long-short-term memory model (DHLSTM) that considers the features of all the nucleotides. We then empirically show that DRCNN and DHLSTM outperform existing state-of-the-art methods under the same datasets, the prediction accuracy of DRCNN model is improved by 5% to 28% for <i>β</i>, <i>δ</i>, <i>ζ</i>, <i>χ</i>, <i>η</i>, and <i>θ</i> angle, and the prediction accuracy of DHLSTM model is improved by 6% to 15% for <i>β</i>, <i>δ</i>, <i>ζ</i>, <i>χ</i>, <i>η</i>, <i>θ</i> angle. The DRCNN model predicts better results than the DHLSTM model and the existing models in the <i>δ</i>, <i>ζ</i>, <i>χ</i>, <i>η</i>, <i>θ</i> angle, and the DHLSTM model predicts better results than the DRCNN model and the existing model in the <i>β</i> and <i>ε</i> angles, and the existing models predicted better results than the DRCNN model and DHLSTM model in the <i>α</i> and <i>γ</i> angles. The DRCNN model and the existing models predict a richer distribution of angles than the DHLSTM model. In terms of model stability, the DHLSTM model is much more stable than the DRCNN model and the existing models, with fewer outliers. The results also show that the <i>α</i> angle and <i>γ</i> angle are the most difficult to predict, the angles of the ring region is more difficult to predict than the angles of the helix region, the model is also not sensitive to the change of the target sequence length, and the deviation of the model prediction angle from the decoys can also be used to evaluate the RNA tertiary structures quality.

An integrated investigation of structural and pathway alteration caused by PIK3CA and TP53 mutations identified in cfDNA of metastatic breast cancer.

In peripheral blood, cell-free DNA (cfDNA) contains circulating tumor DNA (ctDNA), which indicates molecular abnormalities in metastatic breast tumor tissue. The sequencing of cfDNA of Metastatic Breast Cancer (MBC) patients allows assessment of therapy response and noninvasive treatment. In the proposed study, clinically significant alterations in PIK3CA and TP53 genes associated with MBC resulting in a missense substitution of His1047Arg and Arg282Trp from an next-generation sequencing-based multi-gene panel were reported in a cfDNA of a patient with MBC. To investigate the impact of the reported mutation, we used molecular docking, molecular dynamics simulation, network analysis, and pathway analysis. Molecular Docking analysis determined the distinct binding pattern revealing H1047R-ATP complex has a higher number of Hydrogen bonds (H-bonds) and binding affinity with a slight difference compared to the PIK3CA-ATP complex. Following, molecular dynamics simulation for 200 ns, of which H1047R-ATP complex resulted in the instability of PIK3CA. Similarly, for TP53 mutant R282W, the zinc-free state (apo) and zinc-bounded (holo) complexes were investigated for conformational change between apo and holo complexes, of which the holo complex mutant R282W was unstable. To validate the conformational change of PIK3CA and TP53, 80% mutation of H1047R in the kinase domain of p110α expressed ubiquitously in PIK3CA protein that alters PI3K pathway, while R282W mutation in DNA binding helix (H2) region of P53 protein inhibits the transcription factor in P53 pathway causing MBC. According to our findings, the extrinsic (hypoxia, oxidative stress, and acidosis); intrinsic factors (MYC amplification) in PIK3CA and TP53 mutations will provide potential insights for developing novel therapeutic methods for MBC therapy.

Csk αC Helix: A Computational Analysis of an Essential Region for Conformational Transitions.

Conformational changes are an essential feature for the function of some dynamic proteins. Understanding the mechanism of such motions may allow us to identify important properties, which may be directly related to the regulatory function of a protein. Also, this knowledge may be employed for a rational design of drugs that can shift the balance between active and inactive conformations, as well as affect the kinetics of the activation process. Here, the conformational changes in carboxyl-terminal Src kinase, the major catalytic repressor to the Src family of kinases, was investigated, and it was proposed as a functionally related hypothesis. A Cα Structure-Based Model (Cα-SBM) was applied to provide a description of the overall conformational landscape and further analysis complemented by detailed molecular dynamics simulations. As a first approach to Cα-SBM simulations, reversible transitions between active (closed) and inactive (open) forms were modeled as fluctuations between these two energetic basins. It was found that, in addition to the interdomain Carboxyl-terminal SRC Kinase (Csk) correlated motions, a conformational change in the αC helix is required for a complete conformational transition. The result reveals this as an important region of transition control and domain coordination. Restrictions in the αC helix region of the Csk protein were performed, and the analyses showed a direct correlation with the global conformational changes, with this location being propitious for future studies of ligands. Also, the Src Homology 3 (SH3) and SH3 plus Src Homology 2 (SH2) domains were excluded for a direct comparison with experimental results previously published. Simulations where the SH3 was deleted presented a reduction of the transitions during the simulations, while the SH3-SH2 deletion vanishes the Csk transitions, corroborating the experimental results mentioned and linking the conformational changes with the catalytic functionality of Csk. The study was complemented by the introduction of a known kinase inhibitor close to the Csk αC helix region where its consequences for the kinetic behavior and domain displacement of Csk were verified through detailed molecular dynamics. The findings describe the mechanisms involving the Csk αC helix for the transitions and also support the dynamic correlation between SH3 and SH2 domains against the Csk lobes and how local energetic restrictions or interactions in the Csk αC helix can play an important role for long-range motions. The results also allow speculation if the Csk activity is restricted to one specific conformation or a consequence of a state transition, this point being a target for future studies. However, the αC helix is revealed as a potential region for rational drug design.

Computational investigation of functional water molecules in GPCRs bound to G protein or arrestin.

G protein-coupled receptors (GPCRs) are membrane proteins constituting the largest family of drug targets. The activated GPCR binds either the heterotrimeric G proteins or arrestin through its activation cycle. Water molecules have been reported to play a role in GPCR activation. Nevertheless, reported studies are focused on the hydrophobic helical bundle region. How water molecules function in GPCR bound either G protein or arrestin is rarely studied. To address this issue, we carried out computational studies on water molecules in both GPCR/G protein complexes and GPCR/arrestin complexes. Using inhomogeneous fluid theory (IFT), we locate all possible hydration sites in GPCRs binding either to G protein or arrestin. We observe that the number of water molecules on the interaction surface between GPCRs and signal proteins are correlated with the insertion depths of the α5-helix from G-protein or "finger loop" from arrestin in GPCRs. In three out of the four simulation pairs, the interfaces of Rhodopsin, M2R and NTSR1 in the G protein-associated systems show more water-mediated hydrogen-bond networks when compared to these in arrestin-associated systems. This reflects that more functionally relevant water molecules may probably be attracted in G protein-associated structures than that in arrestin-associated structures. Moreover, we find the water-mediated interaction networks throughout the NPxxY region and the orthosteric pocket, which may be a key for GPCR activation. Reported studies show that non-biased agonist, which can trigger both GPCR-G protein and GPCR-arrestin activation signal, can result in pharmacologically toxicities. Our comprehensive studies of the hydration sites in GPCR/G protein complexes and GPCR/arrestin complexes may provide important insights in the design of G-protein biased agonists.

Morphology, molecular phylogeny, and species delimitation within microalgal genera Eubrownia, Spongiococcum, and Chlorococcum (Chlorophyceae, Chlorophyta)

• Our knowledge of their taxonomic and phylogenetic relationships is still deficient. • Some genus and species from Moewusinia clade differ specific feature combinations. • The ITS2 of Eubrownia and Chlorococcum genera is unique within the chlorophytes. Despite that the diversity of chlorophycean microalgae has been studied for a long time, our knowledge of their taxonomic and phylogenetic relationships is still deficient. This study represents one of the first efforts to examine the congruence of the morphological identification of microalgal species from genera Eubrownia, Spongiococcum , and Chlorococcum with DNA-based species delimitation methods. We found unique combinations of morphological, molecular, and ecological features, distinguishing some members of the Moewusinia clade at the genus and species levels. Most genetic groups recovered by GMYC, ABGD, and PTP analyses of 18S rRNA and ITS2 were consistent with their morphological identification. The ITS2 region of genera Eubrownia and Chlorococcum is among the longest within the chlorophytes theretofore studied and has a unique secondary structure with branched helix III and uncommonly long helix IV. Moreover, the molecular signatures were found in the basal stem region of helix I, distinguishing two genera. The CBC analysis is able to effectively discriminate between E. isobilateralis, E. aggregata , and all Chlorococcum species, except for C. infusionum and C. echinozygotum . Thus, we confirmed tentative hypotheses of species, obtained by DNA-based delimitation methods, with morphological data, molecular signatures, ITS2 secondary structure, and CBC criterion.

Frequency, underdiagnosis, and heterogeneity of epidermal growth factor receptor exon 20 insertion mutations using real-world genomic datasets.

Epidermal growth factor receptor (EGFR) exon 20 insertion mutations (ex20ins) account for ≤ 12% of all EGFR-mutant nonsmall cell lung cancers. We analysed real-world datasets to determine the frequency of ex20ins variants, and the ability of polymerase chain reaction (PCR) and next-generation sequencing (NGS) to identify them. Three real-world United States NGS databases were used: GENIE, FoundationInsights, and GuardantINFORM. Mutation profiles consistent with in-frame EGFR ex20ins were summarized. GENIE, FoundationInsights, and GuardantINFORM datasets identified 180, 627, and 627 patients with EGFR ex20ins respectively. The most frequent insertion region of exon 20 was the near loop (~ 70%), followed by the far loop (~ 30%) and the helical (~ 3-6%) regions. GENIE, FoundationInsights, and GuardantINFORM datasets identified 41, 102, and 96 unique variants respectively. An analysis of variants projected that ~ 50% of EGFR ex20ins identified by NGS would have been missed by PCR-based assays. Given the breadth of EGFR ex20ins identified in the real-world US datasets, the ability of PCR to identify these mutations is limited. NGS platforms are more appropriate to identify patients likely to benefit from EGFR ex20ins-targeted therapies.

Dynamics and Allosteric Information Pathways of Unphosphorylated c-Cbl.

Human c-Cbl is a RING-type ligase and plays a central role in the protein degradation cascade. To elucidate its conformational changes related to substrate binding, we performed molecular dynamics simulations of different variants/states of c-Cbl for a cumulative time of 68 μs. Our simulations demonstrate that before the substrate binds, the RING domain samples a broad set of conformational states at a biologically relevant salt concentration, including the closed, partially open, and fully open states, whereas substrate binding leads to a restricted conformational sampling. Phe378 and the C-terminal region play an essential role in stabilizing the partially open state. To visualize the allosteric signal transmission pathways from the substrate-binding site to the 40 Å apart RING domain and identify the critical residues for allostery, we have created a subgraph from the optimal and suboptimal paths. Redundant paths are seen in the SH2 domain where the substrate binds, while the major bottlenecks are found at the junction between the SH2 domain and the linker helix region as well as that between the SH2 domain and the 4H bundle. These bottlenecks separate the paths into two overall routes. The nodes/residues at the bottlenecks on the subgraph are considered allosteric hot spots. This subgraph approach provides a general tool for network visualization and determination of critical residues for allostery. The structurally and allosterically critical residues identified in our work are testable and would provide valuable insights into the emerging strategies for drug discovery, such as targeted protein degradation.

Conformational Change of Transcription Factors from Search to Specific Binding: A lac Repressor Case Study.

In a process known as facilitated diffusion, DNA-binding proteins find their target sites by combining three-dimensional diffusion and one-dimensional scanning of the DNA. Following the trade-off between speed and stability, agile exploration of DNA requires loose binding, whereas, at the DNA target site, the searching protein needs to establish tight interactions with the DNA. To enable both efficient search and stable binding, DNA-binding proteins and DNA often switch conformations upon recognition. Here, we study the one-dimensional diffusion and DNA binding of the dimeric lac repressor (LacI), which was reported to adopt two different conformations when binding different conformations of DNA. Using coarse-grained molecular dynamic simulations, we studied the diffusion and the sequence-specific binding of these conformations of LacI, as well as their truncated or monomeric variants, with two DNA conformations: straight and bent. The simulations were compared to experimental observables. This study supports that linear diffusion along DNA combines tight rotation-coupled groove tracking and rotation-decoupled hopping, where the protein briefly dissociates and reassociates just a few base pairs away. Tight groove tracking is crucial for target-site recognition, while hopping speeds up the overall search process. We investigated the diffusion of different LacI conformations on DNA and show how the flexibility of LacI's hinge regions ensures agility on DNA as well as faithful groove tracking. If the hinge regions instead form α-helices at the protein-DNA interface, tight groove tracking is not possible. On the contrary, the helical hinge region is essential for tight binding to bent, specific DNA, for the formation of the specific complex. Based on our study of different encounter complexes, we argue that the conformational change in LacI and DNA bending are somewhat coupled. Our findings underline the importance of two distinct protein conformations for facilitated diffusion and specific binding, respectively.

A BAK subdomain that binds mitochondrial lipids selectively and releases cytochrome C.

How BAK and BAX induce mitochondrial outer membrane (MOM) permeabilization (MOMP) during apoptosis is incompletely understood. Here we have used molecular dynamics simulations, surface plasmon resonance, and assays for membrane permeabilization in vitro and in vivo to assess the structure and function of selected BAK subdomains and their derivatives. Results of these studies demonstrate that BAK helical regions α5 and α6 bind the MOM lipid cardiolipin. While individual peptides corresponding to these helical regions lack the full biological activity of BAK, tandem peptides corresponding to α4-α5, α5-α6, or α6-α7/8 can localize exogenous proteins to mitochondria, permeabilize liposomes composed of MOM lipids, and cause MOMP in the absence of the remainder of the BAK protein. Importantly, the ability of these tandem helices to induce MOMP under cell-free conditions is diminished by mutations that disrupt the U-shaped helix-turn-helix structure of the tandem peptides or decrease their lipid binding. Likewise, BAK-induced apoptosis in intact cells is diminished by CLS1 gene interruption, which decreases mitochondrial cardiolipin content, or by BAK mutations that disrupt the U-shaped tandem peptide structure or diminish lipid binding. Collectively, these results suggest that BAK structural rearrangements during apoptosis might mobilize helices involved in specific protein-lipid interactions that are critical for MOMP.

Small molecules targeting the disordered transactivation domain of the androgen receptor induce the formation of collapsed helical states

Intrinsically disordered proteins, which do not adopt well-defined structures under physiological conditions, are implicated in many human diseases. Small molecules that target the disordered transactivation domain of the androgen receptor have entered human trials for the treatment of castration-resistant prostate cancer (CRPC), but no structural or mechanistic rationale exists to explain their inhibition mechanisms or relative potencies. Here, we utilize all-atom molecular dynamics computer simulations to elucidate atomically detailed binding mechanisms of the compounds EPI-002 and EPI-7170 to the androgen receptor. Our simulations reveal that both compounds bind at the interface of two transiently helical regions and induce the formation of partially folded collapsed helical states. We find that EPI-7170 binds androgen receptor more tightly than EPI-002 and we identify a network of intermolecular interactions that drives higher affinity binding. Our results suggest strategies for developing more potent androgen receptor inhibitors and general strategies for disordered protein drug design.

Structural and thermodynamic analyses of the β-to-α transformation in RfaH reveal principles of fold-switching proteins.

The two-domain protein RfaH, a paralog of the universally conserved NusG/Spt5 transcription factors, is regulated by autoinhibition coupled to the reversible conformational switch of its 60-residue C-terminal Kyrpides, Ouzounis, Woese (KOW) domain between an α-hairpin and a β-barrel. In contrast, NusG/Spt5-KOW domains only occur in the β-barrel state. To understand the principles underlying the drastic fold switch in RfaH, we elucidated the thermodynamic stability and the structural dynamics of two RfaH- and four NusG/Spt5-KOW domains by combining biophysical and structural biology methods. We find that the RfaH-KOW β-barrel is thermodynamically less stable than that of most NusG/Spt5-KOWs and we show that it is in equilibrium with a globally unfolded species, which, strikingly, contains two helical regions that prime the transition toward the α-hairpin. Our results suggest that transiently structured elements in the unfolded conformation might drive the global folding transition in metamorphic proteins in general.

Computational Design and Experimental Validation of ACE2-Derived Peptides as SARS-CoV-2 Receptor Binding Domain Inhibitors.

The COVID-19 pandemic has caused significant social and economic disruption across the globe. Cellular entry of SARS-CoV-2 into the human body is mediated via binding of the Receptor Binding Domain (RBD) on the viral Spike protein (SARS-CoV-2 RBD) to Angiotensin-Converting Enzyme 2 (ACE2) expressed on host cells. Molecules that can disrupt ACE2:RBD interactions are attractive therapeutic candidates to prevent virus entry into human cells. A computational strategy that combines our Peptide Binding Design (PepBD) algorithm with atomistic molecular dynamics simulations was used to design new inhibitory peptide candidates via sequence iteration starting with a 23-mer peptide, referred to as SBP1. SBP1 is derived from a region of the ACE2 Peptidase Domain α1 helix that binds to the SARS-CoV-2 RBD of the initial Wuhan-Hu-1 strain. Three peptides demonstrated a solution-phase RBD-binding dissociation constant in the micromolar range during tryptophan fluorescence quenching experiments, one peptide did not bind, and one was insoluble at micromolar concentrations. However, in competitive ELISA assays, none of these peptides could outcompete ACE2 binding to SARS-CoV-2-RBD up to concentrations of 50 μM, similar to the parent SBP1 peptide which also failed to outcompete ACE2:RBD binding. Molecular dynamics simulations suggest that P4 would have a good binding affinity for the RBD domain of Beta-B.1.351, Gamma-P.1, Kappa-B.1.617.1, Delta-B.1.617.2, and Omicron-B.1.1.529 variants, but not the Alpha variant. Consistent with this, P4 bound Kappa-B.1.617.1 and Delta-B.1.617.2 RBD with micromolar affinity in tryptophan fluorescence quenching experiments. Collectively, these data show that while relatively short unstructured peptides can bind to SARS-CoV-2 RBD with moderate affinity, they are incapable of outcompeting the strong interactions between RBD and ACE2.

Computational analysis of prodomain cysteines in human TGF-β proteins reveals frequent loss of disulfide-dependent regulation in tumors.

The functionally diverse members of the human Transforming Growth Factor-β (TGF-β) family are tightly regulated. TGF-β regulation includes 2 disulfide-dependent mechanisms-dimerization and partner protein binding. The specific cysteines participating in these regulatory mechanisms are known in just 3 of the 33 human TGF-β proteins. Human prodomain alignments revealed that 24 TGF-β prodomains contain conserved cysteines in 2 highly exposed locations. There are 3 in the region of the β8 helix that mediates dimerization near the prodomain carboxy terminus. There are 2 in the Association region that mediates partner protein binding near the prodomain amino terminus. The alignments predict the specific cysteines contributing to disulfide-dependent regulation of 72% of human TGF-β proteins. Database mining then identified 9 conserved prodomain cysteine mutations and their disease phenotypes in 7 TGF-β proteins. Three common adenoma phenotypes for prodomain cysteine mutations suggested 7 new regulatory heterodimer pairs. Two common adenoma phenotypes for prodomain and binding partner cysteine mutations revealed 17 new regulatory interactions. Overall, the analysis of human TGF-β prodomains suggests a significantly expanded scope of disulfide-dependent regulation by heterodimerization and partner protein binding; regulation that is often lost in tumors.

Broadly neutralizing antibodies to SARS-CoV-2 and other human coronaviruses.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a recently emerged pathogenic human coronavirus that belongs to the sarbecovirus lineage of the genus Betacoronavirus. The ancestor strain has evolved into a number of variants of concern, with the Omicron variant of concern now having many distinct sublineages. The ongoing COVID-19 pandemic caused by SARS-CoV-2 has caused serious damage to public health and the global economy, and one strategy to combat COVID-19 has been the development of broadly neutralizing antibodies for prophylactic and therapeutic use. Many are in preclinical and clinical development, and a few have been approved for emergency use. Here we summarize neutralizing antibodies that target four key regions within the SARS-CoV-2 spike (S) protein, namely the N-terminal domain and the receptor-binding domain in the S1 subunit, and the stem helix region and the fusion peptide region in the S2 subunit. Understanding the characteristics of these broadly neutralizing antibodies will accelerate the development of new antibody therapeutics and provide guidance for the rational design of next-generation vaccines.