Molecular modeling of ARF6 dysregulation caused by mutations in IQSEC2.

IQSEC2 gene mutations are associated with epilepsy, autism, and intellectual disability. The primary function IQSEC2, mediated via its Sec 7 domain, is to act as a guanine nucleotide exchange factor for ARF6. We sought to develop a molecular model, which may explain the aberrant Sec 7 activity on ARF6 of different human IQSEC2 mutations. We integrated experimental data of IQSEC2 mutants with protein structure prediction by the RaptorX server combined with molecular modeling and molecular dynamics simulations. Normally, apocalmodulin (apoCM) binds to IQSEC2 resulting in its N-terminal fragment inhibiting access of its Sec 7 domain to ARF6. An increase in Ca2+ concentration destabilizes the interaction of IQSEC2 with apoCM and removes steric hindrance of Sec 7 binding with ARF6. Mutations at amino acid residue 350 of IQSEC2 result in loss of steric hindrance of Sec 7 binding with ARF6 leading to constitutive activation of ARF6 by Sec 7. On the other hand, a mutation at amino acid residue 359 of IQSEC2 results in constitutive hindrance of Sec 7 binding to ARF6 leading to the loss of the ability of IQSEC2 to activate ARF6. These studies provide a model for dysregulation of IQSEC2 Sec 7 activity by mutant IQSEC2 proteins.Communicated by Ramaswamy H. Sarma.

Phage-Displayed Mimotopes of SARS-CoV-2 Spike Protein Targeted to Authentic and Alternative Cellular Receptors.

The evolution of the SARS-CoV-2 virus during the COVID-19 pandemic was accompanied by the emergence of new heavily mutated viral variants with increased infectivity and/or resistance to detection by the human immune system. To respond to the urgent need for advanced methods and materials to empower a better understanding of the mechanisms of virus's adaptation to human host cells and to the immuno-resistant human population, we suggested using recombinant filamentous bacteriophages, displaying on their surface foreign peptides termed "mimotopes", which mimic the structure of viral receptor-binding sites on the viral spike protein and can serve as molecular probes in the evaluation of molecular mechanisms of virus infectivity. In opposition to spike-binding antibodies that are commonly used in studying the interaction of the ACE2 receptor with SARS-CoV-2 variants in vitro, phage spike mimotopes targeted to other cellular receptors would allow discovery of their role in viral infection in vivo using cell culture, tissue, organs, or the whole organism. Phage mimotopes of the SARS-CoV-2 Spike S1 protein have been developed using a combination of phage display and molecular mimicry concepts, termed here "phage mimicry", supported by bioinformatics methods. The key elements of the phage mimicry concept include: (1) preparation of a collection of p8-type (landscape) phages, which interact with authentic active receptors of live human cells, presumably mimicking the binding interactions of human coronaviruses such as SARS-CoV-2 and its variants; (2) discovery of closely related amino acid clusters with similar 3D structural motifs on the surface of natural ligands (FGF1 and NRP1), of the model receptor of interest FGFR and the S1 spike protein; and (3) an ELISA analysis of the interaction between candidate phage mimotopes with FGFR3 (a potential alternative receptor) in comparison with ACE2 (the authentic receptor).

Development and characterisation of SMURF2-targeting modifiers.

The C2-WW-HECT-domain E3 ubiquitin ligase SMURF2 emerges as an important regulator of diverse cellular processes. To date, SMURF2-specific modulators were not developed. Here, we generated and investigated a set of SMURF2-targeting synthetic peptides and peptidomimetics designed to stimulate SMURF2's autoubiquitination and turnover via a disruption of the inhibitory intramolecular interaction between its C2 and HECT domains. The results revealed the effects of these molecules both in vitro and in cellulo at the nanomolar concentration range. Moreover, the data showed that targeting of SMURF2 with either these modifiers or SMURF2-specific shRNAs could accelerate cell growth in a cell-context-dependent manner. Intriguingly, a concomitant cell treatment with a selected SMURF2-targeting compound and the DNA-damaging drug etoposide markedly increased the cytotoxicity produced by this drug in growing cells. Altogether, these findings demonstrate that SMURF2 can be druggable through its self-destructive autoubiquitination, and inactivation of SMURF2 might be used to affect cell sensitivity to certain anticancer drugs.

How does the exosite of rhomboid protease affect substrate processing and inhibition?

Rhomboid proteases constitute a family of intramembrane serine proteases ubiquitous in all forms of life. They differ in many aspects from their soluble counterparts. We applied molecular dynamics (MD) computational approach to address several challenging issues regarding their catalytic mechanism: How does the exosite of GlpG rhomboid protease control the kinetics efficiency of substrate hydrolysis? What is the mechanism of inhibition by the non-competitive peptidyl aldehyde inhibitors bound to the GlpG rhomboid active site (AS)? What is the underlying mechanism that explains the hypothesis that GlpG rhomboid protease is not adopted for the hydrolysis of short peptides that do not contain a transmembrane domain (TMD)? Two fundamental features of rhomboid catalysis, the enzyme recognition and discrimination of substrates by TMD interactions in the exosite, and the concerted mechanism of non-covalent pre-catalytic complex to covalent tetrahedral complex (TC) conversion, provide answers to these mechanistic questions.

Stepwise Versus Concerted Mechanisms in General-Base Catalysis by Serine Proteases.

General-base catalysis in serine proteases still poses mechanistic challenges despite decades of research. Whether proton transfer from the catalytic Ser to His and nucleophilic attack on the substrate are concerted or stepwise is still under debate, even for the classical Asp-His-Ser catalytic triad. To address these key catalytic steps, the transformation of the Michaelis complex to tetrahedral complex in the covalent inhibition of two prototype serine proteases was studied: chymotrypsin (with the catalytic triad) inhibition by a peptidyl trifluoromethane and GlpG rhomboid (with Ser-His dyad) inhibition by an isocoumarin derivative. The sampled MD trajectories of averaged pKa  values of catalytic residues were QM calculated by the MD-QM/SCRF(VS) method on molecular clusters simulating the active site. Differences between concerted and stepwise mechanisms are controlled by the dynamically changing pKa  values of the catalytic residues as a function of their progressively reduced water exposure, caused by the incoming ligand.

A new method for filtering of reactive "warheads" of transition-state analog protease inhibitors.

In light of the major contribution of the reactive warhead to the binding energy trend in reversible covalent transition-state analog inhibitors of serine and cysteine hydrolases, would it be possible to rationally design and quickly filter such warheads, especially for large-scale screening? The previously defined W1 and W2 covalent descriptors quantitatively account for the energetic effect of the covalent bonds reorganization, accompanying enzyme-inhibitor covalent binding. The quantum mechanically calculated W1 and W2 reflect the warhead binding energy by modeling of the enzyme-inhibitor reaction core. Here, we demonstrate the use of these descriptors for warhead filtering, and examine its scope and limitations. The W1 and W2 descriptors provide a tool for rational design of various warheads as universal building blocks of real inhibitors without the requirement of 3D structural information about the target enzyme or QSAR studies. These warheads could then be used as hit structural templates in the subsequent optimization of inhibitors recognition sites.

Application of EMBM to Structure-Based Design of Warheads for Protease Inhibitors.

Most CADD tools handle non-covalent enzyme inhibitors, despite the growing interest of the pharma industry in covalent inhibitors. We have recently introduced an enzyme mechanism-based method, EMBM, as a computational tool for binding trend analysis and prediction of chemical sites (CS) of reversible covalent enzyme inhibitors. In the current study we demonstrate the utility of EMBM to structure-based applications. In this mode, the energy of the enzyme-inhibitor covalent bond is accounted for by the W1 and W2 covalent descriptors we have developed, whereas the non-covalent interactions between the inhibitor CS and the enzyme active site can be estimated directly on the 3D structure of the enzyme-inhibitor complex.

The Catalytic Machinery of Rhomboid Proteases: Combined MD and QM Simulations.

Rhomboid proteases are a ubiquitous family of intramembrane serine proteases in prokaryotic and eukaryotic organisms that cleave membrane proteins in their transmembrane region. Their catalytic activity is centered at a His-Ser catalytic dyad. We applied molecular dynamics and quantum mechanics calculations in order to clarify the protonation state of the catalytic residues of E. coli GlpG rhomboid protease and how it is affected by the immersion of the enzyme in the membrane. We identified (Nε)H150(dpr)_H254(dpr)_S201(pr) as the protonation (and H150 tautomeric) state of free GlpG in both lipid-solubilized and membrane environments. We used our MD-QM/SCRF(VS) computational protocol to rationalize and predict the trend of pKa change caused by the decrease of water exposure of the active site of GlpG due to ligand binding. The catalytic diad of lipid-solubilized GlpG exists as an H254(+)_S201(-) ion pair at the Michaelis complex stage, with Ser201 ready for nucleophilic attack on the substrate. Therefore, deprotonation of S201 does not contribute to the activation barrier of covalent tetrahedral complex formation. In contrast, both catalytic residues, H254 and S201, are neutral in the Michaelis complex of GlpG in the membrane. Therefore, S201 deprotonation by H254 general base catalysis should contribute to the activation barrier of the covalent tetrahedral complex formation.

The mechanism of papain inhibition by peptidyl aldehydes.

Various mechanisms for the reversible formation of a covalent tetrahedral complex (TC) between papain and peptidyl aldehyde inhibitors were simulated by DFT calculations, applying the quantum mechanical/self consistent reaction field (virtual solvent) [QM/SCRF(VS)] approach. Only one mechanism correlates with the experimental kinetic data. The His-Cys catalytic diad is in an N/SH protonation state in the noncovalent papain-aldehyde Michaelis complex. His159 functions as a general base catalyst, abstracting a proton from the Cys25, whereas the activated thiolate synchronously attacks the inhibitor's carbonyl group. The final product of papain inhibition is the protonated neutral form of the hemithioacetal TC(OH), in agreement with experimental data. The predicted activation barrier g enz≠ = 5.2 kcal mol⁻¹ is close to the experimental value of 6.9 kcal mol⁻¹. An interpretation of the experimentally observed slow binding effect for peptidyl aldehyde inhibitors is presented. The calculated g cat≠ is much lower than the rate determining activation barrier of hemithioacetal formation in water, g w≠, in agreement with the concept that the preorganized electrostatic environment in the enzyme active site is the driving force of enzyme catalysis. We have rationalized the origin of the acidic and basic pK(a)'s on the k2/K(S) versus pH bell-shaped profile of papain inhibition by peptidyl aldehydes.

EMBM - a new enzyme mechanism-based method for rational design of chemical sites of covalent inhibitors.

We introduce an enzyme mechanism-based method (EMBM) aimed at rational design of chemical sites (CS) of reaction coordinate analog inhibitors. The energy of valence reorganization of CS, caused by the formation of the enzyme-inhibitor covalent complex, is accounted for by new covalent descriptors W1 and W2. We considered CS fragments with a carbonyl reactivity center, like in native protease substrates. The W1 and W2 descriptors are calculated quantum mechanically on small molecular clusters simulating the reaction core of the formed covalent tetrahedral complex, anionic TC(O-) or neutral TC(OH). The modeling on a reaction core allows generation of various CS and corresponding TC(O-) and TC(OH) as universal building blocks of real inhibitors and their covalent complexes with serine or cysteine hydrolases. Moreover, the approach avoids the need for 3D structure of the target enzyme, so EMBM may be used for ligand-based design. We have built a chemical site of inhibitors (CSI) databank with pairs of W1 and W2 descriptors precalculated for both CH3O(-) and CH3S(-) nucleophiles for every collected CS fragment. We demonstrated that contribution of a CS fragment to the binding affinity of an inhibitor depends on both its covalent reorganization during the chemical transformation and its noncovalent interactions in the enzyme active site. Consequently, prediction of inhibitors binding trend can be done only by accounting for all of these factors, using W1 and W2 in combination with noncovalent QSAR descriptors.

Challenging a paradigm: theoretical calculations of the protonation state of the Cys25-His159 catalytic diad in free papain.

A central mechanistic paradigm of cysteine proteases is that the His-Cys catalytic diad forms an ion-pair NH(+)/S(-) already in the catalytically active free enzyme. Most molecular modeling studies of cysteine proteases refer to this paradigm as their starting point. Nevertheless, several recent kinetics and X-ray crystallography studies of viral and bacterial cysteine proteases depart from the ion-pair mechanism, suggesting general base catalysis. We challenge the postulate of the ion-pair formation in free papain. Applying our QM/SCRF(VS) molecular modeling approach, we analyzed all protonation states of the catalytic diad in free papain and its SMe derivative, comparing the predicted and experimental pK(a) data. We conclude that the His-Cys catalytic diad in free papain is fully protonated, NH(+)/SH. The experimental pK(a) = 8.62 of His159 imidazole in free papain, obtained by NMR-controlled titration and originally interpreted as the NH(+)/S(-) <==> N/S(-) NH(+)/S(-) <==> N/S(-) equilibrium, is now assigned to the NH(+)/SH <==> N/SH NH(+)/SH <==> N/SH equilibrium.

Peptidyl epoxides extended in the P' direction as cysteine protease inhibitors: effect on affinity and mechanism of inhibition.

Endo peptidyl epoxides, in which the central epoxidic moiety replaces the scissile amide bond of a P(3)-P(3)' peptide, were designed as cysteine proteases inhibitors. The additional P'-S' interactions, relative to those of an exo peptidyl epoxide of the same P(3)-P(1) sequence, significantly improved affinity to the enzymes papain and cathepsin B, but also changed the mode of inhibition from active-site directed inactivation to reversible competitive inhibition. Computational models rationalize the binding affinity and the inhibition mechanism.

Screening of the active site from water by the incoming ligand triggers catalysis and inhibition in serine proteases.

The pKa of the catalytic His57 N(epsilon)H in the tetrahedral complex (TC) of chymotrypsin with trifluoromethyl ketone inhibitors is 4-5 units higher relative to the free enzyme (FE). Such stable TC's, formed with transition state (TS) analog inhibitors, are topologically similar to the catalytic TS. Thus, analysis of this pKa shift may shed light on the role of water solvation in the general base catalysis by histidine. We applied our QM/SCRF(VS) approach to study this shift. The method enables explicit quantum mechanical DFT calculations of large molecular clusters that simulate chemical reactions at the active site (AS) of water solvated enzymes. We derived an analytical expression for the pKa dependence on the degree of water exposure of the ionizable group, and on the total charge in the enzyme AS, Q(A) and Q(B), when the target ionizable functional group (His57 in this study) is in the acidic (A) and basic (B) forms, respectively. Q2(B) > Q2(A) both in the FE and in the TC of chymotrypsin. Therefore, water solvation decreases the relative stability of the protonated histidine in both. Ligand binding reduces the degree of water solvation of the imidazole ring, and consequently elevates the histidine pKa. Thus, the binding of the ligand plays a triggering role that switches on the cascade of catalytic reactions in serine proteases.

The cooperative effect between active site ionized groups and water desolvation controls the alteration of acid/base catalysis in serine proteases.

What is the driving force that alters the catalytic function of His57 in serine proteases between general base and general acid in each step along the enzymatic reaction? The stable tetrahedral complexes (TC) of chymotrypsin with trifluoromethyl ketone transition state analogue inhibitors are topologically similar to the catalytic transition state. Therefore, they can serve as a good model to study the enzyme catalytic reaction. We used DFT quantum mechanical calculations to analyze the effect of solvation and of polar factors in the active site of chymotrypsin on the pKa of the catalytic histidine in FE (the free enzyme), EI (the noncovalent enzyme inhibitor complex), and TC. We demonstrated that the acid/base alteration is controlled by the charged groups in the active site--the catalytic Asp102 carboxylate and the oxyanion. The effect of these groups on the catalytic His is modulated by water solvation of the active site.

Enzyme isoselective inhibitors: a tool for binding-trend analysis.

A transition-state analogue inhibitor that covalently reversibly binds to an enzyme formally consists of two parts: the chemical site, CS and the recognition site, RS. We have experimentally and theoretically demonstrated that the trend of binding affinity in a series of isoselective inhibitors (with identical RS and different CS fragments) depends mainly on their CS fragments. Isoselective inhibitors have the same affinity trend toward different enzymes of the same family with a common catalytic mechanism. Thus, very good correlation between experimentally determined and theoretically calculated Ki values was demonstrated. A practical outcome is the application of the described method as a tool for an expert analysis in virtual screening of inhibitor libraries and in the design of new enzyme inhibitors.

Enzyme isoselective inhibitors: application to drug design.

Common methodologies of computer-assisted drug design focus on noncovalent enzyme-ligand interactions. We introduced enzyme isoselective inhibition trend analysis as a tool for the expert analysis of covalent reversible inhibitors. The methodology is applied to predict the binding affinities of a series of transition-state analogue inhibitors of medicinally important serine and cysteine hydrolases. These inhibitors are isoselective: they have identical noncovalent recognition fragments (RS) and different reactive chemical fragments (CS). Furthermore, it is possible to predict the binding affinities of a series of isoselective inhibitors toward a prototype enzyme and to extrapolate the data to a target medicinally important enzyme of the same family. Rational design of CS fragments followed by conventional RS optimization could be used as a novel approach to drug design.

Identification of protons position in acid-base enzyme catalyzed reactions: the hepatitis C viral NS3 protease.

General acid-base catalysis is a key element of the catalytic activity of most enzymes. Therefore, any explicit molecular modeling of enzyme-catalyzed chemical reactions requires correct identification of protons location on the catalytic groups. In this work, we apply our quantum mechanical/self-consistent reaction field in virtual solvent [QM/SCRF(VS)] method for identification of the position of protons shared by the enzyme catalytic groups and the polar groups of the inhibitor in a covalent tetrahedral complex (TC) of the hepatitis C virus NS3 protease with a peptidyl alpha-ketoacid inhibitor. To identify the relevant protonation states, we have analyzed relative stabilities of R and S configurations of the TC that depend on the specific proton distribution over the polar groups and correlated it with experimental NMR and X-ray crystallography data, both at low and neutral pH ranges. The tentative assignment of the single resonance in the (13)C NMR spectrum of the hemiketal carbon at physiological pH to the S configuration of TC is confirmed. Both R and S configurations are equally stable at acidic pH in our modeling, in good agreement with the (13)C NMR observation.

Is there a weak H-bond --> LBHB transition on tetrahedral complex formation in serine proteases?

The transformation of a weak hydrogen bond in the free enzyme into a low-barrier hydrogen bond (LBHB) in the tetrahedral intermediate has been suggested as an important factor facilitating catalysis in serine proteases. In this work, we examine the structure of the H-bond in the Asp102-His57 diad of serine proteases in the free enzyme and in a covalent tetrahedral complex (TC) with a trifluoromethylketone inhibitor. We apply ab initio quantum mechanical calculations to models consisting of a large molecular fragment of the enzyme active site, and the combined effect of the rest of the protein body and the solvation by surrounding bulk water was simulated by a self-consistent reaction field method in our novel QM/SCRF(VS) approach. Potential profiles of adiabatic proton transfer in the Asp102-His57 diad in these model systems were calculated. We conclude that the hydrogen bond in both the free enzyme and in the enzyme-inhibitor TC is a strong ionic asymmetric one-well hydrogen bond, in contrast to a previous suggestion that it is a weak H-bond in the former and a double-well LBHB in the latter.

Nitrogen inversion in cyclic amines and the bicyclic effect.

The hitherto unsolved problem of the origin of the unusually high nitrogen inversion-rotation (NIR) barriers in 7-azabicyclo[2.2.1]heptanes (the bicyclic effect) was examined using the natural bond orbital (NBO) approach. Reinvestigating the NIR barrier for tropane by DNMR, we found that NIR barriers increase smoothly on going from nitrogen-bridged bicyclic systems of a larger ring size to the smaller ring homologous systems. The experimental NIR barriers are reproduced with good accuracy using the MP2/6-31G level of theory. The NBO analysis for these and other azabicycles led to the conclusion that the height of these barriers is mostly determined by the energy of the sigma-orbitals of the C(alpha)(-)C(beta) bonds as well as the nitrogen lone pair. Thus, the bicyclic effect is actually an extreme case of a common C(alpha-)N-C(alpha) tripyramid geometry-NIR barrier dependence for N-bridged bicyclic amines. By establishing the rate-determining role of the C(alpha-)N-C(alpha) tripyramid fragment for NIR, we have derived the first sufficiently accurate quantitative correlations amine geometry-NIR barrier for monocyclic as well as bicyclic N-H and N-Me amines (i.e., for an amine set which also includes the bicyclic effect systems).