No abstract

No abstract

No abstract

Cold-adapted enzymes feature a lower thermostability and higher catalytic activity compared to their warm-active homologues, which are considered as a consequence of increased flexibility of their molecular structures. The complexity of the (thermo)stability-flexibility-activity relationship makes it difficult to define the strategies and formulate a general theory for enzyme cold adaptation. Here, the psychrophilic serine hydroxymethyltransferase (pSHMT) from Psychromonas ingrahamii and its mesophilic counterpart, mSHMT from Escherichia coli, were subjected to μs-scale multiple-replica molecular dynamics (MD) simulations to explore the cold-adaptation mechanism of the dimeric SHMT. The comparative analyses of MD trajectories reveal that pSHMT exhibits larger structural fluctuations and inter-monomer positional movements, a higher global flexibility, and considerably enhanced local flexibility involving the surface loops and active sites. The largest-amplitude motion mode of pSHMT describes the trends of inter-monomer dissociation and enlargement of the active-site cavity, whereas that of mSHMT characterizes the opposite trends. Based on the comparison of the calculated structural parameters and constructed free energy landscapes (FELs) between the two enzymes, we discuss in-depth the physicochemical principles underlying the stability-flexibility-activity relationships and conclude that (i) pSHMT adopts the global-flexibility mechanism to adapt to the cold environment and, (ii) optimizing the protein-solvent interactions and loosening the inter-monomer association are the main strategies for pSHMT to enhance its flexibility.

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein (RBDCoV2) has a higher binding affinity to the human receptor angiotensin-converting enzyme 2 (ACE2) than the SARS-CoV RBD (RBDCoV). Here, we performed molecular dynamics (MD) simulations, binding free energy (BFE) calculations, and interface residue contact network (IRCN) analysis to explore the mechanistic origin of different ACE2-binding affinities of the two RBDs. The results demonstrate that, when compared to the RBDCoV2-ACE2 complex, RBDCoV-ACE2 features enhanced dynamicsand inter-protein positional movements and increased conformational entropy and conformational diversity. Although the inter-protein electrostatic attractive interactions are the primary determinant for the high ACE2-binding affinities of both RBDs, the significantly enhanced electrostatic attractive interactions between ACE2 and RBDCoV2 determine the higher ACE2-binding affinity of RBDCoV2 than of RBDCoV. Comprehensive comparative analyses of the residue BFE components and IRCNs between the two complexes reveal that it is the residue changes at the RBD interface that lead to the overall stronger inter-protein electrostatic attractive force in RBDCoV2-ACE2, which not only tightens the interface packing and suppresses the dynamics of RBDCoV2-ACE2, but also enhances the ACE2-binding affinity of RBDCoV2. Since the RBD residue changes involving gain/loss of the positively/negatively charged residues can greatly enhance the binding affinity, special attention should be paid to the SARS-CoV-2 variants carrying such mutations, particularly those near or at the binding interfaces with the potential to form hydrogen bonds and/or salt bridges with ACE2.

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein mediates viral entry into host cells through binding to the cell-surface receptor angiotensin-converting enzyme 2 (ACE2). It has been shown that SARS-CoV-2 RBD (RBDCoV2) has a higher binding affinity to human ACE2 than its highly homologous SARS-CoV RBD (RBDCoV), for which the mechanistic reasons still remain to be elucidated. Here, we used the multiple-replica molecular dynamics (MD) simulations, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding free energy calculations, and interface residue contact network (IRCN) analysis approach to explore the mechanistic origin of different ACE2 binding affinities of these two RBDs. The results demonstrate that, when compared to the RBDCoV2-ACE2 complex, the RBDCoV-ACE2 complex features the enhanced overall structural fluctuations and inter-protein positional movements and increased conformational entropy and diversity. The inter-protein electrostatic attractive interactions are a dominant force in determining the high ACE2 affinities of both RBDs, while the significantly strengthened electrostatic forces of attraction of ACE2 to RBDCoV2 determine the higher ACE2 binding affinity of RBDCoV2 than of RBDCoV. Comprehensive comparative analyses of the residue binding free energy components and IRCNs reveal that, although any RBD residue substitution involved in the charge change can significantly impact the inter-protein electrostatic interaction strength, it is the substitutions at the RBD interface that lead to the overall stronger electrostatic attractive force of RBDCoV2-ACE2, which in turn not only tightens the interface packing and suppresses the dynamics of RBDCoV2-ACE2, but also enhances the ACE2 binding affinity of RBDCoV2 compared to that of RBDCoV. Since the RBD residue substitutions involving gain/loss of the positively/negatively charged residues, in particular those near/at the binding interfaces with the potential to form hydrogen bonds and/or salt bridges with ACE2, can greatly enhance the ACE2 binding affinity, the SARS-CoV-2 variants carrying such mutations should be paid special attention to.

BackgroundA number of studies have demonstrated the significant and rapid antidepressant effects of ketamine, which is also known as a neurotoxic and illicit drug. Ketamine and alcohol are increasingly used together in clubs by teenagers and young adults. Previous studies have proven that chronic ketamine consumption induces a delayed and persistent activation of the dopamine (DA) system. However, the rewarding properties of recreational ketamine abuse remain unclear, and the underlying mechanisms of the effects on the DA system after administration of ketamine with ethanol are yet to be explored.MethodsHere, we evaluated the effects of two different doses of ketamine (30 mg/kg and 60 mg/kg) with and without ethanol (0.3156 g/kg) on DA concentration in the rat’s ventral tegmental area (VTA), a vital region in the reward and motivation system. We explored the effects of the combined drug treatment on the expression profiling of the DA metabolism genes, tyrosine hydroxylase, dopa decarboxylase, vesicular monoamine transporter 2, and synaptosomal-associated protein 25, as well as protein expression level of brain-derived neurotrophic factor in the rat’s VTA.ResultsWe found that administration of ketamine with ethanol led to a significant increase of DA in the VTA associated with differential regulation of mRNA levels of the four DA metabolism genes and protein levels of brain-derived neurotrophic factor. Moreover, the rewarding properties of coadministration of ketamine and ethanol were related to dopaminergic neuron activation in the VTA.ConclusionThese results indicated the possibility that combined drug treatment might positively affect the mesencephalic DA reward system.

The indispensable role of the SARS-CoV-2 main protease (Mpro) in the viral replication cycle and its dissimilarity to human proteases make Mpro a promising drug target. In order to identify the non-covalent Mpro inhibitors, we performed a comprehensive study using a combined computational strategy. We first screened the ZINC purchasable compound database using the pharmacophore model generated from the reference crystal structure of Mpro complexed with the inhibitor ML188. The hit compounds were then filtered by molecular docking and predicted parameters of drug-likeness and pharmacokinetics. The final molecular dynamics (MD) simulations identified three effective candidate inhibitors (ECIs) capable of maintaining binding within the substrate-binding cavity of Mpro. We further performed comparative analyses of the reference and effective complexes in terms of dynamics, thermodynamics, binding free energy (BFE), and interaction energies and modes. The results reveal that, when compared to the inter-molecular electrostatic forces/interactions, the inter-molecular van der Waals (vdW) forces/interactions are far more important in maintaining the association and determining the high affinity. Given the un-favorable effects of the inter-molecular electrostatic interactions—association destabilization by the competitive hydrogen bond (HB) interactions and the reduced binding affinity arising from the un-compensable increase in the electrostatic desolvation penalty—we suggest that enhancing the inter-molecular vdW interactions while avoiding introducing the deeply buried HBs may be a promising strategy in future inhibitor optimization.