Activation mechanisms of αVβ3 integrin by binding to fibronectin: A computational study

Mechanism of RGD-conjugated nanodevice binding to its target protein integrin αVβ3 by atomistic molecular dynamics and machine learning

Active targeting strategies have been proposed to enhance the selective uptake of nanoparticles (NPs) by diseased cells, and recent experimental findings have proven the effectiveness of this approach. However, no mechanistic studies have yet revealed the atomistic details of the interactions between ligand-activated NPs and integrins. As a case study, here we investigate, by means of advanced molecular dynamics simulations (MD) and machine learning methods (namely equilibrium MD, binding free energy calculations and training of self-organized maps), the interaction of a cyclic-RGD-conjugated PEGylated TiO2 NP (the nanodevice) with the extracellular segment of integrin αVβ3 (the target), the latter experimentally well-known to be over-expressed in several solid tumors. Firstly, we proved that the cyclic-RGD ligand binding to the integrin pocket is established and kept stable even in the presence of the cumbersome realistic model of the nanodevice. In this respect, the unsupervised machine learning analysis allowed a detailed comparison of the ligand/integrin binding in the presence and in the absence of the nanodevice, which unveiled differences in the chemical features. Then, we discovered that unbound cyclic RGDs conjugated to the NP largely contribute to the interactions between the nanodevice and the integrin. Finally, by increasing the density of cyclic RGDs on the PEGylated TiO2 NP, we observed a proportional enhancement of the nanodevice/target binding. All these findings can be exploited to achieve an improved targeting selectivity and cellular uptake, and thus a more successful clinical outcome.

ECM proteins and cationic polymers coating promote dedifferentiation of patient-derived mature adipocytes to stem cells

Reprogramming of mature adipocytes is an attractive research area due to the plasticity of these cells. Mature adipocytes can be reprogrammed in vitro, transforming them into dedifferentiated fat cells (DFATs), which are considered a new type of stem cell, and thereby have a high potential for use in tissue engineering and regenerative medicine. However, there are still no reports or findings on in vitro controlling the dedifferentiation. Although ceiling culture performed in related studies is a relatively simple method, its yield is low and does not allow manipulation of mature adipocytes to increase or decrease the dedifferentiation. In this study, to understand the role of physicochemical surface effects on the dedifferentiation of patient-derived mature adipocytes, the surfaces of cell culture flasks were coated with extracellular matrix, basement membrane proteins, and cationic/anionic polymers. Extracellular matrix such as fibronectin and collagen type I, and basement membrane proteins such as collagen type IV and laminin strongly promoted dedifferentiation of mature adipocytes, with laminin showing the highest effect with a DFAT ratio of 2.98 (±0.84). Interestingly, cationic polymers also showed a high dedifferentiation effect, but anionic polymers did not, and poly(diallyl dimethylammonium chloride) showed the highest DFAT ratio of 2.27 (±2.8) among the cationic polymers. Protein assay results revealed that serum proteins were strongly adsorbed on the surfaces of the cationic polymer coating, including inducing high mature adipocyte adhesion. This study demonstrates for the first time the possibility of regulating the transformation of mature adipocytes to DFAT stem cells by controlling the physicochemical properties of the surface of conventional cell culture flasks.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Molecular and structural basis of interactions of vitamin D3 hydroxyderivatives with aryl hydrocarbon receptor (AhR): An integrated experimental and computational study

To better understand the molecular and structural basis underlying the interaction of vitamin D3 hydroxyderivatives with AhR, molecular simulation was used to probe the binding of 1,20(OH)₂D3, 1,25(OH)₂D3, 20,23(OH)₂D3 and 20(OH)D3 to AhR. qPCR showed that vitamin D3 derivatives stimulate expression of cyp1A1 and cyp1B1 genes that are downstream targets of AhR signaling. These secosteroids stimulated the translocation of the AhR to the nucleus, as measured by flow cytometry and western blotting. Molecular dynamics simulations were used to model the binding of vitamin D3 derivatives to AhR to examine their influence on the structure, conformation and dynamics of the AhR ligand binding domain (LBD). Binding thermodynamics, conformation, secondary structure, dynamical motion and electrostatic potential of AhR were analyzed. The molecular docking scores and binding free energy were all favorable for the binding of D3 derivatives to the AhR. These established ligands and the D3 derivatives are predicted to have different patterns of hydrogen bond formation with the AhR, and varied residue conformational fluctuations and dynamical motion for the LBD. These changes could alter the shape, size and electrostatic potential distribution of the ligand binding pocket, contributing to the different binding affinities of AhR for the natural ligands and D3 derivatives.

Using ELP Repeats as a Scaffold for De Novo Construction of Gadolinium-Binding Domains within Multifunctional Recombinant Proteins for Targeted Delivery of Gadolinium to Tumour Cells

Three artificial proteins that bind the gadolinium ion (Gd3+) with tumour-specific ligands were de novo engineered and tested as candidate drugs for binary radiotherapy (BRT) and contrast agents for magnetic resonance imaging (MRI). Gd3+-binding modules were derived from calmodulin. They were joined with elastin-like polypeptide (ELP) repeats from human elastin to form the four-centre Gd3+-binding domain (4MBS-domain) that further was combined with F3 peptide (a ligand of nucleolin, a tumour marker) to form the F3-W4 block. The F3-W4 block was taken alone (E2-13W4 protein), as two repeats (E1-W8) and as three repeats (E1-W12). Each protein was supplemented with three copies of the RGD motif (a ligand of integrin αvβ3) and green fluorescent protein (GFP). In contrast to Magnevist (a Gd-containing contrast agent), the proteins exhibited three to four times higher accumulation in U87MG glioma and A375 melanoma cell lines than in normal fibroblasts. The proteins remained for >24 h in tumours induced by Ca755 adenocarcinoma in C57BL/6 mice. They exhibited stability towards blood proteases and only accumulated in the liver and kidney. The technological advantages of using the engineered proteins as a basis for developing efficient and non-toxic agents for early diagnosis of tumours by MRI as well as part of BRT were demonstrated.

Vitamin D3 and its hydroxyderivatives as promising drugs against COVID-19: a computational study

The epidemiologic correlation between the poor prognosis of SARS-CoV-2 infection and vitamin D deficiency has been observed worldwide, however, their molecular mechanisms are not fully understood. In this study, we used combined molecular docking, molecular dynamics simulations and binding free energy analyses to investigate the potentials of vitamin D3 and its hydroxyderivatives as TMPRSS2 inhibitor and to inhibit the SARS-CoV-2 receptor binding domain (RBD) binding to angiotensin-converting enzyme 2 (ACE2), as well as to unveil molecular and structural basis of 1,25(OH)₂D3 capability to inhibit ACE2 and SARS-CoV-2 RBD interactions. The results show that vitamin D3 and its hydroxyderivatives are favorable to bind active site of TMPRSS2 and the binding site(s) between ACE2 and SARS-CoV2-RBD, which indicate that vitamin D3 and its biologically active hydroxyderivatives can serve as TMPRSS2 inhibitor and can inhibit ACE2 binding of SARS-CoV-2 RBD to prevent SARS-CoV-2 entry. Interaction of 1,25(OH)₂D3 with SARS-CoV-2 RBD and ACE2 resulted in the conformation and dynamical motion changes of the binding surfaces between SARS-CoV-2 RBD and ACE2 to interrupt the binding of SARS-CoV-2 RBD with ACE2. The interaction of 1,25(OH)₂D3 with TMPRSS2 also caused the conformational and dynamical motion changes of TMPRSS2, which could affect TMPRSS2 to prime SARS-CoV-2 spike proteins. Our results propose that vitamin D3 and its biologically active hydroxyderivatives are promising drugs or adjuvants in the treatment of COVID-19. Communicated by Ramaswamy H. Sarma.

Comparison of Linear vs. Cyclic RGD Pentapeptide Interactions with Integrin αvβ3 by Molecular Dynamics Simulations

Simple Summary

The integrin α_vβ₃-RGD motif interaction plays a key role in the progression of malignant tumor. Although two typical cyclic and linear RGD short peptides have been widely used in tumor diagnosis and therapy, little is known about the internal dynamic mechanism for different configurations of RGD peptides with different affinities interacting with the integrin α_vβ₃. Our results showed that the cyclic RGD peptide had a more stable configuration in binding to integrins α_vβ₃, which depended on the higher binding energy and higher static electrical energy, especially in the interaction between Asp^RGD-MIDAS. The steered molecular dynamics simulation showed a stronger interaction for the cyclic RGD-integrin α_vβ₃ system than the linear one, with a larger dissociation force (average peak force) and more time to dissociate. Our findings provide insights into the dynamics of integrin α_vβ₃ interactions with linear and cyclic RGD ligands and offer some new therapeutic approaches for the design and development of novel antitumor drugs.

Abstract

Integrin α_vβ₃ interacting with the short Arg-Gly-Asp (RGD) motif plays a critical role in the progression of several types of tumors. However, the effects of the RGD structure (cyclic or linear) with integrin α_vβ₃ at the atomic level remain poorly understood. Here, we performed association and dissociation dynamic simulations for integrin α_vβ₃ in complex with a linear or cyclic pentapeptide by steered molecular dynamics simulations. Compared with cyclic RGD, the linear RGD peptide triggers instability of the configurational changes, mainly resting with the RGD domain due to its flexibility. The main interaction energy between Mg²⁺ and cyclic RGD is much stronger than that of the linear RGD system by the well shield to lessen attacks by free water molecules. The force-dependent dissociation results show that it is easier for linear RGD peptides to leave the active site and much quicker than the cyclic RGD ligand, whereas it is harder to enter the appropriate active binding site in linear RGD. The Ser¹²³-Asp^RGD bond may play a critical role in the allosteric pathway. Our findings provide insights into the dynamics of α_vβ₃ interactions with linear and cyclic RGD ligands and contribute to the application of RGD-based strategies in preclinical therapy.

Engineering of Bio-Adhesive Ligand Containing Recombinant RGD and PHSRN Fibronectin Cell-Binding Domains in Fusion with a Colored Multi Affinity Tag: Simple Approach for Fragment Study from Expression to Adsorption

Engineering of biomimetic motives have emerged as promising approaches to improving cells’ binding properties of biomaterials for tissue engineering and regenerative medicine. In this study, a bio-adhesive ligand including cell-binding domains of human fibronectin (FN) was engineered using recombinant protein technology, a major extracellular matrix (ECM) protein that interacts with a variety of integrins cell-surface’s receptors and other ECM proteins through specific binding domains. 9th and 10th fibronectin type III repeat containing Arginine-Glycine-Aspartic acid (RGD) and Pro-His-Ser-Arg-Asn (PHSRN) synergic site (FNIII9-10) were expressed in fusion with a Colored Multi Affinity Tag (CMAT) to develop a simplified production and characterization process. A recombinant fragment was produced in the bacterial system using E. coli with high yield purified protein by double affinity chromatography. Bio-adhesive surfaces were developed by passive coating of produced fragment onto non adhesive surfaces model. The recombinant fusion protein (CMAT-FNIII9/10) demonstrated an accurate monitoring capability during expression purification and adsorption assay. Finally, biological activity of recombinant FNIII9/10 was validated by cellular adhesion assay. Binding to α5β1 integrins were successfully validated using a produced fragment as a ligand. These results are robust supports to the rational development of bioactivation strategies for biomedical and biotechnological applications.

Fuzzy binding model of molecular interactions between tropoelastin and integrin alphaVbeta3

Tropoelastin is the highly flexible monomer subunit of elastin, required for the resilience of the extracellular matrix in elastic tissues. To elicit biological signaling, multiple sites on tropoelastin bind to cell surface integrins in a poorly understood multifactorial process. We constructed a full atomistic molecular model of the interactions between tropoelastin and integrin α_vβ₃ using ensemble-based computational methodologies. Conformational changes of integrin α_vβ₃ associated with outside-in signaling were more frequently facilitated in an ensemble in which tropoelastin bound the integrin's α1 helix rather than the upstream canonical binding site. Our findings support a model of fuzzy binding, whereby many tropoelastin conformations and defined sites cooperatively interact with multiple α_vβ₃ regions. This model explains prior experimental binding to distinct tropoelastin regions, domains 17 and 36, and points to the cooperative participation of domain 20. Our study highlights the utility of ensemble-based approaches in helping to understand the interactive mechanisms of functionally significant flexible proteins.

Expression and purification of a recombinant ELRL-MAP30 with dual-targeting anti-tumor bioactivity

MAP30 (Momordica antiviral protein 30kD) is a single-chain Ⅰ-type ribosome inactivating protein with a variety of biological activities, including anti-tumor ability. It was reported that MAP30 would serve as a novel and relatively safe agent for prophylaxis and treatment of liver cancer. To determine whether adding two tumor targeting peptides could improve the antitumor activities of MAP30, we genetically modified MAP30 with an RGD motif and a EGFRi motif, which is a ligand with high affinity for αvβ3 integrins and with high affinity for EGFR. The recombinant protein ELRL-MAP30 (rELRL-MAP30) containing a GST-tag was expressed in E. coli. The rELRL-MAP30 was highly expressed in the soluble fraction after induction with 0.15 mM IPTG for 20 h at 16 °C. The purified rELRL-MAP30 appeared as a band on SDS-PAGE. It was identified by western blotting. Cytotoxicity of recombinant protein to HepG2, MDA-MB-231, HUVEC and MCF-7 cells was detected by MTT analysis. Half maximal inhibitory concentration (IC₅₀) values were 54.64 μg/mL, 70.13 μg/mL, 146 μg/mL, 466.4 μg/mL, respectively. Proliferation inhibition assays indicated that rELRL-MAP30 could inhibit the growth of Human liver cancer cell HepG2 effectively. We found that rELRL-MAP30 significantly induced apoptosis in liver cancer cells, as evidenced by nuclear staining of DAPI. In addition, rELRL-MAP30 induced apoptosis in human liver cancer HepG2 cells by up-regulation of Bax as well as down-regulation of Bcl-2. Migration of cell line were markedly inhibited by rELRL-MAP30 in a dose-dependent manner compared to the recombinant MAP30 (rMAP30). In summary, the fusion protein displaying extremely potent cytotoxicity might be highly effective for tumor therapy.

Vitamin D and lumisterol derivatives can act on liver X receptors (LXRs)

The interactions of derivatives of lumisterol (L3) and vitamin D3 (D3) with liver X receptors (LXRs) were investigated. Molecular docking using crystal structures of the ligand binding domains (LBDs) of LXRα and β revealed high docking scores for L3 and D3 hydroxymetabolites, similar to those of the natural ligands, predicting good binding to the receptor. RNA sequencing of murine dermal fibroblasts stimulated with D3-hydroxyderivatives revealed LXR as the second nuclear receptor pathway for several D3-hydroxyderivatives, including 1,25(OH)2D3. This was validated by their induction of genes downstream of LXR. L3 and D3-derivatives activated an LXR-response element (LXRE)-driven reporter in CHO cells and human keratinocytes, and by enhanced expression of LXR target genes. L3 and D3 derivatives showed high affinity binding to the LBD of the LXRα and β in LanthaScreen TR-FRET LXRα and β coactivator assays. The majority of metabolites functioned as LXRα/β agonists; however, 1,20,25(OH)3D3, 1,25(OH)2D3, 1,20(OH)2D3 and 25(OH)D3 acted as inverse agonists of LXRα, but as agonists of LXRβ. Molecular dynamics simulations for the selected compounds, including 1,25(OH)2D3, 1,20(OH)2D3, 25(OH)D3, 20(OH)D3, 20(OH)L3 and 20,22(OH)2L3, showed different but overlapping interactions with LXRs. Identification of D3 and L3 derivatives as ligands for LXRs suggests a new mechanism of action for these compounds.

Heterogeneous role of integrins in fibroblast response to small cyclic mechanical stimulus generated by a nanoporous gold actuator

It is important to understand the effects of mechanical stimulation on cell behaviors for homeostasis. Many studies have been performed on cell responses to mechanical stimuli, but the mechanosensing mechanism is still under debate. In the present study, experiments employing molecular dynamics (MD) simulations concerning the effects of cyclic mechanical stimulus on cell proliferation were performed based on the hypothesis that mechanosensing depends on integrin types. We used a nanoporous gold (NPG) actuator to prevent transfer of a mechanical stimulus via molecules other than integrins. Surprisingly, a small cyclic strain of only 0.5% enhanced the proliferation of fibroblasts. α₅β₁ and α_vβ₃ integrins showed high sensitivity to the mechanical stimulus, whereas α₁β₁ and α₂β₁ integrins exhibited low mechanosensitivity. The MD simulations showed that different conformational changes of the integrin headpiece induced by binding to the ECM led to a difference in mechanosensitivity between αI and αI-less integrin types. Thus, the present study provides evidence to support the hypothesis and suggests the mechanism for the heterogeneous roles of integrins in mechanosensing.

Molecular dynamics simulations to the bidirectional adhesion signaling pathway of integrin αV β3

The bidirectional force transmission process of integrin through the cell membrane is still not well understood. Several possible mechanisms have been discussed in literature on the basis of experimental data, and in this study, we investigate these mechanisms by free and steered molecular dynamics simulations. For the first time, constant velocity pulling on the complete integrin molecule inside a dipalmitoyl-phosphatidylcholine membrane is conducted. From the results, the most likely mechanism for inside-out and outside-in signaling is the switchblade model with further separation of the transmembrane helices.

Fibronectin Regulation of Integrin B1 and SLUG in Circulating Tumor Cells

Metastasis is the leading cause of cancer death worldwide. Circulating tumor cells (CTCs) are a critical step in the metastatic cascade and a good tool to study this process. We isolated CTCs from a syngeneic mouse model of hepatocellular carcinoma (HCC) and a human xenograft mouse model of castration-resistant prostate cancer (CRPC). From these models, novel primary tumor and CTC cell lines were established. CTCs exhibited greater migration than primary tumor-derived cells, as well as epithelial-to-mesenchymal transition (EMT), as observed from decreased E-cadherin and increased SLUG and fibronectin expression. Additionally, when fibronectin was knocked down in CTCs, integrin B1 and SLUG were decreased, indicating regulation of these molecules by fibronectin. Investigation of cell surface molecules and secreted cytokines conferring immunomodulatory advantage to CTCs revealed decreased major histocompatibility complex class I (MHCI) expression and decreased endostatin, C-X-C motif chemokine 5 (CXCL5), and proliferin secretion by CTCs. Taken together, these findings indicate that CTCs exhibit distinct characteristics from primary tumor-derived cells. Furthermore, CTCs demonstrate enhanced migration in part through fibronectin regulation of integrin B1 and SLUG. Further study of CTC biology will likely uncover additional important mechanisms of cancer metastasis.

Design Optimization of Tumor Vasculature-Bound Nanoparticles

Nanotherapy may constitute a promising approach to target tumors with anticancer drugs while minimizing systemic toxicity. Computational modeling can enable rapid evaluation of nanoparticle (NP) designs and numerical optimization. Here, an optimization study was performed using an existing tumor model to find NP size and ligand density that maximize tumoral NP accumulation while minimizing tumor size. Optimal NP avidity lies at lower bound of feasible values, suggesting reduced ligand density to prolong NP circulation. For the given set of tumor parameters, optimal NP diameters were 288 nm to maximize NP accumulation and 334 nm to minimize tumor diameter, leading to uniform NP distribution and adequate drug load. Results further show higher dependence of NP biodistribution on the NP design than on tumor morphological parameters. A parametric study with respect to drug potency was performed. The lower the potency of the drug, the bigger the difference is between the maximizer of NP accumulation and the minimizer of tumor size, indicating the existence of a specific drug potency that minimizes the differential between the two optimal solutions. This study shows the feasibility of applying optimization to NP designs to achieve efficacious cancer nanotherapy, and offers a first step towards a quantitative tool to support clinical decision making.

Modeling the structural and dynamical changes of the epithelial calcium channel TRPV5 caused by the A563T variation based on the structure of TRPV6

TRPV5, transient receptor potential cation channel vanilloid subfamily member 5, is an epithelial Ca²⁺ channel that plays a key role in the active Ca²⁺ reabsorption process in the kidney. A single nucleotide polymorphism (SNP) rs4252499 in the TRPV5 gene results in an A563T variation in the sixth transmembrane (TM) domain of TRPV5. Our previous study indicated that this variation increases the Ca²⁺ transport function of TRPV5. To understand the molecular mechanism, a model of TRPV5 was established based on the newly deposited structure of TRPV6 that has 83.1% amino acid identity with TRPV5 in the modeled region. Computational simulations were performed to study the structural and dynamical differences between the TRPV5 variants with A563 and T563. Consistent with the TRPV1-based simulation, the results indicate that the A563T variation increases the contacts between residues 563 and V540, which is one residue away from the key residue D542 in the Ca²⁺-selective filter. The variation enhanced the stability of the secondary structure of the pore region, decreased the fluctuation of residues around residue 563, and reduced correlated and anti-correlated motion between monomers. Furthermore, the variation increases the pore radius at the selective filter. These findings were confirmed using simulations based on the recently determined structure of rabbit TRPV5. The simulation results provide an explanation for the observation of enhanced Ca²⁺ influx in TRPV5 caused by the A563T variation. The A563T variation is an interesting example of how a residue distant from the Ca²⁺-selective filter influences the Ca²⁺ transport function of the TRPV5 channel. Communicated by Ramaswamy H. Sarma.

Multiscale simulation of the interaction of calreticulin-thrombospondin-1 complex with a model membrane microdomain

Cell surface calreticulin (CRT) binding to thrombospondin-1 (TSP1), regulates cell adhesion, migration, anoikis resistance, and collagen production. Due to the essential role of membrane microdomains in CRT-mediated focal adhesion disassembly, we previously studied the effect of raft-like bilayers on TSP1-CRT interactions with all-atom molecular dynamics (AAMD) simulations. However, the simulated systems of protein on the surface of the bilayer(s) in the explicit solvent are too large for long timescale AAMD simulations due to computational expense. In this study, we adopted a multiscale modeling approach of combining AAMD, coarse-grained molecule dynamics (CGMD), and reversed AAMD (REV AAMD) simulations to investigate the interactions of single CRT or of the TSP1-CRT complex with a membrane microdomain at microsecond timescale. Results showed that CRT conformational stabilization by binding of TSP1 in AAMD simulation was undetectable in CGMD simulation, but it was recovered in REV AAMD simulation. Similarly, interactions of the CRT N-domain and TSP1 with the membrane microdomain were lost in CGMD simulations but they were re-gained in the REV AAMD simulations. There was the higher coordination of the CRT P-domain in the TSP1-CRT complex with the lipid components of membrane microdomain compared to that of single CRT, which could directly affect the conformation of CRT and further mediate CRT recruitment of LDL receptor-related protein for signaling events. This study provides structural and molecular insights into TSP1-CRT interactions in a membrane microdomain environment and demonstrates the feasibility of using multiscale simulations to investigate the interactions between protein and membrane microdomains at a long timescale.

Molecular insights into the specific recognition between the RNA binding domain qRRM2 of hnRNP F and G-tract RNA: A molecular dynamics study

Heterogeneous nuclear ribonucleoprotein F (hnRNP F) controls the expression of various genes through regulating the alternative splicing of pre-mRNAs in the nucleus. It uses three quasi-RNA recognition motifs (qRRMs) to recognize G-tract RNA which contains at least three consecutive guanines. The structures containing qRRMs of hnRNP F in complex with G-tract RNA have been determined by nuclear magnetic resonance (NMR) spectroscopy, shedding light on the recognition mechanism of qRRMs with G-tract RNA. However, knowledge of the recognition details is still lacking. To investigate how qRRMs specifically bind with G-tract RNA and how the mutations of any guanine to an adenine in the G-tract affect the binding, molecular dynamics simulations with binding free energy analysis were performed based on the NMR structure of qRRM2 in complex with G-tract RNA. Simulation results demonstrate that qRRM2 binds strongly with G-tract RNA, but any mutation of the G-tract leads to a drastic reduction of the binding free energy. Further comparisons of the energetic components reveal that van der Waals and non-polar interactions play essential roles in the binding between qRRM2 and G-tract RNA, but the interactions are weakened by the effect of RNA mutations. Structural and dynamical analyses indicate that when qRRM2 binds with G-tract RNA, both qRRM2 and G-tract maintain stabilized structures and dynamics; however, the stability is disrupted by the mutations of the G-tract. These results provide novel insights into the recognition mechanism of qRRM2 with G-tract RNA that are not elucidated by the NMR technique.

The L530R variation associated with recurrent kidney stones impairs the structure and function of TRPV5

TRPV5 is a Ca²⁺-selective channel that plays a key role in the reabsorption of Ca²⁺ ions in the kidney. Recently, a rare L530R variation (rs757494578) of TRPV5 was found to be associated with recurrent kidney stones in a founder population. However, it was unclear to what extent this variation alters the structure and function of TRPV5. To evaluate the function and expression of the TRPV5 variant, Ca²⁺ uptake in Xenopus oocytes and western blot analysis were performed. The L530R variation abolished the Ca²⁺ uptake activity of TRPV5 in Xenopus oocytes. The variant protein was expressed with drastic reduction in complex glycosylation. To assess the structural effects of this L530R variation, TRPV5 was modeled based on the crystal structure of TRPV6 and molecular dynamics simulations were carried out. Simulation results showed that the L530R variation disrupts the hydrophobic interaction between L530 and L502, damaging the secondary structure of transmembrane domain 5. The variation also alters its interaction with membrane lipid molecules. Compared to the electroneutral L530, the positively charged R530 residue shifts the surface electrostatic potential towards positive. R530 is attracted to the negatively charged phosphate group rather than the hydrophobic carbon atoms of membrane lipids. This shifts the pore helix where R530 is located and the D542 residue in the Ca²⁺-selective filter towards the surface of the membrane. These alterations may lead to misfolding of TRPV5, reduction in translocation of the channel to the plasma membrane and/or impaired Ca²⁺ transport function of the channel, and ultimately disrupt TRPV5-mediated Ca²⁺ reabsorption.

Deprotonation states of the two active site water molecules regulate the binding of protein phosphatase 5 with its substrate: A molecular dynamics study

Protein phosphatase 5 (PP5), mainly localized in human brain, can dephosphorylate tau protein whose high level of phosphorylation is related to Alzheimer's disease. Similar to other protein phosphatases, PP5 has a conserved motif in the catalytic domain that contains two binding sites for manganese (Mn²⁺ ) ions. Structural data indicate that two active site water molecules, one bridging the two Mn²⁺ ions and the other terminally coordinated with one of the Mn²⁺ ions (Mn1), are involved in catalysis. Recently, a density functional theory study revealed that the two water molecules can be both deprotonated to keep a neutral active site for catalysis. The theoretical study gives us an insight into the catalytic mechanism of PP5, but the knowledge of how the deprotonation states of the two water molecules affect the binding of PP5 with its substrate is still lacking. To approach this problem, molecular dynamics simulations were performed to model the four possible deprotonation states. Through structural, dynamical and energetic analyses, the results demonstrate that the deprotonation states of the two water molecules affect the structure of the active site including the distance between the two Mn²⁺ ions and their coordination, impact the interaction energy of residues R275, R400 and H304 which directly interact with the substrate phosphoserine, and mediate the dynamics of helix αJ which is involved in regulation of the enzyme's activity. Furthermore, the deprotonation state that is preferable for PP5 binding of its substrate has been identified. These findings could provide new design strategy for PP5 inhibitor.