Mechanistic insights of key host proteins and potential repurposed inhibitors regulating SARS-CoV-2 pathway.

The emergence of pandemic situations originated from severe acute respiratory syndrome (SARS)-CoV-2 and its new variants created worldwide medical emergencies. Due to the non-availability of efficient drugs and vaccines at these emergency hours, repurposing existing drugs can effectively treat patients critically infected by SARS-CoV-2. Finding a suitable repurposing drug with inhibitory efficacy to a host-protein is challenging. A detailed mechanistic understanding of the kinetics, (dis)association pathways, key protein residues facilitating the entry-exit of the drugs with targets are fundamental in selecting these repurposed drugs. Keeping this target as the goal of the paper, the potential repurposing drugs, Nafamostat, Camostat, Silmitasertib, Valproic acid, and Zotatifin with host-proteins HDAC2, CSK22, eIF4E2 are studied to elucidate energetics, kinetics, and dissociation pathways. From an ensemble of independent simulations, we observed the presence of single or multiple dissociation pathways with varying host-proteins-drug systems and quantitatively estimated the probability of unbinding through these specific pathways. We also explored the crucial gateway residues facilitating these dissociation mechanisms. Interestingly, the residues we obtained for HDAC2 and CSK22 are also involved in the catalytic activity. Our results demonstrate how these potential drugs interact with the host machinery and the specific target residues, showing involvement in the mechanism. Most of these drugs are in the preclinical phase, and some are already being used to treat severe COVID-19 patients. Hence, the mechanistic insight presented in this study is envisaged to support further findings of clinical studies and eventually develop efficient inhibitors to treat SARS-CoV-2.

Theoretical investigations of a platinum-water interface using quantum-mechanics-molecular-mechanics based molecular dynamics simulations.

Pt-water interfaces have been of immense interest in the field of energy storage and conversion. Studying this interface using both experimental and theoretical tools is challenging. On the theoretical front, typically one uses classical molecular dynamics (MD) simulations to handle large system sizes or time scales while for a more accurate quantum mechanical description Born Oppenheimer MD (BOMD) is typically used. The latter is limited to smaller system sizes and time-scales. In this study using quantum-mechanics-molecular-mechanics (QMMM), we have performed atomistic MD simulations to have a microscopic understanding of the structure of the Pt-water interface using a system size that is much larger than that accessible when using BOMD simulations. In contrast to recent reports using BOMD simulations, our study reveals that the water molecules typically form two distinct layers above the Pt-surface before they form bulk like structures. Further, we also find that a significant fraction of the water molecules at the interface are pointed towards the surface thereby disrupting the H-bond network. Consistent with this observation, the layer resolved oxygen-oxygen radial distribution function for the water molecules belonging to the solvating water layer shows a high density liquid like behaviour even though the overall water behaves like a low density liquid. A charge transfer analysis reveals that this solvating water layer donates electrons to the Pt atoms in contact with it thereby resulting in the formation of an interface dipole that is pointing towards the surface. Our results suggest that, using QMMM-MD, on one hand it is possible to study more realistic models of solid-liquid interfaces that are inaccessible with BOMD, while on the other hand one also has access to information about such systems that are not obtained from conventional classical MD simulations.

Molecular View of CO2 Capture by Polyethylenimine: Role of Structural and Dynamical Heterogeneity.

The molecular thermodynamics and kinetics of CO2 sorption in Polyethylenimine (PEI) melt have been investigated systematically using GCMC and MD simulations. We elucidate presence of significant structural and dynamic heterogeneity associated with the overall absorption process. CO2 adsorption in a PEI membrane shows a distinct two-stage process of a rapid CO2 adsorption at the interfaces (hundreds of picoseconds) followed by a significantly slower diffusion limited release toward the interior bulk regions of PEI melt (hundreds of nanoseconds to microseconds). The spatial heterogeneity of local structural features of the PEI chains lead to significantly heterogeneous absorption characterized by clustering and trapping of CO2 molecules that then lead to subdiffusive motion of CO2. In the complex interplay of interaction and entropy, the latter emerges out to be the major determining factor with significantly higher solubility of CO2 near the interfaces despite having lower density of binding amine groups. Regions having higher free-volume (entropically favorable) viz. interfaces, pores and loops demonstrate higher CO2 capture ability. Various local structural features of PEI conformations, for example, inter- and intrachain loops, pores of different radii, and di- or tricoordinated pores are explored for their effects on the varying CO2 adsorption abilities.

Mechanism of the formation of microphase separated water clusters in a water-mediated physical network of perfluoropolyether tetraol.

Perfluoropolyether tetraol (PFPE tetraol) possesses a hydrophobic perfluoropolyether chain in the backbone and two hydroxyl groups at each chain terminal, which facilitates the formation of hydrogen bonds with water molecules resulting in the formation an extended physical network. About 3 wt% water was required for the formation of the microphase separated physical network of PFPE tetraol. The mechanism responsible for the microphase separation of water clusters in the physical network was studied using a combination of techniques such as NMR spectroscopy, molecular dynamics (MD) simulations and DSC. MD simulation studies provided evidence for the formation of clusters in the PFPE tetraol physical network and the size of these clusters increased gradually with an increase in the extent of hydration. Both MD simulations and NMR spectroscopy studies revealed that these clusters position themselves away from the hydrophobic backbone or vice versa. The presence of intra- and inter-chain aggregation possibility among hydrophilic groups was evident. DSC results demonstrated the presence of tightly and loosely bound water molecules to the terminal hydroxyl groups of PFPE tetraol through hydrogen bonding. The data from all the three techniques established the formation of a physical network driven by hydrogen bonding between the hydrophilic end groups of PFPE tetraol and water molecules. The flexible nature of the PFPE tetraol backbone and its low solubility parameter favour clustering of water molecules at the terminal groups and result in the formation of a gel.

Self-Exfoliated Guanidinium-Based Ionic Covalent Organic Nanosheets (iCONs).

Covalent organic nanosheets (CONs) have emerged as functional two-dimensional materials for versatile applications. Although π-π stacking between layers, hydrolytic instability, possible restacking prevents their exfoliation on to few thin layered CONs from crystalline porous polymers. We anticipated rational designing of a structure by intrinsic ionic linker could be the solution to produce self-exfoliated CONs without external stimuli. In an attempt to address this issue, we have synthesized three self-exfoliated guanidinium halide based ionic covalent organic nanosheets (iCONs) with antimicrobial property. Self-exfoliation phenomenon has been supported by molecular dynamics (MD) simulation as well. Intrinsic ionic guanidinium unit plays the pivotal role for both self-exfoliation and antibacterial property against both Gram-positive and Gram-negative bacteria. Using such iCONs, we have devised a mixed matrix membrane which could be useful for antimicrobial coatings with plausible medical benefits.

Free Energy of Bare and Capped Gold Nanoparticles Permeating through a Lipid Bilayer.

Herein, we study the permeation free energy of bare and octane-thiol-capped gold nanoparticles (AuNPs) translocating through a lipid membrane. To investigate this, we have pulled the bare and capped AuNPs from bulk water to the membrane interior and estimated the free energy cost. The adsorption of the bare AuNP on the bilayer surface is energetically favorable but further loading inside it requires energy. However, the estimated free-energy barrier for loading the capped AuNP into the lipid membrane is much higher compared to bare AuNP. We also demonstrate the details of the permeation process of bare and capped AuNPs. Bare AuNP induces the curvature in the lipid membrane whereas capped AuNP creates an opening in the interacting monolayer and get inserted into the membrane. The insertion of capped AuNP induces a partial unzipping of the lipid bilayer, which results in the ordering of the local lipids interacting with the nanoparticle. However, bare AuNP disrupts the lipid membrane by pushing the lipid molecules inside the membrane. We also analyze pore formation due to the insertion of capped AuNP into the membrane, which results in water molecules penetrating the hydrophobic region.

Transferability of different classical force fields for right and left handed α-helices constructed from enantiomeric amino acids.

Amino acids can form d and l enantiomers, of which the l enantiomer is abundant in nature. The naturally occurring l enantiomer has a greater preference for a right handed helical conformation, and the d enantiomer for a left handed helical conformation. The other conformations, that is, left handed helical conformations of the l enantiomers and right handed helical conformations of the d enantiomers, are not common. The energetic differences between left and right handed alpha helical peptide chains constructed from enantiomeric amino acids are investigated using quantum chemical calculations (using the M06/6-311g(d,p) level of theory). Further, the performances of commonly used biomolecular force fields (OPLS/AA, CHARMM27/CMAP and AMBER) to represent the different helical conformations (left and right handed) constructed from enantiomeric (D and L) amino acids are evaluated. 5- and 10-mer chains from d and l enantiomers of alanine, leucine, lysine, and glutamic acid, in right and left handed helical conformations, are considered in the study. Thus, in total, 32 α-helical polypeptides (4 amino acids × 4 conformations of 5-mer and 10-mer) are studied. Conclusions, with regards to the performance of the force fields, are derived keeping the quantum optimized geometry as the benchmark, and on the basis of phi and psi angle calculations, hydrogen bond analysis, and different long range helical order parameters.

Compact polar moieties induce lipid-water systems to form discontinuous reverse micellar phase.

The role of molecular interactions in governing lipid mesophase organization is of fundamental interest and has technological implications. Herein, we describe an unusual pathway for monoolein/water reorganization from a bicontinuous mesophase to a discontinuous reverse micellar assembly, directed by the inclusion of polar macromolecules. This pathway is very different from those reported earlier, wherein the Fd3m phase formed only upon addition of apolar oils. Experiments and molecular dynamics simulations indicate that hydrophilic ternary additives capable of inducing discontinuous phase formation must (i) interact strongly with the monoolein head group and (ii) have a compact molecular architecture. We present a detailed investigation that contrasts a monoolein-water system containing polyamidoamine (PAMAM) dendrons with one containing their linear analogs. The Fd3m phase forms only on the addition of PAMAM dendrons but not their linear analogs. Thus, the dendritic architecture of PAMAM plays an important role in determining lipid mesophase behavior. Both dendrons and their linear analogs interact strongly with monoolein through their amine groups. However, while linear polymers adsorb and spread on monoolein, dendrons form aggregates that interact with the lipid. Dendrons induce formation of an intermediate reverse hexagonal phase, which subsequently restructures into the Fd3m phase. Finally, we demonstrate that other additives with compact structures that are known to interact with monoolein, such as branched polyethylenimine and polyhedral silsesquioxane cages, also induce the formation of the Fd3m phase.

Proton transport mechanism of imidazole, triazole and phosphoric acid mixtures from ab initio molecular dynamics simulations.

We performed first principles molecular dynamics simulations to elucidate the mechanism and role of 1,2,3-triazole in proton transport while it is mixed with phosphoric acid (PA) and a phosphoric acid imidazole mixture. PA doped imidazole based polymer acts as an efficient polyelectrolyte membrane for fuel cells. The conductivity of this membrane increases when triazole is added to the system. For the first time we performed ab initio molecular dynamics simulations of complex mixtures of PA, imidazole and triazole. We have quantitatively estimated the structural diffusion and vehicular motion of protons. We found that upon the addition of triazole in PA and the PA imidazole mixture, the structural diffusion of protons increases significantly. The mechanism of proton transport is different when triazole is added to the mixture. We have also identified two different paths for structural diffusion (constructive and non-constructive) that contribute to long and short range proton transport.

Self-assembly of phospholipids on flat supports.

The current study deals with the self-assembly of phospholipids on flat supports using the Martini coarse grain model. We reported here the effect of the hydrophilic and hydrophobic nature of the solid supports on the lipid self-assembly. The hydrophilic and hydrophobic supports were modeled on the basis of water droplet simulations. The present work addresses the self-assembly mechanism of lipids on eight different supports with different strengths of hydrophilicity and hydrophobicity. We demonstrated how interplay between the interactions of lipid and water with the support can guide the lipid self-assembly process. Thereafter, we calculated the energetics of the components of the system to quantify the competitions between water and a lipid head-group with hydrophilic supports. Finally, the properties of the self-assembled bilayers were also analyzed and reported here.

A molecular dynamics study of model SI clathrate hydrates: the effect of guest size and guest-water interaction on decomposition kinetics.

One of the options suggested for methane recovery from natural gas hydrates is molecular replacement of methane by suitable guests like CO2 and N2. This approach has been found to be feasible through many experimental and molecular dynamics simulation studies. However, the long term stability of the resultant hydrate needs to be evaluated; the decomposition rate of these hydrates is expected to depend on the interaction between these guest and water molecules. In this work, molecular dynamics simulation has been performed to illustrate the effect of guest molecules with different sizes and interaction strengths with water on structure I (SI) hydrate decomposition and hence the stability. The van der Waals interaction between water of hydrate cages and guest molecules is defined by Lennard Jones potential parameters. A wide range of parameter spaces has been scanned by changing the guest molecules in the SI hydrate, which acts as a model gas for occupying the small and large cages of the SI hydrate. All atomistic simulation results show that the stability of the hydrate is sensitive to the size and interaction of the guest molecules with hydrate water. The increase in the interaction of guest molecules with water stabilizes the hydrate, which in turn shows a slower rate of hydrate decomposition. Similarly guest molecules with a reasonably small (similar to Helium) or large size increase the decomposition rate. The results were also analyzed by calculating the structural order parameter to understand the dynamics of crystal structure and correlated with the release rate of guest molecules from the solid hydrate phase. The results have been explained based on the calculation of potential energies felt by guest molecules in amorphous water, hydrate bulk and hydrate-water interface regions.

Food restriction increases glutamate receptor-mediated burst firing of dopamine neurons.

Restriction of food intake increases the acquisition of drug abuse behavior and enhances the reinforcing efficacy of those drugs. However, the neurophysiological mechanisms responsible for the interactions between feeding state and drug use are largely unknown. Here we show that chronic mild food restriction increases the burst firing of dopamine neurons in the substantia nigra. Dopamine neurons from food-restricted mice exhibited increased burst firing in vivo, an effect that was enhanced by an injection of the psychomotor stimulant cocaine (10 mg/kg, i.p.). Food restriction also enhanced aspartic acid-induced burst firing of dopamine neurons in an ex vivo brain slice preparation, consistent with an adaptation occurring in the somatodendritic compartment and independent of a circuit mechanism. Enhanced burst firing persisted after 10 d of free feeding following chronic food restriction but was not observed following a single overnight fast. Whole-cell patch-clamp recordings indicated that food restriction also increased electrically evoked AMPAR/NMDAR ratios and increased D2 autoreceptor-mediated desensitization in dopamine neurons. These results identify dopamine neurons in the substantia nigra as a convergence point for the interactions between feeding state and drugs of abuse. Furthermore, increased glutamate transmission combined with decreased autoreceptor inhibition could work in concert to enhance drug efficacy in response to food restriction.

Vesicle structures from bolaamphiphilic biosurfactants: experimental and molecular dynamics simulation studies on the effect of unsaturation on sophorolipid self-assemblies.

The formation of giant-vesicle-like structures by self-assembling linolenic acid sophorolipid (LNSL) molecules is revealed. Sophorolipids belong to the class of bolaamphiphilic glycolipid biosurfactants. Interestingly, the number of double bonds present in the hydrophobic core of sophorolipids is seen to have a great influence on the type of self-assembled structures formed. Dye encapsulation results establish the presence of an aqueous compartment inside the LNSL vesicles. Molecular dynamics simulation (MD) studies suggest the existence of two possible conformations of LNSLs inside the self-assembled structures and that LNSL molecules arrange in layered structures.

Model atomistic protrusions favouring the ordering and retention of water.

The ordering of water molecules near model linear atomistic protrusions is studied using classical molecular dynamics simulations. The protrusions are made up of Lennard-Jones particles of hydrophobic and hydrophilic blocks. Simulations are performed at a range of temperatures and pressures, keeping the position of the protrusions fixed. At different temperatures and pressures, the ordering and residence time of water molecules is enhanced on the surface of the hydrophilic block. Detailed analysis of the systems shows that the surface region is potentially the most energetically favorable for water molecules, which is consistent with the tetrahedral ordering of water molecules. A competition between energetics and structuring is observed from residence time calculations.

A computational study to probe the binding aspects of potent polyphenolic inhibitors of pancreatic lipase.

Pancreatic lipase (PL) is a keen target for anti-obesity therapy that reduces dietary fat absorption. Here, we investigated the binding patterns of 220 PL inhibitors having experimental IC50 values, using molecular docking and binding energy calculations. Screening of these compounds illustrated most of them bound at the catalytic site (S1-S2 channel) and a few compounds are at the non-catalytic site (S2-S3 channel/S1-S3 channel) of PL. This binding pattern could be due to structural uniqueness or bias in conformational search. A strong correlation of pIC50 values with SP/XP docking scores, binding energies (ΔGMMGBSA) assured the binding poses are more true positives. Further, understanding of each class and subclasses of polyphenols indicated tannins preferred non-catalytic site wherein binding energies are underestimated due to huge desolvation energy. In contrast, most of the flavonoids and furan-flavonoids have good binding energies due to strong interactions with catalytic residues. While scoring functions limited the understanding of sub-classes of flavonoids. Hence, focused on 55 potent PL inhibitors of IC50 < 5 µM for better in vivo efficacy. The prediction of bioactivity, drug-likeness properties, led to 14 bioactive compounds. The low root mean square deviation (0.1-0.2 nm) of these potent flavonoids and non-flavonoid/non-polyphenols PL-inhibitor complexes during 100 ns molecular dynamics runs (MD) as well as binding energies obtained from both MD and well-tempered metadynamics, support strong binding to catalytic site. Based on the bioactivity, ADMET properties, and binding affinity data of MD and wt-metaD of potent PL-inhibitors suggests Epiafzelechin 3-O-gallate, Sanggenon C, and Sanggenofuran A shall be promising inhibitors at in vivo conditions.Communicated by Ramaswamy H. Sarma.

Characterization of conformation and interaction of gene delivery vector polyethylenimine with phospholipid bilayer at different protonation state.

Polyethylenimine (PEI) is a pH sensitive polymer possessing stretched and coiled conformation at low and high pH, respectively. It is an efficient gene delivery agent. Thus, the interaction of PEI with the biomembrane is very crucial to understand the gene delivery mechanism. In this report, we have investigated the structural properties of PEI and bilayer due to the interaction of PEI with lipid molecules. PEI has coil structure at high pH while at low pH it is elongated. The neutral PEI chain predominately settles itself at the bilayer water interface. We do not find any disruption or pore formation on the bilayer due to interaction of neutral PEI chain. PEI at low pH gets elongated due to electrostatic interaction between charges of the protonated sites. This protonated PEI chain interacts with bilayer membrane, which leads to formation of water/ion channel through the membrane. We have analyzed the structure of the channel and water dynamics along the channel.

Evidence and characterization of dynamic heterogeneity in binary mixtures of phosphoric acid and benzimidazole.

We report here anomalous diffusions of components in mixtures of monomer of polybenzimidazole, i.e., 2-phenyl-1H,1'H-5,5'-bibenzo[d]imidazole (BI) and phosphoric acid (PA) from molecular dynamics simulations. We have observed initial drop and further increase in self-diffusion constant for both monomer molecule (BI) and PA with gradual increase in PA concentration. The origin of such anomalous diffusion is identified in this work, which happens to be the presence of dynamic heterogeneity in each component of the binary mixture. We characterized microscopic picture of dynamical heterogeneity by finding correlation between dynamical heterogeneity and structural arrangement among the components of the binary system. Different types of H-bonding arrangements in the BI-PA systems at different concentration of PA are observed. The stability of the H-bonded network consisting of different types of H-bonds between BI and PA in the system has been studied by calculating the lifetime of various H-bonds. The results indicate that there are fast and slow moving PA molecules in the mixtures because of coexistence of different types of hydrogen bonds among the components of the mixture.

High-Risk Patient Experiences Associated With an Intensive Primary Care Management Program in the Veterans Health Administration.

Intensive management programs may improve health care experiences among high-risk and complex patients. We assessed patient experience among (1) prior enrollees (n = 59) of an intensive management program (2014-2018); (2) nonenrollees (n = 356) at program sites; and (3) nonprogram site patients (n = 728), using a patient survey based on the Consumer Assessment of Healthcare Providers and Systems in 2019. Outcomes included patient ratings of patient-centered care; overall health care experience; and satisfaction with their usual outpatient care provider. In multivariate models, enrollees were more satisfied with their current provider versus nonenrollees within program sites (adjusted odds ratio 2.36; 95% confidence interval 1.15-4.85).

Clinical Evaluation of Efficacy of Triamcinolone Acetonide with Tacrolimus in the Management of Oral Lichen Planus: A Pilot Prospective Observational Study.

Introduction: Lichen planus (LP) is a relatively common chronic, mucocutaneous disease of autoimmune origin, involves oral mucosa, skin, scalp, nails, and genital mucosa. The prevalence of oral LP (OLP) varies worldwide, commonly seen in middle-aged and elderly people. It usually presents as symmetrical and bilateral or multiple lesions with burning sensation (BS) sometimes accompanied by pain. Corticosteroids and calcineurin inhibitors have shown promising results in the treatment of OLP, but its chronic course and unpredictable exacerbations/remission continues to result in a high degree of morbidity. The study aimed to evaluate the efficacy of intralesional triamcinolone acetonide (injection TA) combined with topical application of TA orabase and Tacrolimus (TAC) ointment for symptomatic cases of OLP. Materials and Methods: The prospective study included 52 symptomatic OLP patients to receive (0.5 ml) intralesional injection of TA once a week for the first 4 weeks followed by one injection in the 6th week along with TA mucosal paste (0.1%.) and TAC ointment (0.03%) in tapering dose till 8th week. The subjective symptoms including BS and pain were assessed on a 10 cm visual analog scale (VAS) and objective signs like size and site of the lesion were scored according to criterion scale modified by Thongprasom et al. Differences were compared after 8 weeks treatment course and follow-up observations were performed at 20th week to record any recurrent lesion. Results: 41 patients (78.8%) had complete remission of disease and 11 (21%) had shown partial improvement. The VAS scores for BS and pain improved significantly. Improvement was also noted with decrease in the average size of active lesions and the number of sites with treatment. The relapse was seen in 17 patients (41%) in the 20th week. Conclusion: TA combined with TAC is a valuable therapeutic option for the treatment of symptomatic OLP. Our findings suggest that patients have shown statistically significant improvement.

Integrated docking and enhanced sampling-based selection of repurposing drugs for SARS-CoV-2 by targeting host dependent factors.

Since the onset of global pandemic, the most focused research currently in progress is the development of potential drug candidates and clinical trials of existing FDA approved drugs for other relevant diseases, in order to repurpose them for the COVID-19. At the same time, several high throughput screenings of drugs have been reported to inhibit the viral components during the early course of infection but with little proven efficacies. Here, we investigate the drug repurposing strategies to counteract the coronavirus infection which involves several potential targetable host proteins involved in viral replication and disease progression. We report the high throughput analysis of literature-derived repurposing drug candidates that can be used to target the genetic regulators known to interact with viral proteins based on experimental and interactome studies. In this work we have performed integrated molecular docking followed by molecular dynamics (MD) simulations and free energy calculations through an expedite in silico process where the number of screened candidates reduces sequentially at every step based on physicochemical interactions. We elucidate that in addition to the pre-clinical and FDA approved drugs that targets specific regulatory proteins, a range of chemical compounds (Nafamostat, Chloramphenicol, Ponatinib) binds to the other gene transcription and translation regulatory proteins with higher affinity and may harbour potential for therapeutic uses. There is a rapid growing interest in the development of combination therapy for COVID-19 to target multiple enzymes/pathways. Our in silico approach would be useful in generating leads for experimental screening for rapid drug repurposing against SARS-CoV-2 interacting host proteins.Communicated by Ramaswamy H. Sarma.