Improved prediction of protein side-chain conformations with SCWRL4

CSAlign and CSAlign-Dock: Structure alignment of ligands considering full flexibility and application to protein-ligand docking

Structure prediction of protein-ligand complexes, called protein-ligand docking, is a critical computational technique that can be used to understand the underlying principle behind the protein functions at the atomic level and to design new molecules regulating the functions. Protein-ligand docking methods have been employed in structure-based drug discovery for hit discovery and lead optimization. One of the important technical challenges in protein-ligand docking is to account for protein conformational changes induced by ligand binding. A small change such as a single side-chain rotation upon ligand binding can hinder accurate docking. Here we report an increase in docking performance achieved by structure alignment to known complex structures. First, a fully flexible compound-to-compound alignment method CSAlign is developed by global optimization of a shape score. Next, the alignment method is combined with a docking algorithm to dock a new ligand to a target protein when a reference protein-ligand complex structure is available. This alignment-based docking method, called CSAlign-Dock, showed superior performance to ab initio docking methods in cross-docking benchmark tests. Both CSAlign and CSAlign-Dock are freely available as a web server at https://galaxy.seoklab.org/csalign.

In silico study predicts a key role of RNA-binding domains 3 and 4 in nucleolin-miRNA interactions

RNA binding proteins (RBPs) regulate many important cellular processes through their interactions with RNA molecules. RBPs are critical for posttranscriptional mechanisms keeping gene regulation in a fine equilibrium. Conversely, dysregulation of RBPs and RNA metabolism pathways is an established hallmark of tumorigenesis. Human nucleolin (NCL) is a multifunctional RBP that interacts with different types of RNA molecules, in part through its four RNA binding domains (RBDs). Particularly, NCL interacts directly with microRNAs (miRNAs) and is involved in their aberrant processing linked with many cancers, including breast cancer. Nonetheless, molecular details of the NCL-miRNA interaction remain obscure. In this study, we used an in silico approach to characterize how NCL targets miRNAs and whether this specificity is imposed by a definite RBD-interface. Here, we present structural models of NCL-RBDs and miRNAs, as well as predict scenarios of NCL-miRNA interactions generated using docking algorithms. Our study suggests a predominant role of NCL RBDs 3 and 4 (RBD3-4) in miRNA binding. We provide detailed analyses of specific motifs/residues at the NCL-substrate interface in both these RBDs and miRNAs. Finally, we propose that the evolutionary emergence of more than two RBDs in NCL in higher organisms coincides with its additional role/s in miRNA processing. Our study shows that RBD3-4 display sequence/structural determinants to specifically recognize miRNA precursor molecules. Moreover, the insights from this study can ultimately support the design of novel antineoplastic drugs aimed at regulating NCL-dependent biological pathways with a causal role in tumorigenesis.

Excitatory and inhibitory D-serine binding to the NMDA receptor

N-methyl-D-aspartate receptors (NMDARs) uniquely require binding of two different neurotransmitter agonists for synaptic transmission. D-serine and glycine bind to one subunit, GluN1, while glutamate binds to the other, GluN2. These agonists bind to the receptor's bi-lobed ligand-binding domains (LBDs), which close around the agonist during receptor activation. To better understand the unexplored mechanisms by which D-serine contributes to receptor activation, we performed multi-microsecond molecular dynamics simulations of the GluN1/GluN2A LBD dimer with free D-serine and glutamate agonists. Surprisingly, we observed D-serine binding to both GluN1 and GluN2A LBDs, suggesting that D-serine competes with glutamate for binding to GluN2A. This mechanism is confirmed by our electrophysiology experiments, which show that D-serine is indeed inhibitory at high concentrations. Although free energy calculations indicate that D-serine stabilizes the closed GluN2A LBD, its inhibitory behavior suggests that it either does not remain bound long enough or does not generate sufficient force for ion channel gating. We developed a workflow using pathway similarity analysis to identify groups of residues working together to promote binding. These conformation-dependent pathways were not significantly impacted by the presence of N-linked glycans, which act primarily by interacting with the LBD bottom lobe to stabilize the closed LBD.

Borna Disease Virus 1 Phosphoprotein Forms a Tetramer and Interacts with Host Factors Involved in DNA Double-Strand Break Repair and mRNA Processing

Determining the structural organisation of viral replication complexes and unravelling the impact of infection on cellular homeostasis represent important challenges in virology. This may prove particularly useful when confronted with viruses that pose a significant threat to human health, that appear unique within their family, or for which knowledge is scarce. Among Mononegavirales, bornaviruses (family Bornaviridae) stand out due to their compact genomes and their nuclear localisation for replication. The recent recognition of the zoonotic potential of several orthobornaviruses has sparked a surge of interest in improving our knowledge on this viral family. In this work, we provide a complete analysis of the structural organisation of Borna disease virus 1 (BoDV-1) phosphoprotein (P), an important cofactor for polymerase activity. Using X-ray diffusion and diffraction experiments, we revealed that BoDV-1 P adopts a long coiled-coil α-helical structure split into two parts by an original β-strand twist motif, which is highly conserved across the members of whole Orthobornavirus genus and may regulate viral replication. In parallel, we used BioID to determine the proximal interactome of P in living cells. We confirmed previously known interactors and identified novel proteins linked to several biological processes such as DNA repair or mRNA metabolism. Altogether, our study provides important structure/function cues, which may improve our understanding of BoDV-1 pathogenesis.

Proteins from Thermophilic Thermus thermophilus Often Do Not Fold Correctly in a Mesophilic Expression System Such as Escherichia coli

Majority of protein structure studies use Escherichia coli (E. coli) and other model organisms as expression systems for other species' genes. However, protein folding depends on cellular environment factors, such as chaperone proteins, cytoplasmic pH, temperature, and ionic concentrations. Because of differences in these factors, especially temperature and chaperones, native proteins in organisms such as extremophiles may fold improperly when they are expressed in mesophilic model organisms. Here we present a methodology of assessing the effects of using E. coli as the expression system on protein structures. We compare these effects between eight mesophilic bacteria and Thermus thermophilus (T. thermophilus), a thermophile, and found that differences are significantly larger for T. thermophilus. More specifically, helical secondary structures in T. thermophilus proteins are often replaced by coil structures in E. coli. Our results show unique directionality in misfolding when proteins in thermophiles are expressed in mesophiles. This indicates that extremophiles, such as thermophiles, require unique protein expression systems in protein folding studies.

Cross-linking of the endolysosomal system reveals potential flotillin structures and cargo

Lysosomes are well-established as the main cellular organelles for the degradation of macromolecules and emerging as regulatory centers of metabolism. They are of crucial importance for cellular homeostasis, which is exemplified by a plethora of disorders related to alterations in lysosomal function. In this context, protein complexes play a decisive role, regulating not only metabolic lysosomal processes but also lysosome biogenesis, transport, and interaction with other organelles. Using cross-linking mass spectrometry, we analyze lysosomes and early endosomes. Based on the identification of 5376 cross-links, we investigate protein-protein interactions and structures of lysosome- and endosome-related proteins. In particular, we present evidence for a tetrameric assembly of the lysosomal hydrolase PPT1 and a heterodimeric structure of FLOT1/FLOT2 at lysosomes and early endosomes. For FLOT1-/FLOT2-positive early endosomes, we identify >300 putative cargo proteins and confirm eleven substrates for flotillin-dependent endocytosis, including the latrophilin family of adhesion G protein-coupled receptors.

Expanding the eukaryotic genetic code with a biosynthesized 21st amino acid

Genetic code expansion technology allows for the use of noncanonical amino acids (ncAAs) to create semisynthetic organisms for both biochemical and biomedical applications. However, exogenous feeding of chemically synthesized ncAAs at high concentrations is required to compensate for the inefficient cellular uptake and incorporation of these components into proteins, especially in the case of eukaryotic cells and multicellular organisms. To generate organisms capable of autonomously biosynthesizing an ncAA and incorporating it into proteins, we have engineered a metabolic pathway for the synthesis of O-methyltyrosine (OMeY). Specifically, we endowed organisms with a marformycins biosynthetic pathway-derived methyltransferase that efficiently converts tyrosine to OMeY in the presence of the co-factor S-adenosylmethionine. The resulting cells can produce and site-specifically incorporate OMeY into proteins at much higher levels than cells exogenously fed OMeY. To understand the structural basis for the substrate selectivity of the transferase, we solved the X-ray crystal structures of the ligand-free and tyrosine-bound enzymes. Most importantly, we have extended this OMeY biosynthetic system to both mammalian cells and the zebrafish model to enhance the utility of genetic code expansion. The creation of autonomous eukaryotes using a 21st amino acid will make genetic code expansion technology more applicable to multicellular organisms, providing valuable vertebrate models for biological and biomedical research.

Changes in antibody binding and functionality after humanizing a murine scFv anti-IFN-α2: From in silico studies to experimental analysis

The structural and dynamic changes introduced during antibody humanization continue to be a topic open to new contributions. For this reason, the study of structural and functional changes of a murine scFv (mu.scFv) anti-rhIFN-α2b after humanization was carried out. As it was shown by long molecular dynamics simulations and circular dichroism analysis, changes in primary sequence affected the tertiary structure of the humanized scFv (hz.scFv): the position of the variable domain of light chain (VL) respective to the variable domain of heavy chain (VH) in each scFv molecule was different. This change mainly impacted on conformation and dynamics of the complementarity-determining region 3 of VH (CDR-H3) which led to changes in the specificity and affinity of humanized scFv (hz.scFv). These observations agree with experimental results that showed a decrease in the antigen-binding strength of hz.scFv, and different capacities of these molecules to neutralize the in vitro rhIFN-α2b biological activity. Besides, experimental studies to characterize antigen-antibody binding showed that mu.scFv and hz.scFv bind to the same antigen area and recognize a conformational epitope, which is evidence of docking results. Finally, the differences between these molecules to neutralize the in vitro rhIFN-α2b biological activity were described as a consequence of the blockade of certain functionally relevant amino acids of the cytokine, after scFv binding. All these observations confirmed that humanization affected the affinity and specificity of hz.scFv and pointed out that two specific changes in the frameworks would be responsible.

A computational study of the structure and function of human Zrt and Irt-like proteins metal transporters: An elevator-type transport mechanism predicted by AlphaFold2

The ZIP (Zrt and Irt-like proteins) protein family includes transporters responsible for the translocation of zinc and other transition metals, such as iron and cadmium, between the extracellular space (or the lumen of organelles) and the cytoplasm. This protein family is present at all the phylogenetic levels, including bacteria, fungi, plants, insects, and mammals. ZIP proteins are responsible for the homeostasis of metals essential for the cell physiology. The human ZIP family consists of fourteen members (hZIP1-hZIP14), divided into four subfamilies: LIV-1, containing nine hZIPs, the subfamily I, with only one member, the subfamily II, which includes three members and the subfamily gufA, which has only one member. Apart from the extracellular domain, typical of the LIV-1 subfamily, the highly conserved transmembrane domain, containing the binuclear metal center (BMC), and the histidine-rich intracellular loop are the common features characterizing the ZIP family. Here is presented a computational study of the structure and function of human ZIP family members. Multiple sequence alignment and structural models were obtained for the 14 hZIP members. Moreover, a full-length three-dimensional model of the hZIP4-homodimer complex was also produced. Different conformations of the representative hZIP transporters were obtained through a modified version of the AlphaFold2 algorithm. The inward and outward-facing conformations obtained suggest that the hZIP proteins function with an “elevator-type” mechanism.

Comparative evaluation of spin-label modeling methods for protein structural studies

Site-directed spin-labeling electron paramagnetic resonance spectroscopy is a powerful technique for the investigation of protein structure and dynamics. Accurate spin-label modeling methods are essential to make full quantitative use of site-directed spin-labeling electron paramagnetic resonance data for protein modeling and model validation. Using a set of double electron-electron resonance data from seven different site pairs on maltodextrin/maltose-binding protein under two different conditions using five different spin labels, we compare the ability of two widely used spin-label modeling methods, based on accessible volume sampling and rotamer libraries, to predict experimental distance distributions. We present a spin-label modeling approach inspired by canonical side-chain modeling methods and compare modeling accuracy with the established methods.

Conformation-Dependent Influences of Hydrophobic Amino Acids in Two In-Register Parallel β-Sheet Amyloids, an α-Synuclein Amyloid and a Local Structural Model of PrPSc

Prions are unconventional pathogens that encode the pathogenic information in conformations of the constituent abnormal isoform of prion protein (PrP^Sc), independently of the nucleotide genome. Therefore, conformational diversity of PrP^Sc underlies the existence of many prion strains and species barriers of prions, although the conformational information is extremely limited. Interestingly, differences between polymorphic or species-specific residues responsible for the species/strain barriers are often caused by conservative replacements between hydrophobic amino acids. This implies that subtle differences among hydrophobic amino acids are significant for PrP^Sc structures. Here we analyzed the influence of different hydrophobic residues on the structures of an in-register parallel β-sheet amyloid of α-synuclein (αSyn) using molecular dynamics (MD) simulation and applied the knowledge from the αSyn amyloid to modeling a local structure of human PrP^Sc encompassing residues 107-143. We found that mutations equivalent to polymorphisms that cause transmission barriers substantially affect the stabilities of the local structures; for example, the G127V mutation, which makes the host resistant to various human prion diseases, greatly destabilized the local structure of the model amyloid. Our study indicates that subtle differences among hydrophobic side chains can considerably affect the interaction network, including hydrogen bonds, and demonstrates specifically how and in what structures hydrophobic residues can exert unique effects on in-register parallel β-sheet amyloids.

Efficient Strategy to Design Protease Inhibitors: Application to Enterovirus 71 2A Protease

One strategy to counter viruses that persistently cause outbreaks is to design molecules that can specifically inhibit an essential multifunctional viral protease. Herein, we present such a strategy using well-established methods to first identify a region present only in viral (but not human) proteases and find peptides that can bind specifically to this "unique" region by maximizing the protease-peptide binding free energy iteratively using single-point mutations starting with the substrate peptide. We applied this strategy to discover pseudosubstrate peptide inhibitors for the multifunctional 2A protease of enterovirus 71 (EV71), a key causative pathogen for hand-foot-and-mouth disease affecting young children, along with coxsackievirus A16. Four peptide candidates predicted to bind EV71 2A protease more tightly than the natural substrate were experimentally validated and found to inhibit protease activity. Furthermore, the crystal structure of the best pseudosubstrate peptide bound to the EV71 2A protease was determined to provide a molecular basis for the observed inhibition. Since the 2A proteases of EV71 and coxsackievirus A16 share nearly identical sequences and structures, our pseudosubstrate peptide inhibitor may prove useful in inhibiting the two key pathogens of hand-foot-and-mouth disease.

Studying protein-protein interaction through side-chain modeling method OPUS-Mut

Protein side chains are vitally important to many biological processes such as protein-protein interaction. In this study, we evaluate the performance of our previous released side-chain modeling method OPUS-Mut, together with some other methods, on three oligomer datasets, CASP14 (11), CAMEO-Homo (65) and CAMEO-Hetero (21). The results show that OPUS-Mut outperforms other methods measured by all residues or by the interfacial residues. We also demonstrate our method on evaluating protein-protein docking pose on a dataset Oligomer-Dock (75) created using the top 10 predictions from ZDOCK 3.0.2. Our scoring function correctly identifies the native pose as the top-1 in 45 out of 75 targets. Different from traditional scoring functions, our method is based on the overall side-chain packing favorableness in accordance with the local packing environment. It emphasizes the significance of side chains and provides a new and effective scoring term for studying protein-protein interaction.

NanoNet: Rapid and accurate end-to-end nanobody modeling by deep learning

Antibodies are a rapidly growing class of therapeutics. Recently, single domain camelid VHH antibodies, and their recognition nanobody domain (Nb) appeared as a cost-effective highly stable alternative to full-length antibodies. There is a growing need for high-throughput epitope mapping based on accurate structural modeling of the variable domains that share a common fold and differ in the Complementarity Determining Regions (CDRs). We develop a deep learning end-to-end model, NanoNet, that given a sequence directly produces the 3D coordinates of the backbone and Cβ atoms of the entire VH domain. For the Nb test set, NanoNet achieves 3.16Å average RMSD for the most variable CDR3 loops and 2.65Å, 1.73Å for the CDR1, CDR2 loops, respectively. The accuracy for antibody VH domains is even higher: 2.38Å RMSD for CDR3 and 0.89Å, 0.96Å for the CDR1, CDR2 loops, respectively. NanoNet run times allow generation of ∼1M nanobody structures in less than 4 hours on a standard CPU computer enabling high-throughput structure modeling. NanoNet is available at GitHub: https://github.com/dina-lab3D/NanoNet.

Dynamics and structural changes of calmodulin upon interaction with the antagonist calmidazolium

BACKGROUND

Calmodulin (CaM) is an evolutionarily conserved eukaryotic multifunctional protein that functions as the major sensor of intracellular calcium signaling. Its calcium-modulated function regulates the activity of numerous effector proteins involved in a variety of physiological processes in diverse organs, from proliferation and apoptosis, to memory and immune responses. Due to the pleiotropic roles of CaM in normal and pathological cell functions, CaM antagonists are needed for fundamental studies as well as for potential therapeutic applications. Calmidazolium (CDZ) is a potent small molecule antagonist of CaM and one the most widely used inhibitors of CaM in cell biology. Yet, CDZ, as all other CaM antagonists described thus far, also affects additional cellular targets and its lack of selectivity hinders its application for dissecting calcium/CaM signaling. A better understanding of CaM:CDZ interaction is key to design analogs with improved selectivity. Here, we report a molecular characterization of CaM:CDZ complexes using an integrative structural biology approach combining SEC-SAXS, X-ray crystallography, HDX-MS, and NMR.

RESULTS

We provide evidence that binding of a single molecule of CDZ induces an open-to-closed conformational reorientation of the two domains of CaM and results in a strong stabilization of its structural elements associated with a reduction of protein dynamics over a large time range. These CDZ-triggered CaM changes mimic those induced by CaM-binding peptides derived from physiological protein targets, despite their distinct chemical natures. CaM residues in close contact with CDZ and involved in the stabilization of the CaM:CDZ complex have been identified.

CONCLUSION

Our results provide molecular insights into CDZ-induced dynamics and structural changes of CaM leading to its inhibition and open the way to the rational design of more selective CaM antagonists. Calmidazolium is a potent and widely used inhibitor of calmodulin, a major mediator of calcium-signaling in eukaryotic cells. Structural characterization of calmidazolium-binding to calmodulin reveals that it triggers open-to-closed conformational changes similar to those induced by calmodulin-binding peptides derived from enzyme targets. These results provide molecular insights into CDZ-induced dynamics and structural changes of CaM leading to its inhibition and open the way to the rational design of more selective CaM antagonists.

G protein coupling and activation of the metabotropic GABAB heterodimer

Metabotropic γ-aminobutyric acid receptor (GABABR), a class C G protein-coupled receptor (GPCR) heterodimer, plays a crucial role in the central nervous system. Cryo-electron microscopy studies revealed a drastic conformational change upon activation and a unique G protein (GP) binding mode. However, little is known about the mechanism for GP coupling and activation for class C GPCRs. Here, we use molecular metadynamics computations to predict the mechanism by which the inactive GP induces conformational changes in the GABABR transmembrane domain (TMD) to form an intermediate pre-activated state. We find that the inactive GP first interacts with TM3, which further leads to the TMD rearrangement and deeper insertion of the α5 helix that causes the Gα subunit to open, releasing GDP, and forming the experimentally observed activated structure. This mechanism provides fresh insights into the mechanistic details of class C GPCRs activation expected to be useful for designing selective agonists and antagonists.

The Single Residue K12 Governs the Exceptional Voltage Sensitivity of Mitochondrial Voltage-Dependent Anion Channel Gating

The voltage-dependent anion channel (VDAC) is a β-barrel channel of the mitochondrial outer membrane (MOM) that passively transports ions, metabolites, polypeptides, and single-stranded DNA. VDAC responds to a transmembrane potential by "gating," i.e. transitioning to one of a variety of low-conducting states of unknown structure. The gated state results in nearly complete suppression of multivalent mitochondrial metabolite (such as ATP and ADP) transport, while enhancing calcium transport. Voltage gating is a universal property of β-barrel channels, but VDAC gating is anomalously sensitive to transmembrane potential. Here, we show that a single residue in the pore interior, K12, is responsible for most of VDAC's voltage sensitivity. Using the analysis of over 40 μs of atomistic molecular dynamics (MD) simulations, we explore correlations between motions of charged residues inside the VDAC pore and geometric deformations of the β-barrel. Residue K12 is bistable; its motions between two widely separated positions along the pore axis enhance the fluctuations of the β-barrel and augment the likelihood of gating. Single channel electrophysiology of various K12 mutants reveals a dramatic reduction of the voltage-induced gating transitions. The crystal structure of the K12E mutant at a resolution of 2.6 Å indicates a similar architecture of the K12E mutant to the wild type; however, 60 μs of atomistic MD simulations using the K12E mutant show restricted motion of residue 12, due to enhanced connectivity with neighboring residues, and diminished amplitude of barrel motions. We conclude that β-barrel fluctuations, governed particularly by residue K12, drive VDAC gating transitions.

Frustration-driven allosteric regulation and signal transmission in the SARS-CoV-2 spike omicron trimer structures: a crosstalk of the omicron mutation sites allosterically regulates tradeoffs of protein stability and conformational adaptability

Dissecting the regulatory principles underlying function and activity of the SARS-CoV-2 spike protein at the atomic level is of paramount importance for understanding the mechanisms of virus transmissibility and immune escape. In this work, we introduce a hierarchical computational approach for atomistic modeling of allosteric mechanisms in the SARS-CoV-2 Omicron spike proteins and present evidence of a frustration-based allostery as an important energetic driver of the conformational changes and spike activation. By examining conformational landscapes and the residue interaction networks in the SARS-CoV-2 Omicron spike protein structures, we have shown that the Omicron mutational sites are dynamically coupled and form a central engine of the allosterically regulated spike machinery that regulates the balance and tradeoffs between conformational plasticity, protein stability, and functional adaptability. We have found that the Omicron mutational sites at the inter-protomer regions form regulatory hotspot clusters that control functional transitions between the closed and open states. Through perturbation-based modeling of allosteric interaction networks and diffusion analysis of communications in the closed and open spike states, we have quantified the allosterically regulated activation mechanism and uncover specific regulatory roles of the Omicron mutations. Atomistic reconstruction of allosteric communication pathways and kinetic modeling using Markov transient analysis reveal that the Omicron mutations form the inter-protomer electrostatic bridges that operate as a network of coupled regulatory switches that could control global conformational changes and signal transmission in the spike protein. The results of this study have revealed distinct and yet complementary roles of the Omicron mutation sites as a network of hotspots that enable allosteric modulation of structural stability and conformational changes which are central for spike activation and virus transmissibility.

Conformational Dynamics and Mechanisms of Client Protein Integration into the Hsp90 Chaperone Controlled by Allosteric Interactions of Regulatory Switches: Perturbation-Based Network Approach for Mutational Profiling of the Hsp90 Binding and Allostery

Understanding the allosteric mechanisms of the Hsp90 chaperone interactions with cochaperones and client protein clientele is fundamental to dissect activation and regulation of many proteins. In this work, atomistic simulations are combined with perturbation-based approaches and dynamic network modeling for a comparative mutational profiling of the Hsp90 binding and allosteric interaction networks in the three Hsp90 maturation complexes with FKBP51 and P23 cochaperones and the glucocorticoid receptor (GR) client. The conformational dynamics signatures of the Hsp90 complexes and dynamics fluctuation analysis revealed how the intrinsic plasticity of the Hsp90 dimer can be modulated by cochaperones and client proteins to stabilize the closed dimer state required at the maturation stage of the ATPase cycle. In silico deep mutational scanning of the protein residues characterized the hot spots of protein stability and binding affinity in the Hsp90 complexes, showing that binding hot spots may often coincide with the regulatory centers that modulate dynamic allostery in the Hsp90 dimer. We introduce a perturbation-based network approach for mutational scanning of allosteric residue potentials and characterize allosteric switch clusters that control mechanism of cochaperone-dependent client recognition and remodeling by the Hsp90 chaperone. The results revealed a conserved network of allosteric switches in the Hsp90 complexes that allow cochaperones and GR protein to become integrated into the Hsp90 system by anchoring to the conformational switch points in the functional Hsp90 regions. This study suggests that the Hsp90 binding and allostery may operate under a regulatory mechanism in which activation or repression of the Hsp90 activity can be pre-encoded in the allosterically regulated Hsp90 dimer motions. By binding directly to the conformational switch centers on the Hsp90, cochaperones and interacting proteins can efficiently modulate the allosteric interactions and long-range communications required for client remodeling and activation.

Optimizing Epitope Conformational Ensembles Using α-Synuclein Cyclic Peptide "Glycindel" Scaffolds: A Customized Immunogen Method for Generating Oligomer-Selective Antibodies for Parkinson's Disease

Effectively presenting epitopes on immunogens, in order to raise conformationally selective antibodies through active immunization, is a central problem in treating protein misfolding diseases, particularly neurodegenerative diseases such as Alzheimer's disease or Parkinson's disease. We seek to selectively target conformations enriched in toxic, oligomeric propagating species while sparing the healthy forms of the protein that are often more abundant. To this end, we computationally modeled scaffolded epitopes in cyclic peptides by inserting/deleting a variable number of flanking glycines ("glycindels") to best mimic a misfolding-specific conformation of an epitope of α-synuclein enriched in the oligomer ensemble, as characterized by a region most readily disordered and solvent-exposed in a stressed, partially denatured protofibril. We screen and rank the cyclic peptide scaffolds of α-synuclein in silico based on their ensemble overlap properties with the fibril, oligomer-model and isolated monomer ensembles. We present experimental data of seeded aggregation that support nucleation rates consistent with computationally predicted cyclic peptide conformational similarity. We also introduce a method for screening against structured off-pathway targets in the human proteome by selecting scaffolds with minimal conformational similarity between their epitope and the same solvent-exposed primary sequence in structured human proteins. Different cyclic peptide scaffolds with variable numbers of glycines are predicted computationally to have markedly different conformational ensembles. Ensemble comparison and overlap were quantified by the Jensen-Shannon divergence and a new measure introduced here, the embedding depth, which determines the extent to which a given ensemble is subsumed by another ensemble and which may be a more useful measure in developing immunogens that confer conformational selectivity to an antibody.

Integrating Conformational Dynamics and Perturbation-Based Network Modeling for Mutational Profiling of Binding and Allostery in the SARS-CoV-2 Spike Variant Complexes with Antibodies: Balancing Local and Global Determinants of Mutational Escape Mechanisms

In this study, we combined all-atom MD simulations, the ensemble-based mutational scanning of protein stability and binding, and perturbation-based network profiling of allosteric interactions in the SARS-CoV-2 spike complexes with a panel of cross-reactive and ultra-potent single antibodies (B1-182.1 and A23-58.1) as well as antibody combinations (A19-61.1/B1-182.1 and A19-46.1/B1-182.1). Using this approach, we quantify the local and global effects of mutations in the complexes, identify protein stability centers, characterize binding energy hotspots, and predict the allosteric control points of long-range interactions and communications. Conformational dynamics and distance fluctuation analysis revealed the antibody-specific signatures of protein stability and flexibility of the spike complexes that can affect the pattern of mutational escape. A network-based perturbation approach for mutational profiling of allosteric residue potentials revealed how antibody binding can modulate allosteric interactions and identified allosteric control points that can form vulnerable sites for mutational escape. The results show that the protein stability and binding energetics of the SARS-CoV-2 spike complexes with the panel of ultrapotent antibodies are tolerant to the effect of Omicron mutations, which may be related to their neutralization efficiency. By employing an integrated analysis of conformational dynamics, binding energetics, and allosteric interactions, we found that the antibodies that neutralize the Omicron spike variant mediate the dominant binding energy hotpots in the conserved stability centers and allosteric control points in which mutations may be restricted by the requirements of the protein folding stability and binding to the host receptor. This study suggested a mechanism in which the patterns of escape mutants for the ultrapotent antibodies may not be solely determined by the binding interaction changes but are associated with the balance and tradeoffs of multiple local and global factors, including protein stability, binding affinity, and long-range interactions.

Extensive Sampling of Molecular Dynamics Simulations to Identify Reliable Protein Structures for Optimized Virtual Screening Studies: The Case of the hTRPM8 Channel

(1) Background: Virtual screening campaigns require target structures in which the pockets are properly arranged for binding. Without these, MD simulations can be used to relax the available target structures, optimizing the fine architecture of their binding sites. Among the generated frames, the best structures can be selected based on available experimental data. Without experimental templates, the MD trajectories can be filtered by energy-based criteria or sampled by systematic analyses. (2) Methods: A blind and methodical analysis was performed on the already reported MD run of the hTRPM8 tetrameric structures; a total of 50 frames underwent docking simulations by using a set of 1000 ligands including 20 known hTRPM8 modulators. Docking runs were performed by LiGen program and involved the frames as they are and after optimization by SCRWL4.0. For each frame, all four monomers were considered. Predictive models were developed by the EFO algorithm based on the sole primary LiGen scores. (3) Results: On average, the MD simulation progressively enhances the performance of the extracted frames, and the optimized structures perform better than the non-optimized frames (EF1% mean: 21.38 vs. 23.29). There is an overall correlation between performances and volumes of the explored pockets and the combination of the best performing frames allows to develop highly performing consensus models (EF1% = 49.83). (4) Conclusions: The systematic sampling of the entire MD run provides performances roughly comparable with those previously reached by using rationally selected frames. The proposed strategy appears to be helpful when the lack of experimental data does not allow an easy selection of the optimal structures for docking simulations. Overall, the reported docking results confirm the relevance of simulating all the monomers of an oligomer structure and emphasize the efficacy of the SCRWL4.0 method to optimize the protein structures for docking calculations.

Rotamer-free protein sequence design based on deep learning and self-consistency

Several previously proposed deep learning methods to design amino acid sequences that autonomously fold into a given protein backbone yielded promising results in computational tests but did not outperform conventional energy function-based methods in wet experiments. Here we present the ABACUS-R method, which uses an encoder-decoder network trained using a multitask learning strategy to predict the sidechain type of a central residue from its three-dimensional local environment, which includes, besides other features, the types but not the conformations of the surrounding sidechains. This eliminates the need to reconstruct and optimize sidechain structures, and drastically simplifies the sequence design process. Thus iteratively applying the encoder-decoder to different central residues is able to produce self-consistent overall sequences for a target backbone. Results of wet experiments, including five structures solved by X-ray crystallography, show that ABACUS-R outperforms state-of-the-art energy function-based methods in success rate and design precision.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

A Method for Assessing the Robustness of Protein Structures by Randomizing Packing Interactions

Many single-domain proteins are not only stable and water-soluble, but they also populate few to no intermediates during folding. This reduces interactions between partially folded proteins, misfolding, and aggregation, and makes the proteins tractable in biotechnological applications. Natural proteins fold thus, not necessarily only because their structures are well-suited for folding, but because their sequences optimize packing and fit their structures well. In contrast, folding experiments on the de novo designed Top7 suggest that it populates several intermediates. Additionally, in de novo protein design, where sequences are designed for natural and new non-natural structures, tens of sequences still need to be tested before success is achieved. Both these issues may be caused by the specific scaffolds used in design, i.e., some protein scaffolds may be more tolerant to packing perturbations and varied sequences. Here, we report a computational method for assessing the response of protein structures to packing perturbations. We then benchmark this method using designed proteins and find that it can identify scaffolds whose folding gets disrupted upon perturbing packing, leading to the population of intermediates. The method can also isolate regions of both natural and designed scaffolds that are sensitive to such perturbations and identify contacts which when present can rescue folding. Overall, this method can be used to identify protein scaffolds that are more amenable to whole protein design as well as to identify protein regions which are sensitive to perturbations and where further mutations should be avoided during protein engineering.

Simultaneous prediction of antibody backbone and side-chain conformations with deep learning

Antibody engineering is becoming increasingly popular in medicine for the development of diagnostics and immunotherapies. Antibody function relies largely on the recognition and binding of antigenic epitopes via the loops in the complementarity determining regions. Hence, accurate high-resolution modeling of these loops is essential for effective antibody engineering and design. Deep learning methods have previously been shown to effectively predict antibody backbone structures described as a set of inter-residue distances and orientations. However, antigen binding is also dependent on the specific conformations of surface side-chains. To address this shortcoming, we created DeepSCAb: a deep learning method that predicts inter-residue geometries as well as side-chain dihedrals of the antibody variable fragment. The network requires only sequence as input, rendering it particularly useful for antibodies without any known backbone conformations. Rotamer predictions use an interpretable self-attention layer, which learns to identify structurally conserved anchor positions across several species. We evaluate the performance of the model for discriminating near-native structures from sets of decoys and find that DeepSCAb outperforms similar methods lacking side-chain context. When compared to alternative rotamer repacking methods, which require an input backbone structure, DeepSCAb predicts side-chain conformations competitively. Our findings suggest that DeepSCAb improves antibody structure prediction with accurate side-chain modeling and is adaptable to applications in docking of antibody-antigen complexes and design of new therapeutic antibody sequences.

ColabFold: making protein folding accessible to all

ColabFold offers accelerated prediction of protein structures and complexes by combining the fast homology search of MMseqs2 with AlphaFold2 or RoseTTAFold. ColabFold's 40-60-fold faster search and optimized model utilization enables prediction of close to 1,000 structures per day on a server with one graphics processing unit. Coupled with Google Colaboratory, ColabFold becomes a free and accessible platform for protein folding. ColabFold is open-source software available at https://github.com/sokrypton/ColabFold and its novel environmental databases are available at https://colabfold.mmseqs.com .

ColabFold: making protein folding accessible to all

ColabFold offers accelerated prediction of protein structures and complexes by combining the fast homology search of MMseqs2 with AlphaFold2 or RoseTTAFold. ColabFold's 40-60-fold faster search and optimized model utilization enables prediction of close to 1,000 structures per day on a server with one graphics processing unit. Coupled with Google Colaboratory, ColabFold becomes a free and accessible platform for protein folding. ColabFold is open-source software available at https://github.com/sokrypton/ColabFold and its novel environmental databases are available at https://colabfold.mmseqs.com .

DLPacker: Deep learning for prediction of amino acid side chain conformations in proteins

Prediction of side chain conformations of amino acids in proteins (also termed "packing") is an important and challenging part of protein structure prediction with many interesting applications in protein design. A variety of methods for packing have been developed but more accurate ones are still needed. Machine learning (ML) methods have recently become a powerful tool for solving various problems in diverse areas of science, including structural biology. In this study, we evaluate the potential of deep neural networks (DNNs) for prediction of amino acid side chain conformations. We formulate the problem as image-to-image transformation and train a U-net style DNN to solve the problem. We show that our method outperforms other physics-based methods by a significant margin: reconstruction RMSDs for most amino acids are about 20% smaller compared to SCWRL4 and Rosetta Packer with RMSDs for bulky hydrophobic amino acids Phe, Tyr, and Trp being up to 50% smaller.

Insights into the client protein release mechanism of the ATP-independent chaperone Spy

Molecular chaperones play a central role in regulating protein homeostasis, and their active forms often contain intrinsically disordered regions (IDRs). However, how IDRs impact chaperone action remains poorly understood. Here, we discover that the disordered N terminus of the prototype chaperone Spy facilitates client release. With NMR spectroscopy and molecular dynamics simulations, we find that the N terminus can bind transiently to the client-binding cavity of Spy primarily through electrostatic interactions mediated by the N-terminal D26 residue. This intramolecular interaction results in a dynamic competition of the N terminus with the client for binding to Spy, which promotes client discharge. Our results reveal the mechanism by which Spy releases clients independent of energy input, thus enriching the current knowledge on how ATP-independent chaperones release their clients and highlighting the importance of synergy between IDRs and structural domains in regulating protein function.

Shape Complementarity Optimization of Antibody–Antigen Interfaces: The Application to SARS-CoV-2 Spike Protein

Many factors influence biomolecule binding, and its assessment constitutes an elusive challenge in computational structural biology. In this aspect, the evaluation of shape complementarity at molecular interfaces is one of the main factors to be considered. We focus on the particular case of antibody–antigen complexes to quantify the complementarities occurring at molecular interfaces. We relied on a method we recently developed, which employs the 2D Zernike descriptors, to characterize the investigated regions with an ordered set of numbers summarizing the local shape properties. Collecting a structural dataset of antibody–antigen complexes, we applied this method and we statistically distinguished, in terms of shape complementarity, pairs of the interacting regions from the non-interacting ones. Thus, we set up a novel computational strategy based on in silico mutagenesis of antibody-binding site residues. We developed a Monte Carlo procedure to increase the shape complementarity between the antibody paratope and a given epitope on a target protein surface. We applied our protocol against several molecular targets in SARS-CoV-2 spike protein, known to be indispensable for viral cell invasion. We, therefore, optimized the shape of template antibodies for the interaction with such regions. As the last step of our procedure, we performed an independent molecular docking validation of the results of our Monte Carlo simulations.

Consecutive Aromatic Residues Are Required for Improved Efficacy of β-Sheet Breakers

Alzheimer’s disease is a fatal neurodegenerative malady which up to very recently did not have approved therapy modifying its course. After controversial approval of aducanumab (monoclonal antibody clearing β-amyloid plaques) by FDA for use in very early stages of disease, possibly new avenue opened for the treatment of patients. In line with this approach is search for compounds blocking aggregation into amyloid oligomers subsequently forming fibrils or compounds helping in getting rid of plaques formed by β-amyloid fibrils. Here we present in silico work on 627 sixtapeptide β-sheet breakers (BSBs) containing consecutive three aromatic residues. Three of these BSBs caused dissociation of one or two β-amyloid chains from U-shaped β-amyloid protofibril model 2BEG after docking and subsequent molecular dynamics simulations. Thorough analysis of our results let us postulate that the first steps of binding these successful BSBs involve π–π interactions with stacked chains of F19 and later also with F20 (F3 and F4 in 2BEG model of protofibril). The consecutive location of aromatic residues in BSBs makes them more attractive for chains of stacked F3 and F4 within the 2BEG model. Spotted by us, BSBs may be prospective lead compounds for an anti-Alzheimer’s therapy.

Mimetic Neural Networks: A Unified Framework for Protein Design and Folding

Recent advancements in machine learning techniques for protein structure prediction motivate better results in its inverse problem-protein design. In this work we introduce a new graph mimetic neural network, MimNet, and show that it is possible to build a reversible architecture that solves the structure and design problems in tandem, allowing to improve protein backbone design when the structure is better estimated. We use the ProteinNet data set and show that the state of the art results in protein design can be met and even improved, given recent architectures for protein folding.

Protposer: the web server that readily proposes protein stabilizing mutations with high PPV

Protein stability is a requisite for most biotechnological and medical applications of proteins. As natural proteins tend to suffer from a low conformational stability ex vivo, great efforts have been devoted toward increasing their stability through rational design and engineering of appropriate mutations. Unfortunately, even the best currently used predictors fail to compute the stability of protein variants with sufficient accuracy and their usefulness as tools to guide the rational stabilisation of proteins is limited. We present here Protposer, a protein stabilising tool based on a different approach. Instead of quantifying changes in stability, Protposer uses structure- and sequence-based screening modules to nominate candidate mutations for subsequent evaluation by a logistic regression model, carefully trained to avoid overfitting. Thus, Protposer analyses PDB files in search for stabilization opportunities and provides a ranked list of promising mutations with their estimated success rates (eSR), their probabilities of being stabilising by at least 0.5 kcal/mol. The agreement between eSRs and actual positive predictive values (PPV) on external datasets of mutations is excellent. When Protposer is used with its Optimal kappa selection threshold, its PPV is above 0.7. Even with less stringent thresholds, Protposer largely outperforms FoldX, Rosetta and PoPMusiC. Indicating the PDB file of the protein suffices to obtain a ranked list of mutations, their eSRs and hints on the likely source of the stabilization expected. Protposer is a distinct, straightforward and highly successful tool to design protein stabilising mutations, and it is freely available for academic use at http://webapps.bifi.es/the-protposer.

Rate-limiting transport of positively charged arginine residues through the Sec-machinery is integral to the mechanism of protein secretion

Transport of proteins across and into membranes is a fundamental biological process with the vast majority being conducted by the ubiquitous Sec machinery. In bacteria, this is usually achieved when the SecY-complex engages the cytosolic ATPase SecA (secretion) or translating ribosomes (insertion). Great strides have been made towards understanding the mechanism of protein translocation. Yet, important questions remain - notably, the nature of the individual steps that constitute transport, and how the proton-motive force (PMF) across the plasma membrane contributes. Here, we apply a recently developed high-resolution protein transport assay to explore these questions. We find that pre-protein transport is limited primarily by the diffusion of arginine residues across the membrane, particularly in the context of bulky hydrophobic sequences. This specific effect of arginine, caused by its positive charge, is mitigated for lysine which can be deprotonated and transported across the membrane in its neutral form. These observations have interesting implications for the mechanism of protein secretion, suggesting a simple mechanism through which the PMF can aid transport by enabling a 'proton ratchet', wherein re-protonation of exiting lysine residues prevents channel re-entry, biasing transport in the outward direction.

Inferring Conformational State of Myosin Motor in an Atomic Force Microscopy Image via Flexible Fitting Molecular Simulations

High-speed atomic force microscopy (HS-AFM) is a powerful technique to image the structural dynamics of biomolecules. We can obtain atomic-resolution structural information from the measured AFM image by superimposing a structural model on the image. We previously developed a flexible fitting molecular dynamics (MD) simulation method that allows for modest conformational changes when superimposed on an AFM image. In this study, for a molecular motor, myosin V (which changes its chemical state), we examined whether the conformationally distinct state in each HS-AFM image could be inferred via flexible fitting MD simulation. We first built models of myosin V bound to the actin filament in two conformational states, the “down-up” and “down-down” states. Then, for the previously obtained HS-AFM image of myosin bound to the actin filament, we performed flexible-fitting MD simulations using the two states. By comparing the fitting results, we inferred the conformational and chemical states from the AFM image.

Integrative ensemble modeling of proteins and their complexes with distance distribution restraints

Many proteins and protein complexes exhibit regions that are intrinsically disordered. Whereas an arsenal of techniques exists to characterize structured proteins or protein regions, characterization of the vast conformational space occupied by intrinsically disordered regions remains a challenging task due the ensemble-averaging nature of many techniques that provide mean value restraints. More representative information can be gained in the form of distribution restraints, such as EPR-derived distance distributions. Previously we developed the ensemble modeling tool MMM, where we partition the macromolecule into structured and unstructured domains and utilize an integrative structural approach with a focus on EPR-derived distance restraints. Here we present the successor program of MMM: MMMx. All the modeling functionality was ported to MMMx and is now accessed by a uniform script format, allowing to combine the different modules at will to modeling pipelines. During the conception of MMMx many of the tools were improved or updated. We discuss the general functionality of MMMx and its modules, and illustrate some of the modeling tools by application examples.

Characterization of Opposing Responses to Phenol by Bacillus subtilis Chemoreceptors

Bacillus subtilis employs 10 chemoreceptors to move in response to chemicals in its environment. While the sensing mechanisms have been determined for many attractants, little is known about the sensing mechanisms for repellents. In this work, we investigated phenol chemotaxis in B. subtilis. Phenol is an attractant at low, micromolar concentrations and a repellent at high, millimolar concentrations. McpA was found to be the principal chemoreceptor governing the repellent response to phenol and other related aromatic compounds. In addition, the chemoreceptors McpC and HemAT were found to govern the attractant response to phenol and related compounds. Using chemoreceptor chimeras, McpA was found to sense phenol using its signaling domain rather than its sensing domain. These observations were substantiated in vitro, where direct binding of phenol to the signaling domain of McpA was observed using saturation transfer difference nuclear magnetic resonance. These results further advance our understanding of B. subtilis chemotaxis and further demonstrate that the signaling domain of B. subtilis chemoreceptors can directly sense chemoeffectors. IMPORTANCE Bacterial chemotaxis is commonly thought to employ a sensing mechanism involving the extracellular sensing domain of chemoreceptors. Some ligands, however, appear to be sensed by the signaling domain. Phenolic compounds, commonly found in soil and root exudates, provide environmental cues for soil microbes like Bacillus subtilis. We show that phenol is sensed as both an attractant and a repellent. While the mechanism for sensing phenol as an attractant is still unknown, we found that phenol is sensed as a repellent by the signaling domain of the chemoreceptor McpA. This study furthers our understanding of the unconventional sensing mechanisms employed by the B. subtilis chemotaxis pathway.

Computer Simulations and Network-Based Profiling of Binding and Allosteric Interactions of SARS-CoV-2 Spike Variant Complexes and the Host Receptor: Dissecting the Mechanistic Effects of the Delta and Omicron Mutations

In this study, we combine all-atom MD simulations and comprehensive mutational scanning of S-RBD complexes with the angiotensin-converting enzyme 2 (ACE2) host receptor in the native form as well as the S-RBD Delta and Omicron variants to (a) examine the differences in the dynamic signatures of the S-RBD complexes and (b) identify the critical binding hotspots and sensitivity of the mutational positions. We also examined the differences in allosteric interactions and communications in the S-RBD complexes for the Delta and Omicron variants. Through the perturbation-based scanning of the allosteric propensities of the SARS-CoV-2 S-RBD residues and dynamics-based network centrality and community analyses, we characterize the global mediating centers in the complexes and the nature of local stabilizing communities. We show that a constellation of mutational sites (G496S, Q498R, N501Y and Y505H) correspond to key binding energy hotspots and also contribute decisively to the key interfacial communities that mediate allosteric communications between S-RBD and ACE2. These Omicron mutations are responsible for both favorable local binding interactions and long-range allosteric interactions, providing key functional centers that mediate the high transmissibility of the virus. At the same time, our results show that other mutational sites could provide a "flexible shield" surrounding the stable community network, thereby allowing the Omicron virus to modulate immune evasion at different epitopes, while protecting the integrity of binding and allosteric interactions in the RBD-ACE2 complexes. This study suggests that the SARS-CoV-2 S protein may exploit the plasticity of the RBD to generate escape mutants, while engaging a small group of functional hotspots to mediate efficient local binding interactions and long-range allosteric communications with ACE2.

Landscape-Based Protein Stability Analysis and Network Modeling of Multiple Conformational States of the SARS-CoV-2 Spike D614G Mutant: Conformational Plasticity and Frustration-Induced Allostery as Energetic Drivers of Highly Transmissible Spike Variants

The structural and functional studies of the SARS-CoV-2 spike protein variants revealed an important role of the D614G mutation that is shared across many variants of concern (VOCs), suggesting the effect of this mutation on the enhanced virus infectivity and transmissibility. The recent structural and biophysical studies provided important evidence about multiple conformational substates of the D614G spike protein. The development of a plausible mechanistic model that can explain the experimental observations from a more unified thermodynamic perspective is an important objective of the current work. In this study, we employed efficient and accurate coarse-grained simulations of multiple structural substates of the D614G spike trimers together with the ensemble-based mutational frustration analysis to characterize the dynamics signatures of the conformational landscapes. By combining the local frustration profiling of the conformational states with residue-based mutational scanning of protein stability and network analysis of allosteric interactions and communications, we determine the patterns of mutational sensitivity in the functional regions and sites of variants. We found that the D614G mutation may induce a considerable conformational adaptability of the open states in the SARS-CoV-2 spike protein without compromising the folding stability and integrity of the spike protein. The results suggest that the D614G mutant may employ a hinge-shift mechanism in which the dynamic couplings between the site of mutation and the interprotomer hinge modulate the interdomain interactions, global mobility change, and the increased stability of the open form. This study proposes that mutation-induced modulation of the conformational flexibility and energetic frustration at the interprotomer interfaces may serve as an efficient mechanism for allosteric regulation of the SARS-CoV-2 spike proteins.

Structure-based dual affinity optimization of a SARS-CoV-1/2 cross-reactive single-domain antibody

The SARS coronavirus 2 (SARS-CoV-2) spike (S) protein binding to the human ACE2 receptor is the molecular event that initiates viral entry into host cells and leads to infection and virus replication. There is a need for agents blocking viral entry into host cells that are cross-reactive with emerging virus variants. VHH-72 is an anti-SARS-CoV-1 single-domain antibody that also exhibits cross-specificity with SARS-CoV-2 but with decreased binding affinity. Here we applied a structure-based approach to affinity-mature VHH-72 for the SARS-CoV-2 spike protein while retaining the original affinity for SARS-CoV-1. This was achieved by employing the computational platform ADAPT in a constrained dual-affinity optimization mode as a means of broadening specificity. Select mutants designed by ADAPT were formatted as fusions with a human IgG1-Fc fragment. These mutants demonstrated improved binding to the SARS-CoV-2 spike protein due to decreased dissociation rates. Functional testing for virus neutralization revealed improvements relative to the parental VHH72-Fc up to 10-fold using a SARS-CoV-2 pseudotyped lentivirus and 20-fold against the SARS-CoV-2 authentic live virus (Wuhan variant). Binding and neutralization improvements were maintained for some other SARS-CoV-2 variants currently in circulation. These improved VHH-72 mutants are predicted to establish novel interactions with the S antigen. They will be useful, alone or as fusions with other functional modules, in the global quest for treatments of COVID-19 infections.

Integrated unsupervised-supervised modeling and prediction of protein-peptide affinities at structural level

Cell signal networks are orchestrated directly or indirectly by various peptide-mediated protein-protein interactions, which are normally weak and transient and thus ideal for biological regulation and medicinal intervention. Here, we develop a general-purpose method for modeling and predicting the binding affinities of protein-peptide interactions (PpIs) at the structural level. The method is a hybrid strategy that employs an unsupervised approach to derive a layered PpI atom-residue interaction (ulPpI[a-r]) potential between different protein atom types and peptide residue types from thousands of solved PpI complex structures and then statistically correlates the potential descriptors with experimental affinities (KD values) over hundreds of known PpI samples in a supervised manner to create an integrated unsupervised-supervised PpI affinity (usPpIA) predictor. Although both the ulPpI[a-r] potential and usPpIA predictor can be used to calculate PpI affinities from their complex structures, the latter seems to perform much better than the former, suggesting that the unsupervised potential can be improved substantially with a further correction by supervised statistical learning. We examine the robustness and fault-tolerance of usPpIA predictor when applied to treat the coarse-grained PpI complex structures modeled computationally by sophisticated peptide docking and dynamics simulation. It is revealed that, despite developed solely based on solved structures, the integrated unsupervised-supervised method is also applicable for locally docked structures to reach a quantitative prediction but can only give a qualitative prediction on globally docked structures. The dynamics refinement seems not to change (or improve) the predictive results essentially, although it is computationally expensive and time-consuming relative to peptide docking. We also perform extrapolation of usPpIA predictor to the indirect affinity quantities of HLA-A*0201 binding epitope peptides and NHERF PDZ binding scaffold peptides, consequently resulting in a good and moderate correlation of the predicted KD with experimental IC50 and BLU on the two peptide sets, with Pearson's correlation coefficients Rp = 0.635 and 0.406, respectively.

Epitope-based peptide vaccine design and elucidation of novel compounds against 3C like protein of SARS-CoV-2

Coronaviruses (CoVs) are positive-stranded RNA viruses with short clubs on their edges. CoVs are pathogenic viruses that infect several animals and plant organisms, as well as humans (lethal respiratory dysfunctions). A noval strain of CoV has been reported and named as SARS-CoV-2. Numerous COVID-19 cases were being reported all over the World. COVID-19 and has a high mortality rate. In the present study, immunoinformatics techniques were utilized to predict the antigenic epitopes against 3C like protein. B-cell epitopes and Cytotoxic T-lymphocyte (CTL) were designed computationally against SARS-CoV-2. Multiple Sequence Alignment (MSA) of seven complete strains (HCoV-229E, HCoV-NL63, HCoV-OC43, HCoV-HKU1, SARS-CoV, MERS-CoV, and SARS-CoV-2) was performed to elucidate the binding domain and interacting residues. MHC-I binding epitopes were evaluated by analyzing the binding affinity of the top-ranked peptides having HLA molecule. By utilizing the docked complexes of CTL epitopes with antigenic sites, the binding relationship and affinity of top-ranked predicted peptides with the MHC-I HLA protein were investigated. The molecular docking analyses were conducted on the ZINC database library and twelve compounds having least binding energy were scrutinized. In conclusion, twelve CTL epitopes (GTDLEGNFY, TVNVLAWLY, GSVGFNIDY, SEDMLNPNY, LSQTGIAV, VLDMCASLK, LTQDHVDIL, TTLNDFNLV, CTSEDMLNP, TTITVNVLA, YNGSPSGVY, and SMQNCVLKL) were identified against SARS-CoV-2.

Enhanced Conformational Sampling with an Adaptive Coarse-Grained Elastic Network Model Using Short-Time All-Atom Molecular Dynamics

Compared to all-atom molecular dynamics (AA-MD) simulations, coarse-grained (CG) MD simulations can significantly reduce calculation costs. However, existing CG-MD methods are unsuitable for sampling structures that depart significantly from the initial structure without any biased force. In this study, we developed a new adaptive CG elastic network model (ENM), in which the dynamic cross-correlation coefficient based on short-time AA-MD of at most ns order is considered. By applying Bayesian optimization to search for a suitable parameter among the vast parameter space of adaptive CG-ENM, we succeeded in reducing the searching cost to approximately 10% of those for random sampling and exhaustive sampling. To evaluate the performance of adaptive CG-ENM, we applied the new methodology to adenylate kinase (ADK) and glutamine binding protein (GBP) in the apo state. The results showed that the structural ensembles explored by adaptive CG-ENM could be considerably more diverse than those by conventional ENMs with enhanced sampling such as temperature replica exchange MD and long-time AA-MD of 1 μs. In particular, some of the structures sampled by adaptive ENM are relatively close to the holo-type structures of ADK and GBP. Furthermore, as a challenging task, to demonstrate the advantages of the CG model with lower calculation cost, we applied our new methodology to a larger biomolecule, integrin (αV) in the inactive state. Then, we sampled various structural ensembles, including extended structures that are apparently different from inactive ones.