Reaction mechanism of the dengue virus serine protease: a QM/MM study

Understanding the Enzyme (S)-Norcoclaurine Synthase Promiscuity to Aldehydes and Ketones

The (S)-norcoclaurine synthase from Thalictrum flavum (TfNCS) stereoselectively catalyzes the Pictet-Spengler reaction between dopamine and 4-hydroxyphenylacetaldehyde to give (S)-norcoclaurine. TfNCS can catalyze the Pictet-Spengler reaction with various aldehydes and ketones, leading to diverse tetrahydroisoquinolines. This substrate promiscuity positions TfNCS as a highly promising enzyme for synthesizing fine chemicals. Understanding carbonyl-containing substrates' structural and electronic signatures that influence TfNCS activity can help expand its applications in the synthesis of different compounds and aid in protein optimization strategies. In this study, we investigated the influence of the molecular properties of aldehydes and ketones on their reactivity in the TfNCS-catalyzed Pictet-Spengler reaction. Initially, we compiled a library of reactive and unreactive compounds from previous publications. We also performed enzymatic assays using nuclear magnetic resonance to identify some reactive and unreactive carbonyl compounds, which were then included in the library. Subsequently, we employed QSAR and DFT calculations to establish correlations between substrate-candidate structures and reactivity. Our findings highlight correlations of structural and stereoelectronic features, including the electrophilicity of the carbonyl group, to the reactivity of aldehydes and ketones toward the TfNCS-catalyzed Pictet-Spengler reaction. Interestingly, experimental data of seven compounds out of fifty-three did not correlate with the electrophilicity of the carbonyl group. For these seven compounds, we identified unfavorable interactions between them and the TfNCS. Our results demonstrate the applications of in silico techniques in understanding enzyme promiscuity and specificity, with a particular emphasis on machine learning methodologies, DFT electronic structure calculations, and molecular dynamic (MD) simulations.

A Protocol to Depict the Proteolytic Processes Using a Combination of Metal-Organic Materials (MOMs), Electron Paramagnetic Resonance (EPR), and Mass Spectrometry (MS)

Proteolysis is a critical biochemical process yet a challenging field to study experimentally due to the self-degradation of a protease and the complex, dynamic degradation steps of a substrate. Mass spectrometry (MS) is the traditional way for proteolytic studies, yet it is challenging when time-resolved, step-by-step details of the degradation process are needed. We recently found a way to resolve the cleavage site, preference/selectivity of cleavage regions, and proteolytic kinetics by combining site-directed spin labeling (SDSL) of protein substrate, time-resolved two-dimensional (2D) electron paramagnetic resonance (EPR) spectroscopy, protease immobilization via metal-organic materials (MOMs), and MS. The method has been demonstrated on a model substrate and protease, yet there is a lack of details on the practical operations to carry out our strategy. Thus, this protocol summarizes the key steps and considerations when carrying out the EPR/MS study on proteolytic processes, which can be generalized to study other protein/polypeptide substrates in proteolysis. Details for the experimental operation and cautions of each step are reported with figures illustrating the concepts. This protocol provides an effective approach to understanding the proteolytic process with the advantages of offering time-resolved, residue-level resolution of structural basis underlying the process. Such information is important for revealing the cleavage site and proteolytic mechanisms of unknown proteases. The advantage of EPR, probing the target substrate regardless of the complexities caused by the proteases and their self-degradation, offers a practically effective, rapid, and easy-to-operate approach to studying proteolysis. Key features • Combining protease immobilization, EPR, spin labeling, and MS experimental methods allows for the analysis of proteolysis process in real time. • Reveals cleavage site, kinetics of product generation, and preference of cleavage regions via time-resolved SDSL-EPR. • MS confirms EPR findings and helps depict the sequences and populations of the cleaved segments in real time. • The demonstrated method can be generalized to other proteins or polypeptide substrates upon proteolysis by other proteases.

Transmembrane Protease Serine 2 Proteolytic Cleavage of the SARS-CoV-2 Spike Protein: A Mechanistic Quantum Mechanics/Molecular Mechanics Study to Inspire the Design of New Drugs To Fight the COVID-19 Pandemic

Despite the development of vaccines against the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, there is an urgent need for efficient drugs to treat infected patients. An attractive drug target is the human transmembrane protease serine 2 (TMPRSS2) because of its vital role in the viral infection mechanism of SARS-CoV-2 by activation of the virus spike protein (S protein). Having in mind that the information derived from quantum mechanics/molecular mechanics (QM/MM) studies could be an important tool in the design of transition-state (TS) analogue inhibitors, we resorted to adiabatic QM/MM calculations to determine the mechanism of the first step (acylation) of proteolytic cleavage of the S protein with atomistic details. Acylation occurred in two stages: (i) proton transfer from Ser441 to His296 concerted with the nucleophilic attack of Ser441 to the substrate's P1-Arg and (ii) proton transfer from His296 to the P1'-Ser residue concerted with the cleavage of the ArgP1-SerP1' peptide bond, with a Gibbs activation energy of 17.1 and 15.8 kcal mol^-1, relative to the reactant. An oxyanion hole composed of two hydrogen bonds stabilized the rate-limiting TS by 8 kcal mol^-1. An analysis of the TMPRSS2 interactions with the high-energy, short-lived tetrahedral intermediate highlighted the limitations of current clinical inhibitors and pointed out specific ways to develop higher-affinity TS analogue inhibitors. The results support the development of more efficient drugs against SARS-CoV-2 using a human target, free from resistance development.

The Q41R mutation in the HCV-protease enhances the reactivity towards MAVS by suppressing non-reactive pathways

Recent experimental findings pointed out a new mutation in the HCV protease, Q41R, responsible for a significant enhancement of the enzyme's reactivity towards the mitochondrial antiviral-signaling protein (MAVS). The Q41R mutation is located rather far from the active site, and its involvement in the overall reaction mechanism is thus unclear. We used classical molecular dynamics and QM/MM to study the acylation reaction of HCV NS3/4A protease variants bound to MAVS and the NS4A/4B substrate and uncovered the indirect mechanism by which the Q41R mutation plays a critical role in the efficient cleavage of the substrate. Our simulations reveal that there are two major conformations of the MAVS H1'(p) residue for the wild type protease and only one conformation for the Q41R mutant. The conformational space of H1'(p) is restricted by the Q41R mutation due to a π-π stacking between H1'(p) and R41 as well as a strong hydrogen bond between the backbone of H57 and the side chain of R41. Further QM/MM calculations indicate that the complex with the conformation ruled out by the Q41R substitution is a non-reactive species due to its higher free energy barrier for the acylation reaction. Based on our calculations, we propose a kinetic mechanism that explains experimental data showing an increase of apparent rate constants for MAVS cleavage in Q41R mutants. Our model predicts that the non-reactive conformation of the enzyme-substrate complex modulates reaction kinetics like an uncompetitive inhibitor.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

An automated protocol for modelling peptide substrates to proteases

BACKGROUND

Proteases are key drivers in many biological processes, in part due to their specificity towards their substrates. However, depending on the family and molecular function, they can also display substrate promiscuity which can also be essential. Databases compiling specificity matrices derived from experimental assays have provided valuable insights into protease substrate recognition. Despite this, there are still gaps in our knowledge of the structural determinants. Here, we compile a set of protease crystal structures with bound peptide-like ligands to create a protocol for modelling substrates bound to protease structures, and for studying observables associated to the binding recognition.

RESULTS

As an application, we modelled a subset of protease-peptide complexes for which experimental cleavage data are available to compare with informational entropies obtained from protease-specificity matrices. The modelled complexes were subjected to conformational sampling using the Backrub method in Rosetta, and multiple observables from the simulations were calculated and compared per peptide position. We found that some of the calculated structural observables, such as the relative accessible surface area and the interaction energy, can help characterize a protease's substrate recognition, giving insights for the potential prediction of novel substrates by combining additional approaches.

CONCLUSION

Overall, our approach provides a repository of protease structures with annotated data, and an open source computational protocol to reproduce the modelling and dynamic analysis of the protease-peptide complexes.

A De Novo Designed Esterase with p-Nitrophenyl Acetate Hydrolysis Activity

Esterases are a large family of enzymes with wide applications in the industry. However, all esterases originated from natural sources, limiting their use in harsh environments or newly- emerged reactions. In this study, we designed a new esterase to develop a new protocol to satisfy the needs for better biocatalysts. The ideal spatial conformation of the serine catalytic triad and the oxygen anion hole at the substrate-binding site was constructed by quantum mechanical calculation. The catalytic triad and oxygen anion holes were then embedded in the protein scaffold using the new enzyme protocol in Rosetta 3. The design results were subsequently evaluated, and optimized designs were used for expression and purification. The designed esterase had significant lytic activities towards p-nitrophenyl acetate, which was confirmed by point mutations. Thus, this study developed a new protocol to obtain novel enzymes that may be useful in unforgiving environments or novel reactions.

Zika virus proteins at an atomic scale: how does structural biology help us to understand and develop vaccines and drugs against Zika virus infection?

In Brazil and in other tropical areas Zika virus infection was directly associated with clinical complications as microcephaly in newborn children whose mothers were infected during pregnancy and the Guillain-Barré syndrome in adults. Recently, research has been focused on developing new vaccines and drug candidates against Zika virus infection since none of those are available. In order to contribute to vaccine and drug development efforts, it becomes important the understanding of the molecular basis of the Zika virus recognition, infection and blockade. To this purpose, it is essential the structural determination of the Zika virus proteins. The genome sequencing of the Zika virus identified ten proteins, being three structural (protein E, protein C and protein prM) and seven non-structural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5). Together, these proteins are the main targets for drugs and antibody recognition. Here we examine new discoveries on high-resolution structural biology of Zika virus, observing the interactions and functions of its proteins identified via state-of-art structural methodologies as X-ray crystallography, nuclear magnetic resonance spectroscopy and cryogenic electronic microscopy. The aim of the present study is to contribute to the understanding of the structural basis of Zika virus infection at an atomic level and to point out similarities and differences to others flaviviruses.

The reaction mechanism of Zika virus NS2B/NS3 serine protease inhibition by dipeptidyl aldehyde: a QM/MM study

Zika virus (ZIKV) infection has become a global public health problem, associated with microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there are no commercially available anti-ZIKV drugs. The viral protease NS2B/NS3, which is involved in viral replication and maturation, is a potential drug target. Peptidomimetic aldehyde inhibitors bind covalently to the catalytic S135 of the NS3 protease. Here, we apply hybrid quantum mechanics/molecular mechanics (QM/MM) free-energy simulations at the PDDG-PM3/ff14SB level to investigate the inhibition mechanism of the ZIKV protease by a dipeptidyl aldehyde inhibitor (acyl-KR-aldehyde). The results show that proton transfer from the catalytic S135 to H51 occurs in concert with nucleophilic addition on the aldehyde warhead by S135. The anionic covalent complex between the dipeptidyl aldehyde and the ZIKV protease is analogous to the tetrahedral intermediate for substrate hydrolysis. Spontaneous protonation by H51 forms the hemiacetal. In addition, we use correlated ab initio QM/MM potential energy path calculations at levels up to LCCSD(T)/(aug)-cc-pVTZ to obtain accurate potential energy profiles of the reaction, which also support a concerted mechanism. These results provide detailed insight into the mechanism of ZIKV protease inhibition by a peptidyl aldehyde inhibitor, which will guide in the design of inhibitors.

Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate

Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol^-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.