Elucidation of cryptic and allosteric pockets within the SARS-CoV-2 protease

Structure and dynamics of pteridine reductase 1: the key phenomena relevant to enzyme function and drug design

Pteridine reductase 1 (PTR1) is a folate and pterin pathway enzyme unique for pathogenic trypanosomatids. As a validated drug target, PTR1 has been the focus of recent research efforts aimed at finding more effective treatments against human parasitic diseases such as leishmaniasis or sleeping sickness. Previous PTR1-centered structural studies highlighted the enzyme characteristics, such as flexible regions around the active site, highly conserved structural waters, and species-specific differences in pocket properties and dynamics, which likely impacts the binding of natural substrates and inhibitors. Furthermore, several aspects of the PTR1 function, such as the substrate inhibition phenomenon and the level of ligand binding cooperativity in the enzyme homotetramer, likely related to the global enzyme dynamics, are poorly known at the molecular level. We postulate that future drug design efforts could greatly benefit from a better understanding of these phenomena through studying both the local and global PTR1 dynamics. This review highlights the key aspects of the PTR1 structure and dynamics relevant to structure-based drug design that could be effectively investigated by modeling approaches. Particular emphasis is given to the perspective of molecular dynamics, what has been accomplished in this area to date, and how modeling could impact the PTR1-targeted drug design in the future.

Molecular docking and molecular dynamics study Lianhua Qingwen granules (LHQW) treats COVID-19 by inhibiting inflammatory response and regulating cell survival

Purpose

2019 Coronavirus disease (COVID-19) is endangering health of populations worldwide. Latest research has proved that Lianhua Qingwen granules (LHQW) can reduce tissue damage caused by inflammatory reactions and relieve patients' clinical symptoms. However, the mechanism of LHQW treats COVID-19 is currently lacking. Therefore, we employed computer simulations to investigate the mechanism of LHQW treats COVID-19 by modulating inflammatory response.

Methods

We employed bioinformatics to screen active ingredients in LHQW and intersection gene targets. PPI, GO and KEGG was used to analyze relationship of intersection gene targets. Molecular dynamics simulations validated the binding stability of active ingredients and target proteins. Binding free energy, radius of gyration and the solvent accessible surface area were analyzed by supercomputer platform.

Results

COVID-19 had 4628 gene targets, LHQW had 1409 gene targets, intersection gene targets were 415. Bioinformatics analysis showed that intersection targets were closely related to inflammation and immunomodulatory. Molecular docking suggested that active ingredients (including: licopyranocoumarin, Glycyrol and 3-3-Oxopropanoic acid) in LHQW played a role in treating COVID-19 by acting on CSF2, CXCL8, CCR5, NLRP3, IFNG and TNF. Molecular dynamics was used to prove the binding stability of active ingredients and protein targets.

Conclusion

The mechanism of active ingredients in LHQW treats COVID-19 was investigated by computer simulations. We found that active ingredients in LHQW not only reduce cell damage and tissue destruction by inhibiting the inflammatory response through CSF2, CXCL8, CCR5 and IFNG, but also regulate cell survival and growth through NLRP3 and TNF thereby reducing apoptosis.

CROMATIC: Cross-Relationship Map of Cavities from Coronaviruses

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiologic agent of COVID-19 disease, has rapidly imposed an urgent need to identify effective drug candidates. In this context, the high resolution and non-redundant beta-Coronavirus protein cavities database is pivotal to help virtual screening protocols. Furthermore, the cross-relationship among cavities can lead to highlighting multitarget therapy chances. Here, we first collect all protein cavities on SARS-CoV-2, SARS-CoV, and MERS-CoV X-ray structures, and then, we compute a similarity map by using molecular interaction fields (MIFs). All the results come together in CROMATIC (CROss-relationship MAp of CaviTIes from Coronaviruses). CROMATIC encloses both a comprehensive and a non-redundant version of the cavities collection and a similarity map revealing, on the one hand, cavities that are conserved among the three Coronaviruses and, on the other hand, unexpected similarities among cavities that can represent a key starting point for multitarget therapy strategies. Similarity analysis was also performed for the available structures of SARS-CoV-2 spike variants, linking sequence mutations to three-dimensional interaction alterations. The CROMATIC repository is freely available to the scientific community at https://github.com/moldiscovery/sars-cromatic.

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Dynamic Profiling of β-Coronavirus 3CL Mpro Protease Ligand-Binding Sites

β-coronavirus (CoVs) alone has been responsible for three major global outbreaks in the 21st century. The current crisis has led to an urgent requirement to develop therapeutics. Even though a number of vaccines are available, alternative strategies targeting essential viral components are required as a backup against the emergence of lethal viral variants. One such target is the main protease (Mpro) that plays an indispensable role in viral replication. The availability of over 270 Mpro X-ray structures in complex with inhibitors provides unique insights into ligand-protein interactions. Herein, we provide a comprehensive comparison of all nonredundant ligand-binding sites available for SARS-CoV2, SARS-CoV, and MERS-CoV Mpro. Extensive adaptive sampling has been used to investigate structural conservation of ligand-binding sites using Markov state models (MSMs) and compare conformational dynamics employing convolutional variational auto-encoder-based deep learning. Our results indicate that not all ligand-binding sites are dynamically conserved despite high sequence and structural conservation across β-CoV homologs. This highlights the complexity in targeting all three Mpro enzymes with a single pan inhibitor.

Finding Druggable Sites in Proteins Using TACTICS

Structure-based drug discovery efforts require knowledge of where drug-binding sites are located on target proteins. To address the challenge of finding druggable sites, we developed a machine-learning algorithm called TACTICS (trajectory-based analysis of conformations to identify cryptic sites), which uses an ensemble of molecular structures (such as molecular dynamics simulation data) as input. First, TACTICS uses k-means clustering to select a small number of conformations that represent the overall conformational heterogeneity of the data. Then, TACTICS uses a random forest model to identify potentially bindable residues in each selected conformation, based on protein motion and geometry. Lastly, residues in possible binding pockets are scored using fragment docking. As proof-of-principle, TACTICS was applied to the analysis of simulations of the SARS-CoV-2 main protease and methyltransferase and the Yersinia pestis aryl carrier protein. Our approach recapitulates known small-molecule binding sites and predicts the locations of sites not previously observed in experimentally determined structures. The TACTICS code is available at https://github.com/Albert-Lau-Lab/tactics_protein_analysis.

Precision therapeutic targets for COVID-19

Beginning in late 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged as a novel pathogen that causes coronavirus disease 2019 (COVID-19). SARS-CoV-2 has infected more than 111 million people worldwide and caused over 2.47 million deaths. Individuals infected with SARS-CoV-2 show symptoms of fever, cough, dyspnea, and fatigue with severe cases that can develop into pneumonia, myocarditis, acute respiratory distress syndrome, hypercoagulability, and even multi-organ failure. Current clinical management consists largely of supportive care as commonly administered treatments, including convalescent plasma, remdesivir, and high-dose glucocorticoids. These have demonstrated modest benefits in a small subset of hospitalized patients, with only dexamethasone showing demonstrable efficacy in reducing mortality and length of hospitalization. At this time, no SARS-CoV-2-specific antiviral drugs are available, although several vaccines have been approved for use in recent months. In this review, we will evaluate the efficacy of preclinical and clinical drugs that precisely target three different, essential steps of the SARS-CoV-2 replication cycle: the spike protein during entry, main protease (MPro) during proteolytic activation, and RNA-dependent RNA polymerase (RdRp) during transcription. We will assess the advantages and limitations of drugs that precisely target evolutionarily well-conserved domains, which are less likely to mutate, and therefore less likely to escape the effects of these drugs. We propose that a multi-drug cocktail targeting precise proteins, critical to the viral replication cycle, such as spike protein, MPro, and RdRp, will be the most effective strategy of inhibiting SARS-CoV-2 replication and limiting its spread in the general population.

High-resolution mining of the SARS-CoV-2 main protease conformational space: supercomputer-driven unsupervised adaptive sampling

We provide an unsupervised adaptive sampling strategy capable of producing μs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs). The global exploration problem is decomposed into a set of separate MD trajectories that can be restarted within a selective process to achieve sufficient phase-space sampling. Accurate statistical properties can be obtained through reweighting. Within this highly parallel setup, the Tinker-HP package can be powered by an arbitrary large number of GPUs on supercomputers, reducing exploration time from years to days. This approach is used to tackle the urgent modeling problem of the SARS-CoV-2 Main Protease (Mpro) producing more than 38 μs of all-atom simulations of its apo (ligand-free) dimer using the high-resolution AMOEBA PFF. The first 15.14 μs simulation (physiological pH) is compared to available non-PFF long-timescale simulation data. A detailed clustering analysis exhibits striking differences between FFs, with AMOEBA showing a richer conformational space. Focusing on key structural markers related to the oxyanion hole stability, we observe an asymmetry between protomers. One of them appears less structured resembling the experimentally inactive monomer for which a 6 μs simulation was performed as a basis for comparison. Results highlight the plasticity of the Mpro active site. The C-terminal end of its less structured protomer is shown to oscillate between several states, being able to interact with the other protomer, potentially modulating its activity. Active and distal site volumes are found to be larger in the most active protomer within our AMOEBA simulations compared to non-PFFs as additional cryptic pockets are uncovered. A second 17 μs AMOEBA simulation is performed with protonated His172 residues mimicking lower pH. Data show the protonation impact on the destructuring of the oxyanion loop. We finally analyze the solvation patterns around key histidine residues. The confined AMOEBA polarizable water molecules are able to explore a wide range of dipole moments, going beyond bulk values, leading to a water molecule count consistent with experimental data. Results suggest that the use of PFFs could be critical in drug discovery to accurately model the complexity of the molecular interactions structuring Mpro.

Candidate Binding Sites for Allosteric Inhibition of the SARS-CoV-2 Main Protease from the Analysis of Large-Scale Molecular Dynamics Simulations

We analyzed a 100 μs MD trajectory of the SARS-CoV-2 main protease by a non-parametric data analysis approach which allows characterizing a free energy landscape as a simultaneous function of hundreds of variables. We identified several conformations that, when visited by the dynamics, are stable for several hundred nanoseconds. We explicitly characterize and describe these metastable states. In some of these configurations, the catalytic dyad is less accessible. Stabilizing them by a suitable binder could lead to an inhibition of the enzymatic activity. In our analysis we keep track of relevant contacts between residues which are selectively broken or formed in the states. Some of these contacts are formed by residues which are far from the catalytic dyad and are accessible to the solvent. Based on this analysis we propose some relevant contact patterns and three possible binding sites which could be targeted to achieve allosteric inhibition.

A Multiscale Coarse-grained Model of the SARS-CoV-2 Virion

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations, however, are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and on-going development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Structural data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data becomes publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.

Significance Statement

This study reports the construction of a molecular model for the SARS-CoV-2 virion and details our multiscale approach towards model refinement. The resulting model and methods can be applied to and enable the simulation of SARS-CoV-2 virions.