Peptide-based allosteric inhibitor targets TNFR1 conformationally active region and disables receptor-ligand signaling complex.

Tumor necrosis factor (TNF) receptor 1 (TNFR1) plays a pivotal role in mediating TNF induced downstream signaling and regulating inflammatory response. Recent studies have suggested that TNFR1 activation involves conformational rearrangements of preligand assembled receptor dimers and targeting receptor conformational dynamics is a viable strategy to modulate TNFR1 signaling. Here, we used a combination of biophysical, biochemical, and cellular assays, as well as molecular dynamics simulation to show that an anti-inflammatory peptide (FKCRRWQWRMKK), which we termed FKC, inhibits TNFR1 activation allosterically by altering the conformational states of the receptor dimer without blocking receptor-ligand interaction or disrupting receptor dimerization. We also demonstrated the efficacy of FKC by showing that the peptide inhibits TNFR1 signaling in HEK293 cells and attenuates inflammation in mice with intraperitoneal TNF injection. Mechanistically, we found that FKC binds to TNFR1 cysteine-rich domains (CRD2/3) and perturbs the conformational dynamics required for receptor activation. Importantly, FKC increases the frequency in the opening of both CRD2/3 and CRD4 in the receptor dimer, as well as induces a conformational opening in the cytosolic regions of the receptor. This results in an inhibitory conformational state that impedes the recruitment of downstream signaling molecules. Together, these data provide evidence on the feasibility of targeting TNFR1 conformationally active region and open new avenues for receptor-specific inhibition of TNFR1 signaling.

Robust Mechanical Destruction to the Cell Membrane of Carbon Nitride Polyaniline (C3N): A Molecular Dynamics Simulation Study.

The drug-resistant bacteria, particularly multidrug-resistant bacteria, has emerged as a major global public health concern posing serious threats to human life and survival. Nanomaterials, including graphene, have shown promise as effective antibacterial agents owing to their unique antibacterial mechanism compared with traditional drugs. Despite the structural similarity to graphene, the potential antibacterial activity of carbon nitride polyaniline (C3N) remains unexplored. In this study, we employed molecular dynamics simulations to investigate the effects of the interaction between the C3N nanomaterial and the bacterial membrane to evaluate the potential antibacterial activity of C3N. Our results suggest that C3N is capable of inserting deep into the bacterial membrane interior, regardless of the presence or absence of positional restraints in the C3N. The insertion process also resulted in local lipid extraction by the C3N sheet. Additional structural analyses revealed that C3N induced significant changes in membrane parameters, including mean square displacement, deuterium order parameters, membrane thickness, and area per lipid. Docking simulations, where all the C3N are restraint to a specific positions, confirmed that C3N can extract lipids from the membrane, indicating the strong interaction between the C3N material and the membrane. Free-energy calculations further revealed that the insertion of the C3N sheet is energetically favorable and that C3N exhibits membrane insertion capacity comparable to that observed for graphene, suggesting their potential for similar antibacterial activity. This study provides the first evidence of the potential antibacterial properties of C3N nanomaterials via bacterial membrane damage and underscores the potential for its use as antibacterial agents in the future applications.

Single-molecule analysis of ligand efficacy in ß2AR-G-protein activation.

G-protein-coupled receptor (GPCR)-mediated signal transduction is central to human physiology and disease intervention, yet the molecular mechanisms responsible for ligand-dependent signalling responses remain poorly understood. In class A GPCRs, receptor activation and G-protein coupling entail outward movements of transmembrane helix 6 (TM6). Here, using single-molecule fluorescence resonance energy transfer imaging, we examine TM6 movements in the ß2 adrenergic receptor (ß2AR) upon exposure to orthosteric ligands with different efficacies, in the absence and presence of the Gs heterotrimer. We show that partial and full agonists differentially affect TM6 motions to regulate the rate at which GDP-bound ß2AR-Gs complexes are formed and the efficiency of nucleotide exchange leading to Gs activation. These data also reveal transient nucleotide-bound ß2AR-Gs species that are distinct from known structures, and provide single-molecule perspectives on the allosteric link between ligand- and nucleotide-binding pockets that shed new light on the G-protein activation mechanism.

Structural Basis of the Binding Mode of the Antineoplastic Compound Motixafortide (BL-8040) in the CXCR4 Chemokine Receptor.

Modulation of the CXCL12-CXCR4 signaling axis is of the utmost importance due to its central involvement in several pathological disorders, including inflammatory diseases and cancer. Among the different currently available drugs that inhibit CXCR4 activation, motixafortide-a best-in-class antagonist of this GPCR receptor-has exhibited promising results in preclinical studies of pancreatic, breast, and lung cancers. However, detailed information on the interaction mechanism of motixafortide is still lacking. Here, we characterize the motixafortide/CXCR4 and CXCL12/CXCR4 protein complexes by using computational techniques including unbiased all-atom molecular dynamics simulations. Our microsecond-long simulations of the protein systems indicate that the agonist triggers changes associated with active-like GPCR conformations, while the antagonist favors inactive conformations of CXCR4. Detailed ligand-protein analysis indicates the importance of motixafortide's six cationic residues, all of which established charge-charge interactions with acidic CXCR4 residues. Furthermore, two synthetic bulky chemical moieties of motixafortide work in tandem to restrict the conformations of important residues associated with CXCR4 activation. Our results not only elucidate the molecular mechanism by which motixafortide interacts with the CXCR4 receptor and stabilizes its inactive states, but also provide essential information to rationally design CXCR4 inhibitors that preserve the outstanding pharmacological features of motixafortide.

Identification of Molecular Determinants in iRhoms1 and 2 That Contribute to the Substrate Selectivity of Stimulated ADAM17.

The metalloprotease ADAM17 is a key regulator of the TNFα, IL-6R and EGFR signaling pathways. The maturation and function of ADAM17 is controlled by the seven-membrane-spanning proteins iRhoms1 and 2. The functional properties of the ADAM17/iRhom1 and ADAM17/iRhom2 complexes differ, in that stimulated shedding of most ADAM17 substrates tested to date can be supported by iRhom2, whereas iRhom1 can only support stimulated shedding of very few ADAM17 substrates, such as TGFα. The first transmembrane domain (TMD1) of iRhom2 and the sole TMD of ADAM17 are important for the stimulated shedding of ADAM17 substrates by iRhom2. However, little is currently known about how the iRhoms interact with different substrates to control their stimulated shedding by ADAM17. To provide new insights into this topic, we tested how various chimeras between iRhom1 and iRhom2 affect the stimulated processing of the EGFR-ligands TGFα (iRhom1- or 2-dependent) and EREG (iRhom2-selective) by ADAM17. This uncovered an important role for the TMD7 of the iRhoms in determining their substrate selectivity. Computational methods utilized to characterize the iRhom1/2/substrate interactions suggest that the substrate selectivity is determined, at least in part, by a distinct accessibility of the substrate cleavage site to stimulated ADAM17. These studies not only provide new insights into why the substrate selectivity of stimulated iRhom2/ADAM17 differs from that of iRhom1/ADAM17, but also suggest new approaches for targeting the release of specific ADAM17 substrates.

Theoretical Study of the Structural Stability, Chemical Reactivity, and Protein Interaction for NMP Compounds as Modulators of the Endocannabinoid System.

The cannabinoid receptors (CB1/CB2) and the T-type calcium channels are involved in disorders associated with both physiological pain and depressive behaviors. Valuable pharmacological species carbazole derivatives such as the NMP-4, NMP-7, and NMP-181 (Neuro Molecular Production) regulate both biological entities. In this work, DFT calculations were performed to characterize theoretically their structural and chemical reactivity properties using the BP86/cc-pVTZ level of theory. The molecular orbital contributions and the chemical reactivity analysis reveal that a major participation of the carbazole group is in the donor-acceptor interactions of the NMP compounds. The DFT analysis on the NMP compounds provides insights into the relevant functional groups involved during the ligand-receptor interactions. Molecular docking analysis is used to reveal possible sites of interaction of the NMP compounds with the Cav3.2 calcium channel. The interaction energy values and reported experimental evidence indicate that the site denominated as "Pore-blocking", which is formed mainly by hydrophobic residues and the T586 residue, is a probable binding site for the NMP compounds.

Novel variants of engineered water soluble mu opioid receptors with extensive mutations and removal of cysteines.

We have shown that water-soluble variants of the human mu opioid receptor (wsMOR) containing a reduced number of hydrophobic residues at the lipid-facing residues of the transmembrane (TM) helices can be expressed in E. coli. In this study, we tested the consequences of increasing the number of mutations on the surface of the transmembrane domain on the receptor's aqueous solubility and ligand binding properties, along with mutation of 11 cysteine residues regardless of their solvent exposure value and location in the protein. We computationally engineered 10 different variants of MOR, and tested four of them for expression in E. coli. We found that all four variants were successfully expressed and could be purified in high quantities. The variants have alpha helical structural content similar to that of the native MOR, and they also display binding affinities for the MOR antagonist (naltrexone) similar to the wsMOR variants we engineered previously that contained many fewer mutations. Furthermore, for these full-length variants, the helical content remains unchanged over a wide range of pH values (pH 6 ~ 9). This study demonstrates the flexibility and robustness of the water-soluble MOR variants with respect to additional designed mutations in the TM domain and changes in pH, whereupon the protein's structural integrity and its ligand binding affinity are maintained. These variants of the full-length MOR with less hydrophobic surface residues and less cysteines can be obtained in large amounts from expression in E. coli and can serve as novel tools to investigate structure-function relationships of the receptor.

Substrate-selective protein ectodomain shedding by ADAM17 and iRhom2 depends on their juxtamembrane and transmembrane domains.

The metalloprotease ADAM17 (a disintegrin and metalloprotease 17) regulates EGF-receptor and TNFα signaling, thereby not only protecting the skin and intestinal barrier, but also contributing to autoimmunity. ADAM17 can be rapidly activated by many stimuli through its transmembrane domain (TMD), with the seven membrane-spanning inactive Rhomboids (iRhom) 1 and 2 implicated as candidate regulatory partners. However, several alternative models of ADAM17 regulation exist that do not involve the iRhoms, such as regulation through disulfide bond exchange or through interaction with charged phospholipids. Here, we report that a non-activatable mutant of ADAM17 with the TMD of betacellulin (BTC) can be rescued by restoring residues from the ADAM17 TMD, but only in Adam17-/- cells, which contain iRhoms, not in iRhom1/2-/- cells. We also provide the first evidence that the extracellular juxtamembrane domains (JMDs) of ADAM17 and iRhom2 regulate the stimulation and substrate selectivity of ADAM17. Interestingly, a point mutation in the ADAM17 JMD identified in a patient with Tetralogy of Fallot, a serious heart valve defect, affects the substrate selectivity of ADAM17 toward Heparin-binding epidermal growth factor like growth factor (HB-EGF), a crucial regulator of heart valve development in mice. These findings provide new insights into the regulation of ADAM17 through an essential interaction with the TMD1 and JMD1 of iRhom2.

The C-terminus of the mu opioid receptor is critical in G-protein interaction as demonstrated by a novel graphene biosensor.

Several water-soluble variants of the human mu opioid receptor (wsMORs) have been designed and expressed, which enables the detection of opioids in the nM to pM range using biosensing platforms. The tools previously developed allowed us to investigate MOR and G-protein interactions in a lipid free system to demonstrate that the lipid bilayer might not be essential for the G-protein recognition and binding. In this study, we are able to investigate G-protein interactions with MOR by using graphene enabled technology, in a lipid free system, with a high sensitivity in a real time manner. A new wsMOR with the native C-terminus was designed, expressed and then immobilized on the surfaces of scalable graphene field effect transistor (GFET)-based biosensors, enabling the recording of wsMOR/G-protein interaction with an electronic readout. G-protein only interacts with the wsMOR in the presence of the native MOR C-terminus with a KA of 32.3±11.1 pM. The electronic readout of such interaction is highly reproducible with little variance across 50 devices in one biosensor array. For devices with receptors that do not have the native C-terminus, no significant electronic response was observed in the presence of G-protein, indicating an absence of interaction. These findings reveal that lipid environment is not essential for the G-protein interaction with MOR, however, the C-terminus of MOR is essential for G-protein recognition and high affinity binding. A system to detect MOR-G protein interaction is developed. wsMOR-G2_Cter provides a novel tool to investigate the role of C terminus in the signaling pathway.

Structural modeling defines transmembrane residues in ADAM17 that are crucial for Rhbdf2-ADAM17-dependent proteolysis.

A disintegrin and metalloproteinase 17 (ADAM17) controls the release of the pro-inflammatory cytokine tumor necrosis factor α (TNFα, also known as TNF) and is crucial for protecting the skin and intestinal barrier by proteolytic activation of epidermal growth factor receptor (EGFR) ligands. The seven-membrane-spanning protein called inactive rhomboid 2 (Rhbdf2; also known as iRhom2) is required for ADAM17-dependent TNFα shedding and crosstalk with the EGFR, and a point mutation (known as sinecure, sin) in the first transmembrane domain (TMD) of Rhbdf2 (Rhbdf2sin) blocks TNFα shedding, yet little is known about the underlying mechanism. Here, we used a structure-function analysis informed by structural modeling to evaluate the interaction between the TMD of ADAM17 and the first TMD of Rhbdf2, and the role of this interaction in Rhbdf2-ADAM17-dependent shedding. Moreover, we show that double mutant mice that are homozygous for Rhbdf2sin/sin and lack Rhbdf1 closely resemble Rhbdf1/2-/- double knockout mice, highlighting the severe functional impact of the Rhbdf2sin/sin mutation on ADAM17 during mouse development. Taken together, these findings provide new mechanistic and conceptual insights into the critical role of the TMDs of ADAM17 and Rhbdf2 in the regulation of the ADAM17 and EGFR, and ADAM17 and TNFα signaling pathways.

Different protonated states at the C-terminal of the amyloid-ß peptide modulate the stability of S-shaped protofibril.

Studies have found strong correlations between polymorphism and structural variations in amyloid-ß (Aß) fibrils and the diverse clinical subtypes of Alzheimer's disease (AD). Thus, a detailed understanding of the conformational behavior of Aß fibrils may be an aid to elucidate the pathological mechanisms involved in AD. However, a key point that has been inadvertently underestimated or dismissed is the role of the protonated state at the C-terminal residue of amyloid-ß peptides, which can give rise to intrinsic differences in the morphology and stability of the fibrils. For instance, the effects of the salt bridge formed between the C-terminal residue A42 and the residue K28 on the S-shaped Aß protofibril structure remain unknown and may be different from those in the U-shaped Aß protofibril structures. To address this effect, we explore the stability of the S-shaped protofibrils capped with different C-terminal modifications, including carboxyl group in its deprotonated (COO-) and protonated (COOH) states, by using molecular dynamics simulations. Our findings indicated that the C-terminal deprotonated protofibril is significantly more stable than its C-terminal protonated counterpart due to a well-defined and highly stable zipper-like salt-bridge-chain formed by the ε-NH3 + groups on the sidechain of residue K28 and the C-terminal COO- group at the A42 residue. The revealed underlying molecular mechanism for the different stability of the protofibrils provides insights into the diversity of polymorphism in Aß fibrils.

Mild Binding of Protein to C2 N Monolayer Reveals Its Suitable Biocompatibility.

The development of biocompatible nanomaterials for smart drug delivery and bioimaging has attracted great interest in recent years in biomedical fields. Here, the interaction between the recently reported nitrogenated graphene (C2 N) and a prototypical protein (villin headpiece HP35) utilizing atomistic molecular dynamics simulations is studied. The simulations reveal that HP35 can form a stable binding with the C2 N monolayer. Although the C2 N-HP35 attractive interactions are constantly preserved, the binding strength between C2 N and the protein is mild and does not cause significant distortion in the protein's structural integrity. This intrinsic biofriendly property of native C2 N is distinct from several widely studied nanomaterials, such as graphene, carbon nanotubes, and MoS2 , which can induce severe protein denaturation. Interestingly, once the protein is adsorbed onto C2 N surface, its transverse migration is highly restricted at the binding sites. This restriction is orchestrated by C2 N's periodic porous structure with negatively charged "holes," where the basic residues-such as lysine-can form stable interactions, thus functioning as "anchor points" in confining the protein displacement. It is suggested that the mild, immobilized protein attraction and biofriendly aspects of C2 N would make it a prospective candidate in bio- and medical-related applications.

Photoinduced Electron Transfer Elicits a Change in the Static Dielectric Constant of a de Novo Designed Protein.

We provide a direct measure of the change in effective dielectric constant (ε(S)) within a protein matrix after a photoinduced electron transfer (ET) reaction. A linked donor-bridge-acceptor molecule, PZn-Ph-NDI, consisting of a (porphinato)Zn donor (PZn), a phenyl bridge (Ph), and a naphthalene diimide acceptor (NDI), is shown to be a "meter" to indicate protein dielectric environment. We calibrated PZn-Ph-NDI ET dynamics as a function of solvent dielectric, and computationally de novo designed a protein SCPZnI3 to bind PZn-Ph-NDI in its interior. Mapping the protein ET dynamics onto the calibrated ET catalogue shows that SCPZnI3 undergoes a switch in the effective dielectric constant following photoinduced ET, from ε(S) ≈ 8 to ε(S) ≈ 3.

Scalable production of highly sensitive nanosensors based on graphene functionalized with a designed G protein-coupled receptor.

We have developed a novel, all-electronic biosensor for opioids that consists of an engineered µ-opioid receptor protein, with high binding affinity for opioids, chemically bonded to a graphene field-effect transistor to read out ligand binding. A variant of the receptor protein that provided chemical recognition was computationally redesigned to enhance its solubility and stability in an aqueous environment. A shadow mask process was developed to fabricate arrays of hundreds of graphene transistors with average mobility of ∼1500 cm(2) V(-1) s(-1) and yield exceeding 98%. The biosensor exhibits high sensitivity and selectivity for the target naltrexone, an opioid receptor antagonist, with a detection limit of 10 pg/mL.

A functional selectivity mechanism at the serotonin-2A GPCR involves ligand-dependent conformations of intracellular loop 2.

With recent progress in determination of G protein-coupled receptor (GPCR) structure with crystallography, a variety of other experimental approaches (e.g., NMR spectroscopy, fluorescent-based assays, mass spectrometry techniques) are also being used to characterize state-specific and ligand-specific conformational states. MD simulations offer a powerful complementary approach to elucidate the dynamic features associated with ligand-specific GPCR conformations. To shed light on the conformational elements and dynamics of the important aspect of GPCR functional selectivity, we carried out unbiased microsecond-length MD simulations of the human serotonin 2A receptor (5-HT(2A)R) in the absence of ligand and bound to four distinct serotonergic agonists. The 5-HT(2A)R is a suitable system to study the structural features involved in the ligand-dependent conformational heterogeneity of GPCRs because it is well-characterized experimentally and exhibits a strong agonist-specific phenotype in that some 5-HT(2A)R agonists induce LSD-like hallucinations, while others lack this psychoactive property entirely. Here we report evidence for structural and dynamic differences in 5-HT(2A)R interacting with such pharmacologically distinct ligands, hallucinogens, and nonhallucinogens obtained from all-atom MD simulations. Differential ligand binding contacts were identified for structurally similar hallucinogens and nonhallucinogens and found to correspond to different conformations in the intracellular loop 2 (ICL2). From the different ICL2 conformations, functional selective phenotypes are suggested through effects on dimerization and/or distinct direct interaction with effector proteins. The findings are presented in the context of currently proposed hallucinogenesis mechanisms, and ICL2 is proposed as a fine-tuning selective switch that can differentiates modes of 5-HT(2A)R activation.

Characterization of a computationally designed water-soluble human µ-opioid receptor variant using available structural information.

BACKGROUND: The recent X-ray crystal structure of the murine µ-opioid receptor (MUR) allowed the authors to reengineer a previously designed water-soluble variant of the transmembrane portion of the human MUR (wsMUR-TM). METHODS: The new variant of water-soluble MUR (wsMUR-TM_v2) was engineered based on the murine MUR crystal structure. This novel variant was expressed in Escherichia coli and purified. The properties of the receptor were characterized and compared with those of wsMUR-TM. RESULTS: Seven residues originally included for mutation in the design of the wsMUR-TM were reverted to their native identities. wsMUR-TM_v2 contains 16% mutations of the total sequence. It was overexpressed and purified with high yield. Although dimers and higher oligomers were observed to form over time, the wsMUR-TM_v2 stayed predominantly monomeric at concentrations as high as 7.5 mg/ml in buffer within a 2-month period. Its secondary structure was predominantly helical and comparable with those of both the original wsMUR-TM variant and the native MUR. The binding affinity of wsMUR-TM_v2 for naltrexone (K(d) approximately 70 nM) was in close agreement with that for wsMUR-TM. The helical content of wsMUR-TM_v2 decreased cooperatively with increasing temperature, and the introduction of sucrose was able to stabilize the protein. CONCLUSIONS: A novel functional wsMUR-TM_v2 with only 16% mutations was successfully engineered, expressed in E. coli, and purified based on information from the crystal structure of murine MUR. This not only provides a novel alternative tool for MUR studies in solution conditions but also offers valuable information for protein engineering and structure-function relations.

Molecular dynamics simulations suggest the potential toxicity of fluorinated graphene to HP35 protein via unfolding the α-helix structure.

Fluorinated graphene, a two-dimensional nanomaterial composed of three atomic layers, a central carbon layer sandwiched between two layers of fluorine atoms, has attracted considerable attention across various fields, particularly for its potential use in biomedical applications. Nonetheless, scant effort has been devoted to assessing the potential toxicological implications of this nanomaterial. In this study, we scrutinize the potential impact of fluorinated graphene on a protein model, HP35 by utilizing extensive molecular dynamics (MD) simulation methods. Our MD results elucidate that upon adsorption to the nanomaterial, HP35 undergoes a denaturation process initiated by the unraveling of the second helix of the protein and the loss of the proteins hydrophobic core. In detail, substantial alterations in various structural features of HP35 ensue, including alterations in hydrogen bonding, Q value, and RMSD. Subsequent analyses underscore that hydrophobic and van der Waals interactions (predominant), alongside electrostatic energy (subordinate), exert influence over the adsorption of HP35 on the fluorinated graphene surface. Mechanistic scrutiny attests that the unrestrained lateral mobility of HP35 on the fluorinated graphene nanomaterial primarily causes the exposure of HP35's hydrophobic core, resulting in the eventual structural denaturation of HP35. A trend in the features of 2D nanostructures is proposed that may facilitate the denaturation process. Our findings not only substantiate the potential toxicity of fluorinated graphene but also unveil the underlying molecular mechanism, which thereby holds significance for the prospective utilization of such nanomaterials in the field of biomedicine.

Villin headpiece unfolding upon binding to boridene mediated by the "anchoring-perturbation" mechanism.

We employ molecular dynamics (MD) simulations to investigate the influence of boridene on the behavior of a protein model, HP35, with the aim of assessing the potential biotoxicity of boridene. Our MD results reveal that HP35 can undergo unfolding via an "anchoring-perturbation" mechanism upon adsorption onto the boridene surface. Specifically, the third helix of HP35 becomes tightly anchored to the boridene surface through strong electrostatic interactions between the abundant molybdenum atoms on the boridene surface and the oxygen atoms on the HP35 backbone. Meanwhile, the first helix, experiencing continuous perturbation from the surrounding water solution over an extended period, suffers from potential breakage of hydrogen bonds, ultimately resulting in its unfolding. Our findings not only propose, for the first time to our knowledge, the "anchoring-perturbation" mechanism as a guiding principle for protein unfolding but also reveal the potential toxicity of boridene on protein structures.

NMR structure and dynamics of a designed water-soluble transmembrane domain of nicotinic acetylcholine receptor.

The nicotinic acetylcholine receptor (nAChR) is an important therapeutic target for a wide range of pathophysiological conditions, for which rational drug designs often require receptor structures at atomic resolution. Recent proof-of-concept studies demonstrated a water-solubilization approach to structure determination of membrane proteins by NMR (Slovic et al., PNAS, 101: 1828-1833, 2004; Ma et al., PNAS, 105: 16537-42, 2008). We report here the computational design and experimental characterization of WSA, a water-soluble protein with ~83% sequence identity to the transmembrane (TM) domain of the nAChR α1 subunit. Although the design was based on a low-resolution structural template, the resulting high-resolution NMR structure agrees remarkably well with the recent crystal structure of the TM domains of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), demonstrating the robustness and general applicability of the approach. NMR T(2) dispersion measurements showed that the TM2 domain of the designed protein was dynamic, undergoing conformational exchange on the NMR timescale. Photoaffinity labeling with isoflurane and propofol photolabels identified a common binding site in the immediate proximity of the anesthetic binding site found in the crystal structure of the anesthetic-GLIC complex. Our results illustrate the usefulness of high-resolution NMR analyses of water-solubilized channel proteins for the discovery of potential drug binding sites.

The single-atom catalytic activity of the hydrogen evolution reaction of the experimentally synthesized boridene 2D material: a density functional theory study.

CONTEXT: Previous theoretical studies have suggested that two-dimensional (2D) MBene materials might display adequate monatomic catalytic activity for the hydrogen evolution reaction (HER). Recently, a study reported the experimental synthesis of a 2D MBene (Mo4/3B2), re-defined as boridene, albeit no effort has been devoted to explore the single-atom catalytic activity for HER of this experimentally synthesized 2D material. Therefore, we herein investigate the single-atom HER performance of the boridene. Interestingly, with Mo defects mixed with single Au and Zn atoms shows excellent hydrogen evolution performance, and the change in the Gibbs free energy ([Formula: see text]) value is close to 0 eV, which can even match the performance of Pt-based materials. Through analysis of the charge density difference and density of states, the mechanism affecting the HER performance is explained at the electronic level. This work provides a new direction for the use of the Mo4/3B2 monolayer two-dimensional materials in the field of single-atom catalysis for HER. METHODS: This study used the DFT calculations in Vienna ab initio simulation package. The GGA-Perdew-Burke-Ernzerhof functional with DFT-D2 correction is used to describe the exchange-correlation interactions. The projection augmented wave is used with the plane wave cutoff of 500 eV. The convergences of energy and force are 10-5 eV and 0.01 eV/Å, respectively. A vacuum layer with a height of 20 Å is set in the Z direction. For geometry optimization, self-consistent, and DOS calculations, the k-point grids sampled in Brillouin zones are 3 × 3 × 1, 9 × 9 × 1, and 9 × 9 × 1, respectively. The AIMD simulation is performed in the canonical ensemble (NVT), and the temperature was maintained at 300 K by Nosé-Hoover thermostats with a time step of 2.0 fs.