Crucial roles of the subnanosecond local dynamics of the flap tips in the global conformational changes of HIV-1 protease

Energetics of Flap Opening in HIV-1 Protease: String Method Calculations

HIV-1 protease (PR) is the viral protein responsible for virion maturation, and its mechanisms of action remain incompletely understood. PR is dimeric and contains two flexible, symmetry-related flaps, which act as a gate to inhibit access to the binding pocket and hold the polypeptide substrate in the binding pocket once bound. Wide flap opening, a conformational change assumed to be necessary for substrate binding, is a rare event in the closed and bound form. In this study, we use molecular dynamics (MD) simulations and advanced MD techniques including temperature acceleration and string method in collective variables to study the conformational changes associated with substrate unbinding of both wild-type and F99Y mutant PR. The F99Y mutation is shown via MD to decouple the closing of previously unrecognized distal pockets from substrate unbinding. To determine whether or not the F99Y mutation affects the energetic cost of wide flap opening, we use string method in collective variables to determine the minimum free-energy mechanism for wide flap opening in concert with distal pocket closing. The results indicate that the major energetic cost in flap opening is disengagement of the two flap-tip Ile50 residues from each other and is not affected by the F99Y mutation.

Non-covalent loading of ionic liquid-functionalized nanoparticles for bovine serum albumin: experiments and theoretical analysis

Biomacromolecule-based nanomaterials have attracted much attention due to their excellent function in sensing, catalysis, medicine, biology and recognition. In this work, a silane-coupling ionic liquid, 1-(3-trimethoxysilylpropyl)-3-methylimidazolium chloride ([TMIM]Cl), was synthesized and applied to prepare ionic liquid-functionalized nanoparticles (SiO₂@IL) using surface grafting technology. By employing multiple non-covalent interactions, including electrostatic interactions, hydrogen bonding and π–π stacking, the obtained functional nanoparticles were able to bind bovine serum albumin (BSA) with strong binding affinity, which has been illustrated through experiments and theoretical calculations. Moreover, the stability of SiO₂@IL further demonstrated that it is promising in applications for biomacromolecule immobilization.

Non-covalent binding between nanosilica and bovine serum albumin has been illustrated by experiments and theoretical calculations.

N-Terminal Acetylation Preserves α-Synuclein from Oligomerization by Blocking Intermolecular Hydrogen Bonds

The abnormal aggregation of α-synuclein (α-Syn) is closely associated with Parkinson's disease. Different post-translational modifications of α-Syn have been identified and contribute distinctly in α-Syn aggregation and cytotoxicity. Recently, α-Syn was reported to be N-terminally acetylated in cells, yet the functional implication of this modification, especially in α-Syn oligomerization, remains unclear. By using a solid-state nanopore system, we found that N-terminal acetylation can significantly decrease α-Syn oligomerization. Replica-exchange molecular dynamics simulations further revealed that addition of an acetyl group at the N-terminus disrupts intermolecular hydrogen bonds, which slows down the initial α-Syn oligomerization. Our finding highlights the essential role of N-terminal acetylation of α-Syn in preserving its native conformation against pathological aggregation.

From two-dimensional nano-sheets to roll-up structures: expanding the family of nanoscroll

The carbon nanoscroll (CNS), a spiral structure rolled up from a flat graphene sheet, has attracted significant attention due to its large surface area and broad electrochemical window. Despite the rapid emergence of many new two-dimensional (2D) materials, studies on nanoscrolls beyond CNS, however, remain limited. In this work, we combine first-principles calculations and theoretical analyses to determine the stable configuration of nanoscrolls constructed from several 2D materials, including graphene, h-BN, MoS₂, graphyne, phosphorene, and graphene oxide. Our study shows that the inner radius of these nanoscrolls ranges from 0.36-2.3 nm, significantly expanding that of CNS (about 1.2 nm). In addition, it is shown that the difference between the inner and outer pressures of nanoscrolls can be used to tune their core sizes. Furthermore, the effect of interlayer friction is found to be significant and may play a dominant role in determining the core sizes. The present study not only expands the nanoscroll family, but also enriches the understanding of nanoscrolls in terms of the effects of pressure and friction, thus widening the application range of nanoscrolls as nano-actuators, ion/water channels, nano-sensors, and hydrogen/energy storage devices, etc.

High Transmembrane Voltage Raised by Close Contact Initiates Fusion Pore

Membrane fusion lies at the heart of neuronal communication but the detailed mechanism of a critical step, fusion pore initiation, remains poorly understood. Here, through atomistic molecular dynamics simulations, a transient pore formation induced by a close contact of two apposed bilayers is firstly reported. Such a close contact gives rise to a high local transmembrane voltage that induces the transient pore formation. Through simulations on two apposed bilayers fixed at a series of given distances, the process in which two bilayers approaching to each other under the pulling force from fusion proteins for membrane fusion was mimicked. Of note, this close contact induced fusion pore formation is contrasted with previous reported electroporation under ad hoc applied external electric field or ionic charge in-balance. We show that the transmembrane voltage increases with the decrease of the distance between the bilayers. Below a critical distance, depending on the lipid composition, the local transmembrane voltage can be sufficiently high to induce the transient pores. The size of these pores is approximately 1~2 nm in diameter, which is large enough to allow passing of neurotransmitters. A resealing of the membrane pores resulting from the neutralization of the transmembrane voltage by ions through the pores was then observed. We also found that the membrane tension can either prolong the lifetime of transient pores or cause them to dilate for full collapse. This result provides a possible mechanism for fusion pore formation and regulation of pathway of fusion process.

Mapping the Dynamics Landscape of Conformational Transitions in Enzyme: The Adenylate Kinase Case

Conformational transition describes the essential dynamics and mechanism of enzymes in pursuing their various functions. The fundamental and practical challenge to researchers is to quantitatively describe the roles of large-scale dynamic transitions for regulating the catalytic processes. In this study, we tackled this challenge by exploring the pathways and free energy landscape of conformational changes in adenylate kinase (AdK), a key ubiquitous enzyme for cellular energy homeostasis. Using explicit long-timescale (up to microseconds) molecular dynamics and bias-exchange metadynamics simulations, we determined at the atomistic level the intermediate conformational states and mapped the transition pathways of AdK in the presence and absence of ligands. There is clearly chronological operation of the functional domains of AdK. Specifically in the ligand-free AdK, there is no significant energy barrier in the free energy landscape separating the open and closed states. Instead there are multiple intermediate conformational states, which facilitate the rapid transitions of AdK. In the ligand-bound AdK, the closed conformation is energetically most favored with a large energy barrier to open it up, and the conformational population prefers to shift to the closed form coupled with transitions. The results suggest a perspective for a hybrid of conformational selection and induced fit operations of ligand binding to AdK. These observations, depicted in the most comprehensive and quantitative way to date, to our knowledge, emphasize the underlying intrinsic dynamics of AdK and reveal the sophisticated conformational transitions of AdK in fulfilling its enzymatic functions. The developed methodology can also apply to other proteins and biomolecular systems.

Pulsed EPR characterization of HIV-1 protease conformational sampling and inhibitor-induced population shifts

The conformational landscape of HIV-1 protease (PR) can be experimentally characterized by pulsed-EPR double electron-electron resonance (DEER). For this characterization, nitroxide spin labels are attached to an engineered cysteine residue in the flap region of HIV-1 PR. DEER distance measurements from spin-labels contained within each flap of the homodimer provide a detailed description of the conformational sampling of apo-enzyme as well as induced conformational shifts as a function of inhibitor binding. The distance distribution profiles are further interpreted in terms of a conformational ensemble scheme that consists of four unique states termed "curled/tucked", "closed", "semi-open" and "wide-open" conformations. Reported here are the DEER results for a drug-resistant variant clinical isolate sequence, V6, in the presence of FDA approved protease inhibitors (PIs) as well as a non-hydrolyzable substrate mimic, CaP2. Results are interpreted in the context of the current understanding of the relationship between conformational sampling, drug resistance, and kinetic efficiency of HIV-1PR as derived from previous DEER and kinetic data for a series of HIV-1PR constructs that contain drug-pressure selected mutations or natural polymorphisms. Specifically, these collective results support the notion that inhibitor-induced closure of the flaps correlates with inhibitor efficiency and drug resistance. This body of work also suggests DEER as a tool for studying conformational sampling in flexible enzymes as it relates to function.

Startling temperature effect on proteins when confined: single molecular level behaviour of human serum albumin in a reverse micelle

The present work reports the effect of confinement, and temperature therein, on the conformational fluctuation dynamics of domain-I of human serum albumin (HSA) by fluorescence correlation spectroscopy (FCS).

Peptide-Graphene Interactions Enhance the Mechanical Properties of Silk Fibroin

Studies reveal that biomolecules can form intriguing molecular structures with fascinating functionalities upon interaction with graphene. Then, interesting questions arise. How does silk fibroin interact with graphene? Does such interaction lead to an enhancement in its mechanical properties? In this study, using large-scale molecular dynamics simulations, we first examine the interaction of graphene with several typical peptide structures of silk fibroin extracted from different domains of silk fibroin, including pure amorphous (P1), pure crystalline (P2), a segment from N-terminal (P3), and a combined amorphous and crystalline segment (P4), aiming to reveal their structural modifications. Our study shows that graphene can have intriguing influences on the structures formed by the peptides with sequences representing different domains of silk fibroin. In general, for protein domains with stable structure and strong intramolecular interaction (e.g., β-sheets), graphene tends to compete with the intramolecular interactions and thus weaken the interchain interaction and reduce the contents of β-sheets. For the silk domains with random or less ordered secondary structures and weak intramolecular interactions, graphene tends to enhance the stability of peptide structures; in particular, it increases the contents of helical structures. Thereafter, tensile simulations were further performed on the representative peptides to investigate how such structure modifications affect their mechanical properties. It was found that the strength and resilience of the peptides are enhanced through their interaction with graphene. The present work reveals interesting insights into the interactions between silk peptides and graphene, and contributes in the efforts to enhance the mechanical properties of silk fibroin.

On the strength of β-sheet crystallites of Bombyx mori silk fibroin

Silk fibroin, a natural multi-domain protein, has attracted great attention due to its superior mechanical properties such as ultra-high strength and stretchability, biocompatibility, as well as its versatile biodegradability and processability. It is mainly composed of β-sheet crystallites and amorphous domains. Although its strength is well known to be controlled by the dissociation of protein chains from β-sheet crystallites, the way that water as the solvent affects its strength and the reason that its theoretically predicted strength is several times higher than experimental measurement remain unclear. We perform all-atom molecular dynamics simulations on a β-sheet crystallite of Bombyx mori silk. We find that water solvent reduces the number and strength of hydrogen bonds between β-chains, and thus greatly weakens the strength of silk fibroin. By dissociating protein chains at different locations from the crystallite, we also find that the pulling strength for the interior chains is several times higher than that for the surface/corner chains, with the former being consistent with the theoretically predicted value, while the latter on par with the experimental value. It is shown that the weakest rupture strength controls the failure strength of silk fibre. Hence, this work sheds light on the role of water in the strength of silk fibroin and also provides clues on the origin of the strength difference between theory and experiment.

Toward in silico biomolecular manipulation through static modes: atomic scale characterization of HIV-1 protease flexibility

Probing biomolecular flexibility with atomic-scale resolution is a challenging task in current computational biology for fundamental understanding and prediction of biomolecular interactions and associated functions. This paper makes use of the static mode method to study HIV-1 protease considered as a model system to investigate the full biomolecular flexibility at the atomic scale, the screening of active site biomechanical properties, the blind prediction of allosteric sites, and the design of multisite strategies to target deformations of interest. Relying on this single calculation run of static modes, we demonstrate that in silico predictive design of an infinite set of complex excitation fields is reachable, thanks to the storage of the static modes in a data bank that can be used to mimic single or multiatom contact and efficiently anticipate conformational changes arising from external stimuli. All along this article, we compare our results to data previously published and propose a guideline for efficient, predictive, and custom in silico experiments.

NMR relaxation in proteins with fast internal motions and slow conformational exchange: model-free framework and Markov state simulations

Calculating NMR relaxation effects for proteins with dynamics on multiple time scales generally requires very long trajectories based on conventional molecular dynamics simulations. In this report, we have built Markov state models from multiple MD trajectories and used the resulting MSM to capture the very fast internal motions of the protein within a free energy basin on a time scale up to hundreds of picoseconds and the more than 3 orders of magnitude slower conformational exchange between macrostates. To interpret the relaxation data, we derive new equations using the model-free framework which includes two slowly exchanging macrostates, each of which also exhibits fast local motions. Using simulations of HIV-1 protease as an example, we show how the populations of slowly exchanging conformational states as well as order parameters for the different states can be determined from the NMR relaxation data.

Identifying the molecular mechanics and binding dynamics characteristics of potent inhibitors to HIV-1 protease

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the primary inhibition targets for chemotherapy of AIDS because of its critical role in the replication cycle of the HIV. In this work, a combinatory coarse-grained and atomistic simulation method was developed for dissecting molecular mechanisms and binding process of inhibitors to the active site of HIV-1 PR, in which 35 typical inhibitors were trialed. We found that the molecular size and stiffness of the inhibitors and the binding energy between the inhibitors and PR play important roles in regulating the binding process. Comparatively, the smaller and more flexible inhibitors have larger binding energy and higher binding rates; they even bind into PR without opening the flaps. In contrast, the larger and stiffer inhibitors have lower binding energy and lower binding rate, and their binding is subjected to the opening and gating of the PR flaps. Furthermore, the components of binding free energy were quantified and analyzed by their dependence on the molecular size, structures, and hydrogen bond networks of inhibitors. Our results also deduce significant dynamics descriptors for determining the quantitative structure and property relationship in potent drug ligands for HIV-1 PR inhibition.

Insights into the dynamics of HIV-1 protease: a kinetic network model constructed from atomistic simulations

The conformational dynamics in the flaps of HIV-1 protease plays a crucial role in the mechanism of substrate binding. We develop a kinetic network model, constructed from detailed atomistic simulations, to determine the kinetic mechanisms of the conformational transitions in HIV-1 PR. To overcome the time scale limitation of conventional molecular dynamics (MD) simulations, our method combines replica exchange MD with transition path theory (TPT) to study the diversity and temperature dependence of the pathways connecting functionally important states of the protease. At low temperatures the large-scale flap opening is dominated by a small number of paths; at elevated temperatures the transition occurs through many structurally heterogeneous routes. The expanded conformation in the crystal structure 1TW7 is found to closely mimic a key intermediate in the flap-opening pathways at low temperature. We investigated the different transition mechanisms between the semi-open and closed forms. The calculated relaxation times reveal fast semi-open ↔ closed transitions, and infrequently the flaps fully open. The ligand binding rate predicted from this kinetic model increases by 38-fold from 285 to 309 K, which is in general agreement with experiments. To our knowledge, this is the first application of a network model constructed from atomistic simulations together with TPT to analyze conformational changes between different functional states of a natively folded protein.

Strength of hydrogen bond network takes crucial roles in the dissociation process of inhibitors from the HIV-1 protease binding pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions.

Structure-based design of carbon nanotubes as HIV-1 protease inhibitors: atomistic and coarse-grained simulations

Nanoparticles such as fullerenes and carbon nanotubes have been extensively studied for biomedical applications. In this paper, we report the design of carbon nanotubes as HIV-1 protease inhibitors. Docking and molecular dynamics calculations are performed using an atomistic model to explore the optimal interaction structure and free energy between the nanotube and HIV-1 protease. A coarse-grained model is then developed based on the atomistic model, allowing us to investigate the dynamic behaviors of the protease in the bound and unbound states. The dynamic process reveals that the carbon nanotube is able to bind to the active site of the protease and prevent the active flaps from opening up, thus blocking the function of the protease. This process is strongly influenced by the size of the nanotube. The binding of carbon nanotubes to an alternative binding site other than the active site is also explored. Therefore, carbon nanotube-based inhibitors have great potential for application as HIV-1 protease inhibitors.