Self-Masked Aldehyde Inhibitors: A Novel Strategy for Inhibiting Cysteine Proteases.

Cysteine proteases comprise an important class of drug targets, especially for infectious diseases such as Chagas disease (cruzain) and COVID-19 (3CL protease, cathepsin L). Peptide aldehydes have proven to be potent inhibitors for all of these proteases. However, the intrinsic, high electrophilicity of the aldehyde group is associated with safety concerns and metabolic instability, limiting the use of aldehyde inhibitors as drugs. We have developed a novel class of self-masked aldehyde inhibitors (SMAIs) for cruzain, the major cysteine protease of the causative agent of Chagas disease-Trypanosoma cruzi. These SMAIs exerted potent, reversible inhibition of cruzain (Ki* = 18-350 nM) while apparently protecting the free aldehyde in cell-based assays. We synthesized prodrugs of the SMAIs that could potentially improve their pharmacokinetic properties. We also elucidated the kinetic and chemical mechanism of SMAIs and applied this strategy to the design of anti-SARS-CoV-2 inhibitors.

Using biochemistry and biophysics to extinguish androgen receptor signaling in prostate cancer.

Castration resistant prostate cancer (CRPC) continues to be androgen receptor (AR) driven. Inhibition of AR signaling in CRPC could be advanced using state-of-the-art biophysical and biochemical techniques. Structural characterization of AR and its complexes by cryo-electron microscopy would advance the development of N-terminal domain (NTD) and ligand-binding domain (LBD) antagonists. The structural basis of AR function is unlikely to be determined by any single structure due to the intrinsic disorder of its NTD, which not only interacts with coregulators but likely accounts for the constitutive activity of AR-splice variants (SV), which lack the LBD and emerge in CRPC. Using different AR constructs lacking the LBD, their effects on protein folding, DNA binding, and transcriptional activity could reveal how interdomain coupling explains the activity of AR-SVs. The AR also interacts with coregulators that promote chromatin looping. Elucidating the mechanisms involved can identify vulnerabilities to treat CRPC, which do not involve targeting the AR. Phosphorylation of the AR coactivator MED-1 by CDK7 is one mechanism that can be blocked by the use of CDK7 inhibitors. CRPC gains resistance to AR signaling inhibitors (ARSI). Drug resistance may involve AR-SVs, but their role requires their reliable quantification by SILAC-mass spectrometry during disease progression. ARSI drug resistance also occurs by intratumoral androgen biosynthesis catalyzed by AKR1C3 (type 5 17ß-hydroxysteroid dehydrogenase), which is unique in that its acts as a coactivator of AR. Novel bifunctional inhibitors that competitively inhibit AKR1C3 and block its coactivator function could be developed using reverse-micelle NMR and fragment-based drug discovery.

On the ability of molecular dynamics force fields to recapitulate NMR derived protein side chain order parameters.

Molecular dynamics (MD) simulations have become a central tool for investigating various biophysical questions with atomistic detail. While many different proxies are used to qualify MD force fields, most are based on largely structural parameters such as the root mean square deviation from experimental coordinates or nuclear magnetic resonance (NMR) chemical shifts and residual dipolar couplings. NMR derived Lipari-Szabo squared generalized order parameter (O(2) ) values of amide NH bond vectors of the polypeptide chain were also often employed for refinement and validation. However, with a few exceptions, side chain methyl symmetry axis order parameters have not been incorporated into experimental reference sets. Using a test set of five diverse proteins, the performance of several force fields implemented in the NAMDD simulation package was examined. It was found that simulations employing explicit water implemented using the TIP3 model generally performed significantly better than those using implicit water in reproducing experimental methyl symmetry axis O(2) values. Overall the CHARMM27 force field performs nominally better than two implementations of the Amber force field. It appeared that recent quantum mechanics modifications to side chain torsional angles of leucine and isoleucine in the Amber force field have significantly hindered proper motional modeling for these residues. There remained significant room for improvement as even the best correlations of experimental and simulated methyl group Lipari-Szabo generalized order parameters fall below an R(2) of 0.8.

Electrochemical and structural coupling of the naphthoquinone amino acid.

As a prelude to engineering artificial energy conversion proteins emulating biology, we examine the inclusion of a synthetic naphthoquinone amino acid in a characterized host-guest protein and determine the effects of its quinone and hydroquinone forms on the helix-coil distribution.

NMR analysis of the alphaIIb beta3 cytoplasmic interaction suggests a mechanism for integrin regulation.

The integrin αIIbß3 is a transmembrane (TM) heterodimeric adhesion receptor that exists in equilibrium between resting and active ligand binding conformations. In resting αIIbß3, the TM and cytoplasmic domains of αIIb and ß3 form a heterodimer that constrains αIIbß3 in its resting conformation. To study the structure and dynamics of the cytoplasmic domain heterodimer, we prepared a disulfide-stabilized complex consisting of portions of the TM domains and the full cytoplasmic domains. NMR and hydrogen-deuterium exchange of this complex in micelles showed that the αIIb cytoplasmic domain is largely disordered, but it interacts with and influences the conformation of the ß3 cytoplasmic domain. The ß3 cytoplasmic domain consists of a stable proximal helix contiguous with the TM helix and two distal amphiphilic helices. To confirm the NMR structure in a membrane-like environment, we studied the ß3 cytoplasmic domain tethered to phospholipid bilayers. Hydrogen-deuterium exchange mass spectrometry, as well as circular dichroism spectroscopy, demonstrated that the ß3 cytoplasmic domain becomes more ordered and helical under these conditions, consistent with our NMR results. Further, these experiments suggest that the two distal helices associate with lipid bilayers but undergo fluctuations that would allow rapid binding of cytoplasmic proteins regulating integrin activation, such as talin and kindlin-3. Thus, these results provide a framework for understanding the kinetics and thermodynamics of protein interactions involving integrin cytoplasmic domains and suggest that such interactions act in a concerted fashion to influence integrin stalk separation and exposure of extracellular ligand binding sites.

Reverse micelle encapsulation of membrane-anchored proteins for solution NMR studies.

Perhaps 5%-10% of proteins bind to membranes via a covalently attached lipid. Posttranslational attachment of fatty acids such as myristate occurs on a variety of viral and cellular proteins. High-resolution information about the nature of lipidated proteins is remarkably sparse, often because of solubility problems caused by the exposed fatty acids. Reverse micelle encapsulation is used here to study two myristoylated proteins in their lipid-extruded states: myristoylated recoverin, which is a switch in the Ca(2+) signaling pathway in vision, and the myristoylated HIV-1 matrix protein, which is postulated to be targeted to the plasma membrane through its binding to phosphatidylinositol-4,5-bisphosphate. Both proteins have been successfully encapsulated in the lipid-extruded state and high-resolution NMR spectra obtained. Both proteins bind their activating ligands in the reverse micelle. This approach seems broadly applicable to membrane proteins with exposed fatty acid chains that have eluded structural characterization by conventional approaches.

Nativelike structure in designed four alpha-helix bundles driven by buried polar interactions.

We show that a single internal polar interaction per helix is sufficient to engender structural specificity in that helix in helical bundle proteins. Furthermore, we use histidine-binding cofactors of different shapes which bind directly into the core, demonstrating that this structural specificity is not the result of a prescribed complimentary, "knobs in holes" core packing. We show that we can switch structural specificity of individual helices on and off by ligating cofactors, singly and in pairs, which bind either one or two histidine ligands. To our knowledge, this is the first demonstration of such extensive manipulation of protein structure by ligand binding, an important result of general interest to those working with self-assembled molecular systems. Finally, as these proteins were designed without the use of computational modeling, we not only demonstrate that designing a uniquely structured cofactor binding protein is not as difficult as is generally believed, we have determined why this is so: hydrophobic core complementarity, which is very difficult to design, is not necessary. Instead, a much simpler design process entails the creation of core polar interactions which themselves can drive conformational specificity.

Characterization of the backbone and side chain dynamics of the CaM-CaMKIp complex reveals microscopic contributions to protein conformational entropy.

Calmodulin is a central mediator of calcium-dependent signal transduction pathways and regulates the activity of a large number of diverse targets. Calcium-dependent interactions of calmodulin with regulated proteins are of generally high affinity but of quite variable thermodynamic origins. Here we investigate the influence of the binding of the calmodulin-binding domain of calmodulin kinase I on the fast internal dynamics of calcium-saturated calmodulin. NMR relaxation was used to probe motion on the backbone (viewed through the backbone amide NH group) and the side chains (viewed through methyl groups). The distribution of the amplitudes of side chain dynamics is trimodal. The microscopic details of side chain motion are compared with those of a thermodynamically and structurally similar complex of calmodulin with the calmodulin-binding domain of the smooth muscle myosin light chain kinase. While there are no significant differences in backbone dynamics and no net change in methyl-bearing side chain dynamics, a large redistribution of the amplitude of methyl dynamics is observed between the two complexes. The variation in dynamics was largely localized to the heterogeneously dynamic target-binding interface, suggesting that differential dynamics of the binding surface plays a functional role in the high-affinity binding interactions of calmodulin. These results begin to reveal a fundamental role for residual protein entropy in molecular recognition by calmodulin.

Arylamide derivatives as peptidomimetic inhibitors of calmodulin.

[structure: see text] Many peptides bind to calmodulin (CaM) in a helical conformation. Here we describe a group of synthetic inhibitors of CaM based on an arylamide scaffold that is intended to mimic smMLCK, a CaM-binding helical peptide. Compound 1 showed a K(i) value of 7.10 +/- 1.48 nM in a fluorescence polarization assay that monitors the strong association of CaM and its peptide ligand mastoparan X. ((1)H,(15)N)-HSQC NMR spectroscopy experiments suggested that 1 binds to CaM in an analogous fashion to that of smMLCK.

Backbone and side chain dynamics of mutant calmodulin-peptide complexes.

The mechanism of long-range coupling of allosteric sites in calcium-saturated calmodulin (CaM) has been explored by characterizing structural and dynamics effects of mutants of calmodulin in complex with a peptide corresponding to the smooth muscle myosin light chain kinase calmodulin-binding domain (smMLCKp). Four CaM mutants were examined: D95N and D58N, located in Ca2+-binding loops; and M124L and E84K, located in the target domain-binding site of CaM. Three of these mutants have altered allosteric coupling either between Ca2+-binding sites (D58N and D95N) or between the target- and Ca2+-binding sites (E84K). The structure and dynamics of the mutant calmodulins in complex with smMLCKp were characterized using solution NMR. Analysis of chemical shift perturbations was employed to detect largely structural perturbations. 15N and 2H relaxation was employed to detect perturbations of the dynamics of the backbone and methyl-bearing side chains of calmodulin. The least median squares method was found to be robust in the detection of perturbed sites. The main chain dynamics of calmodulin are found to be largely unresponsive to the mutations. Three mutants show significantly perturbed dynamics of methyl-bearing side chains. Despite the pseudosymmetric location of Ca2+-binding loop mutations D58N and D95N, the dynamic response of CaM is asymmetric, producing long-range perturbation in D58N and almost none in D95N. The mutations located at the target domain-binding site have quite different effects. For M124L, a local perturbation of the methyl dynamics is observed, while the E84K mutation produces a long-range propagation of dynamic perturbations along the target domain-binding site.

Changes in calmodulin main-chain dynamics upon ligand binding revealed by cross-correlated NMR relaxation measurements.

The fast dynamics of protein backbones are often investigated by nuclear magnetic relaxation experiments that report on the degree of spatial restriction of the amide bond vector. By comparing calmodulin in the peptide-bound and peptide-free states with these classical methods, we observe little difference in the dynamics of the polypeptide main chain (average order parameter decrease of 0.01 unit upon binding). However, when using NMR methods that monitor the mobility of the CO-Calpha bond vector, we reveal a significant reduction of dynamics of the protein main chain (average order parameter decrease of 0.048 units). Previous investigations have suggested that the side-chain dynamics is reduced by an average of 0.07 order parameter units upon ligand binding (Lee, A. L.; Kinnear, S. A.; Wand, A. J. Nat. Struct. Biol. 2000, 7, 72-77). The current findings suggest that the change of the CO-Calpha bond vector dynamics is intermediate between the changes in NH and side-chain dynamics and report a previously undetected loss of main-chain entropy. Weak site-to-site correlations between the different motional indicators are also observed.

De novo design of a D2-symmetrical protein that reproduces the diheme four-helix bundle in cytochrome bc1.

An idealized, water-soluble D(2)-symmetric diheme protein is constructed based on a mathematical parametrization of the backbone coordinates of the transmembrane diheme four-helix bundle in cytochrome bc(1). Each heme is coordinated by two His residues from diagonally apposed helices. In the model, the imidazole rings of the His ligands are held in a somewhat unusual perpendicular orientation as found in cytochrome bc(1), which is maintained by a second-shell hydrogen bond to a Thr side chain on a neighboring helix. The resulting peptide is unfolded in the apo state but assembles cooperatively upon binding to heme into a well-folded tetramer. Each tetramer binds two hemes with high affinity at low micromolar concentrations. The equilibrium reduction midpoint potential varies between -76 mV and -124 mV vs SHE in the reducing and oxidizing direction, respectively. The EPR spectrum of the ferric complex indicates the presence of a low-spin species, with a g(max) value of 3.35 comparable to those obtained for hemes b of cytochrome bc(1) (3.79 and 3.44). This provides strong support for the designed perpendicular orientation of the imidazole ligands. Moreover, NMR spectra show that the protein exists in solution in a unique conformation and is amenable to structural studies. This protein may provide a useful scaffold for determining how second-shell ligands affect the redox potential of the heme cofactor.

Cooperativity and specificity of association of a designed transmembrane peptide.

Thermodynamics studies aimed at quantitatively characterizing free energy effects of amino acid substitutions are not restricted to two state systems, but do require knowing the number of states involved in the equilibrium under consideration. Using analytical ultracentrifugation and NMR methods, we show here that a membrane-soluble peptide, MS1, designed by modifying the sequence of the water-soluble coiled-coil GCN4-P1, exhibits a reversible monomer-dimer-trimer association in detergent micelles with a greater degree of cooperativity in C14-betaine than in dodecyl phosphocholine detergents.

A direct test of the reductionist approach to structural studies of calmodulin activity: relevance of peptide models of target proteins.

Ca(2+)-saturated calmodulin (CaM) directly associates with and activates CaM-dependent protein kinase I (CaMKI) through interactions with a short sequence in its regulatory domain. Using heteronuclear NMR (13)C-(15)N-(1)H correlation experiments, the backbone assignments were determined for CaM bound to a peptide (CaMKIp) corresponding to the CaM-binding sequence of CaMKI. A comparison of chemical shifts for free CaM with those of the CaM. CaMKIp complex indicate large differences throughout the CaM sequence. Using NMR techniques optimized for large proteins, backbone resonance assignments were also determined for CaM bound to the intact CaMKI enzyme. NMR spectra of CaM bound to either the CaMKI enzyme or peptide are virtually identical, indicating that calmodulin is structurally indistinguishable when complexed to the intact kinase or the peptide CaM-binding domain. Chemical shifts of CaM bound to a peptide (smMLCKp) corresponding to the calmodulin-binding domain of smooth muscle myosin light chain kinase are also compared with the CaM. CaMKI complexes. Chemical shifts can differentiate one complex from another, as well as bound versus free states of CaM. In this context, the observed similarity between CaM. CaMKI enzyme and peptide complexes is striking, indicating that the peptide is an excellent mimetic for interaction of calmodulin with the CaMKI enzyme.