LoopGrafter: a web tool for transplanting dynamical loops for protein engineering.

The transplantation of loops between structurally related proteins is a compelling method to improve the activity, specificity and stability of enzymes. However, despite the interest of loop regions in protein engineering, the available methods of loop-based rational protein design are scarce. One particular difficulty related to loop engineering is the unique dynamism that enables them to exert allosteric control over the catalytic function of enzymes. Thus, when engaging in a transplantation effort, such dynamics in the context of protein structure need consideration. A second practical challenge is identifying successful excision points for the transplantation or grafting. Here, we present LoopGrafter (https://loschmidt.chemi.muni.cz/loopgrafter/), a web server that specifically guides in the loop grafting process between structurally related proteins. The server provides a step-by-step interactive procedure in which the user can successively identify loops in the two input proteins, calculate their geometries, assess their similarities and dynamics, and select a number of loops to be transplanted. All possible different chimeric proteins derived from any existing recombination point are calculated, and 3D models for each of them are constructed and energetically evaluated. The obtained results can be interactively visualized in a user-friendly graphical interface and downloaded for detailed structural analyses.

Insights into the Dynamics and Binding of Two Polyprotein Substrate Cleavage Points in the Context of the SARS-CoV-2 Main and Papain-like Proteases.

It is well known that vital enzymes in the replication process of the coronavirus are the SARS-CoV-2 PLpro and SARS-CoV-2 3CLpro, both of which are important targets in the search for anti-coronavirus agents. These two enzymes are responsible for cleavage at various polyprotein sites in the SARS-CoV-2 lifecycle. Herein, the dynamics of the polyprotein cleavage sequences for the boundary between non-structural proteins Nsp1 and Nsp2 (CS1) and between Nsp2 and Nsp3 (CS2) in complex with both the papain-like protein PLpro and the main protease 3CLpro were explored using computational methods. The post dynamics analysis reveals that CS1 and CS2 both have greater stability when complexed with PLpro. Of these two, greater stability is observed for the CS1-PLpro complex, while destabilization resulting in loss of CS2 from the PLpro active site is observed for CS2-PLpro, suggesting the rate of exchange by the papain-like protease is faster for CS2 compared to CS1. On the other hand, the 3CLpro main protease also reveals stability for CS1 suggesting that the main protease could also play a potential role in the cleavage at point CS1. However, destabilization occurs early in the simulation for the complex CLpro-CS2 suggesting a poor interaction and non-plausible protease cleavage of the polyprotein at CS2 by the main protease. These findings could be used as a guide in the development and design of potent COVID-19 antiviral inhibitors that mimic the CS1 cleavage site.

Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: insight through molecular dynamics and hybrid QM/MM studies.

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease were studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme-substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.

Mental health research in the lower-middle-income countries of Africa and Asia during the COVID-19 pandemic: A scoping review.

BACKGROUND: Lower middle - income countries of Africa and Asia have accommodated a large portion of the world's population, where mental health research has been under-prioritized. This study aimed to review all published research on mental health issues related to the COVID-19 pandemic in lower middle - income countries of the Afro-Asian region. METHODS: A search was conducted in the PubMed and PubMed Central databases using the terms "mental health," "psychiatric disorders," "COVID-19", "coronavirus," "Asia," "Africa," "Lower Middle-income countries." All articles published in the English language until 3rd July 2020 were included. RESULT: A total of 133 papers were found in lower-middle-income countries of Africa (n = 11) and Asia (n = 122), which discussed various aspects of mental health in the context of COVID-19. Most of the studies are cross-sectional studies that addressed mental morbidities, psychological reactions, stress, coping among the general population, and focused groups (healthcare workers, students, elderly, patients). Researchers attempted to develop and validate tools that measure certain psychological constructs (fear, anxiety) concerning COVID-19. Online surveys were the primary modality of researching this COVID-19 pandemic. CONCLUSION: COVID-19 related mental health research is scarce in lower-middle-income countries of Africa and Asia. Available researches suggest that mental health challenges during this COVID-19 pandemic are enormous and need attention. There is a need for policy and recommendations to deal with the mental health challenges in lower middle - income countries of Africa and Asia.

An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant.

The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host-guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.

CoVAMPnet: Comparative Markov State Analysis for Studying Effects of Drug Candidates on Disordered Biomolecules.

Computational study of the effect of drug candidates on intrinsically disordered biomolecules is challenging due to their vast and complex conformational space. Here, we developed a comparative Markov state analysis (CoVAMPnet) framework to quantify changes in the conformational distribution and dynamics of a disordered biomolecule in the presence and absence of small organic drug candidate molecules. First, molecular dynamics trajectories are generated using enhanced sampling, in the presence and absence of small molecule drug candidates, and ensembles of soft Markov state models (MSMs) are learned for each system using unsupervised machine learning. Second, these ensembles of learned MSMs are aligned across different systems based on a solution to an optimal transport problem. Third, the directional importance of inter-residue distances for the assignment to different conformational states is assessed by a discriminative analysis of aggregated neural network gradients. This final step provides interpretability and biophysical context to the learned MSMs. We applied this novel computational framework to assess the effects of ongoing phase 3 therapeutics tramiprosate (TMP) and its metabolite 3-sulfopropanoic acid (SPA) on the disordered Aß42 peptide involved in Alzheimer's disease. Based on adaptive sampling molecular dynamics and CoVAMPnet analysis, we observed that both TMP and SPA preserved more structured conformations of Aß42 by interacting nonspecifically with charged residues. SPA impacted Aß42 more than TMP, protecting α-helices and suppressing the formation of aggregation-prone ß-strands. Experimental biophysical analyses showed only mild effects of TMP/SPA on Aß42 and activity enhancement by the endogenous metabolization of TMP into SPA. Our data suggest that TMP/SPA may also target biomolecules other than Aß peptides. The CoVAMPnet method is broadly applicable to study the effects of drug candidates on the conformational behavior of intrinsically disordered biomolecules.

Exploring the concerted mechanistic pathway for HIV-1 PR-substrate revealed by umbrella sampling simulation.

HIV-1 protease (HIV-1 PR) is an essential enzyme for the replication process of its virus, and therefore considered an important target for the development of drugs against the acquired immunodeficiency syndrome (AIDS). Our previous study shows that the catalytic mechanism of subtype B/C-SA HIV-1 PR follows a one-step concerted acyclic hydrolysis reaction process using a two-layered ONIOM B3LYP/6-31++G(d,p) method. This present work is aimed at exploring the proposed mechanism of the proteolysis catalyzed by HIV-1 PR and to ensure our proposed mechanism is not an artefact of a single theoretical technique. Hence, we present umbrella sampling method that is suitable for calculating potential mean force (PMF) for non-covalent ligand/substrate-enzyme association/dissociation interactions which provide thermodynamic details for molecular recognition. The free activation energy results were computed in terms of PMF analysis within the hybrid QM(DFTB)/MM approach. The theoretical findings suggest that the proposed mechanism corresponds in principle with experimental data. Given our observations, we suggest that the QM/MM MD method can be used as a reliable computational technique to rationalize lead compounds against specific targets such as the HIV-1 protease.

Computational modelling of potent ß-secretase (BACE1) inhibitors towards Alzheimer's disease treatment.

Researchers have identified the ß-amyloid precursor protein cleaving enzyme 1 (BACE1) in the multifactorial pathway of Alzheimer's disease (AD) as a drug target. The design and development of molecules to inhibit BACE1 as a potential cure for AD thus remained significant. Herein, we simulated two potent BACE1 inhibitors (AM-6494 and CNP-520) to understand their binding affinity at the atomistic level. AM-6494 is a newly reported potent BACE1 inhibitor with an IC50 value of 0.4 nM in vivo and now picked for preclinical considerations. Umibecestat (CNP-520), which was discontinued at human trials lately, was considered to enable a reasonable evaluation of our results. Using density functional theory (DFT) and Our Own N-layered Integrated molecular Orbital and Molecular Mechanics (ONIOM), we achieved the aim of this investigation. These computational approaches enabled the prediction of the electronic properties of AM-6494 and CNP-520 plus their binding energies when complexed with BACE1. For AM-6494 and CNP-520 interaction with protonated BACE1, the ONIOM calculation gave binding free energy of -62.849 and -33.463 kcal/mol, respectively. In the unprotonated model, we observed binding free energy of -59.758 kcal/mol in AM-6494. Taken together thermochemistry of the process and molecular interaction plot, AM-6494 is more favourable than CNP-520 towards the inhibition of BACE1. The protonated model gave slightly better binding energy than the unprotonated form. However, both models could sufficiently describe ligand binding to BACE1 at the atomistic level. Understanding the detailed molecular interaction of these inhibitors could serve as a basis for pharmacophore exploration towards improved inhibitor design.

From Recognition to Reaction Mechanism: An Overview on the Interactions between HIV-1 Protease and its Natural Targets.

Current investigations on the Human Immunodeficiency Virus Protease (HIV-1 PR) as a druggable target towards the treatment of AIDS require an update to facilitate further development of promising inhibitors with improved inhibitory activities. For the past two decades, up to 100 scholarly reports appeared annually on the inhibition and catalytic mechanism of HIV-1 PR. A fundamental literature review on the prerequisite of HIV-1 PR action leading to the release of the infectious virion is absent. Herein, recent advances (both computationally and experimentally) on the recognition mode and reaction mechanism of HIV-1 PR involving its natural targets are provided. This review features more than 80 articles from reputable journals. Recognition of the natural Gag and Gag-Pol cleavage junctions by this enzyme and its mutant analogs was first addressed. Thereafter, a comprehensive dissect of the enzymatic mechanism of HIV-1 PR on its natural polypeptide sequences from literature was put together. In addition, we highlighted ongoing research topics in which in silico methods could be harnessed to provide deeper insights into the catalytic mechanism of the HIV-1 protease in the presence of its natural substrates at the molecular level. Understanding the recognition and catalytic mechanism of HIV-1 PR leading to the release of an infective virion, which advertently affects the immune system, will assist in designing mechanismbased inhibitors with improved bioactivity.

Theoretical Model for HIV-1 PR That Accounts for Substrate Recognition and Preferential Cleavage of Natural Substrates.

The Human Immunodeficiency Virus type 1 (HIV-1) protease is a crucial target for HIV/AIDS treatment, and understanding its catalytic mechanism is the basis on which HIV-1 enzyme inhibitors are developed. Several experimental studies have indicated that HIV-1 protease facilitates the cleavage of the Gag and Gag-Pol polyproteins and it is highly selective with regard to the cleaved amino acid precursors and physical parameters. However, the main theoretical principles of substrate specificity and recognition remain poorly understood theoretically. By means of a one-step concerted transition state modeling, the recognition of natural substrates by HIV-1 PR subtypes (B and C-SA) was studied. This was carried out to compare the activation free energies at varying peptide bond regions (scissile and nonscissile) within the polypeptide sequence using ONIOM calculations. We studied both P3-P3' and P5-P5' natural substrate systems. For P3-P3' substrates, excellent recognition was observed for the MA-CA family but not for the RH-IN substrates. Satisfactory recognition for the latter was only observed for the longer sequence (P5-P5') after the substrate was subjected to an MD run to maximize the interaction between the enzyme and the substrate. These results indicate that both sequence and structure are important for correct scissile bond recognition of these natural substrates.