Molecular Dynamics Free Energy Calculations Research Articles

In silico analysis of aptamer-RNA conjugate interactions with human transferrin receptor

The human transmembrane protein Transferrin Receptor-1 is regarded as a promising target for the systemic delivery of therapeutic agents, particularly of nucleic acid therapeutics, such as short double stranded RNAs. This ubiquitous receptor is involved in cellular iron uptake, keeping intracellular homeostasis. It is overexpressed in multiple cancer cell types and is internalized via clathrin-mediated endocytosis. In previous studies, a human transferrin receptor-1 RNA aptamer, identified as TR14 ST1–3, was shown to be capable of effectively internalizing into cells in culture and to deliver small, double stranded RNAs in vitro and in vivo, via systemic administration.To understand, at the molecular level, the aptamer binding to the receptor and the impact of conjugation with the therapeutic RNA, a multi-level in silico protocol was employed, including protein-aptamer docking, molecular dynamics simulations and free energy calculations. The competition for the binding pocket, between the aptamer and the natural ligand human Transferrin, was also evaluated.The results show that the aptamer binds to the same region as Transferrin, with residues from the helical domain showing a critical role. Moreover, the conjugation to the therapeutic RNA, was shown not to affect aptamer binding. Overall, this study provides an atomic-level understanding of aptamer association to human Transferrin Receptor-1 and of its conjugation with a short model-therapeutic RNA, providing also important clues for futures studies aiming to deliver other oligonucleotide-based therapeutics via Transferrin Receptor.

H and CO Co-Induced Roughening of Cu Surface in CO2 Electroreduction Conditions.

The dynamic restructuring of Cu has been observed under electrochemical conditions, and it has been hypothesized to underlie the unique reactivity of Cu toward CO2 electroreduction. Roughening is one of the key surface phenomena for Cu activation, whereby numerous atomic vacancies and adatoms form. However, the atomic structure of such surface motifs in the presence of relevant adsorbates has remained elusive. Here, we explore the chemical space of Cu surface restructuring under coverage of CO and H in realistic electroreduction conditions, by combining grand canonical DFT and global optimization techniques, from which we construct a potential-dependent grand canonical ensemble representation. The regime of intermediate and mixed CO and H coverage─where structures exhibit some elevated surface Cu─is thermodynamically unfavorable yet kinetically inevitable. Therefore, we develop a quasi-kinetic Monte Carlo simulation to track the system's evolution during a simulated cathodic scan. We reveal the evolution path of the system across coverage space and identify the accessible metastable structures formed along the way. Chemical bonding analysis is performed on the metastable structures with elevated Cu*CO species to understand their formation mechanism. By molecular dynamics simulations and free energy calculations, the surface chemistry of the Cu*CO species is explored, and we identify plausible mechanisms via which the Cu*CO species may diffuse or dimerize. This work provides rich atomistic insights into the phenomenon of surface roughening and the structure of involved species. It also features generalizable methods to explore the chemical space of restructuring surfaces with mixed adsorbates and their nonequilibrium evolution.

Direct Interaction of Minocycline to p47phox Contributes to its Attenuation of TNF-α-Mediated Neuronal PC12 Cell Death: Experimental and Simulation Validation.

Minocycline, a repurposed approved medication, shows promise in treating neurodegeneration. However, the specific pathways targeted by minocycline remain unclear despite the identification of molecular targets. This study explores minocycline's potential protective effects against TNF-α-mediated neuronal death in PC12 cells, with a focus on unraveling its interactions with key molecular targets. The study begins by exploring minocycline's protective role against TNF-α-mediated neuronal death in PC12 cells, showcasing a substantial reduction in cleaved caspase-3 expression, DNA fragmentation, and intracellular ROS levels following minocycline pretreatment. Subsequently, a comprehensive analysis utilizing pull-down assays, computational docking, mutation analysis, molecular dynamics simulations, and free energy calculations is conducted to elucidate the direct interaction between minocycline and p47phox-the organizer subunit of NADPH oxidase-2 (NOX2) complex. Computational insights, including a literature survey and analysis of key amino acid residues, reveal a potential binding site for minocycline around Trp193 and Cys196. In silico substitutions of Trp193 and Cys196 further confirm their importance in binding with minocycline. These integrated findings underscore minocycline's protective mechanisms, linking its direct interaction with p47phox to the modulation of NOX2 activity and attenuation of NOX-derived ROS generation. Minocycline demonstrates protective effects against TNF-α-induced PC12 cell death, potentially linked to its direct interaction with p47phox. This interaction leads to a reduction in NOX2 complex assembly, ultimately attenuating NOX-derived ROS generation. These findings hold significance for researchers exploring neuroprotection and the development of p47phox inhibitors.

Exploring the Super-Relaxed State of Human Cardiac β-Myosin by Molecular Dynamics Simulations.

Human β-cardiac myosin plays a critical role in generating the mechanical forces necessary for cardiac muscle contraction. This process relies on a delicate dynamic equilibrium between the disordered relaxed state (DRX) and the super-relaxed state (SRX) of myosin. Disruptions in this equilibrium due to mutations can lead to heart diseases. However, the structural characteristics of SRX and the molecular mechanisms underlying pathogenic mutations have remained elusive. To bridge this gap, we conducted molecular dynamics simulations and free energy calculations to explore the conformational changes in myosin. Our findings indicate that the size of the phosphate-binding pocket can serve as a valuable metric for characterizing the transition from the DRX to SRX state. Importantly, we established a global dynamic coupling network within the myosin motor head at the residue level, elucidating how the pathogenic mutation E483K impacts the equilibrium between SRX and DRX through allosteric effects. Our work illuminates molecular details of SRX and offers valuable insights into disease treatment through the regulation of SRX.

Lomitapide as a Potential Estrogen Receptor Inhibitor: A Computational Drug Repurposing Study

Objective: Estrogen receptor (ER) inhibitors have significant therapeutic potential for hormone-dependent cancers and related disorders. Tamoxifen, a well-known selective estrogen receptor modulator, has been widely used as adjuvant therapy for estrogen receptor-positive breast cancer. However, tamoxifen may exhibit a tendency to develop resistance with prolonged usage and particularly elevate the risk of uterine cancer. Therefore, there is a need for the discovery and development of new ER modulators or inhibitors. In this study, we identified potential estrogen receptor inhibitors through computational drug repositioning.
 
 Methods: A set of 2048 compounds, encompassing FDA-approved drugs and active metabolites, were subjected to molecular docking, molecular dynamics simulations, and free energy calculations to evaluate their interaction with estrogen receptor α (ERα).
 
 Results: Among the compounds evaluated, conivaptan, atogepant, and lomitapide exhibited the highest affinities for ERα. Lomitapide displayed a superior docking score (-12 kcal/mol) compared to the established ER inhibitor, tamoxifen (-10 kcal/mol). Further investigation using molecular dynamics simulations and free energy calculations disclosed lomitapide's heightened binding affinity of -380.727 kJ/mol, surpassing tamoxifen's binding affinity of -352.029 kJ/mol.
 
 Conclusion: This comprehensive computational exploration underscores lomitapide's potential as a compelling candidate with an envisaged stronger estrogen receptor affinity than the acknowledged standard, tamoxifen. To validate lomitapide's promise as a novel ER inhibitor, essential in vitro and in vivo studies are suggested. These investigations will provide essential insights into lomitapide's reposition in addressing the challenges tied to hormone-dependent cancers and associated maladies.

Cyclization Scaffolding for Improved Vaccine Immunogen Stability: Application to Tau Protein in Alzheimer's Disease.

Effective scaffolding of immunogens is crucial for generating conformationally selective antibodies through active immunization, particularly in the treatment of protein misfolding diseases such as Alzheimer's and Parkinson's disease. Previous computational work has revealed that a disorder-prone region of the tau protein, when in a stacked form, is predicted to structurally resemble a small, soluble protofibril, having conformational properties similar to those of experimental in vitro tau oligomers. Such an oligomeric structural mimic has the potential to serve as a vaccine immunogen design for Alzheimer's disease. In this study, we developed a cyclization scaffolding method in Rosetta, in which multiple cyclic peptides are stacked into a protofibril. Cyclization results in significant stabilization of protofibril-like structures by constraining the conformational space. Applying this method to the disorder-prone region of the tau fibril, we evaluated the metastability of the cyclized tau immunogen using molecular dynamics simulations, and we identified sequences of two cyclic constructs having high metastability in the protofibril. We then assessed their thermodynamic stability by computing the free energy required to separate a distal chain from the rest of the stacked structure. Our computational results, based on molecular dynamics simulations and free energy calculations, demonstrate that two cyclized constructs, cyclo-(VKSEKLDFKDRVQSKIFyN) and cyclo-(VKSEKLDFKDRVQSKIYvG) (lowercase letters indicate d-form amino acids), possess significantly increased thermodynamic stability in the protofibril over an uncyclized linear construct VKSEKLDFKDRVQSKI. The cyclization scaffolding approach proposed here holds promise as a means to effectively design immunogens for protein misfolding diseases, particularly those involving liposome-conjugated peptide constructs.

Harnessing Brazilian biodiversity database: identification of flavonoids as potential inhibitors of SARS-CoV-2 main protease using computational approaches and all-atom molecular dynamics simulation.

SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) is the etiological agent responsible for the global outbreak of COVID-19 (Coronavirus Disease 2019). The main protease of SARS-CoV-2, Mpro, is a key enzyme that plays a vital role in mediating viral replication and transcription. In this study, a comprehensive computational approach was employed to investigate the binding affinity, selectivity, and stability of natural product candidates as potential new antivirals acting on the viral polyprotein processing mediated by SARS-CoV-2 Mpro. A library of 288 flavonoids extracted from Brazilian biodiversity was screened to select potential Mpro inhibitors. An initial filter based on Lipinski's rule of five was applied, and 204 compounds that did not violate any of the Lipinski rules were selected. The compounds were then docked into the active site of Mpro using the GOLD program, and the poses were subsequently re-scored using MM-GBSA (Molecular Mechanics Generalized Born Surface Area) binding free energy calculations performed by AmberTools23. The top five flavonoids with the best MM-GBSA binding free energy values were selected for analysis of their interactions with the active site residues of the protein. Next, we conducted a toxicity and drug-likeness analysis, and non-toxic compounds were subjected to molecular dynamics simulation and free energy calculation using the MM-PBSA (Molecular Mechanics Poisson-Boltzmann Surface Area) method. It was observed that the five selected flavonoids had lower MM-GBSA binding free energy with Mpro than the co-crystal ligand. Furthermore, these compounds also formed hydrogen bonds with two important residues, Cys145 and Glu166, in the active site of Mpro. Two compounds that passed the drug-likeness filter showed stable conformations during the molecular dynamics simulations. Among these, NuBBE_867 exhibited the best MM-PBSA binding free energy value compared to the crystallographic inhibitor. Therefore, this study suggests that NuBBE_867 could be a potential inhibitor against the main protease of SARS-CoV-2 and may be further examined to confirm our results.

An Efficient Approach to the Accurate Prediction of Mutational Effects in Antigen Binding to the MHC1.

The major histocompatibility complex (MHC) can recognize and bind to external peptides to generate effective immune responses by presenting the peptides to T cells. Therefore, understanding the binding modes of peptide-MHC complexes (pMHC) and predicting the binding affinity of pMHCs play a crucial role in the rational design of peptide vaccines. In this study, we employed molecular dynamics (MD) simulations and free energy calculations with an Alanine Scanning with Generalized Born and Interaction Entropy (ASGBIE) method to investigate the protein-peptide interaction between HLA-A*02:01 and the G9209 peptide derived from the melanoma antigen gp100. The energy contribution of individual residue was calculated using alanine scanning, and hotspots on both the MHC and the peptides were identified. Our study shows that the pMHC binding is dominated by the van der Waals interactions. Furthermore, we optimized the ASGBIE method, achieving a Pearson correlation coefficient of 0.91 between predicted and experimental binding affinity for mutated antigens. This represents a significant improvement over the conventional MM/GBSA method, which yields a Pearson correlation coefficient of 0.22. The computational protocol developed in this study can be applied to the computational screening of antigens for the MHC1 as well as other protein-peptide binding systems.

PHARMACOINFORMATICS ANALYSIS OF MORUS MACROURA FOR DRUG DISCOVERY AND DEVELOPMENT

Objective: Pharmacoinformatics is an innovative approach rapidly evolving in pharmaceutical research and drug development. This study focuses on analysing Morus macroura, a plant species with untapped pharmacological potential. This investigation aims to leverage pharmacoinformatics techniques to unveil the hidden potential of Morus macroura in drug discovery and development.
 Methods: The study includes analyses of protein-protein interactions, deep learning docking, adsorption tests, distribution, metabolism, excretion, molecular dynamics simulations and free energy calculation using Molecular Mechanics Generalized Born Surface Area (MMGBSA).
 Results: Nine active compounds were identified in Morus macroura, namely Andalasin A, Guangsangon K, Guangsangon L, Guangsangon M, Guangsangon N, Macrourone C, Mulberrofuran G, Mulberrofuran K, and Mulberroside C. These compounds exhibit protein-protein interaction activities against a cytochrome P450 monooxygenase that catalyses the conversion of C19 androgens. These plant compounds influence aromatase excess syndrome, deficiency, and ovarian dysgenesis. Regarding drug-likeness, Mulberroside C and Macrourone C demonstrated good absorption potential by adhering to Lipinski's rule of five. Deep learning docking simulations yielded affinity results of-9.62 kcal/mol for Guangsangon M,-10.44 kcal/mol for Macrourone C, and-10.99 kcal/mol for Guangsangon L. Subsequent molecular dynamics simulations indicated that Guangsangon L and Macrourone C remained stable during a 100 ns simulation.
 Conclusion: Morus macroura interacts with important proteins, particularly CYP19A1, which might influence health conditions like aromatase excess syndrome and ovarian dysgenesis. These findings provide potential paths for addressing specific health issues and advancing drug development. Molecular dynamics simulations indicated that Guangsangon L and Macrourone C remained stable during simulation.

Interaction of Asn297-Linked Glycan Ligands with the Fc Fragment of the Immunoglobulin Class G1: A Molecular Dynamics Simulation Study.

The allosteric modulation of the homodimeric H10-03-6 protein to glycan ligands L1 and L2, and the STAB19 protein to glycan ligands L3 and L4, respectively, has been studied by molecular dynamics simulations and free energy calculations. The results revealed that the STAB19 protein has a significantly higher affinity for L3 (-11.38 ± 2.32 kcal/mol) than that for L4 (-5.51 ± 1.92 kcal/mol). However, the combination of the H10-03-6 protein with glycan L2 (1.23 ± 6.19 kcal/mol) is energetically unfavorable compared with that of L1 (-13.96 ± 0.35 kcal/mol). Further, the binding of glycan ligands L3 and L4 to STAB19 would result in the significant closure of the two CH2 domains of the STAB19 conformation with the decrease of the centroid distances between the two CH2 domains compared with the H10-03-6/L1/L2 complex. The CH2 domain closure of STAB19 relates directly to the formation of new hydrogen bonds and hydrophobic interactions between the residues Ser239, Val240, Asp265, Glu293, Asn297, Thr299, Ser337, Asp376, Thr393, Pro395, and Pro396 in STAB19 and glycan ligands L3 and L4, which suggests that these key residues would contribute to the specific regulation of STAB19 to L3 and L4. In addition, the distance analysis revealed that the EF loop in the H10-03-6/L1/L2 model presents a high flexibility and partial disorder compared with the stabilized STAB19/L3/L4 complex. These results will be helpful in understanding the specific regulation through the asymmetric structural characteristics in the CH2 and CH3 domains of the H10-03-6 and STAB19 proteins.

Virtual screening of flavonoids as potential RIPK1 inhibitors for neurodegeneration therapy.

Global prevalence of neurodegenerative diseases such as Alzheimer's disease and Parkinson's disease is increasing gradually, whereas approvals of successful therapeutics for central nervous system disorders are inadequate. Accumulating evidence suggests pivotal roles of the receptor-interacting serine/threonine-protein kinase 1 (RIPK1) in modulating neuroinflammation and necroptosis. Discoveries of potent small molecule inhibitors for RIPK1 with favorable pharmacokinetic properties could thus address the unmet medical needs in treating neurodegeneration. In a structure-based virtual screening, we performed site-specific molecular docking of 4,858 flavonoids against the kinase domain of RIPK1 using AutoDock Vina. We predicted physicochemical descriptors of the top ligands using the SwissADME webserver. Binding interactions of the best ligands and the reference ligand L8D were validated using replicated 500-ns Gromacs molecular dynamics simulations and free energy calculations. From Vina docking, we shortlisted the top 20 flavonoids with the highest binding affinities, ranging from -11.7 to -10.6 kcal/mol. Pharmacokinetic profiling narrowed down the list to three orally bioavailable and blood-brain-barrier penetrant flavonoids: Nitiducarpin, Pinocembrin 7-O-benzoate, and Paratocarpin J. Next, trajectories of molecular dynamics simulations of the top protein-ligand complexes were analyzed for binding interactions. The root-mean-square deviation (RMSD) was 1.191 Å (±0.498 Å), 1.725 Å (±0.828 Å), 1.923 Å (±0.942 Å), 0.972 Å (±0.155 Å) for Nitiducarpin, Pinocembrin 7-O-benzoate, Paratocarpin J, and L8D, respectively. The radius of gyration (Rg) was 2.034 nm (±0.015 nm), 2.0.39 nm (± 0.025 nm), 2.053 nm (±0.021 nm), 2.037 nm (±0.016 nm) for Nitiducarpin, Pinocembrin 7-O-benzoate, Paratocarpin J, and L8D, respectively. The solvent accessible surface area (SASA) was 159.477 nm2 (±3.021 nm2), 159.661 nm2 (± 3.707 nm2), 160.755 nm2 (±4.252 nm2), 156.630 nm2 (±3.521 nm2), for Nitiducarpin, Pinocembrin 7-O-benzoate, Paratocarpin J, and L8D complexes, respectively. Therefore, lower RMSD, Rg, and SASA values demonstrated that Nitiducarpin formed the most stable complex with the target protein among the best three ligands. Finally, 2D protein-ligand interaction analysis revealed persistent hydrophobic interactions of Nitiducarpin with the critical residues of RIPK1, including the catalytic triads and the activation loop residues, implicated in the kinase activity and ligand binding. Our target-based virtual screening identified three flavonoids as strong RIPK1 inhibitors, with Nitiducarpin exhibiting the most potent inhibitory potential. Future in vitro and in vivo studies with these ligands could offer new hope for developing effective therapeutics and improving the quality of life for individuals affected by neurodegeneration.

(De)carboxylation mechanisms of heteroaromatic substrates catalyzed by prenylated FMN-dependent UbiD decarboxylases: An in-silico study

The UbiD enzymes are proposed to catalyze reversible (de)carboxylation reaction of unsaturated carboxylic acids using prenylated flavin mononucleotide (prFMN) as a cofactor. This positions UbiD enzymes as promising candidates for converting CO2 into valuable chemicals. However, their industrial-scale biotransformation is currently constrained by low conversion rates attributed to thermodynamic limitations. To enhance the carboxylation activity of UbiD enzymes, a molecular-level understanding of the (de)carboxylation mechanisms is necessary. In this study, we investigated the reaction mechanisms of heteroaromatic substrates catalyzed by PtHmfF, PaHudA, and AnlnD enzymes using molecular dynamics (MD) simulations and free energy calculations. Our extensive mechanistic study elucidates the mechanisms involved in the formation of the initial prFMN-substrate intermediate. Specifically, we observed nucleophilic attack during decarboxylation, while carboxylation reactions involving furoic acid, pyrrole, and indole tend to favor a 1,3-dipolar cycloaddition mechanism. Furthermore, we identified proton transfer as the rate-limiting step in the carboxylation reaction. In addition, we considered the perspectives of reaction energies and electron transfer to understand the distinct mechanisms underlying decarboxylation and carboxylation. Our calculated free energies are consistent with available experimental kinetics data. Finally, we explored how different rotamers of catalytic residues influence the efficiency of the initial intermediate formation.

Ion Transport and Inhibitor Binding by Human NHE1: Insights from Molecular Dynamics Simulations and Free Energy Calculations.

The human Na+/H+ exchanger (NHE1) plays a crucial role in maintaining intracellular pH by regulating the electroneutral exchange of a single intracellular H+ for one extracellular Na+ across the plasma membrane. Understanding the molecular mechanisms governing ion transport and the binding of inhibitors is of importance in the development of anticancer therapeutics targeting NHE1. In this context, we performed molecular dynamics (MD) simulations based on the recent cryo-electron microscopy (cryo-EM) structures of outward- and inward-facing conformations of NHE1. These simulations allowed us to explore the dynamics of the protein, examine the ion-translocation pore, and confirm that Asp267 is the ion-binding residue. Our free energy calculations did not show a significant difference between Na+ and K+ binding at the ion-binding site. Consequently, Na+ over K+ selectivity cannot be solely explained by differences in ion binding. Our MD simulations involving NHE1 inhibitors (cariporide and amiloride analogues) maintained stable interactions with Asp267 and Glu346. Our study highlights the importance of the salt bridge between the positively charged acylguanidine moiety and Asp267, which appears to play a role in the competitive inhibitory mechanism for this class of inhibitors. Our computational study provides a detailed mechanistic interpretation of experimental data and serves the basis of future structure-based inhibitor design.

Tuning the stability of a model quasicrystal and its approximants with a periodic substrate.

Quasicrystals and their periodic approximants are complex crystalline phases. They have now been observed in many metallic alloys, soft matter systems, and particle simulations. In recent experiments of thin-film perovskites on solid substrates, the type of complex phase was found to change depending on thermodynamic conditions and the type of substrate used. Here, we investigate the effect of a substrate on the relative thermodynamic stability of a two-dimensional model quasicrystal and its approximants. Our simulation model is particles interacting via the Lennard-Jones-Gauss potential. Our numerical methods are molecular dynamics simulations and free energy calculations that take into account phason flips explicitly. For substrates interacting weakly with the particles, we observe an incommensurate-commensurate transition, in which a continuous series of quasicrystal approximants locks into a small number of approximants. Interestingly, we observe that the 3/2 approximant exhibits phason mode fluctuations in thermodynamic equilibrium. Such fluctuations are reminiscent of random tiling and a phenomenon usually associated only with quasiperiodic order. For stronger substrates, we find an enhancement of the stability of the dodecagonal quasicrystal and variants of square lattices. We explain all observed phenomena by the interplay of the model system with the substrate. Our results demonstrate that designing novel complex periodic and quasiperiodic structures by choice of suitable substrates is a promising strategy.

The effect of polymer coating on nanoparticles' interaction with lipid membranes studied by coarse-grained molecular dynamics simulations.

Nanoparticles' (NPs) permeation through cell membranes, whether it happens via passive or active transport, is an essential initial step for their cellular internalization. The NPs' surface coating impacts the way they translocate through the lipid bilayer and the spontaneity of the process. Understanding the molecular details of NPs' interaction with cell membranes allows the design of nanosystems with optimal characteristics for crossing the lipid bilayer: computer simulations are a powerful tool for this purpose. In this work, we have performed coarse-grained molecular dynamics simulations and free energy calculations on spherical titanium dioxide NPs conjugated with polymer chains of different chemical compositions. We have demonstrated that the hydrophobic/hydrophilic character of the chains, more than the nature of their terminal group, plays a crucial role in determining the NPs' interaction with the lipid bilayer and the thermodynamic spontaneity of NPs' translocation from water to the membrane. We envision that this computational work will be helpful to the experimental community in terms of the rational design of NPs for efficient cell membrane permeation.

Relationships between Water's Structure and Solute Affinity at Polypeptoid Brush Surfaces.

Excellent antifouling surfaces are generally thought to create a tightly bound layer of water that resists solute adsorption, and highly hydrophilic surfaces such as those with zwitterionic functionalities are of significant current interest as antifoulant strategies. However, despite significant proofs-of-concept, we still lack a fundamental understanding of how the nanoscopic structure of this hydration layer translates to reduced fouling, how surface chemistry can be tuned to achieve antifouling through hydration water, and why, in particular, zwitterionic surfaces seem so promising. Here, we use molecular dynamics simulations and free energy calculations to investigate the molecular relationships among surface chemistry, hydration water structure, and surface-solute affinity across a variety of surface-decorated chemistries. Specifically, we consider polypeptoid-decorated surfaces that display well-known experimental antifouling capabilities and that can be synthesized sequence specifically, with precise backbone positioning of, e.g., charged groups. Through simulations, we calculate the affinities of a range of small solutes to polypeptoid brush surfaces of varied side-chain chemistries. We then demonstrate that measures of the structure of surface hydration water in response to a particular surface chemistry signal solute-surface affinity; specifically, we find that zwitterionic chemistries produce solute-surface repulsion through highly coordinated hydration water while suppressing tetrahedral structuring around the solute, in contrast to uncharged surfaces that show solute-surface affinity. Based on the relationship of this structural perturbation to the affinity of small-molecule solutes, we propose a molecular mechanism by which zwitterionic surface chemistries enhance solute repulsion, with broader implications for the design of antifouling surfaces.

Mechanistic Insight into how β-Ketoacyl ACP Synthase I (KasA) Recognizes the Fatty Acid Chain Length of its Substrate.

β-ketoacyl ACP synthase I (KasA) has been considered as a promising drug target against Tuberculosis because it is known to play a pivotal role in the survival of Mycobacterium Tuberculosis, a causative agent of Tuberculosis. KasA catalyzes the reaction elongating only the acyl chain that is 16 carbon atoms in length or longer, but the molecular details of how KasA selectively recognizes only the substrate longer than a certain length still remain unknown. In the present study, this challenging subject is addressed, and to this end, molecular dynamics (MD) simulations and free energy calculations for actual substrate binding process are carried out. The results illustrate that the substrate specificity of KasA is highly linked to its cooperativity and this cooperativity is realized through the activation of catalytic residues. Through these results, the mechanistic details of how KasA can be selectively activated only by the substrate with a proper length are suggested.

Targeting of HPV E6 at the binding sites to the host-cell E6AP, p53, and the endoplasmic reticulum-resident chaperone, GRP78

Background: Human papillomavirus (HPV) represents an etiological factor for many cancer types, especially cervical cancer. Its oncoprotein E6 sheds drug designers who aim to stop its cellular protein associations, such as p53 and E6AP. Recently, it was discovered that the host-cell chaperone glucose-regulated protein 78 (GRP78) plays a crucial function in HPV infectivity by association with the viral E6 and E7 proteins. Therefore, we aimed to test small molecules inhibitor that could contradict the association between E6 and cellular factors E6AP, GRP78, and p53. Methods: In this study, molecular docking protocol was elaborated to test 115 small molecule compounds against the three binding sites of HPV E6 to the host-cell proteins; E6AP, p53, and GRP78. After that, molecular dynamics simulation and free energy calculations were performed on the best three complexes. Results: The results reveal the potency of 18 compounds against the HPV E6 at different binding sites, which give lower free energies than paclitaxel (positive control). The best two compounds, hypericin, and anabsinthin, could bind effectively and stably during the 100 ns MD simulation period to HPV E6. The calculated average free energies for hypericin and anabsinthin are −18.76 and −14.40 kcal/mol, respectively. They formed stable complexes with the three binding sites by forming hydrophobic contacts. The key residues that stabilize the two ligands in HPV E6 binding sites are V31, Y32, V62, and Y70 (E6AP), P13, C16, T22, I23 and A46 (p53), and M1, V31, L50, L67, and Q107 (GRP78). Conclusions: The best two compounds, hypericin, and anabsinthin, are potential candidates against HPV E6 at the host-cell factors binding sites, hence could block the oncoprotein activity of E6 in infected cells. Further experimental validation is yet to be performed and suggested as future work. Communicated by Ramaswamy H. Sarma

Phytoconstituents as potential therapeutic agents against COVID-19: a computational study on inhibition of SARS-CoV-2 main protease

The Coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) has become a global health crisis, and the urgent need for effective treatments is evident. One potential target for COVID-19 therapeutics is the main protease (Mpro) of SARS‑CoV‑2, an essential enzyme for viral replication. Natural compounds have been explored as a source of potential inhibitors for Mpro due to their safety and availability. In this study, we employed a computational approach to screen a library of phytoconstituents and identified potential Mpro inhibitors based on their binding affinities and molecular interactions. The top-ranking compounds were further validated through molecular dynamics simulations (MDS) and free energy calculations. As a result of the above procedures, we identified two phytoconstituents, Khelmarin B and Neogitogenin, with appreciable binding affinity and specificity towards the Mpro binding pocket. Our results suggest that Khelmarin B and Neogitogenin could potentially serve as Mpro inhibitors and have the potential to be developed as COVID-19 therapeutics. Further experimental studies are required to confirm the efficacy and safety of these compounds. Communicated by Ramaswamy H. Sarma

The energetics and ion coupling of cholesterol transport through Patched1.

Patched1 (PTCH1) is a tumor suppressor protein of the mammalian Hedgehog (HH) signaling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating ciliary cholesterol accessibility. Using extensive molecular dynamics simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15 to 20 kilojoule per mole for cholesterol export. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between cation binding sites, transmembrane motions, and PTCH1 activity. Using complementary simulations of Dispatched1, we find that transition between "inward-open" and solvent "occluded" states is accompanied by Na+-induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion coupling stoichiometries of PTCH1 transport mechanisms, whereby one to three Na+ or two to three K+ couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.