Accelerating Membrane Simulations with Hydrogen Mass Repartitioning

Structure and mechanism of a phosphotransferase system glucose transporter

Glucose is the primary source of energy for many organisms and is efficiently taken up by bacteria through a dedicated transport system that exhibits high specificity. In Escherichia coli, the glucose-specific transporter IICBGlc serves as the major glucose transporter and functions as a component of the phosphoenolpyruvate-dependent phosphotransferase system. Here, we report cryo-electron microscopy (cryo-EM) structures of the glucose-bound IICBGlc protein. The dimeric transporter embedded in lipid nanodiscs was captured in the occluded, inward- and occluded, outward-facing conformations. Together with biochemical and biophysical analyses, and molecular dynamics (MD) simulations, we provide insights into the molecular basis and dynamics for substrate recognition and binding, including the gates regulating the binding sites and their accessibility. By combination of these findings, we present a mechanism for glucose transport across the plasma membrane. Overall, this work provides molecular insights into the structure, dynamics, and mechanism of the IICBGlc transporter in a native-like lipid environment.

Structural characterization of the POTRA domains from A. baumannii reveals new conformations in BamA

Recent studies have demonstrated BamA, the central component of the β-barrel assembly machinery (BAM), as an important therapeutic target to combat infections caused by Acinetobacter baumannii and other Gram-negative pathogens. Homology modeling indicates BamA in A. baumannii consists of five polypeptide transport-associated (POTRA) domains and a β-barrel membrane domain. We characterized the POTRA domains of BamA from A. baumannii in solution using size-exclusion chromatography small angle X-ray scattering (SEC-SAXS) analysis and determined crystal structures in two conformational states that are drastically different than those previously observed in BamA from other bacteria, indicating that the POTRA domains are even more conformationally dynamic than has been observed previously. Molecular dynamics simulations of the POTRA domains from A. baumannii and Escherichia coli allowed us to identify key structural features that contribute to the observed novel states. Together, these studies expand on our current understanding of the conformational plasticity within BamA across differing bacterial species.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 in apo. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganizations in key salt bridge and hydrogen bond networks contributing to the CB1 activation/inactivation. For instance, D213-Y224 hydrogen bond and D184-K192 salt bridge showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

Generative β-hairpin design using a residue-based physicochemical property landscape

De novo peptide design is a new frontier that has broad application potential in the biological and biomedical fields. Most existing models for de novo peptide design are largely based on sequence homology that can be restricted based on evolutionarily derived protein sequences and lack the physicochemical context essential in protein folding. Generative machine learning for de novo peptide design is a promising way to synthesize theoretical data that are based on, but unique from, the observable universe. In this study, we created and tested a custom peptide generative adversarial network intended to design peptide sequences that can fold into the β-hairpin secondary structure. This deep neural network model is designed to establish a preliminary foundation of the generative approach based on physicochemical and conformational properties of 20 canonical amino acids, for example, hydrophobicity and residue volume, using extant structure-specific sequence data from the PDB. The beta generative adversarial network model robustly distinguishes secondary structures of β hairpin from α helix and intrinsically disordered peptides with an accuracy of up to 96% and generates artificial β-hairpin peptide sequences with minimum sequence identities around 31% and 50% when compared against the current NCBI PDB and nonredundant databases, respectively. These results highlight the potential of generative models specifically anchored by physicochemical and conformational property features of amino acids to expand the sequence-to-structure landscape of proteins beyond evolutionary limits.

Polarizable atomic multipole-based force field for cholesterol

Cholesterol is one of the essential component of lipid in membrane. We present a polarizable atomic multipole force field (FF) for the molecular dynamic simulation of cholesterol. The FF building process follows the computational framework as the atomic multipole optimized energetics for biomolecular applications (AMOEBA) model. In this framework, the electronics parameters, including atomic monopole moments, dipole moments, and quadrupole moments calculated from ab initio calculations in the gas phase, are applied to represent the charge distribution. Furthermore, the many-body polarization is modeled by following the same pattern of distributed atomic polarizabilities. Then, the bilayers composed of two typical phospholipid molecules, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), in a range of different cholesterol concentrations are built and implemented by molecular dynamics (MD) simulations based on the proposed polarizable FF. The simulation results are statistically analyzed to validate the feasibility of the proposed FF. The structural properties of the bilayers are calculated to compare with the related experimental values. The MD values show the same trend of experimental values changes.Communicated by Ramaswamy H. Sarma.

Allosteric coupling of substrate binding and proton translocation in MmpL3 transporter from Mycobacterium tuberculosis

Infections caused by Mycobacterium spp. are very challenging to treat, and multidrug-resistant strains rapidly spread in human populations. Major contributing factors include the unique physiological features of these bacteria, drug efflux, and the low permeability barrier of their outer membrane. Here, we focus on MmpL3 from Mycobacterium tuberculosis, an essential inner membrane transporter of the resistance-nodulation-division superfamily required for the translocation of mycolic acids in the form of trehalose monomycolates (TMM) from the cytoplasm or plasma membrane to the periplasm or outer membrane. The MmpL3-dependent transport of TMM is essential for the growth of M. tuberculosis in vitro, inside macrophages, and in M. tuberculosis-infected mice. MmpL3 is also a validated target for several recently identified anti-mycobacterial agents. In this study, we reconstituted the lipid transport activity of the purified MmpL3 using a two-lipid vesicle system and established the ability of MmpL3 to actively extract phospholipids from the outer leaflet of a lipid bilayer. In contrast, we found that MmpL3 lacks the ability to translocate the same phospholipid substrate across the plasma membrane indicating that it is not an energy-dependent flippase. The lipid extraction activity was modulated by substitutions in critical charged and polar residues of the periplasmic substrate-binding pocket of MmpL3, coupled to the proton transfer activity of MmpL3 and inhibited by a small molecule inhibitor SQ109. Based on the results, we propose a mechanism of allosteric coupling wherein substrate translocation by MmpL3 is coupled to the energy provided by the downhill transfer of protons. The reconstituted activities will facilitate understanding the mechanism of MmpL3-dependent transport of lipids and the discovery of new therapeutic options for Mycobacterium spp. infections.IMPORTANCEMmpL3 from Mycobacterium tuberculosis is an essential transporter involved in the assembly of the mycobacterial outer membrane. It is also an important target in undergoing efforts to discover new anti-tuberculosis drugs effective against multidrug-resistant strains spreading in human populations. The recent breakthrough structural studies uncovered features of MmpL3 that suggested a possible lipid transport mechanism. In this study, we reconstituted and characterized the lipid transport activity of MmpL3 and demonstrated that this activity is blocked by MmpL3 inhibitors and substrate mimics. We further uncovered the mechanism of how the binding of a substrate in the periplasmic domain is communicated to the transmembrane proton relay of MmpL3. The uncovered mechanism and the developed assays provide new opportunities for mechanistic analyses of MmpL3 function and its inhibition.

Insight into molecular basis and dynamics of full-length CRaf kinase in cellular signaling mechanisms

Raf kinases play key roles in signal transduction in cells for regulating proliferation, differentiation, and survival. Despite decades of research into functions and dynamics of Raf kinases with respect to other cytosolic proteins, understanding Raf kinases is limited by the lack of their full-length structures at the atomic resolution. Here, we present the first model of the full-length CRaf kinase obtained from artificial intelligence/machine learning algorithms with a converging ensemble of structures simulated by large-scale temperature replica exchange simulations. Our model is validated by comparing simulated structures with the latest cryo-EM structure detailing close contacts among three key domains and regions of the CRaf. Our simulations identify potentially new epitopes of intramolecule interactions within the CRaf and reveal a dynamical nature of CRaf kinases, in which the three domains can move back and forth relative to each other for regulatory dynamics. The dynamic conformations are then used in a docking algorithm to shed insight into the paradoxical effect caused by vemurafenib in comparison with a paradox breaker PLX7904. We propose a model of Raf-heterodimer/KRas-dimer as a signalosome based on the dynamics of the full-length CRaf.

Membrane curvature sensing and symmetry breaking of the M2 proton channel from Influenza A

The M2 proton channel aids in the exit of mature influenza viral particles from the host plasma membrane through its ability to stabilize regions of high negative Gaussian curvature (NGC) that occur at the neck of budding virions. The channels are homo-tetramers that contain a cytoplasm-facing amphipathic helix (AH) that is necessary and sufficient for NGC generation; however, constructs containing the transmembrane spanning helix, which facilitates tetramerization, exhibit enhanced curvature generation. Here, we used all-atom molecular dynamics (MD) simulations to explore the conformational dynamics of M2 channels in lipid bilayers revealing that the AH is dynamic, quickly breaking the fourfold symmetry observed in most structures. Next, we carried out MD simulations with the protein restrained in four- and twofold symmetric conformations to determine the impact on the membrane shape. While each pattern was distinct, all configurations induced pronounced curvature in the outer leaflet, while conversely, the inner leaflets showed minimal curvature and significant lipid tilt around the AHs. The MD-generated profiles at the protein-membrane interface were then extracted and used as boundary conditions in a continuum elastic membrane model to calculate the membrane-bending energy of each conformation embedded in different membrane surfaces characteristic of a budding virus. The calculations show that all three M2 conformations are stabilized in inward-budding, concave spherical caps and destabilized in outward-budding, convex spherical caps, the latter reminiscent of a budding virus. One of the C2-broken symmetry conformations is stabilized by 4 kT in NGC surfaces with the minimum energy conformation occurring at a curvature corresponding to 33 nm radii. In total, our work provides atomistic insight into the curvature sensing capabilities of M2 channels and how enrichment in the nascent viral particle depends on protein shape and membrane geometry.

Model Mechanism for Lipid Uptake by the Human STARD2/PC-TP Phosphatidylcholine Transfer Protein

The human StAR-related lipid transfer domain protein 2 (STARD2), also known as phosphatidylcholine (PC) transfer protein, is a single-domain lipid transfer protein thought to transfer PC lipids between intracellular membranes. We performed extensive μs-long molecular dynamics simulations of STARD2 of its apo and holo forms in the presence or absence of complex lipid bilayers. The simulations in water reveal ligand-dependent conformational changes. In the 2 μs-long simulations of apo STARD2 in the presence of a lipid bilayer, we observed spontaneous reproducible PC lipid uptake into the protein hydrophobic cavity. We propose that the lipid extraction mechanism involves one to two metastable states stabilized by choline-tyrosine or choline-tryptophane cation-π interactions. Using free energy perturbation, we evaluate that PC-tyrosine cation-π interactions contribute 1.8 and 2.5 kcal/mol to the affinity of a PC-STARD2 metastable state, thus potentially providing a significant decrease of the energy barrier required for lipid desorption.

Structural insights into the mechanism and dynamics of proteorhodopsin biogenesis and retinal scavenging

Microbial ion-pumping rhodopsins (MRs) are extensively studied retinal-binding membrane proteins. However, their biogenesis, including oligomerisation and retinal incorporation, remains poorly understood. The bacterial green-light absorbing proton pump proteorhodopsin (GPR) has emerged as a model protein for MRs and is used here to address these open questions using cryo-electron microscopy (cryo-EM) and molecular dynamics (MD) simulations. Specifically, conflicting studies regarding GPR stoichiometry reported pentamer and hexamer mixtures without providing possible assembly mechanisms. We report the pentameric and hexameric cryo-EM structures of a GPR mutant, uncovering the role of the unprocessed N-terminal signal peptide in the assembly of hexameric GPR. Furthermore, certain proteorhodopsin-expressing bacteria lack retinal biosynthesis pathways, suggesting that they scavenge the cofactor from their environment. We shed light on this hypothesis by solving the cryo-EM structure of retinal-free proteoopsin, which together with mass spectrometry and MD simulations suggests that decanoate serves as a temporary placeholder for retinal in the chromophore binding pocket. Further MD simulations elucidate possible pathways for the exchange of decanoate and retinal, offering a mechanism for retinal scavenging. Collectively, our findings provide insights into the biogenesis of MRs, including their oligomeric assembly, variations in protomer stoichiometry and retinal incorporation through a potential cofactor scavenging mechanism.

Structural determinants for activity of the antidepressant vortioxetine at human and rodent 5-HT3 receptors

Vortioxetine (VTX) is a recently approved antidepressant that targets a variety of serotonin receptors. Here, we investigate the drug's molecular mechanism of operation at the serotonin 5-HT3 receptor (5-HT3R), which features two properties: VTX acts differently on rodent and human 5-HT3R, and VTX appears to suppress any subsequent response to agonists. Using a combination of cryo-EM, electrophysiology, voltage-clamp fluorometry and molecular dynamics, we show that VTX stabilizes a resting inhibited state of the mouse 5-HT3R and an agonist-bound-like state of human 5-HT3R, in line with the functional profile of the drug. We report four human 5-HT3R structures and show that the human receptor transmembrane domain is intrinsically fragile. We also explain the lack of recovery after VTX administration via a membrane partition mechanism.

In Silico Screening Identification of Fatty Acids and Fatty Acid Derivatives with Antiseizure Activity: In Vitro and In Vivo Validation

High fat diets have been used as complementary treatments for seizure disorders for more than a century. Moreover, many fatty acids and derivatives, including the broad-spectrum antiseizure medication valproic acid, have been explored and used as pharmacological agents to treat epilepsy. In this work, we have explored the anticonvulsant potential of a large library of fatty acids and fatty acid derivatives, the LIPID MAPS Structure Database, using structure-based virtual screening to assess their ability to block the voltage-gated sodium channel 1.2 (NaV1.2), a validated target for antiseizure medications. Four of the resulting in silico hits were submitted for experimental confirmation using in vitro patch clamp experiments, and their protective role was evaluated in an acute mice seizure model, the Maximal Electroshock seizure model. These four compounds were found to protect mice against seizures. Two of them exhibited blocking effects on NaV1.2, CaV2.2, and CaV3.1.

Unraveling the mechanism of small molecule induced activation of Staphylococcus aureus signal peptidase IB

Staphylococcus aureus signal peptidase IB (SpsB) is an essential enzyme for protein secretion. While inhibition of its activity by small molecules is a well-precedented mechanism to kill bacteria, the mode of activation is however less understood. We here investigate the activation mechanism of a recently introduced activator, the antibiotic compound PK150, and demonstrate by combined experimental and Molecular Dynamics (MD) simulation studies a unique principle of enzyme stimulation. Mass spectrometric studies with an affinity-based probe of PK150 unravel the binding site of PK150 in SpsB which is used as a starting point for MD simulations. Our model shows the localization of the molecule in an allosteric pocket next to the active site which shields the catalytic dyad from excess water that destabilizes the catalytic geometry. This mechanism is validated by the placement of mutations aligning the binding pocket of PK150. While the mutants retain turnover of the SpsB substrate, no stimulation of activity is observed upon PK150 addition. Overall, our study elucidates a previously little investigated mechanism of enzyme activation and serves as a starting point for the development of future enzyme activators.

The Influence of Ionic Environment on Nucleosome-Mica Interactions Revealed via Molecular Dynamics Simulations

Mica serves as a crucial substrate in Atomic Force Microscopy (AFM) studies for visualizing and characterizing nucleosomes. Nucleosomes interact with the negatively charged mica surface via adsorbed cations. However, the specific influences of monovalent and divalent cations on nucleosome adsorption to the mica surface remain unclear. In this study, we investigated the binding of nucleosomes to the mica surface in the presence of monovalent potassium ions and divalent magnesium ions using molecular dynamics simulations. We also explored the impact of pre-treated mica surfaces on nucleosome binding and structure. Our findings reveal that nucleosome-mica interactions vary depending on the cations present, resulting in distinct effects on nucleosome structure. Notably, nucleosomes bind effectively to a mica surface in the presence of potassium ions with minimal structural perturbations.

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

Enhanced Surface Accessibility of SARS-CoV-2 Omicron Spike Protein Due to an Altered Glycosylation Profile

SARS-CoV-2 spike (S) proteins undergo extensive glycosylation, aiding in proper folding, enhancing stability, and evading host immune surveillance. In this study, we used mass spectrometric analysis to elucidate the N-glycosylation characteristics and disulfide bonding of recombinant spike proteins derived from the SARS-CoV-2 Omicron variant (B.1.1.529) in comparison with the D614G spike variant. Furthermore, we conducted microsecond-long molecular dynamics simulations on spike proteins to resolve how the different N-glycans impact spike conformational sampling in the two variants. Our findings reveal that the Omicron spike protein maintains an overall resemblance to the D614G spike variant in terms of site-specific glycan processing and disulfide bond formation. Nonetheless, alterations in glycans were observed at certain N-glycosylation sites. These changes, in synergy with mutations within the Omicron spike protein, result in increased surface accessibility of the macromolecule, including the ectodomain, receptor-binding domain, and N-terminal domain. Additionally, mutagenesis and pull-down assays reveal the role of glycosylation of a specific sequon (N149); furthermore, the correlation of MD simulation and HDX-MS identified several high-dynamic areas of the spike proteins. These insights contribute to our understanding of the interplay between structure and function, thereby advancing effective vaccination and therapeutic strategies.

Nanopore Translocation Reveals Electrophoretic Force on Noncanonical RNA:DNA Double Helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on the translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of noncanonical RNA:DNA hybrids in electrophoretic transport to the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a noncanonical helix, with distinct transport properties from B-form DD molecules. We find that RD and DD molecules, with the same contour length, move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and the molecular understanding of electrophoretic transport.

Mg2+-dependent mechanism of environmental versatility in a multidrug efflux pump

Tripartite resistance nodulation and cell division multidrug efflux pumps span the periplasm and are a major driver of multidrug resistance among Gram-negative bacteria. The periplasm provides a distinct environment between the inner and outer membranes of Gram-negative bacteria. Cations, such as Mg2+, become concentrated within the periplasm and, in contrast to the cytoplasm, its pH is sensitive to conditions outside the cell. Here, we reveal an interplay between Mg2+ and pH in modulating the dynamics of the periplasmic adaptor protein, AcrA, and its function within the prototypical AcrAB-TolC multidrug efflux pump from Escherichia coli. In the absence of Mg2+, AcrA becomes increasingly plastic within acidic conditions, but when Mg2+ is bound this is ameliorated, resulting in domain specific organisation in neutral to weakly acidic regimes. We establish a unique histidine residue directs these structural dynamics and is essential for sustaining pump efflux activity across acidic, neutral, and alkaline conditions. Overall, we propose Mg2+ conserves the structural mobility of AcrA to ensure optimal AcrAB-TolC function within rapid changing environments commonly faced by the periplasm during bacterial infection and colonization. This work highlights that Mg2+ is an important mechanistic component in this pump class and possibly across other periplasmic lipoproteins.

Thicket and Mesh: How the Outer Membrane Can Resist Tension Imposed by the Cell Wall

The cell envelope of Gram-negative bacteria is composed of an outer membrane (OM) and an inner membrane (IM) and a peptidoglycan cell wall (CW) between them. Combined with Braun's lipoprotein (Lpp), which connects the OM and the CW, and numerous membrane proteins that exist in both OM and IM, the cell envelope creates a mechanically stable environment that resists various physical and chemical perturbations to the cell, including turgor pressure caused by the solute concentration difference between the cytoplasm of the cell and the extracellular environment. Previous computational studies have explored how individual components (OM, IM, and CW) can resist turgor pressure although combinations of them have been less well studied. To that end, we constructed multiple OM-CW systems, including the Lpp connections with the CW under increasing degrees of strain. The results show that the OM can effectively resist the tension imposed by the CW, shrinking by only 3-5% in area even when the CW is stretched to 2.5× its relaxed area. The area expansion modulus of the system increases with increasing CW strain, although the OM remains a significant contributor to the envelope's mechanical stability. Additionally, we find that when the protein TolC is embedded in the OM, its stiffness increases.

Nanopore translocation reveals electrophoretic force on non-canonical RNA:DNA double helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of non-canonical RNA:DNA hybrids in electrophoretic transport with the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a non-canonical helix with distinct transport properties from B-form DD molecules. We find RD and DD molecules with the same contour length move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same within experimental errors. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and molecular understanding of electrophoretic transport.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

The cannabinoid receptor CB1 is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 states in apo conditions. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganization of key salt bridge and hydrogen bond networks known to control CB1 activation between states. For instance, a conserved D213-Y224 hydrogen bond and D184-K192 salt bridge interactions showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

Expanding the CHARMM36 United Atom Chain Model for the Inclusion of Sphingolipids

The inclusion of accurate yet computationally inexpensive lipid force fields (FF) is pertinent for the study of lipids and lipid-containing systems using molecular dynamics (MD). Within the past decade, the implementation and further expansion of a united atom (UA) FF for lipids have been developed in the CHARMM family of FFs. The most recent version of the UA presented more accurate descriptions of lipid properties for several phospholipids with saturated and monounsaturated chains, termed C36UAr. However, the original C36UAr model lacks parameters for an important class of lipids, such as sphingolipids. The focus of this article is to broaden the scope of the C36UAr chain model to incorporate these lipids. In this study, two common sphingolipids, N-palmitoyl sphingomyelin and N-stearoyl sphingomyelin are converted to a UA-chain representation and simulated to investigate the accuracy and speed over the all-atom FF model for sphingolipids. Improvements were found among multiple parameters, for example, in the surface area per lipid (SA/lip) and hydrogen order parameters, over the all-atom simulations of these sphingomyelins in C36, while as much as halving the simulation time for simulations of the same setup otherwise. Thus, the accuracy and efficiency found in this study are consistent with those found in the C36UAr model for phospholipids and expand the application of C36UAr to a wider array of membrane models to better match that in vivo.

Substrate recognition mechanism of the endoplasmic reticulum-associated ubiquitin ligase Doa10

Doa10 (MARCHF6 in metazoans) is a large polytopic membrane-embedded E3 ubiquitin ligase in the endoplasmic reticulum (ER) that plays an important role in quality control of cytosolic and ER proteins. Although Doa10 is highly conserved across eukaryotes, it is not understood how Doa10 recognizes its substrates. Here, we define the substrate recognition mechanism of Doa10 by structural and functional analyses on Saccharomyces cerevisiae Doa10 and its model substrates. Cryo-EM analysis shows that Doa10 has unusual architecture with a large lipid-filled central cavity, and its conserved middle domain forms an additional water-filled lateral tunnel open to the cytosol. Our biochemical data and molecular dynamics simulations suggest that the entrance of the substrate's degron peptide into the lateral tunnel is required for efficient polyubiquitination. The N- and C-terminal membrane domains of Doa10 seem to form fence-like features to restrict polyubiquitination to those proteins that can access the central cavity and lateral tunnel. Our study reveals how extended hydrophobic sequences at the termini of substrate proteins are recognized by Doa10 as a signal for quality control.

A DNA turbine powered by a transmembrane potential across a nanopore

Rotary motors play key roles in energy transduction, from macroscale windmills to nanoscale turbines such as ATP synthase in cells. Despite our abilities to construct engines at many scales, developing functional synthetic turbines at the nanoscale has remained challenging. Here, we experimentally demonstrate rationally designed nanoscale DNA origami turbines with three chiral blades. These DNA nanoturbines are 24-27 nm in height and diameter and can utilize transmembrane electrochemical potentials across nanopores to drive DNA bundles into sustained unidirectional rotations of up to 10 revolutions s-1. The rotation direction is set by the designed chirality of the turbine. All-atom molecular dynamics simulations show how hydrodynamic flows drive this turbine. At high salt concentrations, the rotation direction of turbines with the same chirality is reversed, which is explained by a change in the anisotropy of the electrophoretic mobility. Our artificial turbines operate autonomously in physiological conditions, converting energy from naturally abundant electrochemical potentials into mechanical work. The results open new possibilities for engineering active robotics at the nanoscale.

Dissecting the Interaction Fingerprints and Binding Affinity of BYL719 Analogs Targeting PI3Kα

Phosphatidylinositol-3-kinase Alpha (PI3Kα) is a lipid kinase which regulates signaling pathways involved in cell proliferation. Dysregulation of these pathways promotes several human cancers, pushing for the development of anticancer drugs to target PI3Kα. One such medicinal chemistry campaign at Novartis led to the discovery of BYL719 (Piqray, Alpelicib), a PI3Kα inhibitor approved by the FDA in 2019 for treatment of HR+/HER2-advanced breast cancer with a PIK3CA mutation. Structure-based drug design played a key role in compound design and optimization throughout the discovery process. However, further characterization of potency drivers via structural dynamics and energetic analyses can be advantageous for ensuing PI3Kα programs. Here, our goal is to employ various in-silico techniques, including molecular simulations and machine learning, to characterize 14 ligands from the BYL719 analogs and predict their binding affinities. The structural insights from molecular simulations suggest that although the ligand-hinge interaction is the primary driver of ligand stability at the pocket, the R group positioning at C2 or C6 of pyridine/pyrimidine also plays a major role. Binding affinities predicted via thermodynamic integration (TI) are in good agreement with previously reported IC50s. Yet, computationally demanding techniques such as TI might not always be the most efficient approach for affinity prediction, as in our case study, fast high-throughput techniques were capable of classifying compounds as active or inactive, and one docking approach showed accuracy comparable to TI.

Substrate recognition mechanism of the endoplasmic reticulum-associated ubiquitin ligase Doa10

Doa10 (MARCH6 in metazoans) is a large polytopic membrane-embedded E3 ubiquitin ligase in the endoplasmic reticulum (ER) that plays an important role in quality control of cytosolic and ER proteins. Although Doa10 is highly conserved across eukaryotes, it is not understood how Doa10 recognizes its substrates. Here, we defined the substrate recognition mechanism of Doa10 by structural and functional analyses on Saccharomyces cerevisiae Doa10 and its well-defined degron Deg1. Cryo-EM analysis shows that Doa10 has unusual architecture with a large lipid-filled central cavity, and its conserved middle domain forms an additional water-filled lateral tunnel open to the cytosol. Our biochemical data and molecular dynamics simulations suggest that the entrance of the substrate's degron peptide into the lateral tunnel is required for efficient polyubiquitination. The N- and C-terminal membrane domains of Doa10 seem to form fence-like features to restrict polyubiquitination to those proteins that can access the central cavity and lateral tunnel.

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

Synthesis of Pyrrolo[3,4-b]pyridin-5-ones via Ugi-Zhu Reaction and In Vitro-In Silico Studies against Breast Carcinoma

An Ugi-Zhu three-component reaction (UZ-3CR) coupled in a one-pot manner to a cascade process (N-acylation/aza Diels-Alder cycloaddition/decarboxylation/dehydration) was performed to synthesize a series of pyrrolo[3,4-b]pyridin-5-ones in 20% to 92% overall yields using ytterbium triflate as a catalyst, toluene as a solvent, and microwaves as a heat source. The synthesized molecules were evaluated in vitro against breast cancer cell lines MDA-MB-231 and MCF-7, finding that compound 1f, at a concentration of 6.25 μM, exhibited a potential cytotoxic effect. Then, to understand the interactions between synthesized compounds and the main proteins related to the cancer cell lines, docking studies were performed on the serine/threonine kinase 1 (AKT1) and Orexetine type 2 receptor (Ox2R), finding moderate to strong binding energies, which matched accurately with the in vitro results. Additionally, molecular dynamics were performed between proteins related to the studied cell lines and the three best ligands.

The inhibitory mechanism of a small protein reveals its role in antimicrobial peptide sensing

A large number of small membrane proteins have been uncovered in bacteria, but their mechanism of action has remained mostly elusive. Here, we investigate the mechanism of a physiologically important small protein, MgrB, which represses the activity of the sensor kinase PhoQ and is widely distributed among enterobacteria. The PhoQ/PhoP two-component system is a master regulator of the bacterial virulence program and interacts with MgrB to modulate bacterial virulence, fitness, and drug resistance. A combination of cross-linking approaches with functional assays and protein dynamic simulations revealed structural rearrangements due to interactions between MgrB and PhoQ near the membrane/periplasm interface and along the transmembrane helices. These interactions induce the movement of the PhoQ catalytic domain and the repression of its activity. Without MgrB, PhoQ appears to be much less sensitive to antimicrobial peptides, including the commonly used C18G. In the presence of MgrB, C18G promotes MgrB to dissociate from PhoQ, thus activating PhoQ via derepression. Our findings reveal the inhibitory mechanism of the small protein MgrB and uncover its importance in antimicrobial peptide sensing.

Long-Range Conductivity in Proteins Mediated by Aromatic Residues

Single-molecule measurements show that many proteins, lacking any redox cofactors, nonetheless exhibit electrical conductance on the order of a nanosiemen over 10 nm distances, implying that electrons can transit an entire protein in less than a nanosecond when subject to a potential difference of less than 1 V. This is puzzling because, for fast transport (i.e., a free energy barrier of zero), the hopping rate is determined by the reorganization energy of approximately 0.8 eV, and this sets the time scale of a single hop to at least 1 μs. Furthermore, the Fermi energies of typical metal electrodes are far removed from the energies required for sequential oxidation and reduction of the aromatic residues of the protein, which should further reduce the hopping current. Here, we combine all-atom molecular dynamics (MD) simulations of non-redox-active proteins (consensus tetratricopeptide repeats) with an electron transfer theory to demonstrate a molecular mechanism that can account for the unexpectedly fast electron transport. According to our MD simulations, the reorganization energy produced by the energy shift on charging (the Stokes shift) is close to the conventional value of 0.8 eV. However, the non-ergodic sampling of molecular configurations by the protein results in reaction-reorganization energies, extracted directly from the distribution of the electrostatic energy fluctuations, that are only ∼0.2 eV, which is small enough to enable long-range conductivity, without invoking quantum coherent transport. Using the MD values of the reorganization energies, we calculate a current decay with distance that is in agreement with experiment.

Structural basis of mitochondrial protein import by the TIM23 complex

Mitochondria import nearly all of their approximately 1,000-2,000 constituent proteins from the cytosol across their double-membrane envelope1-5. Genetic and biochemical studies have shown that the conserved protein translocase, termed the TIM23 complex, mediates import of presequence-containing proteins (preproteins) into the mitochondrial matrix and inner membrane. Among about ten different subunits of the TIM23 complex, the essential multipass membrane protein Tim23, together with the evolutionarily related protein Tim17, has long been postulated to form a protein-conducting channel6-11. However, the mechanism by which these subunits form a translocation path in the membrane and enable the import process remains unclear due to a lack of structural information. Here we determined the cryo-electron microscopy structure of the core TIM23 complex (heterotrimeric Tim17-Tim23-Tim44) from Saccharomyces cerevisiae. Contrary to the prevailing model, Tim23 and Tim17 themselves do not form a water-filled channel, but instead have separate, lipid-exposed concave cavities that face in opposite directions. Our structural and biochemical analyses show that the cavity of Tim17, but not Tim23, forms the protein translocation path, whereas Tim23 probably has a structural role. The results further suggest that, during translocation of substrate polypeptides, the nonessential subunit Mgr2 seals the lateral opening of the Tim17 cavity to facilitate the translocation process. We propose a new model for the TIM23-mediated protein import and sorting mechanism, a central pathway in mitochondrial biogenesis.

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.

Insights into substrate transport and water permeation in the mycobacterial transporter MmpL3

Mycobacteria, such as Mycobacterium tuberculosis, are characterized by a uniquely thick and waxy cell envelope that consists of two membranes, with a variety of mycolates comprising their outer membrane (OM). The protein Mycobacterial membrane protein Large 3 (MmpL3) is responsible for the transport of a primary OM component, trehalose monomycolate (TMM), from the inner (cytoplasmic) membrane (IM) to the periplasmic space, a process driven by the proton gradient. Although multiple structures of MmpL3 with bound substrates have been solved, the exact pathway(s) for TMM or proton transport remains elusive. Here, employing molecular dynamics simulations we investigate putative pathways for either transport species. We hypothesized that MmpL3 will cycle through similar conformational states as the related transporter AcrB, which we used as targets for modeling the conformation of MmpL3. A continuous water pathway through the transmembrane region was found in one of these states, illustrating a putative pathway for protons. Additional equilibrium simulations revealed that TMM can diffuse from the membrane into a binding pocket in MmpL3 spontaneously. We also found that acetylation of TMM, which is required for transport, makes it more stable within MmpL3's periplasmic cavity compared with the unacetylated form.

Missense mutations in the membrane domain of PRRT2 affect its interaction with Nav1.2 voltage-gated sodium channels

PRRT2 is a neuronal protein that controls neuronal excitability and network stability by modulating voltage-gated Na⁺ channel (Nav). PRRT2 pathogenic variants cause pleiotropic syndromes including epilepsy, paroxysmal kinesigenic dyskinesia and episodic ataxia attributable to loss-of-function pathogenetic mechanism. Based on the evidence that the transmembrane domain of PRRT2 interacts with Nav1.2/1.6, we focused on eight missense mutations located within the domain that show expression and membrane localization similar to the wild-type protein. Molecular dynamics simulations showed that the mutants do not alter the structural stability of the PRRT2 membrane domain and preserve its conformation. Using affinity assays, we found that the A320V and V286M mutants displayed respectively decreased and increased binding to Nav1.2. Accordingly, surface biotinylation showed an increased Nav1.2 surface exposure induced by the A320V mutant. Electrophysiological analysis confirmed the lack of modulation of Nav1.2 biophysical properties by the A320V mutant with a loss-of-function phenotype, while the V286M mutant displayed a gain-of-function with respect to wild-type PRRT2 with a more pronounced left-shift of the inactivation kinetics and delayed recovery from inactivation. The data confirm the key role played by the PRRT2-Nav interaction in the pathogenesis of the PRRT2-linked disorders and suggest an involvement of the A320 and V286 residues in the interaction site. Given the similar clinical phenotype caused by the two mutations, we speculate that circuit instability and paroxysmal manifestations may arise when PRRT2 function is outside the physiological range.

Improving Speed and Affordability without Compromising Accuracy: Standard Binding Free-Energy Calculations Using an Enhanced Sampling Algorithm, Multiple-Time Stepping, and Hydrogen Mass Repartitioning

Accurate evaluation of protein-ligand binding free energies in silico is of paramount importance for understanding the mechanisms of biological regulation and providing a theoretical basis for drug design and discovery. Based on a series of atomistic molecular dynamics simulations in an explicit solvent, using well-tempered metadynamics extended adaptive biasing force (WTM-eABF) as an enhanced sampling algorithm, the so-called "geometrical route" offers a rigorous theoretical framework for binding affinity calculations that match experimental values. However, although robust, this strategy remains expensive, requiring substantial computational time to achieve convergence of the simulations. Improving the efficiency of the geometrical route, while preserving its reliability through improved ergodic sampling, is, therefore, highly desirable. In this contribution, having identified the computational bottleneck of the geometrical route, to accelerate the calculations we combine (i) a longer time step for the integration of the equations of motion with hydrogen-mass repartitioning (HMR), and (ii) multiple time-stepping (MTS) for collective-variable and biasing-force evaluation. Altogether, we performed 50 independent WTM-eABF simulations in triplicate for the "physical" separation of the Abl kinase-SH3 domain:p41 complex, following different HMR and MTS schemes, while tuning, in distinct protocols, the parameters of the enhanced-sampling algorithm. To demonstrate the consistency and reliability of the results obtained with the best-performing setups, we carried out quintuple simulations. Furthermore, we demonstrated the transferability of our method to other complexes by triplicating a 200 ns separation simulation of nine chosen protocols for the MDM2-p53:NVP-CGM097 complex. [Holzer et al. J. Med. Chem. 2015, 58, 6348-6358.] Our results, based on an aggregate simulation time of 14.4 μs, allowed an optimal set of parameters to be identified, able to accelerate convergence by a factor of three without any noticeable loss of accuracy.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels

The recent determination of cryo-EM structures of voltage-gated sodium (Na_v) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Na_v pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Na_v channels based on homology modeling of the cryo-EM structure of the human Na_v1.4 channel and, in addition, on the recently resolved configuration for Na_v1.2. In particular, single Na⁺ permeation events during standard MD runs suggest that the ion resides in the inner part of the Na_v selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na⁺ translocation through the SF of the homology-based Na_v1.2 model and the cryo-EM Na_v1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na⁺ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Na_v cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Na_v channels.

The intramembrane COOH-terminal domain of PRRT2 regulates voltage-dependent Na+ channels

Proline-rich transmembrane protein 2 (PRRT2) is the single causative gene for pleiotropic paroxysmal syndromes including epilepsy, kinesigenic dyskinesia, episodic ataxia and migraine. PRRT2 is a neuron-specific type-2 membrane protein with a COOH-terminal intramembrane domain and a long proline-rich NH₂-terminal cytoplasmic region. A large array of experimental data indicates that PRRT2 is a neuron stability gene that negatively controls intrinsic excitability by regulating surface membrane localization and biophysical properties of voltage-dependent Na⁺ channels Nav1.2 and Nav1.6, but not Nav1.1. To further investigate the regulatory role of PRRT2, we studied the structural features of this membrane protein with molecular dynamics simulations, and its structure-function relationships with Nav1.2 channels by biochemical and electrophysiological techniques. We found that the intramembrane COOH-terminal region maintains a stable conformation over time, with the first transmembrane domain forming a helix-loop-helix motif within the bilayer. The unstructured NH₂-terminal cytoplasmic region bound to the Nav1.2 better than the isolated COOH-terminal intramembrane domain, mimicking full-length PRRT2, while the COOH-terminal intramembrane domain was able to modulate Na⁺ current and channel biophysical properties, still maintaining the striking specificity for Nav1.2 vs Nav1.1. channels. The results identify PRRT2 as a dual-domain protein in which the NH₂-terminal cytoplasmic region acts as a binding antenna for Na⁺ channels, while the COOH-terminal membrane domain regulates channel exposure on the membrane and its biophysical properties.

Optimal Bond Constraint Topology for Molecular Dynamics Simulations of Cholesterol

We recently observed artificial temperature gradients in molecular dynamics (MD) simulations of phase-separating ternary lipid mixtures using the Martini 2 force field. We traced this artifact to insufficiently converged bond length constraints with typical time steps and default settings for the linear constraint solver (LINCS). Here, we systematically optimize the constraint scaffold of cholesterol. With massive virtual sites in an equimomental arrangement, we accelerate bond constraint convergence while preserving the original cholesterol force field and dynamics. The optimized model does not induce nonphysical temperature gradients even at relaxed LINCS settings and is at least as fast as the original model at the strict LINCS settings required for proper thermal sampling. We provide a python script to diagnose possible problems with constraint convergence for other molecules and force fields. Equimomental constraint topology optimization can also be used to boost constraint convergence in atomistic MD simulations of molecular systems.

Long-time-step molecular dynamics can retard simulation of protein-ligand recognition process

Molecular dynamics (MD) simulation of biologically relevant processes at realistic time scale and atomistic precision is generally limited by prohibitively large computational cost, due to its restriction of using an ultrashort integration time step (1-2 fs). A popular numerical recipe to reduce the associated computational burden is adopting schemes that would allow relatively longer-time-step for MD propagation. Here, we explore the perceived potential of one of the most frequently used long-time-step protocols, namely the hydrogen mass repartitioning (HMR) approach, in alleviating the computational overhead associated with simulation of the kinetic process of protein-ligand recognition events. By repartitioning the mass of heavier atoms to their linked hydrogen atoms, HMR leverages around twofold longer time step than regular simulation, holding promise of significant performance boost. However, our probe into direct simulation of the protein-ligand recognition event, one of the computationally most challenging processes, shows that long-time-step HMR MD simulations do not necessarily translate to a computationally affordable solution. Our investigations spanning cumulative 176 μs in three independent proteins (T4 lysozyme, sensor domain of MopR, and galectin-3) show that long-time-step HMR-based MD simulations can catch the ligand in its act of recognizing the native cavity. But, as a major caveat, the ligand is found to require significantly longer time to identify buried native protein cavity in an HMR MD simulation than regular simulation, thereby defeating the purpose of its usage for performance upgrade. A molecular analysis shows that the longer time required by a ligand to recognize the protein in HMR is rooted in faster diffusion of the ligand, which reduces the survival probability of decisive on-pathway metastable intermediates, thereby slowing down the eventual recognition process at the native cavity. Together, the investigation stresses careful assessment of pitfalls of long-time-step algorithms before attempting to utilize them for higher performance for biomolecular recognition simulations.

Q-RepEx: A Python pipeline to increase the sampling of empirical valence bond simulations

The exploration of chemical systems occurs on complex energy landscapes. Comprehensively sampling rugged energy landscapes with many local minima is a common problem for molecular dynamics simulations. These multiple local minima trap the dynamic system, preventing efficient sampling. This is a particular challenge for large biochemical systems with many degrees of freedom. Replica exchange molecular dynamics (REMD) is an approach that accelerates the exploration of the conformational space of a system, and thus can be used to enhance the sampling of complex biomolecular processes. In parallel, the empirical valence bond (EVB) approach is a powerful approach for modeling chemical reactivity in biomolecular systems. Here, we present an open-source Python-based tool that interfaces with the Q simulation package, and increases the sampling efficiency of the EVB free energy perturbation/umbrella sampling approach by means of REMD. This approach, Q-RepEx, both decreases the computational cost of the associated REMD-EVB simulations, and opens the door to more efficient studies of biochemical reactivity in systems with significant conformational fluctuations along the chemical reaction coordinate.

Structure and supramolecular organization of the canine distemper virus attachment glycoprotein

Canine distemper virus (CDV) is an enveloped RNA morbillivirus that triggers respiratory, enteric, and high incidence of severe neurological disorders. CDV induces devastating outbreaks in wild and endangered animals as well as in domestic dogs in countries associated with suboptimal vaccination programs. The receptor-binding tetrameric attachment (H)-protein is part of the morbilliviral cell entry machinery. Here, we present the cryo-electron microscopy (cryo-EM) structure and supramolecular organization of the tetrameric CDV H-protein ectodomain. The structure reveals that the morbilliviral H-protein is composed of three main domains: stalk, neck, and heads. The most unexpected feature was the inherent asymmetric architecture of the CDV H-tetramer being shaped by the neck, which folds into an almost 90° bent conformation with respect to the stalk. Consequently, two non-contacting receptor-binding H-head dimers, which are also tilted toward each other, are located on one side of an intertwined four helical bundle stalk domain. Positioning of the four protomer polypeptide chains within the neck domain is guided by a glycine residue (G158), which forms a hinge point exclusively in two protomer polypeptide chains. Molecular dynamics simulations validated the stability of the asymmetric structure under near physiological conditions and molecular docking showed that two receptor-binding sites are fully accessible. Thus, this spatial organization of the CDV H-tetramer would allow for concomitant protein interactions with the stalk and head domains without steric clashes. In summary, the structure of the CDV H-protein ectodomain provides new insights into the morbilliviral cell entry system and offers a blueprint for next-generation structure-based antiviral drug discovery.

A multipoint guidance mechanism for β-barrel folding on the SAM complex

Mitochondrial β-barrel proteins are essential for the transport of metabolites, ions and proteins. The sorting and assembly machinery (SAM) mediates their folding and membrane insertion. We report the cryo-electron microscopy structure of the yeast SAM complex carrying an early eukaryotic β-barrel folding intermediate. The lateral gate of Sam50 is wide open and pairs with the last β-strand (β-signal) of the substrate-the 19-β-stranded Tom40 precursor-to form a hybrid barrel in the membrane plane. The Tom40 barrel grows and curves, guided by an extended bridge with Sam50. Tom40's first β-segment (β1) penetrates into the nascent barrel, interacting with its inner wall. The Tom40 amino-terminal segment then displaces β1 to promote its pairing with Tom40's last β-strand to complete barrel formation with the assistance of Sam37's dynamic α-protrusion. Our study thus reveals a multipoint guidance mechanism for mitochondrial β-barrel folding.

Mechanism of hormone and allosteric agonist mediated activation of follicle stimulating hormone receptor

Follicle stimulating hormone (FSH) is an essential glycoprotein hormone for human reproduction, which functions are mediated by a G protein-coupled receptor, FSHR. Aberrant FSH-FSHR signaling causes infertility and ovarian hyperstimulation syndrome. Here we report cryo-EM structures of FSHR in both inactive and active states, with the active structure bound to FSH and an allosteric agonist compound 21 f. The structures of FSHR are similar to other glycoprotein hormone receptors, highlighting a conserved activation mechanism of hormone-induced receptor activation. Compound 21 f formed extensive interactions with the TMD to directly activate FSHR. Importantly, the unique residue H615^7.42 in FSHR plays an essential role in determining FSHR selectivity for various allosteric agonists. Together, our structures provide a molecular basis of FSH and small allosteric agonist-mediated FSHR activation, which could inspire the design of FSHR-targeted drugs for the treatment of infertility and controlled ovarian stimulation for in vitro fertilization.

Discovery of New Zika Protease and Polymerase Inhibitors through the Open Science Collaboration Project OpenZika

The Zika virus (ZIKV) is a neurotropic arbovirus considered a global threat to public health. Although there have been several efforts in drug discovery projects for ZIKV in recent years, there are still no antiviral drugs approved to date. Here, we describe the results of a global collaborative crowdsourced open science project, the OpenZika project, from IBM's World Community Grid (WCG), which integrates different computational and experimental strategies for advancing a drug candidate for ZIKV. Initially, molecular docking protocols were developed to identify potential inhibitors of ZIKV NS5 RNA-dependent RNA polymerase (NS5 RdRp), NS3 protease (NS2B-NS3pro), and NS3 helicase (NS3hel). Then, a machine learning (ML) model was built to distinguish active vs inactive compounds for the cytoprotective effect against ZIKV infection. We performed three independent target-based virtual screening campaigns (NS5 RdRp, NS2B-NS3pro, and NS3hel), followed by predictions by the ML model and other filters, and prioritized a total of 61 compounds for further testing in enzymatic and phenotypic assays. This yielded five non-nucleoside compounds which showed inhibitory activity against ZIKV NS5 RdRp in enzymatic assays (IC50 range from 0.61 to 17 μM). Two compounds thermally destabilized NS3hel and showed binding affinity in the micromolar range (Kd range from 9 to 35 μM). Moreover, the compounds LabMol-301 inhibited both NS5 RdRp and NS2B-NS3pro (IC50 of 0.8 and 7.4 μM, respectively) and LabMol-212 thermally destabilized the ZIKV NS3hel (Kd of 35 μM). Both also protected cells from death induced by ZIKV infection in in vitro cell-based assays. However, while eight compounds (including LabMol-301 and LabMol-212) showed a cytoprotective effect and prevented ZIKV-induced cell death, agreeing with our ML model for prediction of this cytoprotective effect, no compound showed a direct antiviral effect against ZIKV. Thus, the new scaffolds discovered here are promising hits for future structural optimization and for advancing the discovery of further drug candidates for ZIKV. Furthermore, this work has demonstrated the importance of the integration of computational and experimental approaches, as well as the potential of large-scale collaborative networks to advance drug discovery projects for neglected diseases and emerging viruses, despite the lack of available direct antiviral activity and cytoprotective effect data, that reflects on the assertiveness of the computational predictions. The importance of these efforts rests with the need to be prepared for future viral epidemic and pandemic outbreaks.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

GPU-Accelerated All-Atom Particle-Mesh Ewald Continuous Constant pH Molecular Dynamics in Amber

Constant pH molecular dynamics (MD) simulations sample protonation states on the fly according to the conformational environment and user specified pH conditions; however, the current accuracy is limited due to the use of implicit-solvent models or a hybrid solvent scheme. Here, we report the first GPU-accelerated implementation, parametrization, and validation of the all-atom continuous constant pH MD (CpHMD) method with particle-mesh Ewald (PME) electrostatics in the Amber22 pmemd.cuda engine. The titration parameters for Asp, Glu, His, Cys, and Lys were derived for the CHARMM c22 and Amber ff14sb and ff19sb force fields. We then evaluated the PME-CpHMD method using the asynchronous pH replica-exchange titration simulations with the c22 force field for six benchmark proteins, including BBL, hen egg white lysozyme (HEWL), staphylococcal nuclease (SNase), thioredoxin, ribonuclease A (RNaseA), and human muscle creatine kinase (HMCK). The root-mean-square deviation from the experimental pK_a's of Asp, Glu, His, and Cys is 0.76 pH units, and the Pearson's correlation coefficient for the pK_a shifts with respect to model values is 0.80. We demonstrated that a finite-size correction or much enlarged simulation box size can remove a systematic error of the calculated pK_a's and improve agreement with experiment. Importantly, the simulations captured the relevant biology in several challenging cases, e.g., the titration order of the catalytic dyad Glu35/Asp52 in HEWL and the coupled residues Asp19/Asp21 in SNase, the large pK_a upshift of the deeply buried catalytic Asp26 in thioredoxin, and the large pK_a downshift of the deeply buried catalytic Cys283 in HMCK. We anticipate that PME-CpHMD will offer proper pH control to improve the accuracies of MD simulations and enable mechanistic studies of proton-coupled dynamical processes that are ubiquitous in biology but remain poorly understood due to the lack of experimental tools and limitation of current MD simulations.

SARS-CoV-2 spike opening dynamics and energetics reveal the individual roles of glycans and their collective impact

The trimeric spike (S) glycoprotein, which protrudes from the SARS-CoV-2 viral envelope, binds to human ACE2, initiated by at least one protomer's receptor binding domain (RBD) switching from a "down" (closed) to an "up" (open) state. Here, we used large-scale molecular dynamics simulations and two-dimensional replica exchange umbrella sampling calculations with more than a thousand windows and an aggregate total of 160 μs of simulation to investigate this transition with and without glycans. We find that the glycosylated spike has a higher barrier to opening and also energetically favors the down state over the up state. Analysis of the S-protein opening pathway reveals that glycans at N165 and N122 interfere with hydrogen bonds between the RBD and the N-terminal domain in the up state, while glycans at N165 and N343 can stabilize both the down and up states. Finally, we estimate how epitope exposure for several known antibodies changes along the opening path. We find that the BD-368-2 antibody's epitope is continuously exposed, explaining its high efficacy.

Structure of the human heterodimeric transporter 4F2hc-LAT2 in complex with Anticalin, an alternative binding protein for applications in single-particle cryo-EM

Cryo-EM structure determination of relatively small and flexible membrane proteins at high resolution is challenging. Increasing the size and structural features by binding of high affinity proteins to the biomolecular target allows for better particle alignment and may result in structural models of higher resolution and quality. Anticalins are alternative binding proteins to antibodies, which are based on the lipocalin scaffold and show potential for theranostic applications. The human heterodimeric amino acid transporter 4F2hc-LAT2 is a membrane protein complex that mediates transport of certain amino acids and derivatives thereof across the plasma membrane. Here, we present and discuss the cryo-EM structure of human 4F2hc-LAT2 in complex with the anticalin D11vs at 3.2 Å resolution. Relative high local map resolution (2.8-3.0 Å) in the LAT2 substrate binding site together with molecular dynamics simulations indicated the presence of fixed water molecules potentially involved in shaping and stabilizing this region. Finally, the presented work expands the application portfolio of anticalins and widens the toolset of binding proteins to promote high-resolution structure solution by single-particle cryo-EM.

Ligand-mediated Structural Dynamics of a Mammalian Pancreatic KATP Channel

Regulation of pancreatic KATP channels involves orchestrated interactions of their subunits, Kir6.2 and SUR1, and ligands. Previously we reported KATP channel cryo-EM structures in the presence and absence of pharmacological inhibitors and ATP, focusing on the mechanisms by which inhibitors act as pharmacological chaperones of KATP channels (Martin et al., 2019). Here we analyzed the same cryo-EM datasets with a focus on channel conformational dynamics to elucidate structural correlates pertinent to ligand interactions and channel gating. We found pharmacological inhibitors and ATP enrich a channel conformation in which the Kir6.2 cytoplasmic domain is closely associated with the transmembrane domain, while depleting one where the Kir6.2 cytoplasmic domain is extended away into the cytoplasm. This conformational change remodels a network of intra- and inter-subunit interactions as well as the ATP and PIP2 binding pockets. The structures resolved key contacts between the distal N-terminus of Kir6.2 and SUR1's ABC module involving residues implicated in channel function and showed a SUR1 residue, K134, participates in PIP2 binding. Molecular dynamics simulations revealed two Kir6.2 residues, K39 and R54, that mediate both ATP and PIP2 binding, suggesting a mechanism for competitive gating by ATP and PIP2.

Dual A1/A3 Adenosine Receptor Antagonists: Binding Kinetics and Structure-Activity Relationship Studies Using Mutagenesis and Alchemical Binding Free Energy Calculations

Drugs targeting adenosine receptors (AR) can provide treatment for diseases. We report the identification of 7-(phenylamino)-pyrazolo[3,4-c]pyridines L2-L10, A15, and A17 as low-micromolar to low-nanomolar A₁R/A₃R dual antagonists, with 3-phenyl-5-cyano-7-(trimethoxyphenylamino)-pyrazolo[3,4-c]pyridine (A17) displaying the highest affinity at both receptors with a long residence time of binding, as determined using a NanoBRET-based assay. Two binding orientations of A17 produce stable complexes inside the orthosteric binding area of A₁R in molecular dynamics (MD) simulations, and we selected the most plausible orientation based on the agreement with alanine mutagenesis supported by affinity experiments. Interestingly, for drug design purposes, the mutation of L250^6.51 to alanine increased the binding affinity of A17 at A₁R. We explored the structure-activity relationships against A₁R using alchemical binding free energy calculations with the thermodynamic integration coupled with the MD simulation (TI/MD) method, applied on the whole G-protein-coupled receptor-membrane system, which showed a good agreement (r = 0.73) between calculated and experimental relative binding free energies.