Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations

Structure based computational RNA design towards MafA transcriptional repressor implicated in multiple myeloma

Multiple myeloma is recognized as the second most common hematological cancer. MafA transcriptional repressor is an established mediator of myelomagenesis. While there are multitude of drugs available for targeting various effectors in multiple myeloma, current literature lacks a candidate RNA based MafA modulator. Thus, using the structure of MafA homodimer-consensus target DNA, a computational effort was implemented to design a novel RNA based chemical modulator against MafA. First, available MafA-consensus DNA structure was employed to generate an RNA library. This library was further subjected to global docking to select the most plausible RNA candidates, preferring to bind DNA binding region of MafA. Following global docking, MD-ready complexes that were prepared via local docking program, were subjected to 500 ns of MD simulations. First, each of these MD simulations were analyzed for relative binding free energy through MM-PBSA method, which pointed towards a strong RNA based MafA binder, RNA1. Second, through a detailed MD analysis, RNA1 was shown to prefer binding to a single monomer of the dimeric DNA binding domain of MafA using higher number of hydrophobic interactions compared with positive control MafA-DNA complex. At the final phase, a principal component analyses was conducted, which led us to identify the actual interaction region of RNA1 and MafA monomer. Overall, to our knowledge, this is the first computational study that presents an RNA molecule capable of potentially targeting MafA protein. Furthermore, limitations of our study together with possible future implications of RNA1 in multiple myeloma were also discussed.

Complex state transitions of the bacterial cell division protein FtsZ

The key bacterial cell division protein FtsZ can adopt multiple conformations, and prevailing models suggest that transitions of FtsZ subunits from the closed to open state are necessary for filament formation and stability. Using all-atom molecular dynamics simulations, we analyzed state transitions of Staphylococcus aureus FtsZ as a monomer, dimer, and hexamer. We found that monomers can adopt intermediate states but preferentially adopt a closed state that is robust to forced reopening. Dimer subunits transitioned between open and closed states, and dimers with both subunits in the closed state remained highly stable, suggesting that open-state conformations are not necessary for filament formation. Mg2+ strongly stabilized the conformation of GTP-bound subunits and the dimer filament interface. Our hexamer simulations indicate that the plus end subunit preferentially closes and that other subunits can transition between states without affecting inter-subunit stability. We found that rather than being correlated with subunit opening, inter-subunit stability was strongly correlated with catalytic site interactions. By leveraging deep-learning models, we identified key intrasubunit interactions governing state transitions. Our findings suggest a greater range of possible monomer and filament states than previously considered and offer new insights into the nuanced interplay between subunit states and the critical role of nucleotide hydrolysis and Mg2+ in FtsZ filament dynamics.

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.

Membrane Stabilization of Helical Previtamin D Conformers as Possible Enhancement of Vitamin D Photoproduction

Photoinduced vitamin D formation occurs 10-15-fold faster in phospholipid bilayers (PLB) than in isotropic solution. It has been hypothesized that amphipatic interactions of the PLB with the rotationally flexible previtamin D (Pre) stabilize its helical conformers, enhancing thermal intramolecular [1,7]-hydrogen transfer, forming vitamin D. To test this hypothesis, we carried out molecular dynamics (MD) simulations of Pre in a PLB composed of dipalmitoylphosphatidylcholine (DPPC). We designed a classical force field capable of accurately describing the equilibrium composition of Pre conformers. Using adaptive biasing force MD simulations, we determined the free energy of Pre conformers in isotropic environments (hexane and gas-phase) and in the anisotropic environment of a DPPC PLB. We find a total increase of 25.5% of the population of both helical conformers (+20.5% g+Zg+ and +5% g-Zg-) in DPPC compared to hexane. In view of ab initio simulations, showing that hydrogen transfer occurs in both helical conformers, our study strongly suggests the validity of the initial hypothesis. Regarding the amphipatic interactions of Pre with the PLB, we find that, similar to cholesterol (Chol) and 7-dehydrocholesterol (7-DHC), Pre entertains hydrogen bonds mainly to the carbonyl groups of DPPC and, to a lesser extent, with phosphate oxygen atoms and rarely to water molecules at the interface. We further report order parameters of the Pre/DPPC system, which are slightly smaller than those for Chol/DPPC and 7-DHC/DPPC, but larger than for pure DPPC. This indicates a loss in membrane viscosity upon photochemical ring-opening of 7-DHC to form Pre.

Lipophilic compounds restore function to neurodevelopmental-associated KCNQ3 mutations

A major driver of neuronal hyperexcitability is dysfunction of K+ channels, including voltage-gated KCNQ2/3 channels. Their hyperpolarized midpoint of activation and slow activation and deactivation kinetics produce a current that regulates membrane potential and impedes repetitive firing. Inherited mutations in KCNQ2 and KCNQ3 are linked to a wide spectrum of neurodevelopmental disorders (NDDs), ranging from benign familial neonatal seizures to severe epileptic encephalopathies and autism spectrum disorders. However, the impact of these variants on the molecular mechanisms underlying KCNQ3 channel function remains poorly understood and existing treatments have significant side effects. Here, we use voltage clamp fluorometry, molecular dynamic simulations, and electrophysiology to investigate NDD-associated variants in KCNQ3 channels. We identified two distinctive mechanisms by which loss- and gain-of function NDD-associated mutations in KCNQ3 affect channel gating: one directly affects S4 movement while the other changes S4-to-pore coupling. MD simulations and electrophysiology revealed that polyunsaturated fatty acids (PUFAs) primarily target the voltage-sensing domain in its activated conformation and form a weaker interaction with the channel's pore. Consistently, two such compounds yielded partial and complete functional restoration in R227Q- and R236C-containing channels, respectively. Our results reveal the potential of PUFAs to be developed into therapies for diverse KCNQ3-based channelopathies.

The retinal chromophore environment in an inward light-driven proton pump studied by solid-state NMR and hydrogen-bond network analysis

Inward proton pumping is a relatively new function for microbial rhodopsins, retinal-binding light-driven membrane proteins. So far, it has been demonstrated for two unrelated subgroups of microbial rhodopsins, xenorhodopsins and schizorhodopsins. A number of recent studies suggest unique retinal-protein interactions as being responsible for the reversed direction of proton transport in the latter group. Here, we use solid-state NMR to analyze the retinal chromophore environment and configuration in an inward proton-pumping Antarctic schizorhodopsin. Using fully 13C-labeled retinal, we have assigned chemical shifts for every carbon atom and, assisted by structure modelling and molecular dynamics simulations, made a comparison with well-studied outward proton pumps, identifying locations of the unique protein-chromophore interactions for this functional subclass of microbial rhodopsins. Both the NMR results and molecular dynamics simulations point to the distinctive polar environment in the proximal part of the retinal, which may result in a hydration pattern dramatically different from that of the outward proton pumps, causing the reversed proton transport.

Structural Basis for Antagonist Binding to Vasopressin V1b Receptor Revealed by the Molecular Dynamics Simulations

The human V1b receptor (V1bR) is primarily expressed in the corticotropic cells of the anterior pituitary where it is involved in the regulation of the hypothalamic-pituitary-adrenal (HPA) axis. The activation of V1bR induces the secretion of adrenocorticotropin hormone (ACTH) from the anterior pituitary cells which, in turn, stimulates the production of cortisol via the adrenal cortex. Clinical studies have demonstrated the chronic dysfunction of the HPA axis in patients with several psychiatric disorders. Thus, the inhibition of the V1b receptor and normalizing the HPA axis hyperactivity is a promising approach to the treatment of many stress-related disorders such as anxiety and depression. Nelivaptan is a selective V1bR antagonist that can be used for this purpose and an excellent molecule to study how antagonists interact with V1bR, especially since in recent years the experimental structures of vasopressin V2 and oxytocin receptors were solved, providing high-similarity templates for homology modeling of V1bR. Therefore, in this work, six independent molecular dynamics simulations of a V1bR-nelivaptan complex in a fully hydrated lipid bilayer, yielding a total simulation time of 6.0 μs, have been conducted. In the lowest-energy complexes obtained in this work and proposed to be the most probable structure of the V1bR-nelivaptan complex, the location of the ligand inside the receptor pocket is very similar to that of the other ligands observed in the experimental structures of the vasopressin/oxytocin receptor family. The receptor-ligand interaction has been analyzed and described, revealing the details of the molecular mechanism of this antagonist binding to V1bR and a probable contribution of L2005×40 and T2035×43 to binding selectivity.

Predicting of novel homoserine dehydrogenase inhibitors against Paracoccidioides brasiliensis: integrating in silico and in vitro approaches

Aim: To search for potential inhibitors to homoserine dehydrogenase (HSD) in Paracoccidioides brasiliensis the causative agent of paracoccidioidomycosis, an infection with a high mortality rate in Brazil.Materials & methods: The enzyme was modeled and used in the virtual screening of the compounds. The library was first screened by the Autodock, in which 66 molecules were better ranked than substrate, and then, also evaluated by the Molegro and Gold programs.Results: The HS23 and HS87 molecules were selected in common by the three programs, and ADME/Tox evaluation indicates they are not toxic. The molecular dynamics of PbHSD bonded to ligands showed stable complexes until 50 ns. To validate the results, compounds were purchased for assays of minimum inhibitory concentration (MIC), minimum fungicidal concentration (MFC), synergic profile with Amphotericin B (AmB) and cytotoxicity. The two molecules presented MIC of 32 μg/ml and MFC of 64 μg/ml against the P. brasiliensis (strain Pb18). They also showed synergistic activity with AmB and a lack of toxicity against Hela and Vero cell lines.Conclusion: These results suggest that the HS23 and HS87 are promising candidates as PbHSD inhibitors and may be used as hits for the development of new drugs against paracoccidioidomycosis.

Multi-method analysis revealed the mechanism of substrate selectivity in NHase: A gatekeeper residue at the activity center

Recognizing the critical need to elucidate the molecular determinants of this selectivity offers a pathway to engineer enzymes with broader and more versatile catalytic capabilities. Through integrated methods including phylogenetic analysis, molecular docking, and structural analysis, we identified a pivotal amino acid residue, αTrp116, linking the substrate binding pocket and the active site of a NHase from Pseudonocardia thermophila JCM 3095 (PtNHase). This residue acts as a crucial determinant of substrate specificity within the NHase enzyme. The mutant αW116R modified the substrate specificity of PtNHase, significantly enhancing its catalytic efficiency towards aromatic substrates. The catalytic activity for aromatic compounds such as 3-Cyanopyridine was 14-fold that of the wild-type, whereas its activity for aliphatic substrates diminished to one-sixth. MD simulations revealed that replacing αTrp116 with Arg allowed aromatic nitrile substrates to achieve more favorable conformations within the active site. Based on the mutant αW116R, we further constructed a combinatorial variant Pt-4, tailored for aromatic substrates, which exhibited an enzyme activity 50 times that of the wild-type. These results highlight the critical influence of amino acid residues in the enzyme's active site on substrate specificity and offer fresh perspectives and approaches for the evolution of enzymes.

Potassium dependent structural changes in the selectivity filter of HERG potassium channels

The fine tuning of biological electrical signaling is mediated by variations in the rates of opening and closing of gates that control ion flux through different ion channels. Human ether-a-go-go related gene (HERG) potassium channels have uniquely rapid inactivation kinetics which are critical to the role they play in regulating cardiac electrical activity. Here, we exploit the K+ sensitivity of HERG inactivation to determine structures of both a conductive and non-conductive selectivity filter structure of HERG. The conductive state has a canonical cylindrical shaped selectivity filter. The non-conductive state is characterized by flipping of the selectivity filter valine backbone carbonyls to point away from the central axis. The side chain of S620 on the pore helix plays a central role in this process, by coordinating distinct sets of interactions in the conductive, non-conductive, and transition states. Our model represents a distinct mechanism by which ion channels fine tune their activity and could explain the uniquely rapid inactivation kinetics of HERG.

Exploring the phase behavior of C8-BTBT-C8 at ambient and high temperatures: insights and challenges from molecular dynamics simulations

C8-BTBT-C8 is one promising candidate for the development of high-performance electronic devices based on thin-film technologies. Its monoclinic polymorph has a well-established role in thin-film growth. Yet, quite little information is available about its dynamics on the molecular scale, and the structures of the mesophases which form at high temperature (about 100 K above ambient temperature). The present study is devoted to the analysis of such phases, with the ultimate goal of developing molecular models. Already at ambient temperature, our molecular dynamics simulations reveal a rich conformational behavior of the alkyl side chains, with gauche conformations as leading structural defects. Heating promotes the formation of a stacking faulted mesophase (380 K), and a smectic phase, at 385 K, upon side chain melting. Although more disordered, this phase bears several analogies with the smectic A phase, experimentally observed at 382.5 K. At higher temperatures, the increase in configurational disorder is brought by molecular diffusion and other phenomena, finally leading to an isotropic molten phase. Our in-depth analysis, complemented by hot-stage polarizing microscopy data, provides interesting insights into this material, highlighting the challenges associated with the modeling of soft semiconducting systems.

phosaa14SB and phosaa19SB: Updated Amber Force Field Parameters for Phosphorylated Amino Acids

Phosphorylated amino acids are involved in many cell regulatory networks; proteins containing these post-translational modifications are widely studied both experimentally and computationally. Simulations are used to investigate a wide range of structural and dynamic properties of biomolecules, such as ligand binding, enzyme-reaction mechanisms, and protein folding. However, the development of force field parameters for the simulation of proteins containing phosphorylated amino acids using the Amber program has not kept pace with the development of parameters for standard amino acids, and it is challenging to model these modified amino acids with accuracy comparable to proteins containing only standard amino acids. In particular, the popular ff14SB and ff19SB models do not contain parameters for phosphorylated amino acids. Here, the dihedral parameters for the side chains of the most common phosphorylated amino acids are trained against reference data from QM calculations adopting the ff14SB approach, followed by validation against experimental data. Library files and corresponding parameter files are provided, with versions that are compatible with both ff14SB and ff19SB.

Getting the intermolecular forces correct: introducing the ASTA strategy for a water model

Having a force field for water providing good bulk properties is paramount for modern studies of most biological systems. Some of the most common three-site force fields are TIP3, SPC/ε or OPC3, providing a decent range of bulk properties. That does not mean though, that they have realistic inter-atomic forces. These force fields have been parameterized with a top-down approach, meaning, by fitting the force field parameters to the experimental bulk properties. This approach has been the governing strategy also for many variants of four- and more-site models. We test a bottom-up approach, in which the force field is parameterized by optimizing the non-bonded inter-atomic forces. Our philosophy is that correct inter-atomic forces lead to correct geometrical and dynamical properties. The first system we try to optimize with the accurately system tailored atomic (ASTA) approach is water, but we aim to eventually probe other systems in the future as well. We applied our ASTA strategy to find a good set of parameters providing accurate bulk properties for the simple three-site force field forms, and also for AMOEBA, a more detailed and polarizable force field. Even though our bottom-up approach did not provide satisfactory results for the simple three-site force fields (with fixed charges), for the case of the AMOEBA force field it led to a modification of the original strategy, giving very good intra- and inter-molecular forces, as compared to accurate quantum chemically calculated reference forces. At the same time, important bulk properties, in this study restricted to the density and diffusion, were accurately reproduced with respect to the experimental values.

Structural basis of the obligatory exchange mode of human neutral amino acid transporter ASCT2

ASCT2 is an obligate exchanger of neutral amino acids, contributing to cellular amino acid homeostasis. ASCT2 belongs to the same family (SLC1) as Excitatory Amino Acid Transporters (EAATs) that concentrate glutamate in the cytosol. The mechanism that makes ASCT2 an exchanger rather than a concentrator remains enigmatic. Here, we employ cryo-electron microscopy and molecular dynamics simulations to elucidate the structural basis of the exchange mechanism of ASCT2. We establish that ASCT2 binds three Na+ ions per transported substrate and visits a state that likely acts as checkpoint in preventing Na+ ion leakage, both features shared with EAATs. However, in contrast to EAATs, ASCT2 retains one Na+ ion even under Na+-depleted conditions. We demonstrate that ASCT2 cannot undergo the structural transition in TM7 that is essential for the concentrative transport cycle of EAATs. This structural rigidity and the high-affinity Na+ binding site effectively confine ASCT2 to an exchange mode.

Residues in the fructose-binding pocket are required for ketohexokinase-A activity

Excessive fructose consumption is a primary contributor to the global surges in obesity, cancer, and metabolic syndrome. Fructolysis is not robustly regulated and is initiated by ketohexokinase (KHK). In this study, we determined the crystal structure of KHK-A, one of two human isozymes of KHK, in the apo-state at 1.85 Å resolution, and we investigated the roles of residues in the fructose-binding pocket by mutational analysis. Introducing alanine at D15, N42, or N45 inactivated KHK-A, whereas mutating R141 or K174 reduced activity and thermodynamic stability. Kinetic studies revealed that the R141A and K174A mutations reduced fructose affinity by 2- to 4-fold compared to WT KHK-A, without affecting ATP affinity. Molecular dynamics simulations provided mechanistic insights into the potential roles of the mutated residues in ligand coordination and the maintenance of an open state in one monomer and a closed state in the other. Protein-protein interactome analysis indicated distinct expression patterns and downregulation of partner proteins in different tumor tissues, warranting a reevaluation of KHK's role in cancer development and progression. The connections between different cancer genes and the KHK signaling pathway suggest that KHK is a potential target for preventing cancer metastasis. This study enhances our understanding of KHK-A's structure and function and offers valuable insights into potential targets for developing treatments for obesity, cancer, and metabolic syndrome.

Prediction of the essential intermolecular contacts for side-binding of VASP on F-actin

Vasodilator-stimulated phosphoprotein (VASP) family proteins play a crucial role in mediating the actin network architecture in the cytoskeleton. The Ena/VASP homology 2 (EVH2) domain in each of the four identical arms of the tetrameric VASP consists of a loading poly-Pro region, a G-actin-binding domain (GAB), and an F-actin-binding domain (FAB). Together, the poly-Pro, GAB, and FAB domains allow VASP to bind to sides of actin filaments in a bundle, and recruit profilin-G-actin to processively elongate the filaments. The atomic resolution structure of the ternary complex, consisting of the loading poly-Pro region and GAB domain of VASP with profilin-actin, has been solved over a decade ago; however, a detailed structure of the FAB-F-actin complex has not been resolved to date. Experimental insights, based on homology of the FAB domain with the C region of WASP, have been used to hypothesize that the FAB domain binds to the cleft between subdomains 1 and 3 of F-actin. Here, in order to develop our understanding of the VASP-actin complex, we first augment known structural information about the GAB domain binding to actin with the missing FAB domain-actin structure, which we predict using homology modeling and docking simulations. In earlier work, we used mutagenesis and kinetic modeling to study the role of domain-level binding-unbinding kinetics of Ena/VASP on actin filaments in a bundle, specifically on the side of actin filaments. We further look at the nature of the side-binding of the FAB domain of VASP at the atomistic level using our predicted structure, and tabulate effective mutation sites on the FAB domain that would disrupt the VASP-actin complex. We test the binding affinity of Ena with mutated FAB domain using total internal reflection fluorescence microscopy experiments. The binding affinity of VASP is affected significantly for the mutant, providing additional support for our predicted structure.

Hydrogen-Bonding and Hydrophobic Interaction Networks as Structural Determinants of Microbial Rhodopsin Function

Microbial pump rhodopsins are highly versatile light-driven membrane proteins that couple protein conformational dynamics with ion translocation across the cell membranes. Understanding how microbial pump rhodopsins use specific amino acid residues at key functional sites to control ion selectivity and ion pumping direction is of general interest for membrane transporters, and could guide site-directed mutagenesis for optogenetics applications. To enable direct comparisons between proteins with different sequences we implement, for the first time, a unique numbering scheme for the microbial pump rhodopsin residues, NS-mrho. We use NS-mrho to show that distinct microbial pump rhodopsins typically have hydrogen-bond networks that are less conserved than anticipated from the amino acid residue conservation, whereas their hydrophobic interaction networks are largely conserved. To illustrate the role of the hydrogen-bond networks as structural elements that determine the functionality of microbial pump rhodopsins, we performed experiments, atomic-level simulations, and hydrogen bond network analyses on GR, the outward proton pump from Gloeobacter violaceus, and KR2, the outward sodium pump from Krokinobacter eikastus. The experiments indicate that multiple mutations that recover KR2 amino acid residues in GR not only fail to convert it into a sodium pump, but completely inactivate GR by abolishing photoisomerization of the retinal chromophore. This observation could be attributed to the drastically altered hydrogen-bond interaction network identified with simulations and network analyses. Taken together, our findings suggest that functional specificity could be encoded in the collective hydrogen-bond network of microbial pump rhodopsins.

Covalent inhibition of pro-apoptotic BAX

BCL-2-associated X protein (BAX) is a promising therapeutic target for activating or restraining apoptosis in diseases of pathologic cell survival or cell death, respectively. In response to cellular stress, BAX transforms from a quiescent cytosolic monomer into a toxic oligomer that permeabilizes the mitochondria, releasing key apoptogenic factors. The mitochondrial lipid trans-2-hexadecenal (t-2-hex) sensitizes BAX activation by covalent derivatization of cysteine 126 (C126). In this study, we performed a disulfide tethering screen to discover C126-reactive molecules that modulate BAX activity. We identified covalent BAX inhibitor 1 (CBI1) as a compound that selectively derivatizes BAX at C126 and inhibits BAX activation by triggering ligands or point mutagenesis. Biochemical and structural analyses revealed that CBI1 can inhibit BAX by a dual mechanism of action: conformational constraint and competitive blockade of lipidation. These data inform a pharmacologic strategy for suppressing apoptosis in diseases of unwanted cell death by covalent targeting of BAX C126.

Biochemical and in silico characterization of glycosyltransferases from red sweet cherry (Prunus avium L.) reveals their broad specificity toward phenolic substrates

Polyphenolic compounds are a class of phytonutrients that play important roles in plants and contribute to human health when incorporated into our diet through fruit consumption. A large proportion occur as glycoconjugates but the enzymes responsible for their glycosylation are poorly characterized. Here, we report the biochemical and structural characterization of two glycosyltransferases from sweet cherry named PaUGT1 and PaUGT2. Both are promiscuous glucosyltransferases active on diverse anthocyanidins and flavonols, as well as phenolic acids in the case of PaUGT1. They also exhibit weaker galactosyltransferase activity. The expression of the gene encoding PaUGT1, the most active of the two proteins, follows anthocyanin accumulation during fruit ripening, suggesting that this enzyme is the primary glycosyltransferase involved in flavonoid glycosylation in sweet cherry. It can potentially be used to synthesize diverse glycoconjugates of flavonoids for integration into bioactive formulations, and for generating new fruit cultivars with enhanced health-promoting properties using breeding methods.

Interior pH-sensing residue of human voltage-gated proton channel Hv1 is histidine 168

The molecular mechanisms governing the human voltage-gated proton channel hHv1 remain elusive. Here, we used membrane-enabled hybrid-solvent continuous constant pH molecular dynamics (CpHMD) simulations with pH replica exchange to further evaluate the structural models of hHv1 in the closed (hyperpolarized) and open (depolarized) states recently obtained with MD simulations and explore potential pH-sensing residues. The CpHMD titration at a set of symmetric pH conditions revealed three residues that can gain or lose protons upon channel depolarization. Among them, residue H168 at the intracellular end of the S3 helix switches from the deprotonated to the protonated state and its protonation is correlated with the increased tilting of the S3 helix during the transition from the closed to the open state. Thus, the simulation data suggest H168 as an interior pH sensor, in support of a recent finding based on electrophysiological experiments of Hv1 mutants. We propose that protonation of H168 acts as a key that unlocks the closed channel configuration by increasing the flexibility of the S2-S3 linker, which increases the tilt angle of S3 and enhances the mobility of the S4 helix, thus promoting channel opening. Our work represents an important step toward deciphering the pH-dependent gating mechanism of hHv1.

Unlocking the unfolded structure of ubiquitin: Combining time-resolved x-ray solution scattering and molecular dynamics to generate unfolded ensembles

The unfolding dynamics of ubiquitin were studied using a combination of x-ray solution scattering (XSS) and molecular dynamics (MD) simulations. The kinetic analysis of the XSS ubiquitin signals showed that the protein unfolds through a two-state process, independent of the presence of destabilizing salts. In order to characterize the ensemble of unfolded states in atomic detail, the experimental XSS results were used as a constraint in the MD simulations through the incorporation of x-ray scattering derived potential to drive the folded ubiquitin structure toward sampling unfolded states consistent with the XSS signals. We detail how biased MD simulations provide insight into unfolded states that are otherwise difficult to resolve and underscore how experimental XSS data can be combined with MD to efficiently sample structures away from the native state. Our results indicate that ubiquitin samples unfolded in states with a high degree of loss in secondary structure yet without a collapse to a molten globule or fully solvated extended chain. Finally, we propose how using biased-MD can significantly decrease the computational time and resources required to sample experimentally relevant nonequilibrium states.

Dual mechanisms contribute to enhanced voltage dependence of an electric fish potassium channel

The voltage dependence of different voltage-gated potassium channels, described by the voltage at which half of the channels are open (V1/2), varies over a range of 80 mV and is influenced by factors such as the number of positive gating charges and the identity of the hydrophobic amino acids in the channel's voltage sensor (S4). Here we explore by experimental manipulations and molecular dynamics simulation the contributions of two derived features of an electric fish potassium channel (Kv1.7a) that is among the most voltage-sensitive Shaker family potassium channels known. These are a patch of four contiguous negatively charged glutamates in the S3-S4 extracellular loop and a glutamate in the S3b helix. We find that these negative charges affect V1/2 by separate, complementary mechanisms. In the closed state, the S3-S4 linker negative patch reduces the membrane surface charge biasing the channel to enter the open state while, upon opening, the negative amino acid in the S3b helix faces the second (R2) gating charge of the voltage sensor electrostatically biasing the channel to remain in the open state. This work highlights two evolutionary novelties that illustrate the potential influence of negatively charged amino acids in extracellular loops and adjacent helices to voltage dependence.

Dissecting the neuroprotective interaction between the BH4 domain of BCL-w and the IP3 receptor

BCL-w is a BCL-2 family protein that promotes cell survival in tissue- and disease-specific contexts. The canonical anti-apoptotic functionality of BCL-w is mediated by a surface groove that traps the BCL-2 homology 3 (BH3) α-helices of pro-apoptotic members, blocking cell death. A distinct N-terminal portion of BCL-w, termed the BCL-2 homology 4 (BH4) domain, selectively protects axons from paclitaxel-induced degeneration by modulating IP3 receptors, a noncanonical BCL-2 family target. Given the potential of BCL-w BH4 mimetics to prevent or mitigate chemotherapy-induced peripheral neuropathy, we sought to characterize the interaction between BCL-w BH4 and the IP3 receptor, combining "staple" and alanine scanning approaches with molecular dynamics simulations. We generated and identified stapled BCL-w BH4 peptides with optimized IP3 receptor binding and neuroprotective activities. Point mutagenesis further revealed the sequence determinants for BCL-w BH4 specificity, providing a blueprint for therapeutic targeting of IP3 receptors to achieve neuroprotection.

Computational Insights into the Intramolecular Aromatic C-C Coupling Catalyzed by the Cytochrome P450 Enzyme CYP121 from Mycobacterium tuberculosis

CYP121 is a P450 enzyme that catalyzes the intramolecular C-C coupling of its native substrate, dicyclotyrosine (cYY). According to previous suggestions, when the cosubstrate peracetic acid was used to generate Cpd I, the substrate cYY was suggested to participate in the cleavage of the O-O bond; however, whether cYY is involved in the formation of Cpd I and how two distant aromatic carbon atoms are activated are still unclear. Here, we constructed computational models and performed QM/MM calculations to clarify the reaction mechanism. On the basis of our calculation results, cYY is not involved in the formation of Cpd I, and the C-C coupling reaction starts from hydrogen abstraction. In the second stage, the substrate should first undergo a complex conformational change, leading to two phenolic hydroxyls of cYY close to each other. In the subsequent reaction, the resultant Cpd II again abstracts a hydrogen atom from the proximal tyrosine to generate the diradical intermediate. In addition, the C-C coupling occurs in the active site, but the final aromatization may be a nonenzymatic reaction. In general, the intramolecular C-C coupling requires two basic conditions, including the active site having good flexibility and the substrate itself having a suitable and rotatable skeleton.

Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment

This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid-base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

Ophthalmic drug effects on the amyloidogenesis of a transforming growth factor β-induced protein (TGFBIp) peptide fragment

Drugs that can treat one disease may either be detrimental or beneficial toward another due to possible cross-interactions. Therefore, care in choosing a suitable drug for patients with multiple diseases is crucial in successful patient management. This study explores several currently available ophthalmic drugs used to treat common ocular diseases to understand how they can affect the amyloidogenesis of a transforming growth factor β-induced protein (TGFBIp) peptide fragment found in abundance in the corneal protein aggregation deposits of lattice corneal dystrophy (LCD) patients. Results from this study provided supporting evidence that some drugs intended to treat other diseases can enhance or inhibit fibrillar aggregation of TGFBIp peptide, which may have potential implication of affecting the disease progression of LCD by either worsening or ameliorating it. Comparisons of the different properties of ophthalmic compounds explored in this study may also provide some guidance for future design of drugs geared toward the treatment of LCD.

Structural Evolution of Bacterial Polyphosphate Degradation Enzyme for Phosphorus Cycling

Living organisms ranging from bacteria to animals have developed their own ways to accumulate and store phosphate during evolution, in particular as the polyphosphate (polyP) granules in bacteria. Degradation of polyP into phosphate is involved in phosphorus cycling, and exopolyphosphatase (PPX) is the key enzyme for polyP degradation in bacteria. Thus, understanding the structure basis of PPX is crucial to reveal the polyP degradation mechanism. Here, it is found that PPX structure varies in the length of ɑ-helical interdomain linker (ɑ-linker) across various bacteria, which is negatively correlated with their enzymatic activity and thermostability - those with shorter ɑ-linkers demonstrate higher polyP degradation ability. Moreover, the artificial DrPPX mutants with shorter ɑ-linker tend to have more compact pockets for polyP binding and stronger subunit interactions, as well as higher enzymatic efficiency (kcat/Km) than that of DrPPX wild type. In Deinococcus-Thermus, the PPXs from thermophilic species possess a shorter ɑ-linker and retain their catalytic ability at high temperatures (70 °C), which may facilitate the thermophilic species to utilize polyP in high-temperature environments. These findings provide insights into the interdomain linker length-dependent evolution of PPXs, which shed light on enzymatic adaption for phosphorus cycling during natural evolution and rational design of enzyme.

Development, validation and analysis of a human profurin 3D model using comparative modeling and molecular dynamics simulations

The emergence of new viruses can lead to the outbreak of pandemics as occurred at the end of 2019 with the coronavirus disease (or COVID-19). The fastest way to effectively control viral infections is to develop broad-spectrum antivirals that can fight at least an entire class of viruses. Profurin, the furin precursor propeptide, is responsible for the autoactivation step which is crucial for the maturation of several viral substrates. This role makes the study of furin and profurin interactions interesting for the development of new potential broad-spectrum antivirals for the treatment against several human viral diseases. Since there is no 3D model of profurin published in the literature or deposited in a database, this work reports the development, validation and analysis of a profurin 3D model using comparative modeling and molecular dynamics. The model is available in ModelArchive at https://www.modelarchive.org/doi/10.5452/ma-ct8l7. The usage of this model will make possible further studies of molecular docking and MD simulations of the profurin-furin system, in the design of new potential broad-spectrum antivirals for the treatment against several human viral diseases.Communicated by Ramaswamy H. Sarma.

Structure and slow dynamics of protein hydration water with cryopreserving DMSO and trehalose upon cooling

We study, through molecular dynamics simulations, three aqueous solutions with one lysozyme protein and three different concentrations of trehalose and dimethyl sulfoxide (DMSO). We analyze the structural and dynamical properties of the protein hydration water upon cooling. We find that trehalose plays a major role in modifying the structure of the network of HBs between water molecules in the hydration layer of the protein. The dynamics of hydration water presents, in addition to the α-relaxation, typical of glass formers, a slower long-time relaxation process, which greatly slows down the dynamics of water, particularly in the systems with trehalose, where it becomes dominant at low temperatures. In all the solutions, we observe, from the behavior of the α-relaxation times, a shift of the Mode Coupling Theory crossover temperature and the fragile-to-strong crossover temperature toward higher values with respect to bulk water. We also observe a strong-to-strong crossover from the temperature behavior of the long-relaxation times. In the aqueous solution with only DMSO, the transition shifts to a lower temperature than in the case with only lysozyme reported in the literature. We observe that the addition of trehalose to the mixture has the opposite effect of restoring the original location of the strong-to-strong crossover. In all the solutions analyzed in this work, the observed temperature of the protein dynamical transition is slightly shifted at lower temperatures than that of the strong-to-strong crossover, but their relative order is the same, showing a correlation between the motion of the protein and that of the hydration water.

Applications of the Newly Developed Force-Field Parameters Uncover a Dynamic Nature of Ω-Loop C in the Lys-Ligated Alkaline Form of Cytochrome c

Lys-ligated cytochromes make up an emerging family of heme proteins. Density functional theory calculations on the amine/imidazole-ligated c-type ferric heme were employed to develop force-field parameters for molecular dynamics (MD) simulations of structural and dynamic features of these proteins. The new force-field parameters were applied to the alkaline form of yeast iso-1 cytochrome c to rationalize discrepancies resulting from distinct experimental conditions in prior structural studies and to provide insights into the mechanisms of the alkaline transition. Our simulations have revealed the dynamic nature of Ω-loop C in the Lys-ligated protein and its unfolding in the Lys-ligated conformer having this loop in the same position as in the native Met-ligated protein. The proximity of Tyr67 or Tyr74 to the Lys ligand of ferric heme iron suggests a possible mechanism of the backward alkaline transition where a proton donor Tyr assists in Lys dissociation. The developed force-field parameters will be useful in structural and dynamic characterization of other native or engineered Lys-ligated heme proteins.

A Practical Guide to All-Atom and Coarse-Grained Molecular Dynamics Simulations Using Amber and Gromacs: A Case Study of Disulfide-Bond Impact on the Intrinsically Disordered Amyloid Beta

Intrinsically disordered proteins (IDPs) pose challenges to conventional experimental techniques due to their large-scale conformational fluctuations and transient structural elements. This work presents computational methods for studying IDPs at various resolutions using the Amber and Gromacs packages with both all-atom (Amber ff19SB with the OPC water model) and coarse-grained (Martini 3 and SIRAH) approaches. The effectiveness of these methodologies is demonstrated by examining the monomeric form of amyloid-β (Aβ42), an IDP, with and without disulfide bonds at different resolutions. Our results clearly show that the addition of a disulfide bond decreases the β-content of Aβ42; however, it increases the tendency of the monomeric Aβ42 to form fibril-like conformations, explaining the various aggregation rates observed in experiments. Moreover, analysis of the monomeric Aβ42 compactness, secondary structure content, and comparison between calculated and experimental chemical shifts demonstrates that all three methods provide a reasonable choice to study IDPs; however, coarse-grained approaches may lack some atomistic details, such as secondary structure recognition, due to the simplifications used. In general, this study not only explains the role of disulfide bonds in Aβ42 but also provides a step-by-step protocol for setting up, conducting, and analyzing molecular dynamics (MD) simulations, which is adaptable for studying other biomacromolecules, including folded and disordered proteins and peptides.

Structural Determinants of Peptide Nanopore Formation

We have evolved the nanopore-forming macrolittin peptides from the bee venom peptide melittin using successive generations of synthetic molecular evolution. Despite their sequence similarity to the broadly membrane permeabilizing cytolytic melittin, the macrolittins have potent membrane selectivity. They form nanopores in synthetic bilayers made from 1-palmitoyl, 2-oleoyl-phosphatidylcholine (POPC) at extremely low peptide concentrations and yet have essentially no cytolytic activity against any cell membrane, even at high concentration. Here, we explore the structural determinants of macrolittin nanopore stability in POPC bilayers using atomistic molecular dynamics simulations and experiments on macrolittins and single-site variants. Simulations of macrolittin nanopores in POPC bilayers show that they are stabilized by an extensive, cooperative hydrogen bond network comprised of the many charged and polar side chains interacting with each other via bridges of water molecules and lipid headgroups. Lipid molecules with unusual conformations participate in the H-bond network and are an integral part of the nanopore structure. To explore the role of this H-bond network on membrane selectivity, we swapped three critical polar residues with the nonpolar residues found in melittin. All variants have potency, membrane selectivity, and cytotoxicity that were intermediate between a cytotoxic melittin variant called MelP5 and the macrolittins. Simulations showed that the variants had less organized H-bond networks of waters and lipids with unusual structures. The membrane-spanning, cooperative H-bond network is a critical determinant of macrolittin nanopore stability and membrane selectivity. The results described here will help guide the future design and optimization of peptide nanopore-based applications.

Allosteric competition and inhibition in AMPA receptors

Excitatory neurotransmission is principally mediated by α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA)-subtype ionotropic glutamate receptors (AMPARs). Negative allosteric modulators are therapeutic candidates that inhibit AMPAR activation and can compete with positive modulators to control AMPAR function through unresolved mechanisms. Here we show that allosteric inhibition pushes AMPARs into a distinct state that prevents both activation and positive allosteric modulation. We used cryo-electron microscopy to capture AMPARs bound to glutamate, while a negative allosteric modulator, GYKI-52466, and positive allosteric modulator, cyclothiazide, compete for control of the AMPARs. GYKI-52466 binds in the ion channel collar and inhibits AMPARs by decoupling the ligand-binding domains from the ion channel. The rearrangement of the ligand-binding domains ruptures the cyclothiazide site, preventing positive modulation. Our data provide a framework for understanding allostery of AMPARs and for rational design of therapeutics targeting AMPARs in neurological diseases.

Discovery of novel flavonoid derivatives as potential dual inhibitors against α-glucosidase and α-amylase: virtual screening, synthesis, and biological evaluation

Diabetes mellitus is one of the top ten causes of death worldwide, accounting for 6.7 million deaths in 2021, and is one of the most rapidly growing global health emergencies of this century. Although several classes of therapeutic drugs have been invented and applied in clinical practice, diabetes continues to pose a serious and growing threat to public health and places a tremendous burden on those affected and their families. The strategy of reducing carbohydrate digestibility by inhibiting the activities of α-glucosidase and α-amylase is regarded as a promising preventative treatment for type 2 diabetes. In this study, we investigated the dual inhibitory effect against two polysaccharide hydrolytic enzymes of flavonoid derivatives from an in-house chemical database. By combining molecular docking and structure-activity relationship analysis, twelve compounds with docking energies less than or equal to - 8.0 kcal mol-1 and containing required structural features for dual inhibition of the two enzymes were identified and subjected to chemical synthesis and in vitro evaluation. The obtained results showed that five compounds exhibited dual inhibitory effects on the target enzymes with better IC50 values than the approved positive control acarbose. Molecular dynamics simulations were performed to elucidate the binding of these flavonoids to the enzymes. The predicted pharmacokinetic and toxicological properties suggest that these compounds are viable for further development as type 2 diabetes drugs.

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system's bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

A rare case of uncharacterized autoinflammatory disease: Patient carrying variations in NLRP3 and TNFRSF1A genes

Tumor necrosis factor type 1A receptor-associated periodic syndrome (TRAPS) and cryopyrin-associated autoinflammatory syndrome (CAPS) are rare monogenic autoinflammatory diseases (AIDs) mainly caused by pathogenic variations in the TNFRSF1A and NLRP3 genes, respectively. Here, we describe a unique patient presenting with symptoms overlapping both TRAPS and CAPS, without known pathogenic variants in the respective genes. The patient harbored the p.Val200Met variation in NLRP3 and the p.Ser226Cys variation in TNFRSF1A, prompting us to delve deeper into the functional analysis due to conflicting or inconclusive pathogenicity interpretations of the variants across various databases. Molecular dynamics analysis of the p.Val200Met variation in NLRP3 revealed a rigid conformation in the helical domain 2 subdomain of the NACHT domain. This increased rigidity suggests a potential mechanism by which this variation supports the assembly of the NLRP3 inflammasome. Notably, the patient's peripheral mononuclear blood cells demonstrated an elevated IL-1β response upon lipopolysaccharides (LPS) induction. Subsequent initiation of anti-IL-1β therapy resulted in a significant alleviation of the patient's symptoms, further supporting our hypothesis. We interpret these findings as suggestive of a potential pathophysiological role for the NLPR3 p.Val200Met variation in shaping the patient's clinical phenotype, which was also supported by clinical and genetic analysis of the family. This case underscores the complexity of the genetic landscape in AIDs and highlights the value of combining family genetic and functional data to refine the understanding and management of such challenging cases.

Theoretical Insights into N-Glycoside Bond Cleavage of 5-Carboxycytosine by Thymine DNA Glycosylase: A QM/MM Study

Thymine DNA glycosylase (TDG)-mediated excision of 5-formylcytosine and 5-carboxylcytosine (5-caC) is a critical step in active DNA demethylation. Herein, we employed a combined quantum mechanics/molecular mechanics approach to investigate the reaction mechanism of TDG-catalyzed N-glycosidic bond cleavage of 5-caC. The calculated results show that TDG-catalyzed 5-caC excision follows a concerted (SN2) mechanism in which glycosidic bond dissociation is coupled with nucleophile attack. Protonation of the 5-caC anion contributes to the cleavage of the N-glycoside bond, in which the N3-protonated zwitterion and imino tautomers are more favorable than carboxyl-protonated amino tautomers. This is consistent with the experimental data. Furthermore, our results reveal that the configuration rearrangement process of the protonated 5-caC would lower the stability of the N-glycoside bond and substantially reduce the barrier height for the subsequent C1'-N1 bond cleavage. This should be attributed to the smaller electrostatic repulsion between the leaving base and the negative phosphate group as a result of the structural rearrangement.

Benzamide Trimethoprim Derivatives as Human Dihydrofolate Reductase Inhibitors-Molecular Modeling and In Vitro Activity Study

Human dihydrofolate reductase (hDHFR) is an essential cellular enzyme, and inhibiting its activity is a promising strategy for cancer therapy. We have chosen the trimethoprim molecule (TMP) as a model compound in our search for a new class of hDHFR inhibitors. We incorporated an amide bond, a structural element typical of netropsin, a ligand that binds selectively in the minor groove of DNA, into the molecules of TMP analogs. In this work, we present previously obtained and evaluated eleven benzamides (JW1-JW8; MB1, MB3, MB4). Recently, these compounds were specifically projected as potential inhibitors of the enzymes acetylcholinesterase (AChE) and β-secretase (BACE1). JW8 was most active against AChE, with an inhibitory concentration of AChE IC50 = 0.056 µM, while the IC50 for donepezil was 0.046 µM. This compound was also the most active against the BACE1 enzyme. The IC50 value was 9.01 µM compared to that for quercetin, with IC50 = 4.89 µM. All the benzamides were active against hDHFR, with IC50 values ranging from 4.72 to 20.17 µM, and showed activity greater than TMP (55.26 µM). Quantitative results identified the derivatives JW2 and JW8 as the most promising. A molecular modeling study demonstrates that JW2 interacts strongly with the key residue Gly-117, while JW8 interacts strongly with Asn-64 and Arg-70. Furthermore, JW2 and JW8 demonstrate the ability to stabilize the hDHFR enzyme, despite forming fewer hydrogen bonds with the protein compared to reference ligands. It can be concluded that this class of compounds certainly holds great promise for good active leads in medicinal chemistry.

Membrane binding and lipid-protein interaction of the C2 domain from coagulation factor V

Anchoring of coagulation factors to anionic regions of the membrane involves the C2 domain as a key player. The rate of enzymatic reactions of the coagulation factors is increased by several orders of magnitude upon membrane binding. However, the precise mechanisms behind the rate acceleration remain unclear, primarily because of a lack of understanding of the conformational dynamics of the C2-containing factors and corresponding complexes. We elucidate the membrane-bound form of the C2 domain from human coagulation factor V (FV-C2) by characterizing its membrane binding the specific lipid-protein interactions. Employing all-atom molecular dynamics simulations and leveraging the highly mobile membrane-mimetic (HMMM) model, we observed spontaneous binding of FV-C2 to a phosphatidylserine (PS)-containing membrane within 2-25 ns across twelve independent simulations. FV-C2 interacted with the membrane through three loops (spikes 1-3), achieving a converged, stable orientation. Multiple HMMM trajectories of the spontaneous membrane binding provided extensive sampling and ample data to examine the membrane-induced effects on the conformational dynamics of C2 as well as specific lipid-protein interactions. Despite existing crystal structures representing presumed "open" and "closed" states of FV-C2, our results revealed a continuous distribution of structures between these states, with the most populated structures differing from both "open" and "closed" states observed in crystal environments. Lastly, we characterized a putative PS-specific binding site formed by K23, Q48, and S78 located in the groove enclosed by spikes 1-3 (PS-specificity pocket), suggesting a different orientation of a bound headgroup moiety compared to previous proposals based upon analysis of static crystal structures.

The impact of hydroxyapatite crystal structures and protein interactions on bone's mechanical properties

Hydroxyapatite (HAP) constitutes the primary mineral component of bones, and its crystal structure, along with the surface interaction with proteins, significantly influences the outstanding mechanical properties of bone. This study focuses on natural hydroxyapatite, constructing a surface model with calcium vacancy defects. Employing a representative model of aspartic acid residues, we delve into the adsorption mechanism on the crystal surface and scrutinize the adsorption forms of amino acid residues on HAP and calcium-deficient hydroxyapatite (CDHA) surfaces. The research also explores the impact of different environments on adsorption energy. Furthermore, a simplified sandwich structure of crystal-polypeptide-crystal is presented, analyzing the distribution of amino acid residue adsorption sites on the crystal surface of the polypeptide fragment. This investigation aims to elucidate how the stick-slip mechanism of polypeptide molecules on the crystal surface influences the mechanical properties of the system. By uncovering the interface mechanical behavior between HAP and osteopontin peptides, this article offers valuable theoretical insights for the construction and biomimetic design of biocomposites.

In vitro and in silico docking and molecular dynamic of antimicrobial activities, alpha-glucosidase, and anti-inflammatory activity of compounds from the aerial parts of Mussaenda saigonensis

Twelve compounds were isolated from Mussaenda saigonensis aerial parts through phytochemical analysis and the genus Mussaenda is the first place where the compounds 4-6 and 11-12 have been found. Based on the ability to inhibit NO production in RAW264.7 cells, compound 2 has demonstrated the strongest anti-inflammatory activity in vitro with an IC50 of 7.6 μM, as opposed to L-NMMA's IC50 of 41.3 μM. Compound 12 was found to be the most effective inhibitor of alpha-glucosidase enzyme in vitro, with an IC50 value of 42.4 μM (compared to 168 μM for acarbose). Compounds 1-12 were evaluated in vitro for antimicrobial activity using the paper dish method. Compound 11 demonstrated strong antifungal activity against M. gypseum with a MIC value of 50 μM. In silico docking for antimicrobial activity, pose 90 or compound 11 docked well to the 2VF5 enzyme, PDB, which explains why compound 11 had the highest activity in vitro. Entry 2/pose 280 demonstrated excellent anti-inflammatory activity in silico. The stability of the complex between pose 280 and the 4WCU enzyme for anti-inflammatory activity has been assessed using molecular dynamics over a simulation course ranging from 0 to 100 ns. It has been found to be stable from 60 and 100 ns. The Tyr 159 (95%, H-bond via water bridge), Asp 318 (200%, multiple contacts), Met 273 (75%, hydrophobic interaction via water bridge), and Gln 369 (75%, H-bond via water bridge) interacted well within the time range of 0 to 100 ns. It has more hydrophilic or polar pharmacokinetics.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Using atomistic solution scattering modelling to elucidate the role of the Fc glycans in human IgG4

Human immunoglobulin G (IgG) exists as four subclasses IgG1-4, each of which has two Fab subunits joined by two hinges to a Fc subunit. IgG4 has the shortest hinge with 12 residues. The Fc subunit has two glycan chains, but the importance of glycosylation is not fully understood in IgG4. Here, to evaluate the stability and structure of non-glycosylated IgG4, we performed a multidisciplinary structural study of glycosylated and deglycosylated human IgG4 A33 for comparison with our similar study of human IgG1 A33. After deglycosylation, IgG4 was found to be monomeric by analytical ultracentrifugation; its sedimentation coefficient of 6.52 S was reduced by 0.27 S in reflection of its lower mass. X-ray and neutron solution scattering showed that the overall Guinier radius of gyration RG and its cross-sectional values after deglycosylation were almost unchanged. In the P(r) distance distribution curves, the two M1 and M2 peaks that monitor the two most common distances within IgG4 were unchanged following deglycosylation. Further insight from Monte Carlo simulations for glycosylated and deglycosylated IgG4 came from 111,382 and 117,135 possible structures respectively. Their comparison to the X-ray and neutron scattering curves identified several hundred best-fit models for both forms of IgG4. Principal component analyses showed that glycosylated and deglycosylated IgG4 exhibited different conformations from each other. Within the constraint of unchanged RG and M1-M2 values, the glycosylated IgG4 models showed more restricted Fc conformations compared to deglycosylated IgG4, but no other changes. Kratky plots supported this interpretation of greater disorder upon deglycosylation, also observed in IgG1. Overall, these more variable Fc conformations may demonstrate a generalisable impact of deglycosylation on Fc structures, but with no large conformational changes in IgG4 unlike those seen in IgG1.

Structural basis for partial agonism in 5-HT3A receptors

Hyperactivity of serotonin 3 receptors (5-HT3R) underlies pathologies associated with irritable bowel syndrome and chemotherapy-induced nausea and vomiting. Setrons, a class of high-affinity competitive antagonists, are used in the treatment of these conditions. Although generally effective for chemotherapy-induced nausea and vomiting, the use of setrons for treating irritable bowel syndrome has been impaired by adverse side effects. Partial agonists are now being considered as an alternative strategy, with potentially less severe side effects than full antagonists. However, a structural understanding of how these ligands work is lacking. Here, we present high-resolution cryogenic electron microscopy structures of the mouse 5-HT3AR in complex with partial agonists (SMP-100 and ALB-148471) captured in pre-activated and open-like conformational states. Molecular dynamics simulations were used to assess the stability of drug-binding poses and interactions with the receptor over time. Together, these studies reveal mechanisms for the functional differences between orthosteric partial agonists, full agonists and antagonists of the 5-HT3AR.

Interaction of Solanum tuberosum L. translation initiation factors eIF4E with potato virus Y VPg: Apprehend and avoid

Potato virus Y (PVY) is one of the most dangerous agricultural pathogens that causes substantial harm to vegetative propagated crops, such as potatoes (Solanum tuberosum L.). A necessary condition for PVY infection is an interaction between the plant cap-binding translation initiation factors eIF4E and a viral protein VPg, which mimics the cap-structure. In this study, we identified the point mutations in potato eIF4E1 and eIF4E2 that disrupt VPg binding while preserving the functional activity. For the structural interpretation of the obtained results, molecular models of all the studied forms of eIF4E1 and eIF4E2 were constructed and analyzed via molecular dynamics. The results of molecular dynamics simulations corresponds to the biochemical results and suggests that the β1β2 loop plays a key role in the stabilization of both eIF4E-cap and eIF4E-VPg complexes.

Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields

Molecular dynamics simulations play a pivotal role in elucidating the dynamic behaviors of RNA structures, offering a valuable complement to traditional methods such as nuclear magnetic resonance or X-ray. Despite this, the current precision of RNA force fields lags behind that of protein force fields. In this work, we systematically compared the performance of four RNA force fields (ff99bsc0χOL3, AMBERDES, ff99OL3_CMAP1, AMBERMaxEnt) across diverse RNA structures. Our findings highlight significant challenges in maintaining stability, particularly with regard to cross-strand and cross-loop hydrogen bonds. Furthermore, we observed the limitations in accurately describing the conformations of nonhelical structural motif, terminal nucleotides, and also base pairing and base stacking interactions by the tested RNA force fields. The identified deficiencies in existing RNA force fields provide valuable insights for subsequent force field development. Concurrently, these findings offer recommendations for selecting appropriate force fields in RNA simulations.

Computational study of the mechanism of binding of antifungal icofungipen in the active site of eukaryotic isoleucyl tRNA synthetase from Candida albicans

The eukaryotic fungal species Candida albicans is a critical infective pathogenic agent. The β-amino acid, Icofungipen, is an effective inhibitor of Candida albicans. Icofungipen binds at the active site of the isoleucyl tRNA synthetase (IleRS) from Candida albicans (CaIleRS) and halts protein translation in fungus. In the present work, we have investigated the mechanism of binding of Icogungipen (abbreviated as IFP). Molecular dynamics (MD) simulations show that the carboxylic acid group of IFP in the CaIleRS: IFP complex is more oriented towards the Connective Polypeptide (CP) core loop compared to the carboxylic acid group of Ile in the CaIleRS: Ile complex. The Arg 410 of the CP core loop near the substrate is extended towards the IFP. Due to the difference in the conformation of residues of the CP core loop, the KMSKR loop is more proximal to the CP core loop in CaIleRS: IFP. The editing domain which is covalently linked with the CP core loop in the CaIleRS: IFP complex is also oriented in such a way that the active site cavity is narrow and longer. The metadynamics calculation shows that the IFP is trapped in a deeper potential well compared to Ile which is due to the effective closure of the gateway of the active site by KMSKR and CP core loop. The thin, long shape of the active site and the closed gate of the active site in CaIleRS: IFP complex is responsible for the effective capture of IFP relative to Ile in the active site.Communicated by Ramaswamy H. Sarma.

Multiscale Study on the Intramolecular C-S Bond Formation Catalyzed by P450 Monooxygenase CxnD Involved in the Biosynthesis of Chuangxinmycin: The Critical Roles of Noncrystal Water Molecule and Conformational Change

Cytochrome P450 monooxygenase CxnD catalyzes intramolecular C-S bond formation in the biosynthesis of chuangxinmycin, which is representative of the synthesis of sulfur-containing natural heterocyclic compounds. The intramolecular cyclization usually requires the activation of two reaction sites and a large conformational change; thus, illuminating its detailed reaction mechanism remains challengeable. Here, the reaction pathway of CxnD-catalyzed C-S bond formation was clarified by a series of calculations, including Gaussian accelerated molecular dynamics simulations and quantum mechanical-molecular mechanical calculations. Our results revealed that the C-S formation follows a diradical coupling mechanism. CxnD first employs Cpd I to abstract the hydrogen atom from the imino group of the indole ring, and then, the resulted Cpd II further extracts another hydrogen atom from the thiol group of the side chain to afford a diradical intermediate, in which a noncrystal water molecule entering into the active site after the formation of Cpd I was proved to play an indispensable role. Moreover, the diradical intermediate cannot directly perform the coupling reaction. It should first undergo a series of conformational changes leading to the proximity of two reaction sites. It is the flexibility of the active site of the enzyme and the side chain of the substrate that makes the diradical coupling to be successful.

Hydration numbers of biologically relevant divalent metal cations from ab initio molecular dynamics and continuum solvation methods

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca2+, Cd2+, Cu2+, Fe2+, Hg2+, Mg2+, Ni2+, Pb2+, and Zn2+, we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M(H2O)2]2+ aqua complexes to [M(H2O)9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M(H2O)n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion-water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data.