Within this project, an experimental pdb-structure for a protease-peptide complex was used to build a structural model for a protease-dimer bound to a CA-SP1 site within the hexameric, immature HIV Gag-lattice. This represents one of the very first structural models for how proteases access the Gag-lattice for subsequent cleavage. Performing atomistic simulations starting from this complex structure, we have monitored that protease-binding induces structural changes in the lattice but not significant remodeling. To better understand large-scale remodeling as it occurs during HIV maturation, we have developed a coarse-grained model for this system. One of the main challenges here is to accurately describe affinities between different CA-SP1 monomers as well as between protease and CA-SP1, especially since binding/unbinding-events are difficult to sample for such large systems. We have resolved this issue by developing a new machine-learning method that represents a regularized version of relative entropy minimization. This new approach includes a quadratic bias on an empirical binding-affinity in its loss-function. In this way, it allows for accurate description of the complex structure as well as for the binding affinity. The method has been designed to guarantee maximum practicality, which means that no additional terms in the coarse-grained force fields or implementations in the simulation engine are required. Furthermore, we have demonstrated that hexameric CA-SP1 unit-models are transferable to simulate an entire immature HIV lattice. To simulate large-scale immature lattice remodeling triggered by proteolytic cleavages, we have developed a simulation routine in openmm that includes a CA-SP1 cleavage mechanism in the presence of protease. Large-scale simulations have been performed with this new model, yielding fundamental new insight into how proteolytic cleavages drive remodeling of the immature CA-SP1 lattice. In conclusion, these simulations support a model of delayed outlier formation (including dissociation of cleaved particles) after proteolytic cleavage expanding from lattice edges and gaps. Thus, the simulations align with the hypothesis of a combination of displacive remodeling and partial disassembly rather than a uniform conversion mechanism for HIV maturation.
Development of a new coarse-grained model for RNA has revealed a dominant impact of sugar puckering on the thermodynamic behavior of RNA. Moreover, simulations performed by us show critical misdescription of sugar puckering also in atomistic force fields. For this reason, we have devised a new machine learning workflow to optimize sugar puckering in RNA based on experimental (NMR) data. Optimized atomistic and coarse-grained models will allow for future merging in a hybrid model.