Μs Molecular Dynamics Research Articles

Molecular Dynamics (MD) Simulations Provide Insights into the Activation Mechanisms of 5-HT2A Receptors

Recent breakthroughs in the determination of atomic resolution 3-D cryo-electron microscopy structures of membrane proteins present an unprecedented opportunity for drug discovery. Structure-based drug discovery utilizing in silico methods enables the study of dynamic connectivity of stable conformations induced by the drug in achieving its effect. With the ever-expanding computational power, simulations of this type reveal protein dynamics in the nano-, micro-, and even millisecond time scales. In the present study, aiming to characterize the protein dynamics of the 5HT2A receptor stimulated by ligands (agonist/antagonist), we performed 1 µs MD simulations on 5HT2A/DOI (agonist), 5HT2A/GSK215083 (antagonist), and 5HT2A (APO, no ligand) systems. The crystal structure of 5HT2A/zotepine (antagonist) (PDB: 6A94) was used to set up the simulation systems in a lipid bilayer environment. We found the monitoring of the ionic lock residue pair (R3.50-E6.30) of 5HT2A in MD simulations to be a good approximation of the effects of agonists (ionic lock breakage) or antagonists (ionic lock formation) on receptor activation. We further performed analyses of the MD trajectories, including Principal Component Analysis (PCA), hydrogen bond, salt bridge, and hydrophobic interaction network analyses, and correlation between residues to identify key elements of receptor activation. Our results suggest that in order to trigger receptor activation, DOI must interact with 5HT2A through residues V5.39, G5.42, S5.43, and S5.46 on TM5, inducing significant conformational changes in the backbone angles of G5.42 and S5.43. DOI also interacted with residues W6.48 (toggle switch) and F6.51 on TM6, causing major conformational shifts in the backbone angles of F6.44 and V6.45. These structural changes were transmitted to the intracellular ends of TM5, TM6, and ICL3, resulting in the breaking of the ionic lock and subsequent G protein activation. The studies could be helpful in future design of selective agonists/antagonists for various serotonin receptors (5HT1A, 5HT2A, 5HT2B, 5HT2C, and 5HT7) involved in detrimental disorders, such as addiction and schizophrenia.

ELUCIDATING THE DIFFERENTIAL IMPACTS OF EQUIVALENT GATING-CHARGE MUTATIONS IN VOLTAGE-GATED SODIUM CHANNELS.

Voltage-gated sodium (Nav) channels are pivotal for cellular signaling and mutations in Nav channels can lead to excitability disorders in cardiac, muscular, and neural tissues. A major cluster of pathological mutations localizes in the voltage-sensing domains (VSDs), resulting in either gain-of-function (GoF), loss-of-function (LoF) effects, or both. However, the mechanism behind this functional divergence of mutations at equivalent positions remains elusive. Through hotspot analysis, we identified three gating charges (R1, R2, and R3) as major mutational hotspots in VSDs. The same amino-acid substitutions at equivalent gating-charge positions in VSDI and VSDII of the cardiac sodium channel Nav1.5 show differential gating-property impacts in electrophysiology measurements. We conducted 120 μs molecular dynamics (MD) simulations on wild-type and six mutants to elucidate the structural basis of their differential impacts. Our μs-scale MD simulations with applied external electric fields captured VSD state transitions and revealed the differential structural dynamics between equivalent R-to-Q mutants. Notably, we observed transient leaky conformations in some mutants during structural transitions, offering a detailed structural explanation for gating-pore currents. Our salt-bridge network analysis uncovered VSD-specific and state-dependent interactions among gating charges, countercharges, and lipids. This detailed analysis elucidated how mutations disrupt critical electrostatic interactions, thereby altering VSD permeability and modulating gating properties. By demonstrating the crucial importance of considering the specific structural context of each mutation, our study represents a significant leap forward in understanding structure-function relationships in Nav channels. Our work establishes a robust framework for future investigations into the molecular basis of ion channel-related disorders.

Identifying Key Binding Interactions Between the Cardiac L-Type Calcium Channel and Calmodulin Using Molecular Dynamics Simulations.

Defects in the binding of the calcium sensing protein calmodulin (CaM) to the L-type calcium channel (CaV1.2) or to the ryanodine receptor type 2 (RyR2) can lead to dangerous cardiac arrhythmias with distinct phenotypes, such as long-QT syndrome (LQTS) and catecholaminergic ventricular tachycardia (CPVT). Certain CaM mutations lead to LQTS while other mutations lead to CPVT, but the mechanisms by which a specific mutation can lead to each disease phenotype are not well-understood. In this study, we use long, 2 μs molecular dynamics simulations and a multitrajectory approach to identify the key binding interactions between the IQ domain of CaV1.2 and CaM. Five key interactions are found between CaV1.2 and CaM in the C-lobe, 1 in the central linker, and 2 in the N-lobe. In addition, while 5 key interactions appear between residues 120-149 in the C-lobe of CaM when it interacts with CaV1.2, only 1 key interaction is found within this region of CaM when it interacts with the RyR2. We show that this difference in the distribution of key interactions correlates with the known distribution of CaM mutations that lead to LQTS or CPVT. This correlation suggests that a disruption of key binding interactions is a plausible mechanism that can lead to these two different disease phenotypes.

Pharmacophore Virtual Screening Identifies Riboflavin as an Inhibitor of the Schistosome Cathepsin B1 Protease with Antiparasitic Activity.

Schistosomiasis is a neglected disease of poverty that affects over 200 million people worldwide and relies on a single drug for therapy. The cathepsin B1 cysteine protease (SmCB1) of Schistosoma mansoni has been investigated as a potential target. Here, a structure-based pharmacophore virtual screening (VS) approach was used on a data set of approved drugs to identify potential antischistosomal agents targeting SmCB1. Pharmacophore (PHP) models underwent validation through receiver operating characteristics curves achieving values >0.8. The data highlighted riboflavin (RBF) as a compound of particular interest. A 1 μs molecular dynamics simulation demonstrated that RBF altered the conformation of SmCB1, causing the protease's binding site to close around RBF while maintaining the protease's overall integrity. RBF inhibited the activity of SmCB1 at low micromolar values and killed the parasite in vitro. Finally, in a murine model of S. mansoni infection, oral administration of 100 mg/kg RBF for 7 days significantly decreased worm burdens by ∼20% and had a major impact on intestinal and fecal egg burdens, which were decreased by ∼80%.

Adaptive Workflows of Machine Learning Illuminate the Sequential Operation Mechanism of the TAK1's Allosteric Network.

Allostery is a fundamental mechanism driving biomolecular processes that holds significant therapeutic concern. Our study rigorously investigates how two distinct machine-learning algorithms uniquely classify two already close-to-active DFG-in states of TAK1, differing just by the presence or absence of its allosteric activator TAB1, from an ensemble mixture of conformations (obtained from 2.4 μs molecular dynamics (MD) simulations). The novelty, however, lies in understanding the deeper algorithmic potentials to systematically derive a diverse set of differential residue connectivity features that reconstruct the essential mechanistic architecture for TAK1-TAB1 allostery in such a close-to-active biochemical scenario. While the recursive, random forest-based workflow displays the potential of conducting discretized, hierarchical derivation of allosteric features, a multilayer perceptron-based approach gains considerable efficacy in revealing fluid connected patterns of features when hybridized with mutual information scoring. Interestingly, both pipelines benchmark similar directions of functional conformational changes for TAK1's activation. The findings significantly advance the depth of mechanistic understanding by highlighting crucial activation signatures along a directed C-lobe → activation loop → ATP pocket channel of information flow, including (1) the αF-αE biterminal alignments and (2) the "catalytic" drift of the activation loop toward kinase active site. Besides, some novel allosteric hotspots (K253, Y206, N189, etc.) are further recognized as TAB1 sensors, transducers, and responders, including a benchmark E70 mutation site, precisely mapping the important structural segments for sequential allosteric execution. Hence, our work demonstrates how to navigate through greater structural depths and dimensions of dynamic allosteric machineries just by leveraging standard ML methods in suitable streamlined workflows adaptive to the specific system and objectives.

Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: A Symmetry-Based Perspective.

We report on a new method for the characterization of local structures in proteins based on extensive molecular dynamics (MD) simulations, here, 1 μs in length. The N-H bond of the Rho GTPase binding domain of plexin-B1 (RBD) serves as a probe and the potential, u(MD), which restricts its internal motion, as a qualifier of the local dynamic structure. u(MD) is derived from the MD trajectory as a function of the polar angles, (θ, φ), which specify the N-H orientation in the protein. u(MD) is statistical in character yielding empirical descriptions. To establish more insightful methodical descriptions, we develop a comprehensive method which approximates u(MD) by combinations of analytical Wigner functions that belong to the D2h point group. These combinations, called u(simulated), make it possible to gain a new perspective of local dynamic structures in proteins based on explicit potentials/free energy surfaces and associated probability densities, entropy, and ordering. A simpler method was developed previously using 100 ns MD simulations. In that case, the traditional "perpendicular N-H ordering" setting centered at Cα-Cα with (θ, φ) = (90, 90) and generally, featuring positive φ, prevailed. u(MD) derived from 1 μs MD simulations is considerably more complex requiring substantial model enhancement. The enhanced method applies to the well-structured sections of the RBD. It only applies partly to its loops where u(MD) extends into the negative-φ region where we detect nonperpendicular N-H ordering. This arrangement requires devising new reference structures and making substantial algorithmic changes, to be performed in future work. Here, we focus on developing the comprehensive method and using it to investigate perpendicular ordering settings. We find that secondary structures (loops) exhibit varying (virtually invariant) potentials with Ag, B2u, and B1u (Ag and B2u) D2h symmetry. Application to RBD dimerization and RBD binding to the GTPase Rac1 is described in the subsequent article. Applications to other probes, proteins, and biological functions, based on explicit local potentials, probability densities, entropy, and ordering, are possible.

Understanding the effects of ethanol on the liposome bilayer structure using microfluidic-based time-resolved small-angle X-ray scattering and molecular dynamics simulations.

Lipid nanoparticles (LNPs) are essential carrier particles in drug delivery systems, particularly in ribonucleic acid delivery. In preparing lipid-based nanoparticles, microfluidic-based ethanol injection may produce precisely size-controlled nanoparticles. Ethanol is critical in LNP formation and post-treatment processes and affects liposome size, structure, lamellarity, and drug-loading efficiency. However, the effects of time-dependent changes in the ethanol concentration on the structural dynamics of liposomes are not clearly understood. Herein, we investigated ethanol-induced lipid bilayer changes in liposomes on a time scale from microseconds to tens of seconds using a microfluidic-based small-angle X-ray scattering (SAXS) measurement system coupled with molecular dynamics (MD) simulations. The time-resolved SAXS measurement system revealed that single unilamellar liposomes were converted to multilamellar liposomes within 0.8 s of contact with ethanol, and the d-spacing was decreased from 6.1 (w/o ethanol) to 4.4 nm (80% ethanol) with increasing ethanol concentration. We conducted 1 μs MD simulations to understand the molecular-level structural changes in the liposomes. The MD simulations revealed that the changes in the lamellar structure caused by ethanol at the molecular level could explain the structural changes in the liposomes observed via time-resolved SAXS. Therefore, the post-treatment process to remove residual ethanol is critical in liposome formation.

Development of 2-Aminotetralin-Type Serotonin 5-HT1 Agonists: Molecular Determinants for Selective Binding and Signaling at 5-HT1A, 5-HT1B, 5-HT1D, and 5-HT1F Receptors.

The serotonin (5-hydroxytryptamine, 5-HT) 5-HT1 G-protein coupled receptor subtypes (5-HT1A/1B/1D/1E/1F) share a high sequence homology, confounding development of subtype-specific ligands. This study used a 5-HT1 structure-based ligand design approach to develop subtype-selective ligands using a 5-substituted-2-aminotetralin (5-SAT) chemotype, leveraging results from pharmacological, molecular modeling, and mutagenesis studies to delineate molecular determinants for 5-SAT binding and function at 5-HT1 subtypes. 5-SATs demonstrated high affinity (Ki ≤ 25 nM) and at least 50-fold stereoselective preference ([2S] > [2R]) at 5-HT1A, 5-HT1B, and 5-HT1D receptors but essentially nil affinity (Ki > 1 μM) at 5-HT1F receptors. The 5-SATs tested were agonists with varying degrees of potency and efficacy, depending on chemotype substitution and 5-HT1 receptor subtype. Models were built from the 5-HT1A (cryo-EM), 5-HT1B (crystal), and 5-HT1D (cryo-EM) structures, and 5-SATs underwent docking studies with up to 1 μs molecular dynamics simulations. 5-SAT interactions observed at positions 3.33, 5.38, 5.42, 5.43, and 7.39 of 5-HT1 subtypes were confirmed with point mutation experiments. Additional 5-SATs were designed and synthesized to exploit experimental and computational results, yielding a new full efficacy 5-HT1A agonist with 100-fold selectivity over 5-HT1B/1D receptors. The results presented lay the foundation for the development of additional 5-HT1 subtype selective ligands for drug discovery purposes.

Effect of Guanine Adduct Size, Shape, and Linker Type on the Conformation of Adducted DNA: A DFT and Molecular Dynamics Study.

DNA is damaged through various exogenous sources (e.g., automobile exhaust, tobacco smoke, and processed foods), which can yield diverse C8-dG bulky aryl adducts. Adducts are known to induce structural changes to DNA that can lead to various biological outcomes, ranging from cell death to diseases such as cancer. Unfortunately, the relationship between the chemical composition of the damaged product, the adducted DNA structure, and the biological consequences is not well understood, which limits the development of disease detection and prevention strategies. The present study uses density functional theory (DFT) calculations and quintuplicate 1 μs molecular dynamics (MD) simulations to characterize the structure of DNA containing 21 model C8-dG adducts that systematically differ in size (phenyl to pyrenyl), shape (α (2,3), β (3,4) fusion, or ring substitution), and nucleobase-aryl group linkage (N, O, and C-linked). DFT calculations reveal that the inherent structural features of the G nucleobase adducts are impacted by linker type and bulky moiety shape, but not size, with the conformational flexibility reducing with α-ring fusion and linker composition as N > O > C. These structural properties are maintained in nucleoside models, which also reveal an increased propensity for anti-to-syn rotation about the glycosidic bond with N < O < C linker type. Although these diverse chemical features do not influence the global structure of adducted DNA, the adducts differentially impact the conformation local to the adducted site, including the relative populations of structures with the bulky moiety in the major groove (B conformer) and intercalated (stacked) into the helix (S conformer). Specifically, while the smallest phenyl adducts favor the B conformation and the largest pyrenyl-derived adducts stabilize the S conformation, the B/S ratio decreases with an increase in ring size and N > O > C linker composition. The shape and size (length) of the adduct can further finetune the B/S ratio, with β-fused naphthyl or α-fused phenanthryl N-linked adducts and O or C-linked adducts containing ring substitution increasing the prevalence of the S adducted DNA conformation. Overall, this work uncovers the significant effect of bulky moiety size and linker type, as well as the lesser impact of aryl group shape, on adducted DNA structure, which suggests differential replication and repair outcomes, and thereby represents an important step toward rationalizing connections between the structure and biological consequences of diverse DNA adducts.

Robust AMBER Force Field Parameters for Glutathionylated Cysteines.

S-glutathionylation is an oxidative post-translational modification, which is involved in the regulation of many cell signaling pathways. Increasing amounts of studies show that it is crucial in cell homeostasis and deregulated in several pathologies. However, the effect of S-glutathionylation on proteins' structure and activity is poorly understood, and a drastic lack of structural information at the atomic scale remains. Studies based on the use of molecular dynamics simulations, which can provide important information about modification-induced modulation of proteins' structure and function, are also sparse, and there is no benchmarked force field parameters for this modified cysteine. In this contribution, we provide robust AMBER parameters for S-glutathionylation, which we tested extensively against experimental data through a total of 33 μs molecular dynamics simulations. We show that our parameter set efficiently describes the global and local structural properties of S-glutathionylated proteins. These data provide the community with an important tool to foster new investigations into the effect of S-glutathionylation on protein dynamics and function, in a common effort to unravel the structural mechanisms underlying its critical role in cellular processes.

Microsecond dynamics of H10N7 influenza neuraminidase reveals the plasticity of loop regions and drug resistance due to the R292K mutation.

At the beginning of the last century, multiple pandemics caused by influenza (flu) viruses severely impacted public health. Despite the development of vaccinations and antiviral medications to prevent and control impending flu outbreaks, unforeseen novel strains and continuously evolving old strains continue to represent a serious threat to human life. Therefore, the recently identified H10N7, for which not much data is available for rational structure-based drug design, needs to be further explored. Here, we investigated the structural dynamics of neuraminidase N7 upon binding of inhibitors, and the drug resistance mechanisms against the oseltamivir (OTV) and laninamivir (LNV) antivirals due to the crucial R292K mutation on the N7 using the computational microscope, molecular dynamics (MD) simulations. In this study, each system underwent long 2 × 1 μs MD simulations to answer the conformational changes and drug resistance mechanisms. These long time-scale dynamics simulations and free energy landscapes demonstrated that the mutant systems showed a high degree of conformational variation compared to their wildtype (WT) counterparts, and the LNV-bound mutant exhibited an extended 150-loop conformation. Further, the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculation and MM/GBSA free energy decomposition were used to characterize the binding of OTV and LNV with WT, and R292K mutated N7, revealing the R292K mutation as drug-resistant, facilitated by a decline in binding interaction and a reduction in the dehydration penalty. Due to the broader binding pocket cavity of the smaller K292 mutant residue relative to the wildtype, the drug carboxylate to K292 hydrogen bonding was lost, and the area surrounding the K292 residue was more accessible to water molecules. This implies that drug resistance could be reduced by strengthening the hydrogen bond contacts between N7 inhibitors and altered N7, creating inhibitors that can form a hydrogen bond to the mutant K292, or preserving the closed cavity conformations.

Evaluating the Cysteine-Rich and Catalytic Subdomains of Human Tyrosinase and OCA1-Related Mutants Using 1 μs Molecular Dynamics Simulation.

The inherited disorder oculocutaneous albinism type 1 (OCA1) is caused by mutations in the TYR gene encoding tyrosinase (Tyr), an enzyme essential to producing pigments throughout the human body. The intramelanosomal domain of Tyr consists of the cysteine-rich and tyrosinase catalytic subdomains, which are essential for enzymatic activity. In protein unfolding, the roles of these subdomains are not well established. Here, we performed six molecular dynamics simulations at room temperature for Tyr and OCA1-related mutant variants P406L and R402Q intramelanosomal domains. The proteins were simulated for 1 μs in water and urea to induce unfolding. In urea, we observed increases in surface area, decreases in intramolecular hydrogen bonding, and decreases in hydrophobic interactions, suggesting a 'molten globule' state for each protein. Between all conditions, the cysteine-rich subdomain remains stable, whereas the catalytic subdomain shows increased flexibility. This flexibility is intensified by the P406L mutation, while R402Q increases the catalytic domain's rigidity. The cysteine-rich subdomain is rigid, preventing the protein from unfolding, whereas the flexibility of the catalytic subdomain accommodates mutational changes that could inhibit activity. These findings match the conclusions from our experimental work suggesting the function alteration by the P406L mutation, and the potential role of R402Q as a polymorphism.

Drug Binding to BamA Targets Its Lateral Gate.

BamA, the core component of the β-barrel assembly machinery (BAM) complex, is an outer-membrane protein (OMP) in Gram-negative bacteria. Its function is to insert and fold substrate OMPs into the outer membrane (OM). Evidence suggests that BamA follows the asymmetric hybrid-barrel model where the first and last strands of BamA separate, a process known as lateral gate opening, to allow nascent substrate OMP β-strands to sequentially insert and fold through β-augmentation. Recently, multiple lead compounds that interfere with BamA's function have been identified. We modeled and then docked one of these compounds into either the extracellular loops of BamA or the open lateral gate. With the compound docked in the loops, we found that the lateral gate remains closed during 5 μs molecular dynamics simulations. The same compound when docked in the open lateral gate stays bound to the β16 strand of BamA during the simulation, which would prevent substrate OMP folding. In addition, we simulated mutants of BamA that are resistant to one or more of the identified lead compounds. In these simulations, we observed a differing degree and/or frequency of opening of BamA's lateral gate compared to BamA-apo, suggesting that the mutations grant resistance by altering the dynamics at the gate. We conclude that the compounds act by inhibiting BamA lateral gate opening and/or binding of substrate, thus preventing subsequent OMP folding and insertion.

Computation of Absolute Binding Free Energies for Noncovalent Inhibitors with SARS-CoV-2 Main Protease.

Accurate, routine calculation of absolute binding free energies (ABFEs) for protein-ligand complexes remains a key goal of computer-aided drug design since it can enable screening and optimization of drug candidates. For development and testing of related methods, it is important to have high-quality datasets. To this end, from our own experimental studies, we have selected a set of 16 inhibitors of the SARS-CoV-2 main protease (Mpro) with structural diversity and well-distributed BFEs covering a 5 kcal/mol range. There is also minimal structural uncertainty since X-ray crystal structures have been deposited for 12 of the compounds. For methods testing, we report ABFE results from 2 μs molecular dynamics (MD) simulations using free energy perturbation (FEP) theory. The correlation of experimental and computed results is encouraging, with a Pearson's r2 of 0.58 and a Kendall τ of 0.24. The results indicate that current FEP-based ABFE calculations can be used for identification of active compounds (hits). While their accuracy for lead optimization is not yet sufficient, this activity remains addressable in separate lead series by relative BFE calculations.

Intrinsic proclivity of left-handed conformation in large Nest motif peptides inferred from molecular dynamics

The ‘Nest’ motif plays a functional role in protein owing to its ligand binding potential aided by geometric concavity. The presence of less favored left-handed conformation (L-state) in its structure makes this concavity possible and in shaping the native chemical environment amenable to stable binding interactions. To understand the persistent appearance of L-state torsion in the Nest motif, we analyzed 0.5μs Molecular Dynamics (MD) simulation trajectories of 35 six-residue peptides (out of a total of 50 large Nest sequences of ≥6 residues) identified in our previous study. Analysis of the MD trajectories of the individual peptides reveals initial L-state in 60% of the peptides persists for >40% of the trajectory. Further, Nests with different sequences appear to adopt a specific conformational state driven by the neighboring L-state residues. The sequences also possess short secondary structures and amino acid repeats, suggesting evolutionary conservation and the specific role of amino acids in locally predisposing the torsion angle to the L-state. These findings help us to understand how L-state conformation is an essential prerequisite in stabilizing the Nest motif and shed light on the sequence-structure-function paradigm in the rational design of peptides and peptidomimetics for therapeutics. Communicated by Ramaswamy H. Sarma

Molecular insights in repurposing selective COX-2 inhibitor celecoxib against matrix metalloproteinases in potentiating delayed wound healing: a molecular docking and MMPB/SA based analysis of molecular dynamic simulations

Matrix metalloproteinases (MMPs) are proteolytic enzymes that play a role in healing, including reducing inflammation, promoting fibroblast and keratinocyte migration, and modifying scar tissue. Due to their pleiotropic functions in the wound-healing process in diabetic wounds, MMPs constitute a significant cause of delayed wound closure. COX-2 inhibitors are proven to inhibit inflammation. The present study aims to repurpose celecoxib against MMP-2, MMP-8 and MMP-9 through in silico approaches, such as molecular docking, molecular dynamics, and MMPB/SA analysis. We considered five selective COX-2 inhibitors (celecoxib, etoricoxib, lumiracoxib, rofecoxib and valdecoxib) for our study against MMPs. Based on molecular docking study and hydrogen bonding pattern, celecoxib in complex with three MMPs was further analyzed using 1 µs (1000 ns) molecular dynamics simulation and MMPB/SA techniques. These studies identified that celecoxib exhibited significant binding affinity −8.8, −7.9 and −8.3 kcal/mol, respectively, against MMP-2, MMP-8 and MMP-9. Celecoxib formed hydrogen bonding and hydrophobic (π–π) interactions with crucial substrate pocket amino acids, which may be accountable for their inhibitory nature. The MMPB/SA studies showed that electrostatic and van der Waal energy terms favoured the total free binding energy component, while polar solvation terms were highly disfavored. The in silico analysis of the secondary structures showed that the celecoxib binding conformation maintains relatively stable along the simulation trajectories. These findings provide some key clues regarding the accommodation of celecoxib in the substrate binding S1’ pocket and also provide structural insights and challenges in repurposing drugs as new MMP inhibitors with anti-inflammatory and anti-inflammatory wound-healing properties. Communicated by Ramaswamy H. Sarma

Interaction of the Inhibitory Peptides ShK and HmK with the Voltage-Gated Potassium Channel KV1.3: Role of Conformational Dynamics.

Peptide toxins that adopt the ShK fold can inhibit the voltage-gated potassium channel KV1.3 with IC50 values in the pM range and are therefore potential leads for drugs targeting autoimmune and neuroinflammatory diseases. Nuclear magnetic resonance (NMR) relaxation measurements and pressure-dependent NMR have shown that, despite being cross-linked by disulfide bonds, ShK itself is flexible in solution. This flexibility affects the local structure around the pharmacophore for the KV1.3 channel blockade and, in particular, the relative orientation of the key Lys and Tyr side chains (Lys22 and Tyr23 in ShK) and has implications for the design of KV1.3 inhibitors. In this study, we have performed molecular dynamics (MD) simulations on ShK and a close homologue, HmK, to probe the conformational space occupied by the Lys and Tyr residues, and docked the different conformations with a recently determined cryo-EM structure of the KV1.3 channel. Although ShK and HmK have 60% sequence identity, their dynamic behaviors are quite different, with ShK sampling a broad range of conformations over the course of a 5 μs MD simulation, while HmK is relatively rigid. We also investigated the importance of conformational dynamics, in particular the distance between the side chains of the key dyad Lys22 and Tyr23, for binding to KV1.3. Although these peptides have quite different dynamics, the dyad in both adopts a similar configuration upon binding, revealing a conformational selection upon binding to KV1.3 in the case of ShK. Both peptides bind to KV1.3 with Lys22 occupying the pore of the channel. Intriguingly, the more flexible peptide, ShK, binds with significantly higher affinity than HmK.

Critical Assessment of Self-Consistency Checks in the All-Atom Molecular Dynamics Simulation of Intrinsically Disordered Proteins.

All atom simulations can be used to quantify conformational properties of Intrinsically Disordered Proteins (IDP). However, simulations must satisfy convergence checks to ensure observables computed from simulation are reliable and reproducible. While absolute convergence is purely a theoretical concept requiring infinitely long simulation, a more practical, yet rigorous, approach is to impose Self Consistency Checks (SCCs) to gain confidence in the simulated data. Currently there is no study of SCCs in IDPs, unlike their folded counterparts. In this paper, we introduce different criteria for self-consistency checks for IDPs. Next, we impose these SCCs to critically assess the performance of different simulation protocols using the N terminal domain of HIV Integrase and the linker region of SARS-CoV-2 Nucleoprotein as two model IDPs. All simulation protocols begin with all-atom implicit solvent Monte Carlo (MC) simulation and subsequent clustering of MC generated conformations to create the representative structures of the IDPs. These representative structures serve as the initial structure for subsequent molecular dynamics (MD) runs with explicit solvent. We conclude that generating multiple short (∼3 μs) MD simulation trajectories─all starting from the most representative MC generated conformation─and merging them is the protocol of choice due to (i) its ability to satisfy multiple SCCs, (ii) consistently reproducing experimental data, and (iii) the efficiency of running independent trajectories in parallel by harnessing multiple cores available in modern GPU clusters. Running one long trajectory (greater than 20 μs) can also satisfy the first two criteria but is less desirable due to prohibitive computation time. These findings help resolve the challenge of identifying a usable starting configuration, provide an objective measure of SCC, and establish rigorous criteria to determine the minimum length (for one long simulation) or number of trajectories needed in all-atom simulation of IDPs.

Conformationally Selective 2-Aminotetralin Ligands Targeting the alpha2A- and alpha2C-Adrenergic Receptors.

Many important physiological processes are mediated by alpha2A- and alpha2C-adrenergic receptors (α2Rs), a subtype of class A G protein-coupled receptors (GPCRs). However, α2R signaling is poorly understood, and there are few approved medications targeting these receptors. Drug discovery aimed at α2Rs is complicated by the high degree of binding pocket homology between α2AR and α2CR, which confounds ligand-mediated selective activation or inactivation of signaling associated with a particular subtype. Meanwhile, α2R signaling is complex and it is reported that activating α2AR is beneficial in many clinical contexts, while activating α2CR signaling may be detrimental to these positive effects. Here, we report on a novel 5-substituted-2-aminotetralin (5-SAT) chemotype that, depending on substitution, has diverse pharmacological activities at α2Rs. Certain lead 5-SAT analogues act as partial agonists at α2ARs, while functioning as inverse agonists at α2CRs, a novel pharmacological profile. Leads demonstrate high potency (e.g., EC50 < 2 nM) at the α2AR and α2CRs regarding Gαi-mediated inhibition of adenylyl cyclase and production of cyclic adenosine monophosphate (cAMP). To help understand the molecular basis of 5-SAT α2R multifaceted functional activity, α2AR and α2CR molecular models were built from the crystal structures and 1 μs molecular dynamics (MD) simulations and molecular docking experiments were performed for a lead 5-SAT with α2AR agonist and α2CR inverse agonist activity, i.e., (2S)-5-(2'-fluorophenyl)-N,N-dimethyl-1,2,3,4-tetrahydronaphthalen-2-amine (FPT), in comparison to the FDA-approved (for opioid withdrawal symptoms) α2AR/α2CR agonist lofexidine. Results reveal several interactions between FPT and α2AR and α2CR amino acids that may impact the functional activity. The computational data in conjunction with experimental in vitro affinity and function results provide information to understand ligand stabilization of functionally distinct GPCR conformations regarding α2AR and α2CRs.

Structural insights into thrombolytic activity of destabilase from medicinal leech

Destabilase from the medical leech Hirudo medicinalis belongs to the family of i-type lysozymes. It has two different enzymatic activities: microbial cell walls destruction (muramidase activity), and dissolution of the stabilized fibrin (isopeptidase activity). Both activities are known to be inhibited by sodium chloride at near physiological concentrations, but the structural basis remains unknown. Here we present two crystal structures of destabilase, including a 1.1 Å-resolution structure in complex with sodium ion. Our structures reveal the location of sodium ion between Glu34/Asp46 residues, which were previously recognized as a glycosidase active site. While sodium coordination with these amino acids may explain inhibition of the muramidase activity, its influence on previously suggested Ser49/Lys58 isopeptidase activity dyad is unclear. We revise the Ser49/Lys58 hypothesis and compare sequences of i-type lysozymes with confirmed destabilase activity. We suggest that the general base for the isopeptidase activity is His112 rather than Lys58. pKa calculations of these amino acids, assessed through the 1 μs molecular dynamics simulation, confirm the hypothesis. Our findings highlight the ambiguity of destabilase catalytic residues identification and build foundations for further research of structure–activity relationship of isopeptidase activity as well as structure-based protein design for potential anticoagulant drug development.