Endogenous expression of inactive lysine deacetylases reveals deacetylation-dependent cellular mechanisms.

Acetylation of lysine residues is an important and common post-translational regulatory mechanism occurring on thousands of non-histone proteins. Lysine deacetylases (KDACs or HDACs) are a family of enzymes responsible for removing acetylation. To identify the biological mechanisms regulated by individual KDACs, we created HT1080 cell lines containing chromosomal point mutations, which endogenously express either KDAC6 or KDAC8 having single inactivated catalytic domain. Engineered HT1080 cells expressing inactive KDA6 or KDAC8 domains remained viable and exhibited enhanced acetylation on known substrate proteins. RNA-seq analysis revealed that many changes in gene expression were observed when KDACs were inactivated, and that these gene sets differed significantly from knockdown and knockout cell lines. Using GO ontology, we identified several critical biological processes associated specifically with catalytic activity and others attributable to non-catalytic interactions. Treatment of wild-type cells with KDAC-specific inhibitors Tubastatin A and PCI-34051 resulted in gene expression changes distinct from those of the engineered cell lines, validating this approach as a tool for evaluating in-cell inhibitor specificity and identifying off-target effects of KDAC inhibitors. Probing the functions of specific KDAC domains using these cell lines is not equivalent to doing so using previously existing methods and provides novel insight into the catalytic functions of individual KDACs by investigating the molecular and cellular changes upon genetic inactivation.

Tracing the birth of structural domains from loops during protein evolution.

The structures and functions of proteins are embedded into the loop scaffolds of structural domains. Their origin and evolution remain mysterious. Here, we use a novel graph-theoretical approach to describe how modular and non-modular loop prototypes combine to form folded structures in protein domain evolution. Phylogenomic data-driven chronologies reoriented a bipartite network of loops and domains (and its projections) into 'waterfalls' depicting an evolving 'elementary functionome' (EF). Two primordial waves of functional innovation involving founder 'p-loop' and 'winged-helix' domains were accompanied by an ongoing emergence and reuse of structural and functional novelty. Metabolic pathways expanded before translation functionalities. A dual hourglass recruitment pattern transferred scale-free properties from loop to domain components of the EF network in generative cycles of hierarchical modularity. Modeling the evolutionary emergence of the oldest P-loop and winged-helix domains with AlphFold2 uncovered rapid convergence towards folded structure, suggesting that a folding vocabulary exists in loops for protein fold repurposing and design.

Structural basis for the allosteric modulation of rhodopsin by nanobody binding to its extracellular domain.

Rhodopsin is a prototypical G protein-coupled receptor (GPCR) critical for vertebrate vision. Research on GPCR signaling states has been facilitated using llama-derived nanobodies (Nbs), some of which bind to the intracellular surface to allosterically modulate the receptor. Extracellularly binding allosteric nanobodies have also been investigated, but the structural basis for their activity has not been resolved to date. Here, we report a library of Nbs that bind to the extracellular surface of rhodopsin and allosterically modulate the thermodynamics of its activation process. Crystal structures of Nb2 in complex with native rhodopsin reveal a mechanism of allosteric modulation involving extracellular loop 2 and native glycans. Nb2 binding suppresses Schiff base deprotonation and hydrolysis and prevents intracellular outward movement of helices five and six - a universal activation event for GPCRs. Nb2 also mitigates protein misfolding in a disease-associated mutant rhodopsin. Our data show the power of nanobodies to modulate the photoactivation of rhodopsin and potentially serve as therapeutic agents for disease-associated rhodopsin misfolding.

Engineering gain-of-function mutants of a WW domain by dynamics and structural analysis.

Proteins gain optimal fitness such as foldability and function through evolutionary selection. However, classical studies have found that evolutionarily designed protein sequences alone cannot guarantee foldability, or at least not without considering local contacts associated with the initial folding steps. We previously showed that foldability and function can be restored by removing frustration in the folding energy landscape of a model WW domain protein, CC16, which was designed based on Statistical Coupling Analysis (SCA). Substitutions ensuring the formation of five local contacts identified as "on-path" were selected using the closest homolog native folded sequence, N21. Surprisingly, the resulting sequence, CC16-N21, bound to Group I peptides, while N21 did not. Here, we identified single-point mutations that enable N21 to bind a Group I peptide ligand through structure and dynamic-based computational design. Comparison of the docked position of the CC16-N21/ligand complex with the N21 structure showed that residues at positions 9 and 19 are important for peptide binding, whereas the dynamic profiles identified position 10 as allosterically coupled to the binding site and exhibiting different dynamics between N21 and CC16-N21. We found that swapping these positions in N21 with matched residues from CC16-N21 recovers nature-like binding affinity to N21. This study validates the use of dynamic profiles as guiding principles for affecting the binding affinity of small proteins.

Protoplasting and Fluorescence-Activated Cell Sorting of the Shoot Apical Meristem Cell Types.

The shoot apical meristems (SAMs) are located at the tip of the shoot apex. The SAM harbors stem cells that divide continually to provide cells for developing above-ground organs. Several important developmental events occur in SAMs, such as stem cell maintenance, organ differentiation, and flowering commitment which are under genetic control. The SAM is a collection of specialized cells organized in specific spatial domains. Deciphering the gene regulatory networks, guided by the developmental and environmental signals, in these discrete cell types is essential to decoding the SAM function. Here, I provide updates to the previously published protocols for the protoplasting and subsequent purification through fluorescence-activated cell sorting (FACS) of SAM cell types (Reddy, Fluorescence activated cell sorting of shoot apical meristem cell types. In: Riechmann JL, Wellmer F (eds) Flower development. Methods in molecular biology, vol 1110. Humana, New York, pp 315-321, 2014), which has provided genome-wide gene expression patterns at a single cell-type resolution.

Base Dynamics in the HhaI Protein Binding Site.

Protein-DNA interactions play an important role in numerous biological functions within the living cell. In many of these interactions, the DNA helix is significantly distorted upon protein-DNA complex formation. The HhaI restriction-modification system is one such system, where the methylation target is flipped out of the helix when bound to the methyltransferase. However, the base flipping mechanism is not well understood. The dynamics of the binding site of the HhaI methyltransferase and endonuclease (underlined) within the DNA oligomer [d(G1A2T3A4G5C6G7C8T9A10T11C12)]2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs deuterated on the base of nucleotides within and flanking the [5'-GCGC-3']2 sequence indicate that all of these positions are structurally flexible. Previously, conformational flexibility within the phosphodiester backbone and furanose ring within the target sequence has been observed and hypothesized to play a role in the distortion mechanism. However, whether that distortion was occurring through an active or passive mechanism remained unclear. These NMR data demonstrate that although the [5'-GCGC-3']2 sequence is dynamic, the target cytosine is not passively flipping out of the double-helix on the millisecond-picosecond time scale. Additionally, although previous studies have shown that both the furanose ring and phosphodiester backbone experience a change in dynamics upon methylation, which may play a role in recognition and cleavage by the endonuclease, our observations here indicate that methylation has no effect on the dynamics of the base itself.

Drug Discovery Targeting Nuclear Receptor Binding SET Domain Protein 2 (NSD2).

Nuclear receptor binding SET domain proteins (NSDs) catalyze the mono- or dimethylation of histone 3 lysine 36 (H3K36me1 and H3K36me2), using S-adenosyl-l-methionine (SAM) as a methyl donor. As a key member of the NSD family of proteins, NSD2 plays an important role in the pathogenesis and progression of various diseases such as cancers, inflammations, and infectious diseases, serving as a promising drug target. Developing potent and specific NSD2 inhibitors may provide potential novel therapeutics. Several NSD2 inhibitors and degraders have been discovered while remaining in the early stage of drug development. Excitingly, KTX-1001, a selective NSD2 inhibitor, has entered clinical trials. In this Perspective, the structures and functions of NSD2, its roles in various human diseases, and the recent advances in drug discovery strategies targeting NSD2 have been summarized. The challenges, opportunities, and future directions for developing NSD2 inhibitors and degraders are also discussed.

Quantitative flow cytometric selection of tau conformational nanobodies specific for pathological aggregates.

Single-domain antibodies, also known as nanobodies, are broadly important for studying the structure and conformational states of several classes of proteins, including membrane proteins, enzymes, and amyloidogenic proteins. Conformational nanobodies specific for aggregated conformations of amyloidogenic proteins are particularly needed to better target and study aggregates associated with a growing class of associated diseases, especially neurodegenerative disorders such as Alzheimer's and Parkinson's diseases. However, there are few reported nanobodies with both conformational and sequence specificity for amyloid aggregates, especially for large and complex proteins such as the tau protein associated with Alzheimer's disease, due to difficulties in selecting nanobodies that bind to complex aggregated proteins. Here, we report the selection of conformational nanobodies that selectively recognize aggregated (fibrillar) tau relative to soluble (monomeric) tau. Notably, we demonstrate that these nanobodies can be directly isolated from immune libraries using quantitative flow cytometric sorting of yeast-displayed libraries against tau aggregates conjugated to quantum dots, and this process eliminates the need for secondary nanobody screening. The isolated nanobodies demonstrate conformational specificity for tau aggregates in brain samples from both a transgenic mouse model and human tauopathies. We expect that our facile approach will be broadly useful for isolating conformational nanobodies against diverse amyloid aggregates and other complex antigens.

Metabolic changes and retinal remodeling in Heterozygous CRX mutant cats (CRXRDY/+).

CRX is a transcription factor essential for normal photoreceptor development and survival. The CRXRdy cat has a naturally occurring truncating mutation in CRX and is a large animal model for dominant Leber congenital amaurosis. This study investigated retinal remodeling that occurs as photoreceptors degenerate. CRXRdy/+ cats from 6 weeks to 10 years of age were investigated. In vivo structural changes of retinas were analyzed by fundus examination, confocal scanning laser ophthalmoscopy and spectral domain optical coherence tomography. Histologic analyses included immunohistochemistry for computational molecular phenotyping with macromolecules and small molecules. Affected cats had a cone-led photoreceptor degeneration starting in the area centralis. Initially there was preservation of inner retinal cells such as bipolar, amacrine and horizontal cells but with time migration of the deafferented neurons occurred. Early in the process of degeneration glial activation occurs ultimately resulting in formation of a glial seal. With progression the macula-equivalent area centralis developed severe atrophy including loss of retinal pigmentary epithelium. Microneuroma formation occured in advanced stages as more marked retinal remodeling occurred. This study indicates that retinal degeneration in the CrxRdy/+ cat retina follows the progressive, phased revision of retina that have been previously described for retinal remodeling. These findings suggest that therapy dependent on targeting inner retinal cells may be useful in young adults with preserved inner retinas prior to advanced stages of retinal remodeling and neuronal cell loss.

P62/SQSTM1 binds with claudin-2 to target for selective autophagy in stressed intestinal epithelium.

Impaired autophagy promotes Inflammatory Bowel Disease (IBD). Claudin-2 is upregulated in IBD however its role in the pathobiology remains uncertain due to its complex regulation, including by autophagy. Irrespective, claudin-2 expression protects mice from DSS colitis. This study was undertaken to examine if an interplay between autophagy and claudin-2 protects from colitis and associated epithelial injury. Crypt culture and intestinal epithelial cells (IECs) are subjected to stress, including starvation or DSS, the chemical that induces colitis in-vivo. Autophagy flux, cell survival, co-immunoprecipitation, proximity ligation assay, and gene mutational studies are performed. These studies reveal that under colitis/stress conditions, claudin-2 undergoes polyubiquitination and P62/SQSTM1-assisted degradation through autophagy. Inhibiting autophagy-mediated claudin-2 degradation promotes cell death and thus suggest that claudin-2 degradation promotes autophagy flux to promote cell survival. Overall, these data inform for the previously undescribed role for claudin-2 in facilitating IECs survival under stress conditions, which can be harnessed for therapeutic advantages.

Pan-cancer structurome reveals overrepresentation of beta sandwiches and underrepresentation of alpha helical domains.

The recent progress in the prediction of protein structures marked a historical milestone. AlphaFold predicted 200 million protein models with an accuracy comparable to experimental methods. Protein structures are widely used to understand evolution and to identify potential drug targets for the treatment of various diseases, including cancer. Thus, these recently predicted structures might convey previously unavailable information about cancer biology. Evolutionary classification of protein domains is challenging and different approaches exist. Recently our team presented a classification of domains from human protein models released by AlphaFold. Here we evaluated the pan-cancer structurome, domains from over and under expressed proteins in 21 cancer types, using the broadest levels of the ECOD classification: the architecture (A-groups) and possible homology (X-groups) levels. Our analysis reveals that AlphaFold has greatly increased the three-dimensional structural landscape for proteins that are differentially expressed in these 21 cancer types. We show that beta sandwich domains are significantly overrepresented and alpha helical domains are significantly underrepresented in the majority of cancer types. Our data suggest that the prevalence of the beta sandwiches is due to the high levels of immunoglobulins and immunoglobulin-like domains that arise during tumor development-related inflammation. On the other hand, proteins with exclusively alpha domains are important elements of homeostasis, apoptosis and transmembrane transport. Therefore cancer cells tend to reduce representation of these proteins to promote successful oncogeneses.

Kß X-ray Emission Spectroscopy of Cu(I)-Lytic Polysaccharide Monooxygenase: Direct Observation of the Frontier Molecular Orbital for H2O2 Activation.

Lytic polysaccharide monooxygenases (LPMOs) catalyze the degradation of recalcitrant carbohydrate polysaccharide substrates. These enzymes are characterized by a mononuclear Cu(I) active site with a three-coordinate T-shaped "His-brace" configuration including the N-terminal histidine and its amine group as ligands. This study explicitly investigates the electronic structure of the d10 Cu(I) active site in a LPMO using Kß X-ray emission spectroscopy (XES). The lack of inversion symmetry in the His-brace site enables the 3d/p mixing required for intensity in the Kß valence-to-core (VtC) XES spectrum of Cu(I)-LPMO. These Kß XES data are correlated to density functional theory (DFT) calculations to define the bonding, and in particular, the frontier molecular orbital (FMO) of the Cu(I) site. These experimentally validated DFT calculations are used to evaluate the reaction coordinate for homolytic cleavage of the H2O2 O-O bond and understand the contribution of this FMO to the low barrier of this reaction and how the geometric and electronic structure of the Cu(I)-LPMO site is activated for rapid reactivity with H2O2.

Ancient transcriptional regulators can easily evolve new pair-wise cooperativity.

Cells regulate gene expression by the specific binding of transcription regulators to cis-regulatory sequences. Pair-wise cooperativity between regulators-whereby two different regulators physically interact and bind DNA in a cooperative manner-is common and permits complex modes of gene regulation. Over evolutionary timescales, the formation of new combinations of regulators represents a major source of phenotypic novelty, facilitating new network structures. How functional, pair-wise cooperative interactions arise between regulators is poorly understood, despite the abundance of examples in extant species. Here, we explore a protein-protein interaction between two ancient transcriptional regulators-the homeodomain protein Matα2 and the MADS box protein Mcm1-that was gained approximately 200 million y ago in a clade of ascomycete yeasts that includes Saccharomyces cerevisiae. By combining deep mutational scanning with a functional selection for cooperative gene expression, we tested millions of possible alternative evolutionary solutions to this interaction interface. The artificially evolved, functional solutions are highly degenerate, with diverse amino acid chemistries permitted at all positions but with widespread epistasis limiting success. Nonetheless, approximately ~45% of the random sequences sampled function as well or better in controlling gene expression than the naturally evolved sequence. From these variants (which are unconstrained by historical contingency), we discern structural rules and epistatic constraints governing the emergence of cooperativity between these two transcriptional regulators. This work provides a mechanistic basis for long-standing observations of transcription network plasticity and highlights the importance of epistasis in the evolution of new protein-protein interactions.

Inhibiting NLRP3 signaling in aging podocytes improves their life- and health-span.

The decrease in the podocyte's lifespan and health-span that typify healthy kidney aging cause a decrease in their normal structure, physiology and function. The ability to halt and even reverse these changes becomes clinically relevant when disease is superimposed on an aged kidney. RNA-sequencing of podocytes from middle-aged mice showed an inflammatory phenotype with increases in the NLRP3 inflammasome, signaling for IL2/Stat5, IL6 and TNF, interferon gamma response, allograft rejection and complement, consistent with inflammaging. Furthermore, injury-induced NLRP3 signaling in podocytes was further augmented in aged mice compared to young ones. The NLRP3 inflammasome (NLRP3, Caspase-1, IL1ß IL-18) was also increased in podocytes of middle-aged humans. Higher transcript expression for NLRP3 in human glomeruli was accompanied by reduced podocyte density and increased global glomerulosclerosis and glomerular volume. Pharmacological inhibition of NLRP3 with MCC950, or gene deletion, reduced podocyte senescence and the genes typifying aging in middle-aged mice, which was accompanied by an improved podocyte lifespan and health-span. Moreover, modeling the injury-dependent increase in NLRP3 signaling in human kidney organoids confirmed the anti-senescence effect of MC9950. Finally, NLRP3 also impacted liver aging. Together, these results suggest a critical role for the NLRP3 inflammasome in podocyte and liver aging.

Plasma membrane association and resistosome formation of plant helper immune receptors.

Intracellular plant immune receptors, termed NLRs (Nucleotide-binding Leucine-rich repeat Receptors), confer effector-triggered immunity. Sensor NLRs are responsible for pathogen effector recognition. Helper NLRs function downstream of sensor NLRs to transduce signaling and induce cell death and immunity. Activation of sensor NLRs that contain TIR (Toll/interleukin-1receptor) domains generates small molecules that induce an association between a downstream heterodimer signalosome of EDS1 (EnhancedDisease Susceptibility 1)/SAG101 (Senescence-AssociatedGene 101) and the helper NLR of NRG1 (NRequired Gene 1). Autoactive NRG1s oligomerize and form calcium signaling channels largely localized at the plasma membrane (PM). The molecular mechanisms of helper NLR PM association and effector-induced NRG1 oligomerization are not well characterized. We demonstrate that helper NLRs require positively charged residues in their N-terminal domains for phospholipid binding and PM association before and after activation, despite oligomerization and conformational changes that accompany activation. We demonstrate that effector activation of a TIR-containing sensor NLR induces NRG1 oligomerization at the PM and that the cytoplasmic pool of EDS1/SAG101 is critical for cell death function. EDS1/SAG101 cannot be detected in the oligomerized NRG1 resistosome, suggesting that additional unknown triggers might be required to induce the dissociation of EDS1/SAG101 from the previously described NRG1/EDS1/SAG101 heterotrimer before subsequent NRG1 oligomerization. Alternatively, the conformational changes resulting from NRG1 oligomerization abrogate the interface for EDS1/SAG101 association. Our data provide observations regarding dynamic PM association during helper NLR activation and underpin an updated model for effector-induced NRG1 resistosome formation.

Elevated FBXW10 drives hepatocellular carcinoma tumorigenesis via AR-VRK2 phosphorylation-dependent GAPDH ubiquitination in male transgenic mice.

Hepatocellular carcinoma (HCC), the most common liver cancer, occurs mainly in men, but the underlying mechanism remains to be further explored. Here, we report that ubiquitinated glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is responsible for HCC tumorigenesis in males. Mechanistically, FBXW10 promotes GAPDH polyubiquitination and activation; VRK2-dependent phosphorylation of GAPDH Ser151 residue is critical for GAPDH ubiquitination and activation. Activated GAPDH interacts with TRAF2, leading to upregulation of the canonical and noncanonical NF-κB pathways, and increases PD-L1 and AR-VRK2 expression, followed by induction of immune evasion, HCC tumorigenesis, and metastasis. Notably, the GAPDH inhibitor koningic acid (KA) activates immune response and protects against FBXW10-driven HCC in vivo. In HCC clinical samples, the expression of active GAPDH is positively correlated with that of FBXW10 and VRK2. We propose that the FBXW10/AR/VRK2/GAPDH/NF-κB axis is critical for HCC tumorigenesis in males. Targeting this axis with KA is a potential therapeutic strategy for male HCC patients.

Targeting the Spike Receptor Binding Domain Class V Cryptic Epitope by an Antibody with Pan-Sarbecovirus Activity.

Novel therapeutic monoclonal antibodies (MAbs) must accommodate comprehensive breadth of activity against diverse sarbecoviruses and high neutralization potency to overcome emerging variants. Here, we report the crystal structure of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) receptor binding domain (RBD) in complex with MAb WRAIR-2063, a moderate-potency neutralizing antibody with exceptional sarbecovirus breadth, that targets the highly conserved cryptic class V epitope. This epitope overlaps substantially with the spike protein N-terminal domain (NTD) -interacting region and is exposed only when the spike is in the open conformation, with one or more RBDs accessible. WRAIR-2063 binds the RBD of SARS-CoV-2 WA-1, all variants of concern (VoCs), and clade 1 to 4 sarbecoviruses with high affinity, demonstrating the conservation of this epitope and potential resiliency against variation. We compare structural features of additional class V antibodies with their reported neutralization capacity to further explore the utility of the class V epitope as a pan-sarbecovirus vaccine and therapeutic target. IMPORTANCE Characterization of MAbs against SARS-CoV-2, elicited through vaccination or natural infection, has provided vital immunotherapeutic options for curbing the COVID-19 pandemic and has supplied critical insights into SARS-CoV-2 escape, transmissibility, and mechanisms of viral inactivation. Neutralizing MAbs that target the RBD but do not block ACE2 binding are of particular interest because the epitopes are well conserved within sarbecoviruses and MAbs targeting this area demonstrate cross-reactivity. The class V RBD-targeted MAbs localize to an invariant site of vulnerability, provide a range of neutralization potency, and exhibit considerable breadth against divergent sarbecoviruses, with implications for vaccine and therapeutic development.

Mapping roles of active site residues in the acceptor site of the PA3944 Gcn5-related N-acetyltransferase enzyme.

An increased understanding of how the acceptor site in Gcn5-related N-acetyltransferase (GNAT) enzymes recognizes various substrates provides important clues for GNAT functional annotation and their use as chemical tools. In this study, we explored how the PA3944 enzyme from Pseudomonas aeruginosa recognizes three different acceptor substrates, including aspartame, NANMO, and polymyxin B, and identified acceptor residues that are critical for substrate specificity. To achieve this, we performed a series of molecular docking simulations and tested methods to identify acceptor substrate binding modes that are catalytically relevant. We found that traditional selection of best docking poses by lowest S scores did not reveal acceptor substrate binding modes that were generally close enough to the donor for productive acetylation. Instead, sorting poses based on distance between the acceptor amine nitrogen atom and donor carbonyl carbon atom placed these acceptor substrates near residues that contribute to substrate specificity and catalysis. To assess whether these residues are indeed contributors to substrate specificity, we mutated seven amino acid residues to alanine and determined their kinetic parameters. We identified several residues that improved the apparent affinity and catalytic efficiency of PA3944, especially for NANMO and/or polymyxin B. Additionally, one mutant (R106A) exhibited substrate inhibition toward NANMO, and we propose scenarios for the cause of this inhibition based on additional substrate docking studies with R106A. Ultimately, we propose that this residue is a key gatekeeper between the acceptor and donor sites by restricting and orienting the acceptor substrate within the acceptor site.

Analysis of neutral mutational drift in an allosteric enzyme.

Neutral mutational drift is an important source of biological diversity that remains underexploited in fundamental studies of protein biophysics. This study uses a synthetic transcriptional circuit to study neutral drift in protein tyrosine phosphatase 1B (PTP1B), a mammalian signaling enzyme for which conformational changes are rate limiting. Kinetic assays of purified mutants indicate that catalytic activity, rather than thermodynamic stability, guides enrichment under neutral drift, where neutral or mildly activating mutations can mitigate the effects of deleterious ones. In general, mutants show a moderate activity-stability tradeoff, an indication that minor improvements in the activity of PTP1B do not require concomitant losses in its stability. Multiplexed sequencing of large mutant pools suggests that substitutions at allosterically influential sites are purged under biological selection, which enriches for mutations located outside of the active site. Findings indicate that the positional dependence of neutral mutations within drifting populations can reveal the presence of allosteric networks and illustrate an approach for using synthetic transcriptional systems to explore these mutations in regulatory enzymes.

Coupling of conformational dynamics and inhibitor binding in the phosphodiesterase-5 family.

Phosphodiesterase-5 (PDE5) is responsible for regulating the concentration of the second messenger molecule cGMP by hydrolyzing it into 5'-GMP. PDE5 is implicated in erectile dysfunction and cardiovascular diseases. The substrate binding site in the catalytic domain of PDE5 is surrounded by several dynamic structural motifs (including the α14 helix, M-loop, and H-loop) that are known to switch between inactive and active conformational states via currently unresolved structural intermediates. We evaluated the conformational dynamics of these structural motifs in the apo state and upon binding of an allosteric inhibitor (evodiamine) or avanafil, a competitive inhibitor. We employed enhanced sampling-based replica exchange solute scaling (REST2) method, principal component analysis (PCA), time-lagged independent component analysis (tICA), molecular dynamics (MD) simulations, and well-tempered metadynamics simulations to probe the conformational changes in these structural motifs. Our results support a regulatory mechanism for PDE5, where the α14 helix alternates between an inward (lower activity) conformation and an outward (higher activity) conformation that is accompanied by the folding/unfolding of the α8' and α8″ helices of the H-loop. When the allosteric inhibitor evodiamine is bound to PDE5, the inward (inactive) state of the α14 helix is preferred, thus preventing substrate access to the catalytic site. In contrast, competitive inhibitors of PDE5 block catalysis by occupying the active site accompanied by stabilization of the outward conformation of the α14 helix. Defining the conformational dynamics underlying regulation of PDE5 activation will be helpful in rational design of next-generation small molecules modulators of PDE5 activity.