Biomedical subjects

Roberto D Lins

Publications and source records attributed to Roberto D Lins.

9 recordsLinked to original sources

PubMed2025

Design of nanobody targeting SARS-CoV-2 spike glycoprotein using CDR-grafting assisted by molecular simulation and machine learning.

The design of proteins capable effectively binding to specific protein targets is crucial for developing therapies, diagnostics, and vaccine candidates for viral infections. Here, we introduce a complementarity-determining region (CDR) grafting approach for designing nanobodies (Nbs) that target specific epitopes, with the aid of computer simulation and machine learning. As a proof-of-concept, we designed, evaluated, and characterized a high-affinity Nb against the spike protein of SARS-CoV-2, the causative agent of the COVID-19 pandemic. The designed Nb, referred to as Nb Ab.2, was synthesized and displayed high-affinity for both the purified receptor-binding domain protein and to the virus-like particle, demonstrating affinities of 9 nM and 60 nM, respectively, as measured with microscale thermophoresis. Circular dichroism showed the designed protein's structural integrity and its proper folding, whereas molecular dynamics simulations provided insights into the internal dynamics of Nb Ab.2. This study shows that our computational pipeline can be used to efficiently design high-affinity Nbs with diagnostic and prophylactic potential, which can be tailored to tackle different viral targets.

Spike Glycoprotein, Coronavirus↗

PubMed2006

Protein under pressure: molecular dynamics simulation of the arc repressor.

Experimental nuclear magnetic resonance results for the Arc Repressor have shown that this dimeric protein dissociates into a molten globule at high pressure. This structural change is accompanied by a modification of the hydrogen-bonding pattern of the intermolecular beta-sheet: it changes its character from intermolecular to intramolecular with respect to the two monomers. Molecular dynamics simulations of the Arc Repressor, as a monomer and a dimer, at elevated pressure have been performed with the aim to study this hypothesis and to identify the major structural and dynamical changes of the protein under such conditions. The monomer appears less stable than the dimer. However, the complete dissociation has not been seen because of the long timescale needed to observe this phenomenon. In fact, the protein structure altered very little when increasing the pressure. It became slightly compressed and the dynamics of the side-chains and the unfolding process slowed down. Increasing both, temperature and pressure, a tendency of conversion of intermolecular into intramolecular hydrogen bonds in the beta-sheet region has been detected, supporting the mentioned hypothesis. Also, the onset of denaturation of the separated chains was observed.

Amino Acid Sequence↗

PubMed2005

Variational particle number approach for rational compound design.

Within density functional theory, a variational particle number approach for rational compound design (RCD) is presented. An expression for RCD is obtained in terms of minimization of a suitably defined energy penalty functional whose gradients are the nuclear and the electronic chemical potential. Using combined quantum and molecular mechanics, a nonpeptidic anticancer drug candidate is designed.

Journal Article↗

PubMed2005

An improved nucleic acid parameter set for the GROMOS force field.

Over the past decades, the GROMOS force field for biomolecular simulation has primarily been developed for performing molecular dynamics (MD) simulations of polypeptides and, to a lesser extent, sugars. When applied to DNA, the 43A1 and 45A3 parameter sets of the years 1996 and 2001 produced rather flexible double-helical structures, in which the Watson-Crick hydrogen-bonding content was more limited than expected. To improve on the currently available parameter sets, the nucleotide backbone torsional-angle parameters and the charge distribution of the nucleotide bases are reconsidered based on quantum-chemical data. The new 45A4 parameter set resulting from this refinement appears to perform well in terms of reproducing solution NMR data and canonical hydrogen bonding. The deviation between simulated and experimental observables is now of the same order of magnitude as the uncertainty in the experimental values themselves.

Algorithms↗

PubMed2005

A new GROMOS force field for hexopyranose-based carbohydrates.

A new parameter set (referred to as 45A4) is developed for the explicit-solvent simulation of hexopyranose-based carbohydrates. This set is compatible with the most recent version of the GROMOS force field for proteins, nucleic acids, and lipids, and the SPC water model. The parametrization procedure relies on: (1) reassigning the atomic partial charges based on a fit to the quantum-mechanical electrostatic potential around a trisaccharide; (2) refining the torsional potential parameters associated with the rotations of the hydroxymethyl, hydroxyl, and anomeric alkoxy groups by fitting to corresponding quantum-mechanical profiles for hexopyranosides; (3) adapting the torsional potential parameters determining the ring conformation so as to stabilize the (experimentally predominant) (4)C(1) chair conformation. The other (van der Waals and nontorsional covalent) parameters and the rules for third and excluded neighbors are taken directly from the most recent version of the GROMOS force field (except for one additional exclusion). The new set is general enough to define parameters for any (unbranched) hexopyranose-based mono-, di-, oligo- or polysaccharide. In the present article, this force field is validated for a limited set of monosaccharides (alpha- and beta-D-glucose, alpha- and beta-D-galactose) and disaccharides (trehalose, maltose, and cellobiose) in solution, by comparing the results of simulations to available experimental data. More extensive validation will be the scope of a forthcoming article. (c) 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1400-1412, 2005.

Carbohydrate Conformation↗

PubMed2004

Trehalose-protein interaction in aqueous solution.

A variety of sugars are known to enhance the stability of biomaterials. Trehalose, a nonreducing disaccharide composed of two alpha, alpha(1 --> 1)-linked D-glucopyranose units, appears to be one of the most effective protectants. Both in vivo and in vitro, trehalose protects biostructures such as proteins and membranes from damage due to dehydration, heat, or cold. However, despite the significant amount of experimental data on this disaccharide, no clear picture of the molecular mechanism responsible for its stabilizing properties has emerged yet. Three major hypotheses (water-trehalose hydrogen-bond replacement, coating by a trapped water layer, and mechanical inhibition of the conformational fluctuations) have been proposed to explain the stabilizing effect of trehalose on proteins. To investigate the nature of protein-trehalose-water interactions in solution at the molecular level, two molecular dynamics simulations of the protein lysozyme in solution at room temperature have been carried out, one in the presence (about 0.5 M) and one in the absence of trehalose. The results show that the trehalose molecules cluster and move toward the protein, but neither completely expel water from the protein surface nor form hydrogen bonds with the protein. Furthermore, the coating by trehalose does not significantly reduce the conformational fluctuations of the protein compared to the trehalose-free system. Based on these observations, a model is proposed for the interaction of trehalose molecules with a protein in moderately concentrated solutions, at room temperature and on the nanosecond timescale.

Computer Simulation↗

PubMed2004

Prediction of HIV-1 integrase/viral DNA interactions in the catalytic domain by fast molecular docking.

This study details the separate analyses of binding specificity of HIV-1 integrase (IN) and viral B-DNA forms through ligand-receptor docking studies by means of a fast molecular docking method. The application of solvated electrostatics with the University of Houston Brownian Dynamics Program (UHBD) and configurational sampling by the Daughter of Turnip (DOT) docking program resulted in the computation of energies of more than 113 billion configurations for each ligand-receptor docking study, a procedure considered computationally intractable a few years ago. A specific binding pattern of viral DNA to the IN catalytic domain region has been predicted as a result of these calculations. In a representative docked configuration, we observe the 3'-hydroxyl of the conserved deoxyadenosine to be close to one of the two divalent metal ions that are necessary for catalysis. A superimposition of our energy-minimized docked complex on representative structures from a molecular dynamics (MD) simulation of a crystallographically resolved IN/inhibitor complex revealed an overlap of viral DNA with the inhibitor, indicating that the bound inhibitor might operate by blocking substrate binding. The DOT docking calculation also identified a second, adjacent DNA-binding site, which we believe is the nonspecific host DNA binding site. The binding pattern predicted by DOT complements previous electrostatics, MD simulation, photo-cross-linking, and mutagenesis studies and also provides a further refinement of the IN/viral DNA binding interaction as a basis for new structure-based design efforts.

Catalytic Domain↗

PubMed2004

Interaction of the disaccharide trehalose with a phospholipid bilayer: a molecular dynamics study.

The disaccharide trehalose is well known for its bioprotective properties. Produced in large amounts during stress periods in the life of organisms able to survive potentially damaging conditions, trehalose plays its protective role by stabilizing biostructures such as proteins and lipid membranes. In this study, molecular dynamics simulations are used to investigate the interaction of trehalose with a phospholipid bilayer at atomistic resolution. Simulations of the bilayer in the absence and in the presence of trehalose at two different concentrations (1 or 2 molal) are carried out at 325 K and 475 K. The results show that trehalose is able to minimize the disruptive effect of the elevated temperature and stabilize the bilayer structure. At both temperature, trehalose is found to interact directly with the bilayer through hydrogen bonds. However, the water molecules at the bilayer surface are not completely replaced. At high temperature, the protective effect of trehalose is correlated with a significant increase in the number of trehalose-bilayer hydrogen bonds, predominantly through an increase in the number of trehalose molecules bridging three or more lipid molecules.

Chemical Phenomena↗

PubMed2003

A consistent potential energy parameter set for lipids: dipalmitoylphosphatidylcholine as a benchmark of the GROMOS96 45A3 force field.

The performance of the GROMOS96 parameter set 45A3 developed for aliphatic alkanes is tested on a bilayer of dipalmitoylphosphatidylcholine (DPPC) in water in the liquid-crystalline L(alpha) phase. Variants of the force-field parameter set as well as different sets of simulation conditions or simulation parameter sets are evaluated. In the case of the force-field parameters, the van der Waals constants for the non-bonded interaction of the ester carbonyl carbon and the partial charges and charge group definition of the phosphatidylcholine head group are examined. On the methodological side, different cut-off distances for the non-bonded interactions, use of a reaction-field force due to long-range electrostatic interactions, the frequency of removal of the centre of mass motion and the strength of the coupling of the pressure of the system to the pressure bath are tested. The area per lipid, as a measure of structure, the order parameters of the chain carbons, as a measure of membrane fluidity, and the translational diffusion of the lipids in the plane of the bilayer are calculated and compared with experimental values. An optimal set of simulation parameters for which the GROMOS96 parameter set 45A3 yields a head group area, chain order parameters and a lateral diffusion coefficient in accordance with the experimental data is listed.

1,2-Dipalmitoylphosphatidylcholine↗