CHARMM-GUI Enhanced Sampler for various collective variables and enhanced sampling methods

Nav1.2 channel mutations preventing fast inactivation lead to SCN2A encephalopathy

SCN2A gene-related early-infantile developmental and epileptic encephalopathy (EI-DEE) is a rare and severe disorder that manifests in early infancy. SCN2A mutations affecting the fast inactivation gating mechanism can result in altered voltage dependence and incomplete inactivation of the encoded neuronal Nav1.2 channel and lead to abnormal neuronal excitability. In this study, we evaluated clinical data of seven missense Nav1.2 variants associated with DEE and performed molecular dynamics simulations, patch-clamp electrophysiology and dynamic clamp real-time neuronal modelling to elucidate the molecular and neuron-scale phenotypic consequences of the mutations. The N1662D mutation almost completely prevented fast inactivation without affecting activation. The comparison of wild-type and N1662D channel structures suggested that the ambifunctional hydrogen bond formation between residues N1662 and Q1494 is essential for fast inactivation. Fast inactivation could also be prevented with engineered Q1494A or Q1494L Nav1.2 channel variants, whereas Q1494E or Q149K variants resulted in incomplete inactivation and persistent current. Molecular dynamics simulations revealed a reduced affinity of the hydrophobic IFM-motif to its receptor site with N1662D and Q1494L variants relative to wild-type. These results demonstrate that the interactions between N1662 and Q1494 underpin the stability and the orientation of the inactivation gate and are essential for the development of fast inactivation. Six DEE-associated Nav1.2 variants, with mutations mapped to channel segments known to be implicated in fast inactivation were also evaluated. Remarkably, the L1657P variant also prevented fast inactivation and produced biophysical characteristics that were similar to those of N1662D, whereas the M1501V, M1501T, F1651C, P1658S and A1659V variants resulted in biophysical properties that were consistent with gain-of-function and enhanced action potential firing of hybrid neurons in dynamic action potential clamp experiments. Paradoxically, low densities of N1662D or L1657P currents potentiated action potential firing, whereas increased densities resulted in sustained depolarization. Our results provide novel structural insights into the molecular mechanism of Nav1.2 channel fast inactivation and inform treatment strategies for SCN2A-related EI-DEE. The contribution of non-inactivating Nav1.2 channels to neuronal excitability may constitute a distinct cellular mechanism in the pathogenesis of SCN2A-related DEE.

Elucidating Protein Structures in the Gas Phase: Traversing Configuration Space with Biasing Methods

Achieving accurate characterization of protein structures in the gas phase continues to be a formidable challenge. To tackle this issue, the present study employs Molecular Dynamics (MD) simulations in tandem with enhanced sampling techniques (methods designed to efficiently explore protein conformations). The objective is to identify suitable structures of proteins by contrasting their calculated Collision Cross-Section (CCS) with those observed experimentally. Significant discrepancies were observed between the initial MD-simulated and experimentally measured CCS values through Ion Mobility-Mass Spectrometry (IMS-MS). To bridge this gap, we employed two distinct enhanced sampling methods, Harmonic Biasing Potential and Adaptive Biasing Force, which help the proteins overcome energy barriers to adopt more compact configurations. These techniques leverage the radius of gyration as a reaction coordinate (guiding parameter), guiding the system toward compressed states that potentially match experimental configurations more closely. The guiding forces are only employed to overcome existing barriers and are removed to allow the protein to naturally arrive at a potential gas phase configuration. The results demonstrated close alignment (within ∼4%) between simulated and experimental CCS values despite using different strengths and/or methods, validating their efficacy. This work lays the groundwork for future studies aimed at optimizing biasing methods and expanding the collective variables used for more accurate gas-phase structural predictions.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

The Putative Role of TIM-3 Variants in Polyendocrine Autoimmunity: Insights from a WES Investigation

Autoimmune polyglandular syndrome (APS) comprises a complex association of autoimmune pathological conditions. APS Type 1 originates from loss-of-function mutations in the autoimmune regulator (AIRE) gene. APS2, APS3 and APS4 are linked to specific HLA alleles within the major histocompatibility complex, with single-nucleotide polymorphisms (SNPs) in non-HLA genes also contributing to disease. In general, variability in the AIRE locus and the presence of heterozygous loss-of-function mutations can impact self-antigen presentation in the thymus. In this study, whole-exome sequencing (WES) was performed on a sixteen-year-old female APS3A/B patient to investigate the genetic basis of her complex phenotype. The analysis identified two variants (p.Arg111Trp and p.Thr101Ile) of the hepatitis A virus cell receptor 2 gene (HAVCR2) encoding for the TIM-3 (T cell immunoglobulin and mucin domain 3) protein. These variants were predicted, through in silico analysis, to impact protein structure and stability, potentially influencing the patient's autoimmune phenotype. While confocal microscopy analysis revealed no alteration in TIM-3 fluorescence intensity between the PBMCs isolated from the patient and those of a healthy donor, RT-qPCR showed reduced TIM-3 expression in the patient's unfractionated PBMCs. A screening conducted on a cohort of thirty APS patients indicated that the p.Thr101Ile and p.Arg111Trp mutations were unique to the proband. This study opens the pathway for the search of TIM-3 variants possibly linked to complex autoimmune phenotypes, highlighting the potential of novel variant discovery in contributing to APS classification and diagnosis.

Structural insights into the molecular recognition of integrin αVβ3 by RGD-containing ligands: The role of the specificity-determining loop (SDL)

Integrin αVβ3 is a prominent member of the "RGD-recognizing" integrin family of cell surface receptors. αVβ3 binds to various extracellular matrix (ECM) proteins and oxysterols such as 25-hydroxycholesterol, is implicated in several diseases, including cancer metastasis, lung fibrosis, inflammation, and autoimmune diseases, and is pursued as a valuable therapeutic target. Despite enormous efforts to seek a pure antagonist, to date, no single drug candidate has successfully reached clinics due to associated partial agonism and toxicity issues. Developing effective and safe inhibitors require a thorough understanding of the molecular interactions and structural changes related to the receptor's activation and inhibition mechanisms. This study offers a comprehensive residue-residue contact and network analyses of the ligand-binding β-propeller βI domains (headpiece) based on all available experimental structures of integrin αVβ3 in unliganded, agonist-, antagonist-, and antibody-bound states. The analyses reveal many critical interactions that were not reported before and show that specific orientation and interactions of residues from the specificity-determining loop (SDL) are critical in molecular recognition and regulation. Also, the network analysis reveals that residues from the nearby allosteric site (site II) connect to the primary RGD-binding site via SDL, which likely acts as an interface between the two sites. Our results provide valuable insights into molecular interactions, structural changes, distinct features of the active and inactive headpiece conformations, the role of SDL in ligand recognition, and SDL-mediated allostery. Thus, the insights from this study may facilitate the designing of pure antagonists or site II-mediated allosteric modulators to integrin αVβ3 to treat various diseases.

Antiviral drug recognition and elevator-type transport motions of CNT3

Nucleoside analogs have broad clinical utility as antiviral drugs. Key to their systemic distribution and cellular entry are human nucleoside transporters. Here, we establish that the human concentrative nucleoside transporter 3 (CNT3) interacts with antiviral drugs used in the treatment of coronavirus infections. We report high-resolution single-particle cryo-electron microscopy structures of bovine CNT3 complexed with antiviral nucleosides N4-hydroxycytidine, PSI-6206, GS-441524 and ribavirin, all in inward-facing states. Notably, we found that the orally bioavailable antiviral molnupiravir arrests CNT3 in four distinct conformations, allowing us to capture cryo-electron microscopy structures of drug-loaded outward-facing and drug-loaded intermediate states. Our studies uncover the conformational trajectory of CNT3 during membrane transport of a nucleoside analog antiviral drug, yield new insights into the role of interactions between the transport and the scaffold domains in elevator-like domain movements during drug translocation, and provide insights into the design of nucleoside analog antiviral prodrugs with improved oral bioavailability.

Insights into dynamic changes in ADC-7 and P99 cephalosporinases using small angle x-ray scattering (SAXS)

To counter the emergence of β-lactamase (BL) mediated resistance, design of new β-lactamase inhibitors (BLIs) is critical. Many high-resolution crystallographic structures of BL complexed with BLIs are available. However, their impact on BLI design is struggling to keep pace with novel and emerging variants. Small angle x-ray scattering (SAXS) in combination with molecular modeling is a useful tool to determine dynamic structures of macromolecules in solution. An important application of SAXS is to determine the conformational changes that occur when BLI bind to BL. To probe if conformational dynamics occur in class C cephalosporinases, we studied SAXS profiles of two clinically relevant class C β-lactamases, Acinetobacter baumannii ADC-7 and Enterobacter cloacae P99 in apo format complexed with BLIs. Importantly, SAXS data analysis demonstrated that in solution, these representative class C enzymes remain monomeric and did not show the associated assemblies that were seen in various crystal structures. SAXS data acquired for ADC-7 and P99, in apo and inhibitor bound states, clearly showed that these enzymes undergo detectable conformational changes, and these class C β-lactamases also close upon binding inhibitors as does BlaC. Further analysis revealed that addition of inhibitor led to the compacting of a range of residues around the active site, indicating that the conformational changes that both P99 and ADC-7 undergo are central to inhibitor recognition and efficacy. Our findings support the importance of exploring conformational changes using SAXS analysis in the design of future BLIs.Communicated by Ramaswamy H. Sarma.

CHARMM-GUI PDB Reader and Manipulator: Covalent Ligand Modeling and Simulation

Molecular modeling and simulation serve an important role in exploring biological functions of proteins at the molecular level, which is complementary to experiments. CHARMM-GUI (https://www.charmm-gui.org) is a web-based graphical user interface that generates complex molecular simulation systems and input files, and we have been continuously developing and expanding its functionalities to facilitate various complex molecular modeling and make molecular dynamics simulations more accessible to the scientific community. Currently, covalent drug discovery emerges as a popular and important field. Covalent drug forms a chemical bond with specific residues on the target protein, and it has advantages in potency for its prolonged inhibition effects. Even though there are higher demands in modeling PDB protein structures with various covalent ligand types, proper modeling of covalent ligands remains challenging. This work presents a new functionality in CHARMM-GUI PDB Reader & Manipulator that can handle a diversity of ligand-amino acid linkage types, which is validated by a careful benchmark study using over 1,000 covalent ligand structures in RCSB PDB. We hope that this new functionality can boost the modeling and simulation study of covalent ligands.

Extracellular ISG15 triggers ISGylation via a type-I interferon independent non-canonical mechanism to regulate host response during virus infection

Type-I interferons (IFN) induce cellular proteins with antiviral activity. One such protein is Interferon Stimulated Gene 15 (ISG15). ISG15 is conjugated to proteins during ISGylation to confer antiviral activity and regulate cellular activities associated with inflammatory and neurodegenerative diseases and cancer. Apart from ISGylation, unconjugated free ISG15 is also released from cells during various conditions, including virus infection. The role of extracellular ISG15 during virus infection was unknown. We show that extracellular ISG15 triggers ISGylation and acts as a soluble antiviral factor to restrict virus infection via an IFN-independent mechanism. Specifically, extracellular ISG15 acts post-translationally to markedly enhance the stability of basal intracellular ISG15 protein levels to support ISGylation. Furthermore, extracellular ISG15 interacts with cell surface integrin (α5β1 integrins) molecules via its RGD-like motif to activate the integrin-FAK (Focal Adhesion Kinase) pathway resulting in IFN-independent ISGylation. Thus, our studies have identified extracellular ISG15 protein as a new soluble antiviral factor that confers IFN-independent non-canonical ISGylation via the integrin-FAK pathway by post-translational stabilization of intracellular ISG15 protein.

Exploiting in silico structural analysis to introduce emerging genotype-phenotype correlations in DHCR24-related sterol biosynthesis disorder: a case study

Desmosterolosis is a rare sterol biosynthesis disorder characterized by multiple congenital anomalies, failure to thrive, severe developmental delay, progressive epileptic encephalopathy, and elevated levels of desmosterol caused by biallelic mutations of DHCR24 encoding 3-β-hydroxysterol Δ-24-reductase. DHCR24 is regarded as the key enzyme of cholesterol synthesis in the metabolism of brain cholesterol as it catalyzes the reduction of the Δ-24 double bond of sterol intermediates during cholesterol biosynthesis. To date, 15 DHCR24 variants, detected in 2 related and 14 unrelated patients, have been associated with the desmosterolosis disorder. Here, we describe a proband harboring the never-described DHCR24 homozygous missense variant NM_014762.4:c.506T>C, NP_055577.1:p.M169T, whose functional validation was confirmed through biochemical assay. By using molecular dynamics simulation techniques, we investigated the impact of this variant on the protein stability and interaction network with the flavin adenine dinucleotide cofactor, thereby providing a preliminary assessment of its mechanistic role in comparison to all known pathogenic variants, the wild-type protein, and a known benign DHCR24 variant. This report expands the clinical and molecular spectra of the DHCR24-related disorder, reports on a novel DHCR24 deleterious variant associated with desmosterolosis, and gives new insights into genotype-phenotype correlations.

Alchemical Free Energy Workflows for the Computation of Protein-Ligand Binding Affinities

Alchemical free energy methods can be used for the efficient computation of relative binding free energies during preclinical drug discovery stages. In recent years, this has been facilitated further by the implementation of workflows that enable non-experts to quickly and consistently set up the required simulations. Given the correct input structures, workflows handle the difficult aspects of setting up perturbations, including consistently defining the perturbable molecule, its atom mapping and topology generation, perturbation network generation, running of the simulations via different sampling methods, and analysis of the results. Different academic and commercial workflows are discussed, including FEW, FESetup, FEPrepare, CHARMM-GUI, Transformato, PMX, QLigFEP, TIES, ProFESSA, PyAutoFEP, BioSimSpace, FEP+, Flare, and Orion. These workflows differ in various aspects, such as mapping algorithms or enhanced sampling methods. Some workflows can accommodate more than one molecular dynamics (MD) engine and use external libraries for tasks. Differences between workflows can present advantages for different use cases, however a lack of interoperability of the workflows' components hinders systematic comparisons.

Analogs of α-conotoxin PnIC selectively inhibit α7β2- over α7-only subtype nicotinic acetylcholine receptors via a novel allosteric mechanism

This study was undertaken to identify and characterize the first ligands capable of selectively identifying nicotinic acetylcholine receptors containing α7 and β2 subunits (α7β2-nAChR subtype). Basal forebrain cholinergic neurons express α7β2-nAChR. Here, they appear to mediate neuronal dysfunction induced by the elevated levels of oligomeric amyloid-β associated with early Alzheimer's disease. Additional work indicates that α7β2-nAChR are expressed across several further critically important cholinergic and GABAergic neuronal circuits within the central nervous system. Further studies, however, are significantly hindered by the inability of currently available ligands to distinguish heteromeric α7β2-nAChR from the closely related and more widespread homomeric α7-only-nAChR subtype. Functional screening using two-electrode voltage-clamp electrophysiology identified a family of α7β2-nAChR-selective analogs of α-conotoxin PnIC (α-CtxPnIC). A combined electrophysiology, functional kinetics, site-directed mutagenesis, and molecular dynamics approach was used to further characterize the α7β2-nAChR selectivity and site of action of these α-CtxPnIC analogs. We determined that α7β2-nAChR selectivity of α-CtxPnIC analogs arises from interactions at a site distinct from the orthosteric agonist-binding site shared between α7β2- and α7-only-nAChR. As numerous previously identified α-Ctx ligands are competitive antagonists of orthosteric agonist-binding sites, this study profoundly expands the scope of use of α-Ctx ligands (which have already provided important nAChR research and translational breakthroughs). More immediately, analogs of α-CtxPnIC promise to enable, for the first time, both comprehensive mapping of the distribution of α7β2-nAChR and detailed investigations of their physiological roles.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

CHARMM-GUI Enhanced Sampler for various collective variables and enhanced sampling methods

Enhanced sampling methodologies modifying underlying Hamiltonians can be used for the systems with a rugged potential energy surface that makes it hard to observe convergence using conventional unbiased molecular dynamics (MD) simulations. We present CHARMM-GUI Enhanced Sampler, a web-based tool to prepare various enhanced sampling simulations inputs with user-selected collective variables (CVs). Enhanced Sampler provides inputs for the following nine methods: accelerated MD, Gaussian accelerated MD, conformational flooding, metadynamics, adaptive biasing force, steered MD, temperature replica exchange MD, replica exchange solute tempering 2, and replica exchange umbrella sampling for the method-implemented MD packages including AMBER, CHARMM, GENESIS, GROMACS, NAMD, and OpenMM. Users only need to select a group of atoms via intuitive web-implementation in order to define commonly used nine CVs of interest: center of mass based distance, angle, dihedral, root-mean-square-distance, radius of gyration, distance projected on axis, two types of angles projected on axis, and coordination numbers. The enhanced sampling methods are tested with several biological systems to illustrate their efficiency over conventional MD. Enhanced Sampler with carefully optimized system-dependent parameters will help users to get meaningful results from their enhanced sampling simulations.