Influence of the Water Model on the Structure and Interactions of the GPR40 Protein with the Lipid Membrane and the Solvent: Rigid versus Flexible Water Models.

G protein-coupled receptors (GPCR) are responsible for modulating various physiological functions and are thus related to the pathophysiology of different diseases. Being potential therapeutic targets, multiple computational methodologies have been developed to analyze their behavior and interactions with other species. The solvent, on the other hand, has received much less attention. In this work, we analyzed the effect of four explicit water models on the structure and interactions of the GPR40 receptor in its apo form. We employed the rigid SPC/E and TIP4P models, and their flexible versions, the FBA/ϵ and TIP4P/ϵflex. We explored the structural changes and their correlation with some bulk dynamic properties of water. Our results showed an adverse effect on the conservation of the secondary structure of the receptor with all the models due to the breaking of the intramolecular hydrogen bond network, being more evident for the TIP4P models. Notably, all four models brought the receptor to states similar to the active one, modifying the intracellular part of the TM5 and TM6 domains in a "hinge" type movement, allowing the opening of the structure. Regarding the dynamic properties, the rigid models showed results comparable to those obtained in other studies on membrane systems. However, flexible models exhibit disparities in the molecular representation of systems. Surprisingly, the FBA/ϵ model improves the molecular picture of several properties, even though their agreement with bulk diffusion is poorer. These findings reinforce our idea that exploring other water models or improving the current ones, to better represent the membrane interface, can lead to a positive impact on the description of the signal transduction mechanisms and the search of new drugs by targeting these receptors.

Co-occurring pathogenic variants in 6q27 associated with dementia spectrum disorders in a Peruvian family.

Evidence suggests that there may be racial differences in risk factors associated with the development of Alzheimer's disease and related dementia (ADRD). We used whole-genome sequencing analysis and identified a novel combination of three pathogenic variants in the heterozygous state (UNC93A: rs7739897 and WDR27: rs61740334; rs3800544) in a Peruvian family with a strong clinical history of ADRD. Notably, the combination of these variants was present in two generations of affected individuals but absent in healthy members of the family. In silico and in vitro studies have provided insights into the pathogenicity of these variants. These studies predict that the loss of function of the mutant UNC93A and WDR27 proteins induced dramatic changes in the global transcriptomic signature of brain cells, including neurons, astrocytes, and especially pericytes and vascular smooth muscle cells, indicating that the combination of these three variants may affect the neurovascular unit. In addition, known key molecular pathways associated with dementia spectrum disorders were enriched in brain cells with low levels of UNC93A and WDR27. Our findings have thus identified a genetic risk factor for familial dementia in a Peruvian family with an Amerindian ancestral background.

Study of the Rv1417 and Rv2617c Membrane Proteins and Their Interactions with Nicotine Derivatives as Potential Inhibitors of Erp Virulence-Associated Factor in Mycobacterium tuberculosis: An In Silico Approach.

The increasing emergence of Mycobacterium tuberculosis (Mtb) strains resistant to traditional anti-tuberculosis drugs has alarmed health services worldwide. The search for new therapeutic targets and effective drugs that counteract the virulence and multiplication of Mtb represents a challenge for the scientific community. Several studies have considered the erp gene a possible therapeutic target in the last two decades, since its disruption negatively impacts Mtb multiplication. This gene encodes the exported repetitive protein (Erp), which is located in the cell wall of Mtb. In vitro studies have shown that the Erp protein interacts with two putative membrane proteins, Rv1417 and Rv2617c, and the impairment of their interactions can decrease Mtb replication. In this study, we present five nicotine analogs that can inhibit the formation of heterodimers and trimers between these proteins. Through DFT calculations, molecular dynamics, docking, and other advanced in silico techniques, we have analyzed the molecular complexes, and show the effect these compounds have on protein interactions. The results show that four of these analogs can be possible candidates to counteract the pathogenicity of Mtb. This study aims to combine research on the Erp protein as a therapeutic target in the search for new drugs that serve to create new therapies against tuberculosis disease.

Novel In Silico Insights into Rv1417 and Rv2617c as Potential Protein Targets: The Importance of the Medium on the Structural Interactions with Exported Repetitive Protein (Erp) of Mycobacterium tuberculosis.

Nowadays, tuberculosis is the second leading cause of death from a monopathogenic transmitted disease, only ahead of COVID-19. The role of exported repetitive protein (Erp) in the virulence of Mycobacterium tuberculosis has been extensively demonstrated. In vitro and in vivo assays have identified that Erp interacts with Rv1417 and Rv2617c proteins, forming putative transient molecular complexes prior to localization to the cell envelope. Although new insights into the interactions and functions of Erp have emerged over the years, knowledge about its structure and protein-protein interactions at the atomistic level has not been sufficiently explored. In this work, we have combined several in silico methodologies to gain new insights into the structural relationship between these proteins. Two system conditions were evaluated by MD simulations: Rv1417 and Rv2617c embedded in a lipid membrane and another with a semi-polar solvent to mimic the electrostatic conditions on the membrane surface. The Erp protein was simulated as an unanchored structure. Stabilized structures were docked, and complexes were evaluated to recognize the main residues involved in protein-protein interactions. Our results show the influence of the medium on the structural conformation of proteins. Globular conformations were favored under high polarity conditions and showed a higher energetic affinity in complex formation. Meanwhile, disordered conformations were favored under semi-polar conditions and an increase in the number of contacts between residues was observed. In addition, the electrostatic potential analysis showed remarkable changes in protein interactions due to the polarity of the medium, demonstrating the relevance of Erp protein in heterodimer formation. On the other hand, contact analysis showed that several C-terminal residues of Erp were involved in the protein interactions, which seems to contradict experimental observations; however, these complexes could be transient forms. The findings presented in this work are intended to open new perspectives in the studies of Erp protein molecular interactions and to improve the knowledge about its function and role in the virulence of Mycobacterium tuberculosis.

In Silico Analysis of the Antagonist Effect of Enoxaparin on the ApoE4-Amyloid-Beta (Aß) Complex at Different pH Conditions.

Apolipoprotein E4 (ApoE4) is thought to increase the risk of developing Alzheimer's disease. Several studies have shown that ApoE4-Amyloid ß (Aß) interactions can increment amyloid depositions in the brain and that this can be augmented at low pH values. On the other hand, experimental studies in transgenic mouse models have shown that treatment with enoxaparin significantly reduces cortical Aß levels, as well as decreases the number of activated astrocytes around Aß plaques. However, the interactions between enoxaparin and the ApoE4-Aß proteins have been poorly explored. In this work, we combine molecular dynamics simulations, molecular docking, and binding free energy calculations to elucidate the molecular properties of the ApoE4-Aß interactions and the competitive binding affinity of the enoxaparin on the ApoE4 binding sites. In addition, we investigated the effect of the environmental pH levels on those interactions. Our results showed that under different pH conditions, the closed form of the ApoE4 protein, in which the C-terminal domain folds into the protein, remains stabilized by a network of hydrogen bonds. This closed conformation allowed the generation of six different ApoE4-Aß interaction sites, which were energetically favorable. Systems at pH5 and 6 showed the highest energetic affinity. The enoxaparin molecule was found to have a strong energetic affinity for ApoE4-interacting sites and thus can neutralize or disrupt ApoE4-Aß complex formation.

In Silico Analysis of Metabolites from Peruvian Native Plants as Potential Therapeutics against Alzheimer's Disease.

BACKGROUND: Despite research on the molecular bases of Alzheimer's disease (AD), effective therapies against its progression are still needed. Recent studies have shown direct links between AD progression and neurovascular dysfunction, highlighting it as a potential target for new therapeutics development. In this work, we screened and evaluated the inhibitory effect of natural compounds from native Peruvian plants against tau protein, amyloid beta, and angiotensin II type 1 receptor (AT1R) pathologic AD markers. METHODS: We applied in silico analysis, such as virtual screening, molecular docking, molecular dynamics simulation (MD), and MM/GBSA estimation, to identify metabolites from Peruvian plants with inhibitory properties, and compared them to nicotinamide, telmisartan, and grapeseed extract drugs in clinical trials. RESULTS: Our results demonstrated the increased bioactivity of three plants' metabolites against tau protein, amyloid beta, and AT1R. The MD simulations indicated the stability of the AT1R:floribundic acid, amyloid beta:rutin, and tau:brassicasterol systems. A polypharmaceutical potential was observed for rutin due to its high affinity to AT1R, amyloid beta, and tau. The metabolite floribundic acid showed bioactivity against the AT1R and tau, and the metabolite brassicasterol showed bioactivity against the amyloid beta and tau. CONCLUSIONS: This study has identified molecules from native Peruvian plants that have the potential to bind three pathologic markers of AD.

Structural and Functional Analysis of Female Sex Hormones against SARS-CoV-2 Cell Entry.

Emerging evidence suggests that males are more susceptible to severe infection by the SARS-CoV-2 virus than females. A variety of mechanisms may underlie the observed gender-related disparities including differences in sex hormones. However, the precise mechanisms by which female sex hormones may provide protection against SARS-CoV-2 infectivity remains unknown. Here we report new insights into the molecular basis of the interactions between the SARS-CoV-2 spike (S) protein and the human ACE2 receptor. We further report that glycosylation of the ACE2 receptor enhances SARS-CoV-2 infectivity. Importantly, estrogens can disrupt glycan-glycan interactions and glycan-protein interactions between the human ACE2 and the SARS-CoV-2 thereby blocking its entry into cells. In a mouse model of COVID-19, estrogens reduced ACE2 glycosylation and thereby alveolar uptake of the SARS-CoV-2 spike protein. These results shed light on a putative mechanism whereby female sex hormones may provide protection from developing severe infection and could inform the development of future therapies against COVID-19.

Structural and functional analysis of female sex hormones against SARS-Cov2 cell entry.

Emerging evidence suggests that males are more susceptible to severe infection by the SARS-CoV-2 virus than females. A variety of mechanisms may underlie the observed gender-related disparities including differences in sex hormones. However, the precise mechanisms by which female sex hormones may provide protection against SARS-CoV-2 infectivity remains unknown. Here we report new insights into the molecular basis of the interactions between the SARS-CoV-2 spike (S) protein and the human ACE2 receptor. We further observed that glycosylation of the ACE2 receptor enhances SARS-CoV-2 infectivity. Importantly estrogens can disrupt glycan-glycan interactions and glycan-protein interactions between the human ACE2 and the SARS-CoV2 thereby blocking its entry into cells. In a mouse model, estrogens reduced ACE2 glycosylation and thereby alveolar uptake of the SARS-CoV-2 spike protein. These results shed light on a putative mechanism whereby female sex hormones may provide protection from developing severe infection and could inform the development of future therapies against COVID-19.

Force Field Parametrization from the Hirshfeld Molecular Electronic Density.

The Hirshfeld charges are linearly increased to reproduce the experimental dielectric constant of 10 polar solvents having values between 13 (pyridine) and 182 ( N-methylformamide). The OPLS/AA force field is used to obtain the new parameters. The surface tension and liquid density are also target properties to determine the new nonbonding parameters. The charge scaling factor is between 1.2 and 1.3. In addition, properties that were not used in the parametrization procedure, such as the heat of vaporization, self-diffusion coefficient, shear viscosity, isothermal compressibility, and volumetric expansion coefficient are obtained. Binary mixtures of amide/water and amide/amide are also studied. The original parameters of OPLS/AA, CGenFF, and GAFF force fields are evaluated. The TIP4P/ε force field is used to simulate water. The results from this work with the new parameters, for both pure components and binary mixtures, are in better agreement with experimental data than those obtained with the original values for most of the calculated properties. The maximum density of N-methylformamide in aqueous solutions is correctly predicted only with the new parameters. The high value of the dielectric constant of acetamide, formamide, and N-methylformamide is discussed in terms of the chain formation from the hydrogen bond interactions.

Force Field Benchmark of the TraPPE_UA for Polar Liquids: Density, Heat of Vaporization, Dielectric Constant, Surface Tension, Volumetric Expansion Coefficient, and Isothermal Compressibility.

The transferable potential for a phase equilibria force field in its united-atom version, TraPPE_UA, is evaluated for 41 polar liquids that include alcohols, thiols, ethers, sulfides, aldehydes, ketones, and esters to determine its ability to reproduce experimental properties that were not included in the parametrization procedure. The intermolecular force field parameters for pure components were fit to reproduce experimental boiling temperature, vapor-liquid coexisting densities, and critical point (temperature, density, and pressure) using Monte Carlo simulations in different ensembles. The properties calculated in this work are liquid density, heat of vaporization, dielectric constant, surface tension, volumetric expansion coefficient, and isothermal compressibility. Molecular dynamics simulations were performed in the gas and liquid phases, and also at the liquid-vapor interface. We found that relative error between calculated and experimental data is 1.2% for density, 6% for heat of vaporization, and 6.2% for surface tension, in good agreement with the experimental data. The dielectric constant is systematically underestimated, and the relative error is 37%. Evaluating the performance of the force field to reproduce the volumetric expansion coefficient and isothermal compressibility requires more experimental data.