Full scale structural, mechanical and dynamical properties of HIV-1 liposomes

Structural maturation of the matrix lattice is not required for HIV-1 particle infectivity

During HIV-1 maturation, the matrix (MA) lattice underlying the viral membrane undergoes a structural rearrangement, and the newly released capsid (CA) protein forms a mature CA. While it is well established that CA formation is essential for particle infectivity, the functional role of MA structural maturation remains unclear. Here, we examine maturation of an MA triple mutant, L20K/E73K/A82T, which, despite replicating similarly to wild-type (WT) in some cell lines, exhibits distinct biochemical behaviors that suggest altered MA-MA interactions. Cryo-electron tomography with subtomogram averaging reveals that, although the MA lattice in immature L20K/E73K/A82T virions closely resembles that of the WT, mature L20K/E73K/A82T virions lack a detectable MA lattice. All-atom molecular dynamics simulations suggest that this absence results from destabilized inter-trimer MA interactions in mature L20K/E73K/A82T mutant virions. These findings suggest that an ordered, membrane-associated mature MA lattice is not essential for HIV-1 infectivity, providing insights into the structural requirements for HIV-1 particle maturation and generation of infectious particles.

Structural maturation of the matrix lattice is not required for HIV-1 particle infectivity

HIV-1 assembly is initiated by the binding of Gag polyproteins to the inner leaflet of the plasma membrane, mediated by the myristylated matrix (MA) domain of Gag. Subsequent to membrane binding, Gag oligomerizes and buds as an immature, non-infectious virus particle, which, upon cleavage of the Gag precursor by the viral protease, transforms into a mature, infectious virion. During maturation, the MA lattice underlying the viral membrane undergoes a structural rearrangement and the newly released capsid (CA) protein forms a mature capsid that encloses the viral genome. While it is well established that formation of the mature capsid is essential to particle infectivity, the functional role of MA structural maturation remains unclear. Here, we examine MA maturation of an MA triple mutant, L20K/E73K/A82T, which exhibits distinct biochemical behaviours. The L20K/E73K/A82T mutant is a revertant derived by propagating the L20K mutant, which exhibits reduced infectivity and increased association of the Gag polyprotein with membranes. L20K/E73K/A82T replicates similarly to wild type but retains the increased Gag membrane binding properties of L20K. L20K/E73K/A82T MA also sediments to high-density fractions in sucrose gradients after detergent treatment under conditions that fully solubilize WT MA, suggesting enhanced MA-MA interactions. Cryo-electron tomography with subtomogram averaging reveals that the immature MA lattice of L20K/E73K/A82T closely resembles the wild type. However, mature virions of the triple mutant lack a detectable MA lattice, in stark contrast to both the wild type and L20K mutant. All-atom molecular dynamics simulations suggest that this absence results from destabilized inter-trimer interactions in the mature L20K/E73K/A82T MA. Furthermore, introducing additional mutations designed to disrupt the mature MA lattice does not impair particle infectivity. These findings suggest that an ordered, membrane-associated mature MA lattice is not essential for HIV-1 infectivity, providing new insights into the structural plasticity of the matrix during maturation and its functional role in the viral lifecycle.

A comprehensive insight from molecular docking and dynamics with clinical investigation on the impact of direct oral anticoagulants on atheroprotective protein in atrial fibrillation

BACKGROUND

Direct oral anticoagulants (DOACs) have high potency against their therapeutic target and are widely used in the treatment of atrial fibrillation (AF). Most DOACs are often claimed to have adverse effects due to off-target inhibition of essential proteins. Human serum paraoxonase 1 (PON1), one of the essential proteins, known for its anti-inflammatory and antioxidant properties, could be affected by DOACs. Thus, a comparative evaluation of DOACs and their effect on PON1 protein will aid in recommending the most effective DOACs for AF treatment. This study aimed to assess the impact of DOACs on PON1 through a combination of computational and experimental analyses.

METHODS

We focus on apixaban, dabigatran, and rivaroxaban, the most recommended DOACs in AF treatment, for their impact on PON1 through molecular docking and molecular dynamics (MD) simulation to elucidate the binding affinity and drug-protein structural stability. This investigation revealed the most influential DOACs on the PON1 protein. Then experimental validation was performed in DOAC-treated AF participants (n = 42; 19 treated with dabigatran and 23 treated with rivaroxaban) compared to a healthy control group (n = 22) through gene expression analysis in peripheral blood mononuclear cells (PBMC) and serum enzyme concentration.

RESULTS

Our computational investigation showed rivaroxaban (-4.24 kcal/mol) exhibited a lower affinity against the PON1 protein compared to apixaban (-5.97 kcal/mol) and dabigatran (-9.03 kcal/mol) through molecular docking. Dabigatran holds complex interactions with PON1 at GLU53, TYR197, SER193, and ASP269 by forming hydrogen bonds. Additionally, MD simulation revealed that dabigatran disrupts PON1 stability, which may contribute functional changes. Further experimental validation revealed a significant down-regulation (p < 0.05) of PON1 gene expression in PBMC and decreased serum PON1 enzyme concentration on DOAC treatment. Rivaroxaban as about 48% has inhibitory percentage and dabigatran as about 75% of inhibitory percentage compared to healthy control.

CONCLUSION

Overall, our computational and experimental results clearly show the higher inhibitory effect of dabigatran than rivaroxaban. Hence, rivaroxaban will be a better drug candidate for improving the outcome of AF.

Modeling membranes in situ

Molecular dynamics simulations of cellular membranes have come a long way-from simple model lipid bilayers to multicomponent systems capturing the crowded and complex nature of real cell membranes. In this opinionated minireview, we discuss the current challenge to simulate the dynamics of membranes in their native environment, in situ, with the prospect of reaching the level of whole cells and cell organelles using an integrative modeling framework.

Disruption of Transmembrane Phosphatidylserine Asymmetry by HIV-1 Incorporated SERINC5 Is Not Responsible for Virus Restriction

Host restriction factor SERINC5 (SER5) incorporates into the HIV-1 membrane and inhibits infectivity by a poorly understood mechanism. Recently, SER5 was found to exhibit scramblase-like activity leading to the externalization of phosphatidylserine (PS) on the viral surface, which has been proposed to be responsible for SER5's antiviral activity. This and other reports that document modulation of HIV-1 infectivity by viral lipid composition prompted us to investigate the role of PS in regulating SER5-mediated HIV-1 restriction. First, we show that the level of SER5 incorporation into virions correlates with an increase in PS levels in the outer leaflet of the viral membrane. We developed an assay to estimate the PS distribution across the viral membrane and found that SER5, but not SER2, which lacks antiviral activity, abrogates PS asymmetry by externalizing this lipid. Second, SER5 incorporation diminished the infectivity of pseudoviruses produced from cells lacking a flippase subunit CDC50a and, therefore, exhibited a higher baseline level of surface-accessible PS. Finally, exogenous manipulation of the viral PS levels utilizing methyl-alpha-cyclodextrin revealed a lack of correlation between external PS and virion infectivity. Taken together, our study implies that the increased PS exposure to SER5-containing virions itself is not directly linked to HIV-1 restriction.

The structure and physical properties of a packaged bacteriophage particle

A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.

pyCapsid: identifying dominant dynamics and quasi-rigid mechanical units in protein shells

SUMMARY

pyCapsid is a Python package developed to facilitate the characterization of the dynamics and quasi-rigid mechanical units of protein shells and other protein complexes. The package was developed in response to the rapid increase of high-resolution structures, particularly capsids of viruses, requiring multiscale biophysical analyses. Given a protein shell, pyCapsid generates the collective vibrations of its amino-acid residues, identifies quasi-rigid mechanical regions associated with the disassembly of the structure, and maps the results back to the input proteins for interpretation. pyCapsid summarizes the main results in a report that includes publication-quality figures.

AVAILABILITY AND IMPLEMENTATION

pyCapsid's source code is available under MIT License on GitHub. It is compatible with Python 3.8-3.10 and has been deployed in two leading Python package-management systems, PIP and Conda. Installation instructions and tutorials are available in the online documentation and in the pyCapsid's YouTube playlist. In addition, a cloud-based implementation of pyCapsid is available as a Google Colab notebook. pyCapsid Colab does not require installation and generates the same report and outputs as the installable version. Users can post issues regarding pyCapsid in the repository's issues section.

Correlation between chemical denaturation and the unfolding energetics of Acanthamoeba actophorin

The actin filament network is in part remodeled by the action of a family of filament severing proteins that are responsible for modulating the ratio between monomeric and filamentous actin. Recent work on the protein actophorin from the amoeba Acanthamoeba castellani identified a series of site-directed mutations that increase the thermal stability of the protein by 22°C. Here, we expand this observation by showing that the mutant protein is also significantly stable to both equilibrium and kinetic chemical denaturation, and employ computer simulations to account for the increase in thermal or chemical stability through an accounting of atomic-level interactions. Specifically, the potential of mean force (PMF) can be obtained from steered molecular dynamics (SMD) simulations in which a protein is unfolded. However, SMD can be inefficient for large proteins as they require large solvent boxes, and computationally expensive as they require increasingly many SMD trajectories to converge the PMF. Adaptive steered molecular dynamics (ASMD) overcomes the second of these limitations by steering the particle in stages, which allows for convergence of the PMF using fewer trajectories compared with SMD. Use of the telescoping water scheme within ASMD partially overcomes the first of these limitations by reducing the number of waters at each stage to only those needed to solvate the structure within a given stage. In the PMFs obtained from ASMD, the work of unfolding Acto-2 was found to be higher than the Acto-WT by approximately 120 kCal/mol and reflects the increased stability seen in the chemical denaturation experiments. The evolution of the average number of hydrogen bonds and number of salt bridges during the pulling process provides a mechanistic view of the structural changes of the actophorin protein as it is unfolded, and how it is affected by the mutation in concert with the energetics reported through the PMF.

Understanding Virus Structure and Dynamics through Molecular Simulations

Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.

Rheumatoid arthritis treatment with zoledronic acid, a potentialinhibitorofGWAS-derived pharmacogenetics STAT3 and IL2 targets

Rheumatoid arthritis (RA) is an inflammatory condition that primarily affects the joints and progress to affect other vital organs. Variety of drugs are being recommended to control the disease progression that benefits patients to perform day-to-day activities. Few of these RA drugs have noticeable side effects; therefore, it's crucial to choose the appropriate drug for treating RA with an understanding of the disease's pathophysiology. Herein, we investigated the RA genes from GWAS data to construct protein-protein interaction (PPI) network and to define appropriate drug targets for RA. The predicted drug targets were screened with the known RA drugs based on molecular docking. Further, the molecular dynamics simulations were performed to comprehend the conformational changes and stability of the targets upon binding of the selected top ranked RA drug. As a result, our constructed protein network from GWAS data revealed, STAT3 and IL2 could be potential pharmacogenetics targets that interlink most of the RA genes encoding proteins. These interlinked proteins of both the targets showed involvement in cell signaling, immune response, and TNF signaling pathway. Among the 192 RA drugs investigated, zoledronic acid had the lowest binding energy that inhibit both STAT3 (-6.307 kcal/mol) and IL2 (-6.231 kcal/mol). Additionally, STAT3 and IL2 trajectories on zoledronic acid binding exhibit notable differences in MD simulations as compared to a drug-free environment. Also, the in vitro assessment with the zoledronic acid confirms the outcome of our computational study. Overall, our study identify zoledronic acid could be potential inhibitor against these targets, that will benefits patients with RA. Comparative efficiency assessments between the RA drugs through clinical trials are needed to validate our findings in the treatment of RA.

Implementation of Telescoping Boxes in Adaptive Steered Molecular Dynamics

Long-time dynamical processes, such as those involving protein unfolding and ligand interactions, can be accelerated and realized through steered molecular dynamics (SMD). The challenge has been the extraction of information from such simulations that generalize for complex nonequilibrium processes. The use of Jarzynski's equality opened the possibility of determining the free energy along the steered coordinate, but sampling over the nonequilibrium trajectories is slow to converge. Adaptive steered molecular dynamics (ASMD) and other related techniques have been introduced to overcome this challenge through the use of stages. Here, we take advantage of these stages to address the numerical cost that arises from the required use of very large solvent boxes. We introduce telescoping box schemes within adaptive steered molecular dynamics (ASMD) in which we adjust the solvent box between stages and thereby vary (and optimize) the required number of solvent molecules. We have benchmarked the method on a relatively long α-helical peptide, Ala₃₀, with respect to the potential of mean force and hydrogen bonds. We show that the use of telescoping boxes introduces little numerical error while significantly reducing the computational cost.