Engineering topics

Lopez Bautista, Cesar Augusto

Publications and source records attributed to Lopez Bautista, Cesar Augusto.

6 records

DOE OSTI2024

In Vitro Evolution and Computational Approaches to Predict, Prevent and Control Future Pandemics

The natural tendency of virus to mutate and the under-sampling of the environment makes it difficult to become aware of the emergence of new viral strains with pandemic potential. Being able to predict what mutations make a virus more infective might allow to spot such strains with minimal sampling and potentially allow to predict/prevent the next pandemic. The team attempted to mimic natural viral mutations and recombination through computational and experimental methods producing a variety of mutants of a SARS-COV-2 protein (receptor binding domain, RBD, of spike protein) responsible for viral entry in mammalian cells. The library of mutants was then interrogated for ability and lack-there-of to interact with the host cell receptor ushering viral entry, Angiotensin-converting enzyme 2 (ACE2). The negative and positive data set is intended to “teach the rules” of virus-host receptor interaction. Additionally, the positive clones were used to screen a set of antibody mutants designed to widen the breadth of viral mutants recognition, to demonstrate that this kind of libraries could also be a tool to produce antibody therapeutics impervious to viral mutation, even before a pandemic strain is discovered.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

HIV Env assembling is tightly modulated by lipid composition and cytoplasmic tail integrity

HIV envelope glycoprotein (Env) is a homotrimeric transmembrane protein which binds to receptors CD4, CCR5- CXCR4, leading to a cascade of conformational changes ultimately resulting in the entrance of the virus into the host cell. Its whole structure comprises an extracellular region, a transmembrane domain and a cytoplasmic region. While many efforts are carried in order to develop an effective vaccine specifically targeting the extracellular domain, less efforts are focus into regions close to the transmembrane region; mostly limited by the lack of structural knowledge within this domain. Thus, it is desirable to have a molecular description of the trimeric association of the transmembrane domain, as well as the membrane properties affecting their structure. We have carried molecular dynamic (MD) simulations of HIV transmembrane region in different lipid environments, mimicking the complex viral membrane. Our results indicate that the lipid environment heavily affects the dynamics and association of the transmembrane helices, modulating its structure. In addition, we have investigated the effect of the cytoplasmic tail (CT). Unexpectedly, this region has a tremendous effect on the lateral organization of the transmembrane. These results clearly highlights the importance of this region in the overall structure of the Env protein and point towards its key role for the in-silico design of effective vaccines towards this domain.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. In this work, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (~98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ~1.5 million atoms is discussed in the context of MuMMI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Large Scale MD to Predict Epitope Regions in HIV Env [Slides]

Highly dense carbohydrates located on the surface of the HIV Env protein play a key role in immune evasion. Such evolutionary adaptation hampers any attempt to obtain a full mechanistic understanding of the role played by the glycans in protecting the virus against an effective immune response. Moreover, and due to their chemical variability, an accurate molecular understanding of the so called “glycan shield” is still limited by the lack of effective resolution of state-of-the-art experimental technics. Here, we have used extensive computational modelling in order to fill this gap, addressing the presence of a large glycan variability as observed experimentally. Based on an automated pipeline, we were able to assemble, set-up and simulate via Molecular dynamics hundreds of different glycosylated Env variants at nearly atomic resolution, recapitulating the glycosylation distributions observed experimentally. Results from these simulations were subjected to machine learning and very accurate prediction of simulation derived glycan shielding areas of each glycan as a function of static sequence features. Such predictive models of per-glycan shielding, incorporating both glycan dynamics and heterogeneity, were used to develop a novel sequence-based glycan shield mapping strategy. Parallel to these studies, we also developed an accurate machine learning approach to predict glycan heterogeneity data using sequence features and found good prediction accuracy.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Computational Engineering of Human Apolipoprotein for Enhanced Interaction with Bacterial Toxins

Sepsis is caused by an over-activation of the immune system in response to the presence of bacterial endotoxins in the bloodstream. LPS exemplifies a potent immune stimulator which can lead to an excessive activation leading to an overly pro-inflammatory reaction and if left untreated results in multi-organ failure and death. Human lipoproteins (HLP) have been shown to neutralize LPS and facilitate its clearance from circulation via the liver. However, host-derived lipoproteins are inhomogeneous in content and size, presenting a challenge for therapeutic applications. In this work, we have used IC facilities in order to computationally engineer HLP for an enhanced interaction with LPS. This work is a continuation of our previous efforts to derive an optimal Nanodisc (ND) lipid mixture for enhanced association with the endotoxin. We have performed fully atomistic Molecular dynamic (MD) simulations in order to carefully identify the most optimal residues in HLP, which can potentially lead to a stronger association to LPS.

42 ENGINEERING↗

DOE OSTI2021

Development of Martini 2.2 parameters for N -glycans: a case study of the HIV-1 Env glycoprotein dynamics

N-linked glycans are ubiquitous in nature and play key roles in biology. For example, glycosylation of pathogenic proteins is a common immune evasive mechanism, hampering the development of successful vaccines. Due to their chemical variability and complex dynamics, an accurate molecular understanding of glycans is still limited by the lack of effective resolution of current experimental approaches. Here, we have developed and implemented a reductive model based on the popular Martini 2.2 coarse-grained force field for the computational study of N-glycosylation. We used the HIV-1 Env as a direct applied example of a highly glycosylated protein. Finally, our results indicate that the model not only reproduces many observables in very good agreement with a fully atomistic force field but also can be extended to study large amount of glycosylation variants, a fundamental property that can aid in the development of drugs and vaccines.

59 BASIC BIOLOGICAL SCIENCES↗