Modeling of Disordered Protein Structures Using Monte Carlo Simulations and Knowledge-Based Statistical Force Fields

Structure prediction of linear and cyclic peptides using CABS-flex

The structural modeling of peptides can be a useful aid in the discovery of new drugs and a deeper understanding of the molecular mechanisms of life. Here we present a novel multiscale protocol for the structure prediction of linear and cyclic peptides. The protocol combines two main stages: coarse-grained simulations using the CABS-flex standalone package and an all-atom reconstruction-optimization process using the Modeller program. We evaluated the protocol on a set of linear peptides and two sets of cyclic peptides, with cyclization through the backbone and disulfide bonds. A comparison with other state-of-the-art tools (APPTEST, PEP-FOLD, ESMFold and AlphaFold implementation in ColabFold) shows that for most cases, AlphaFold offers the highest resolution. However, CABS-flex is competitive, particularly when it comes to short linear peptides. As demonstrated, the protocol performance can be further improved by combination with the residue-residue contact prediction method or more efficient scoring. The protocol is included in the CABS-flex standalone package along with online documentation to aid users in predicting the structure of peptides and mini-proteins.

Possibilities of Using De Novo Design for Generating Diverse Functional Food Enzymes

Food enzymes have an important role in the improvement of certain food characteristics, such as texture improvement, elimination of toxins and allergens, production of carbohydrates, enhancing flavor/appearance characteristics. Recently, along with the development of artificial meats, food enzymes have been employed to achieve more diverse functions, especially in converting non-edible biomass to delicious foods. Reported food enzyme modifications for specific applications have highlighted the significance of enzyme engineering. However, using direct evolution or rational design showed inherent limitations due to the mutation rates, which made it difficult to satisfy the stability or specific activity needs for certain applications. Generating functional enzymes using de novo design, which highly assembles naturally existing enzymes, provides potential solutions for screening desired enzymes. Here, we describe the functions and applications of food enzymes to introduce the need for food enzymes engineering. To illustrate the possibilities of using de novo design for generating diverse functional proteins, we reviewed protein modelling and de novo design methods and their implementations. The future directions for adding structural data for de novo design model training, acquiring diversified training data, and investigating the relationship between enzyme-substrate binding and activity were highlighted as challenges to overcome for the de novo design of food enzymes.

Computational analysis of protein stability and allosteric interaction networks in distinct conformational forms of the SARS-CoV-2 spike D614G mutant: reconciling functional mechanisms through allosteric model of spike regulation

In this study, we used an integrative computational approach to examine molecular mechanisms underlying functional effects of the D614G mutation by exploring atomistic modeling of the SARS-CoV-2 spike proteins as allosteric regulatory machines. We combined coarse-grained simulations, protein stability and dynamic fluctuation communication analysis with network-based community analysis to examine structures of the native and mutant SARS-CoV-2 spike proteins in different functional states. Through distance fluctuations communication analysis, we probed stability and allosteric communication propensities of protein residues in the native and mutant SARS-CoV-2 spike proteins, providing evidence that the D614G mutation can enhance long-range signaling of the allosteric spike engine. By combining functional dynamics analysis and ensemble-based alanine scanning of the SARS-CoV-2 spike proteins we found that the D614G mutation can improve stability of the spike protein in both closed and open forms, but shifting thermodynamic preferences towards the open mutant form. Our results revealed that the D614G mutation can promote the increased number of stable communities and allosteric hub centers in the open form by reorganizing and enhancing the stability of the S1-S2 inter-domain interactions and restricting mobility of the S1 regions. This study provides atomistic-based view of allosteric communications in the SARS-CoV-2 spike proteins, suggesting that the D614G mutation can exert its primary effect through allosterically induced changes on stability and communications in the residue interaction networks.Communicated by Ramaswamy H. Sarma.

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.

AlphaFold 2 and NMR Spectroscopy: Partners to Understand Protein Structure, Dynamics and Function

The artificial intelligence program AlphaFold 2 is revolutionizing the field of protein structure determination as it accurately predicts the 3D structure of two thirds of the human proteome. Its predictions can be used directly as structural models or indirectly as aids for experimental structure determination using X-ray crystallography, CryoEM or NMR spectroscopy. Nevertheless, AlphaFold 2 can neither afford insight into how proteins fold, nor can it determine protein stability or dynamics. Rare folds or minor alternative conformations are also not predicted by AlphaFold 2 and the program does not forecast the impact of post translational modifications, mutations or ligand binding. The remaining third of human proteome which is poorly predicted largely corresponds to intrinsically disordered regions of proteins. Key to regulation and signaling networks, these disordered regions often form biomolecular condensates or amyloids. Fortunately, the limitations of AlphaFold 2 are largely complemented by NMR spectroscopy. This experimental approach provides information on protein folding and dynamics as well as biomolecular condensates and amyloids and their modulation by experimental conditions, small molecules, post translational modifications, mutations, flanking sequence, interactions with other proteins, RNA and virus. Together, NMR spectroscopy and AlphaFold 2 can collaborate to advance our comprehension of proteins.

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.

Efficient Monte Carlo Sampling for Molecular Systems Using Continuous Normalizing Flow

Monte Carlo molecular simulation is a powerful computational method for simulating molecular behavior. It generates samples of the possible states of molecular systems. To generate a sample efficiently, it is advantageous to avoid suggesting extremely high-energy states that would never become possible states. In this study, we propose a new sampling method for Monte Carlo molecular simulation, that is, a continuous normalizing molecular flow (CNMF) method, which can create various probabilistic distributions of molecular states from some initial distribution. The CNMF method generates samples by solving a first-order differential equation with two-body intermolecular interaction terms. We also develop specific probabilistic distributions using CNMF called inverse square flow, which yields distributions with zero probability density when molecule pairs are in close proximity, whereas probability densities are compressed uniformly from the initial distribution in all other cases. Using inverse square flow, we demonstrate that Monte Carlo molecular simulation is more efficient than the standard simulation. Although the increased computational costs of the CNMF method are non-negligible, this method is feasible for parallel computation and has the potential for expansion.

Allosteric Determinants of the SARS-CoV-2 Spike Protein Binding with Nanobodies: Examining Mechanisms of Mutational Escape and Sensitivity of the Omicron Variant

Structural and biochemical studies have recently revealed a range of rationally engineered nanobodies with efficient neutralizing capacity against the SARS-CoV-2 virus and resilience against mutational escape. In this study, we performed a comprehensive computational analysis of the SARS-CoV-2 spike trimer complexes with single nanobodies Nb6, VHH E, and complex with VHH E/VHH V nanobody combination. We combined coarse-grained and all-atom molecular simulations and collective dynamics analysis with binding free energy scanning, perturbation-response scanning, and network centrality analysis to examine mechanisms of nanobody-induced allosteric modulation and cooperativity in the SARS-CoV-2 spike trimer complexes with these nanobodies. By quantifying energetic and allosteric determinants of the SARS-CoV-2 spike protein binding with nanobodies, we also examined nanobody-induced modulation of escaping mutations and the effect of the Omicron variant on nanobody binding. The mutational scanning analysis supported the notion that E484A mutation can have a significant detrimental effect on nanobody binding and result in Omicron-induced escape from nanobody neutralization. Our findings showed that SARS-CoV-2 spike protein might exploit the plasticity of specific allosteric hotspots to generate escape mutants that alter response to binding without compromising activity. The network analysis supported these findings showing that VHH E/VHH V nanobody binding can induce long-range couplings between the cryptic binding epitope and ACE2-binding site through a broader ensemble of communication paths that is less dependent on specific mediating centers and therefore may be less sensitive to mutational perturbations of functional residues. The results suggest that binding affinity and long-range communications of the SARS-CoV-2 complexes with nanobodies can be determined by structurally stable regulatory centers and conformationally adaptable hotspots that are allosterically coupled and collectively control resilience to mutational escape.

Conformational Flexibility and Local Frustration in the Functional States of the SARS-CoV-2 Spike B.1.1.7 and B.1.351 Variants: Mutation-Induced Allosteric Modulation Mechanism of Functional Dynamics and Protein Stability

Structural and functional studies of the SARS-CoV-2 spike proteins have recently determined distinct functional states of the B.1.1.7 and B.1.351 spike variants, providing a molecular framework for understanding the mechanisms that link the effect of mutations with the enhanced virus infectivity and transmissibility. A detailed dynamic and energetic analysis of these variants was undertaken in the present work to quantify the effects of different mutations on functional conformational changes and stability of the SARS-CoV-2 spike protein. We employed the efficient and accurate coarse-grained (CG) simulations of multiple functional states of the D614G mutant, B.1.1.7 and B.1.351 spike variants to characterize conformational dynamics of the SARS-CoV-2 spike proteins and identify dynamic signatures of the functional regions that regulate transitions between the closed and open forms. By combining molecular simulations with full atomistic reconstruction of the trajectories and the ensemble-based mutational frustration analysis, we characterized how the intrinsic flexibility of specific spike regions can control functional conformational changes required for binding with the host-cell receptor. Using the residue-based mutational scanning of protein stability, we determined protein stability hotspots and identified potential energetic drivers favoring the receptor-accessible open spike states for the B.1.1.7 and B.1.351 spike variants. The results suggested that modulation of the energetic frustration at the inter-protomer interfaces can serve as a mechanism for allosteric couplings between mutational sites and the inter-protomer hinges of functional motions. The proposed mechanism of mutation-induced energetic frustration may result in greater adaptability and the emergence of multiple conformational states in the open form. This study suggested that SARS-CoV-2 B.1.1.7 and B.1.351 variants may leverage the intrinsic plasticity of functional regions in the spike protein for mutation-induced modulation of protein dynamics and allosteric regulation to control binding with the host cell receptor.

Multiscale Modeling of Amyloid Fibrils Formed by Aggregating Peptides Derived from the Amyloidogenic Fragment of the A-Chain of Insulin

Computational prediction of molecular structures of amyloid fibrils remains an exceedingly challenging task. In this work, we propose a multi-scale modeling procedure for the structure prediction of amyloid fibrils formed by the association of ACC_1-13 aggregation-prone peptides derived from the N-terminal region of insulin’s A-chain. First, a large number of protofilament models composed of five copies of interacting ACC_1-13 peptides were predicted by application of CABS-dock coarse-grained (CG) docking simulations. Next, the models were reconstructed to all-atom (AA) representations and refined during molecular dynamics (MD) simulations in explicit solvent. The top-scored protofilament models, selected using symmetry criteria, were used for the assembly of long fibril structures. Finally, the amyloid fibril models resulting from the AA MD simulations were compared with atomic force microscopy (AFM) imaging experimental data. The obtained results indicate that the proposed multi-scale modeling procedure is capable of predicting protofilaments with high accuracy and may be applied for structure prediction and analysis of other amyloid fibrils.

Protein-Protein Docking with Large-Scale Backbone Flexibility Using Coarse-Grained Monte-Carlo Simulations

Most of the protein–protein docking methods treat proteins as almost rigid objects. Only the side-chains flexibility is usually taken into account. The few approaches enabling docking with a flexible backbone typically work in two steps, in which the search for protein–protein orientations and structure flexibility are simulated separately. In this work, we propose a new straightforward approach for docking sampling. It consists of a single simulation step during which a protein undergoes large-scale backbone rearrangements, rotations, and translations. Simultaneously, the other protein exhibits small backbone fluctuations. Such extensive sampling was possible using the CABS coarse-grained protein model and Replica Exchange Monte Carlo dynamics at a reasonable computational cost. In our proof-of-concept simulations of 62 protein–protein complexes, we obtained acceptable quality models for a significant number of cases.

Molecular Dynamics Scoring of Protein-Peptide Models Derived from Coarse-Grained Docking

One of the major challenges in the computational prediction of protein-peptide complexes is the scoring of predicted models. Usually, it is very difficult to find the most accurate solutions out of the vast number of sometimes very different and potentially plausible predictions. In this work, we tested the protocol for Molecular Dynamics (MD)-based scoring of protein-peptide complex models obtained from coarse-grained (CG) docking simulations. In the first step of the scoring procedure, all models generated by CABS-dock were reconstructed starting from their original C-alpha trace representations to all-atom (AA) structures. The second step included geometry optimization of the reconstructed complexes followed by model scoring based on receptor-ligand interaction energy estimated from short MD simulations in explicit water. We used two well-known AA MD force fields, CHARMM and AMBER, and a CG MARTINI force field. Scoring results for 66 different protein-peptide complexes show that the proposed MD-based scoring approach can be used to identify protein-peptide models of high accuracy. The results also indicate that the scoring accuracy may be significantly affected by the quality of the reconstructed protein receptor structures.

Comparative Perturbation-Based Modeling of the SARS-CoV-2 Spike Protein Binding with Host Receptor and Neutralizing Antibodies: Structurally Adaptable Allosteric Communication Hotspots Define Spike Sites Targeted by Global Circulating Mutations

In this study, we used an integrative computational approach to examine molecular mechanisms and determine functional signatures underlying the role of functional residues in the SARS-CoV-2 spike protein that are targeted by novel mutational variants and antibody-escaping mutations. Atomistic simulations and functional dynamics analysis are combined with alanine scanning and mutational sensitivity profiling of the SARS-CoV-2 spike protein complexes with the ACE2 host receptor and the REGN-COV2 antibody cocktail(REG10987+REG10933). Using alanine scanning and mutational sensitivity analysis, we have shown that K417, E484, and N501 residues correspond to key interacting centers with a significant degree of structural and energetic plasticity that allow mutants in these positions to afford the improved binding affinity with ACE2. Through perturbation-based network modeling and community analysis of the SARS-CoV-2 spike protein complexes with ACE2, we demonstrate that E406, N439, K417, and N501 residues serve as effector centers of allosteric interactions and anchor major intermolecular communities that mediate long-range communication in the complexes. The results provide support to a model according to which mutational variants and antibody-escaping mutations constrained by the requirements for host receptor binding and preservation of stability may preferentially select structurally plastic and energetically adaptable allosteric centers to differentially modulate collective motions and allosteric interactions in the complexes with the ACE2 enzyme and REGN-COV2 antibody combination. This study suggests that the SARS-CoV-2 spike protein may function as a versatile and functionally adaptable allosteric machine that exploits the plasticity of allosteric regulatory centers to fine-tune response to antibody binding without compromising the activity of the spike protein.

Integrated Biophysical Modeling of the SARS-CoV-2 Spike Protein Binding and Allosteric Interactions with Antibodies

Structural and biochemical studies of the severe acute respiratory syndrome (SARS)-CoV-2 spike glycoproteins and complexes with highly potent antibodies have revealed multiple conformation-dependent epitopes highlighting conformational plasticity of spike proteins and capacity for eliciting specific binding and broad neutralization responses. In this study, we used coevolutionary analysis, molecular simulations, and perturbation-based hierarchical network modeling of the SARS-CoV-2 spike protein complexes with a panel of antibodies targeting distinct epitopes to explore molecular mechanisms underlying binding-induced modulation of dynamics and allosteric signaling in the spike proteins. Through coevolutionary analysis of the SARS-CoV-2 spike proteins, we identified highly coevolving hotspots and functional clusters that enable a functional cross-talk between distant allosteric regions in the SARS-CoV-2 spike complexes with antibodies. Coarse-grained and all-atom molecular dynamics simulations combined with mutational sensitivity mapping and perturbation-based profiling of the SARS-CoV-2 receptor-binding domain (RBD) complexes with CR3022 and CB6 antibodies enabled a detailed validation of the proposed approach and an extensive quantitative comparison with the experimental structural and deep mutagenesis scanning data. By combining in silico mutational scanning, perturbation-based modeling, and network analysis of the SARS-CoV-2 spike trimer complexes with H014, S309, S2M11, and S2E12 antibodies, we demonstrated that antibodies can incur specific and functionally relevant changes by modulating allosteric propensities and collective dynamics of the SARS-CoV-2 spike proteins. The results provide a novel insight into regulatory mechanisms of SARS-CoV-2 S proteins showing that antibody-escaping mutations can preferentially target structurally adaptable energy hotspots and allosteric effector centers that control functional movements and allosteric communication in the complexes.

Genetic Algorithm Embedded with a Search Space Dimension Reduction Scheme for Efficient Peptide Structure Predictions

The computational determination of peptide conformations is a challenging task of finding minima in a high dimensional space. By combining the sampling efficiency of the genetic algorithm (GA) and the dimensionality reduction resulted from the backbone dihedral angle correlations, named as the path matrix (PM) method, a new searching algorithm, parallel microgenetic algorithm (PMGA), is proposed. Meanwhile, PMGA employs the density functional theory based energy as the fitness function and performs local geometry optimizations to enhance the reliability of its GA encoding strategy. Tests on peptides with up to eight amino-acid residues show PMGA is quite efficient for providing high-quality conformational coverages. The computational cost of the PMGA search increases slowly with the number of amino-acid residues in a peptide, with no sign of deterioration on the searching results for the increased length of the peptide. The PMGA method should therefore be useful for determining the conformations of oligopeptide, studying the protein-ligand interactions, and designing the peptide-based drugs.

Protocols for Fast Simulations of Protein Structure Flexibility Using CABS-Flex and SURPASS

Conformational flexibility of protein structures can play an important role in protein function. The flexibility is often studied using computational methods since experimental characterization can be difficult. Depending on protein system size, computational tools may require large computational resources or significant simplifications in the modeled systems to speed up calculations. In this work, we present the protocols for efficient simulations of flexibility of folded protein structures that use coarse-grained simulation tools of different resolutions: medium, represented by CABS-flex, and low, represented by SUPRASS. We test the protocols using a set of 140 globular proteins and compare the results with structure fluctuations observed in MD simulations, ENM modeling, and NMR ensembles. As demonstrated, CABS-flex predictions show high correlation to experimental and MD simulation data, while SURPASS is less accurate but promising in terms of future developments.

Dynamic Network Modeling of Allosteric Interactions and Communication Pathways in the SARS-CoV-2 Spike Trimer Mutants: Differential Modulation of Conformational Landscapes and Signal Transmission via Cascades of Regulatory Switches

The rapidly growing body of structural and biochemical studies of the SARS-CoV-2 spike glycoprotein has revealed a variety of distinct functional states with radically different arrangements of the receptor-binding domain, highlighting a remarkable function-driven conformational plasticity and adaptability of the spike proteins. In this study, we examined molecular mechanisms underlying conformational and dynamic changes in the SARS-CoV-2 spike mutant trimers through the lens of dynamic analysis of allosteric interaction networks and atomistic modeling of signal transmission. Using an integrated approach that combined coarse-grained molecular simulations, protein stability analysis, and perturbation-based modeling of residue interaction networks, we examined how mutations in the regulatory regions of the SARS-CoV-2 spike protein can differentially affect dynamics and allosteric signaling in distinct functional states. The results of this study revealed key functional regions and regulatory centers that govern collective dynamics, allosteric interactions, and control signal transmission in the SARS-CoV-2 spike proteins. We found that the experimentally confirmed regulatory hotspots that dictate dynamic switching between conformational states of the SARS-CoV-2 spike protein correspond to the key hinge sites and global mediating centers of the allosteric interaction networks. The results of this study provide a novel insight into allosteric regulatory mechanisms of SARS-CoV-2 spike proteins showing that mutations at the key regulatory positions can differentially modulate distribution of states and determine topography of signal communication pathways operating through state-specific cascades of control switch points. This analysis provides a plausible strategy for allosteric probing of the conformational equilibrium and therapeutic intervention by targeting specific hotspots of allosteric interactions and communications in the SARS-CoV-2 spike proteins.

Molecular Simulations and Network Modeling Reveal an Allosteric Signaling in the SARS-CoV-2 Spike Proteins

The development of computational strategies for the quantitative characterization of the functional mechanisms of SARS-CoV-2 spike proteins is of paramount importance in efforts to accelerate the discovery of novel therapeutic agents and vaccines combating the COVID-19 pandemic. Structural and biophysical studies have recently characterized the conformational landscapes of the SARS-CoV-2 spike glycoproteins in the prefusion form, revealing a spectrum of stable and more dynamic states. By employing molecular simulations and network modeling approaches, this study systematically examined functional dynamics and identified the regulatory centers of allosteric interactions for distinct functional states of the wild-type and mutant variants of the SARS-CoV-2 prefusion spike trimer. This study presents evidence that the SARS-CoV-2 spike protein can function as an allosteric regulatory engine that fluctuates between dynamically distinct functional states. Perturbation-based modeling of the interaction networks revealed a key role of the cross-talk between the effector hotspots in the receptor binding domain and the fusion peptide proximal region of the SARS-CoV-2 spike protein. The results have shown that the allosteric hotspots of the interaction networks in the SARS-CoV-2 spike protein can control the dynamic switching between functional conformational states that are associated with virus entry to the host receptor. This study offers a useful and novel perspective on the underlying mechanisms of the SARS-CoV-2 spike protein through the lens of allosteric signaling as a regulatory apparatus of virus transmission that could open up opportunities for targeted allosteric drug discovery against SARS-CoV-2 proteins and contribute to the rapid response to the current and potential future pandemic scenarios.

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

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

Impact of Early Pandemic Stage Mutations on Molecular Dynamics of SARS-CoV-2 Mpro

A new coronavirus (SARS-CoV-2) is a global threat to world health and economy. Its dimeric main protease (Mpro), which is required for the proteolytic cleavage of viral precursor proteins, is a good candidate for drug development owing to its conservation and the absence of a human homolog. Improving our understanding of Mpro behavior can accelerate the discovery of effective therapies to reduce mortality. All-atom molecular dynamics (MD) simulations (100 ns) of 50 mutant Mpro dimers obtained from filtered sequences from the GISAID database were analyzed using root-mean-square deviation, root-mean-square fluctuation, Rg, averaged betweenness centrality, and geometry calculations. The results showed that SARS-CoV-2 Mpro essentially behaves in a similar manner to its SAR-CoV homolog. However, we report the following new findings from the variants: (1) Residues GLY15, VAL157, and PRO184 have mutated more than once in SARS CoV-2; (2) the D48E variant has lead to a novel "TSEEMLN"" loop at the binding pocket; (3) inactive apo Mpro does not show signs of dissociation in 100 ns MD; (4) a non-canonical pose for PHE140 widens the substrate binding surface; (5) dual allosteric pockets coinciding with various stabilizing and functional components of the substrate binding pocket were found to display correlated compaction dynamics; (6) high betweenness centrality values for residues 17 and 128 in all Mpro samples suggest their high importance in dimer stability-one such consequence has been observed for the M17I mutation whereby one of the N-fingers was highly unstable. (7) Independent coarse-grained Monte Carlo simulations suggest a relationship between the rigidity/mutability and enzymatic function. Our entire approach combining database preparation, variant retrieval, homology modeling, dynamic residue network (DRN), relevant conformation retrieval from 1-D kernel density estimates from reaction coordinates to other existing approaches of structural analysis, and data visualization within the coronaviral Mpro is also novel and is applicable to other coronaviral proteins.

Exploring Mechanisms of Communication Switching in the Hsp90-Cdc37 Regulatory Complexes with Client Kinases through Allosteric Coupling of Phosphorylation Sites: Perturbation-Based Modeling and Hierarchical Community Analysis of Residue Interaction Networks

Understanding molecular principles underlying chaperone-based modulation of kinase client activity is critically important to dissect functions and activation mechanisms of many oncogenic proteins. The recent experimental studies have suggested that phosphorylation sites in the Hsp90 and Cdc37 proteins can serve as conformational communication switches of chaperone regulation and kinase interactions. However, a mechanism of allosteric coupling between phosphorylation sites in the Hsp90 and Cdc37 during client binding is poorly understood, and the molecular signatures underpinning specific roles of phosphorylation sites in the Hsp90 regulation remain unknown. In this work, we employed a combination of evolutionary analysis, coarse-grained molecular simulations together with perturbation-based network modeling and scanning of the unbound and bound Hsp90 and Cdc37 structures to quantify allosteric effects of phosphorylation sites and identify unique signatures that are characteristic for communication switches of kinase-specific client binding. By using network-based metrics of the dynamic intercommunity bridgeness and community centrality, we characterize specific signatures of phosphorylation switches involved in allosteric regulation. Through perturbation-based analysis of the dynamic residue interaction networks, we show that mutations of kinase-specific phosphorylation switches can induce long-range effects and lead to a global rewiring of the allosteric network and signal transmission in the Hsp90-Cdc37-kinase complex. We determine a specific group of phosphorylation sites in the Hsp90 where mutations may have a strong detrimental effect on allosteric interaction network, providing insight into the mechanism of phosphorylation-induced communication switching. The results demonstrate that kinase-specific phosphorylation switches of communications in the Hsp90 may be partly predisposed for their regulatory role based on preexisting allosteric propensities.

Docking of peptides to GPCRs using a combination of CABS-dock with FlexPepDock refinement

The structural description of peptide ligands bound to G protein-coupled receptors (GPCRs) is important for the discovery of new drugs and deeper understanding of the molecular mechanisms of life. Here we describe a three-stage protocol for the molecular docking of peptides to GPCRs using a set of different programs: (1) CABS-dock for docking fully flexible peptides; (2) PD2 method for the reconstruction of atomistic structures from C-alpha traces provided by CABS-dock and (3) Rosetta FlexPepDock for the refinement of protein–peptide complex structures and model scoring. We evaluated the proposed protocol on the set of seven different GPCR–peptide complexes (including one containing a cyclic peptide), for which crystallographic structures are available. We show that CABS-dock produces high resolution models in the sets of top-scored models. These sets of models, after reconstruction to all-atom representation, can be further improved by Rosetta high-resolution refinement and/or minimization, leading in most of the cases to sub-Angstrom accuracy in terms of interface root-mean-square-deviation measure.

Prediction of Protein Tertiary Structure via Regularized Template Classification Techniques

We discuss the use of the regularized linear discriminant analysis (LDA) as a model reduction technique combined with particle swarm optimization (PSO) in protein tertiary structure prediction, followed by structure refinement based on singular value decomposition (SVD) and PSO. The algorithm presented in this paper corresponds to the category of template-based modeling. The algorithm performs a preselection of protein templates before constructing a lower dimensional subspace via a regularized LDA. The protein coordinates in the reduced spaced are sampled using a highly explorative optimization algorithm, regressive–regressive PSO (RR-PSO). The obtained structure is then projected onto a reduced space via singular value decomposition and further optimized via RR-PSO to carry out a structure refinement. The final structures are similar to those predicted by best structure prediction tools, such as Rossetta and Zhang servers. The main advantage of our methodology is that alleviates the ill-posed character of protein structure prediction problems related to high dimensional optimization. It is also capable of sampling a wide range of conformational space due to the application of a regularized linear discriminant analysis, which allows us to expand the differences over a reduced basis set.

Computational Investigation on the p53-MDM2 Interaction Using the Potential of Mean Force Study

Murine double minute 2 (MDM2) proteins are found to be overproduced by many human tumors in order to inhibit the functioning of p53 molecules, a tumor suppressor protein. Thus, reactivating p53 functioning in cancer cells by disrupting p53-MDM2 interactions may offer a significant approach in cancer treatment. However, the structural characterization of the p53-MDM2 complex at the atomistic level and the mechanism of binding/unbinding of the p53-MDM2 complex still remain unclear. Therefore, we demonstrate here the probable binding (unbinding) pathway of transactivation domain 1 of p53 during the formation (dissociation) of the p53-MDM2 complex in terms of free energy as a function of reaction coordinate from the potential of mean force (PMF) study using two different force fields: ff99SB and ff99SB-ILDN. From the PMF plot, we noticed the PMF to have a minimum value at a p53-MDM2 separation of 12 Å, with a dissociation energy of 30 kcal mol^-1. We also analyzed the conformational dynamics and stability of p53 as a function of its distance of separation from MDM2. The secondary structure content (helix and turns) in p53 was found to vary with its distance of separation from MDM2. The p53-MDM2 complex structure with lowest potential energy was isolated from the ensemble at the reaction coordinate corresponding to the minimum PMF value and subjected to molecular dynamics simulation to identify the interface surface area, interacting residues at the interface, and the stability of the complex. The simulation results highlight the importance of hydrogen bonds and the salt bridge between Lys94 of MDM2 and Glu17 of p53 in the stability of the p53-MDM2 complex. We also carried out the binding free energy calculations and the per residue energy decomposition analyses of the interface residues of the p53-MDM2 complex. We found that the binding affinity between MDM2 and p53 is indeed high [ΔG _bind = -7.29 kcal mol^-1 from molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) and ΔG _bind = -53.29 kcal mol^-1 from molecular mechanics/generalized borne surface area]. The total binding energy obtained using the MM/PBSA method was noticed to be closer to the experimental values (-6.4 to -9.0 kcal mol^-1). The p53-MDM2 complex binding profile was observed to follow the same trend even in the duplicate simulation run and also in the simulation carried out with different force fields. We found that Lys51, Leu54, Tyr100, and Tyr104 from MDM2 and the residues Phe19, Trp23, and Leu26 from p53 provide the highest energy contributions for the p53-MDM2 interaction. Our findings highlight the prominent structural and binding characteristics of the p53-MDM2 complex that may be useful in designing potential inhibitors to disrupt the p53-MDM2 interactions.

Integrated Computational Approaches and Tools forAllosteric Drug Discovery

Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.

Flexible docking of peptides to proteins using CABS-dock

Molecular docking of peptides to proteins can be a useful tool in the exploration of the possible peptide binding sites and poses. CABS-dock is a method for protein-peptide docking that features significant conformational flexibility of both the peptide and the protein molecules during the peptide search for a binding site. The CABS-dock has been made available as a web server and a standalone package. The web server is an easy to use tool with a simple web interface. The standalone package is a command-line program dedicated to professional users. It offers a number of advanced features, analysis tools and support for large-sized systems. In this article, we outline the current status of the CABS-dock method, its recent developments, applications, and challenges ahead.

Using Machine Learning in Accuracy Assessment of Knowledge-Based Energy and Frequency Base Likelihood in Protein Structures

Many aspects of the study of protein folding and dynamics have been affected by the accumulation of data about native protein structures and recent advances in machine learning. Computational methods for predicting protein structures from their sequences are now heavily based on machine learning tools and on approaches that extract knowledge and rules from data using probabilistic models. Many of these methods use scoring functions to determine which structure best fits a native protein sequence. Using computational approaches, we obtained two scoring functions: knowledge-based energy and likelihood of base frequency, and we compared their accuracy in measuring the sequence structure fit. We compared the machine learning models’ accuracy of predictions for knowledge-based energy and likelihood values to validate our results, showing that likelihood is a more accurate scoring function than knowledge-based energy.

Computational reconstruction of atomistic protein structures from coarse-grained models

Three-dimensional protein structures, whether determined experimentally or theoretically, are often too low resolution. In this mini-review, we outline the computational methods for protein structure reconstruction from incomplete coarse-grained to all atomistic models. Typical reconstruction schemes can be divided into four major steps. Usually, the first step is reconstruction of the protein backbone chain starting from the C-alpha trace. This is followed by side-chains rebuilding based on protein backbone geometry. Subsequently, hydrogen atoms can be reconstructed. Finally, the resulting all-atom models may require structure optimization. Many methods are available to perform each of these tasks. We discuss the available tools and their potential applications in integrative modeling pipelines that can transfer coarse-grained information from computational predictions, or experiment, to all atomistic structures.

Establishing Computational Approaches Towards Identifying Malarial Allosteric Modulators: A Case Study of Plasmodium falciparum Hsp70s

Combating malaria is almost a never-ending battle, as Plasmodium parasites develop resistance to the drugs used against them, as observed recently in artemisinin-based combination therapies. The main concern now is if the resistant parasite strains spread from Southeast Asia to Africa, the continent hosting most malaria cases. To prevent catastrophic results, we need to find non-conventional approaches. Allosteric drug targeting sites and modulators might be a new hope for malarial treatments. Heat shock proteins (HSPs) are potential malarial drug targets and have complex allosteric control mechanisms. Yet, studies on designing allosteric modulators against them are limited. Here, we identified allosteric modulators (SANC190 and SANC651) against P. falciparum Hsp70-1 and Hsp70-x, affecting the conformational dynamics of the proteins, delicately balanced by the endogenous ligands. Previously, we established a pipeline to identify allosteric sites and modulators. This study also further investigated alternative approaches to speed up the process by comparing all atom molecular dynamics simulations and dynamic residue network analysis with the coarse-grained (CG) versions of the calculations. Betweenness centrality (BC) profiles for PfHsp70-1 and PfHsp70-x derived from CG simulations not only revealed similar trends but also pointed to the same functional regions and specific residues corresponding to BC profile peaks.

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence-structure-function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well-defined three-dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.